Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 167,842:1 +T 71.1% coding (359/2244 nt) fhuA → ferrichrome outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913167,8421.T71.1% 238.8 / 100.1 97coding (359/2244 nt)fhuAferrichrome outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base . (17/11);  new base T (39/30);  total (56/41)
Fisher's exact test for biased strand distribution p-value = 8.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.70e-01

Only 100 of 101 total aligned reads displayed.

ACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGG  >  NC_000913/167697‑167988
                                                                                                                                                  |                                                                                                                                                  
aCCGATACGCCGATTCAAAAAGTGCCACCGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTAtt                                                                                                                                                  >  5:154526/1‑147 (MQ=255)
   gATACGCCTTTTCAAAAAGTGCCACAGTCTATTTCTGTTGTGACCGCCGACGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAATAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTAttgacc                                                                                                                                              <  2:136353/149‑6 (MQ=255)
        gCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCt                                                                                                                                           >  4:154575/1‑148 (MQ=255)
             ttCAAAAAGTGCCACAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATCa                                                                                                                                      >  2:372429/1‑148 (MQ=255)
             ttCAAAAAGTGCCACAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATCa                                                                                                                                      <  1:372429/148‑1 (MQ=255)
                aaaaaGTGCCACAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATCATTcg                                                                                                                                  >  8:398209/1‑149 (MQ=255)
                aaaaaGTGCCACAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTAt                                                                                                                                                   <  8:178711/131‑1 (MQ=255)
                aaaaaGTGCCACAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTAt                                                                                                                                                   >  7:178711/1‑131 (MQ=255)
                aaaaaGTGCCACAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTcgc                                                                                                                                 <  7:62416/149‑1 (MQ=255)
                 aaaaGTGCCACAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTcgcg                                                                                                                                >  6:70563/1‑149 (MQ=255)
                        ccACAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATCATTcgc                                                                                                                                 >  1:24892/1‑142 (MQ=255)
                        ccACAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATCATTcgc                                                                                                                                 <  2:24892/142‑1 (MQ=255)
                        ccACAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATCATTCGCGGCTTTg                                                                                                                          >  1:144327/1‑149 (MQ=255)
                        ccACAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATCATTCGCGGCTTTg                                                                                                                          >  7:332749/1‑149 (MQ=255)
                         cacaGTCTATTTCTTTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCAGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGTGTGTCTCTGTTGGTACGCGTGTCTCATCCAACACCTATTGACCACCTGATCATTCTCGGCTTTGc                                                                                                                         <  1:238499/149‑1 (MQ=255)
                         cacaGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATCATTCGCGGCTTTGc                                                                                                                         >  6:154218/1‑149 (MQ=255)
                         cacaGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATCATTCGCGGCTTTGc                                                                                                                         >  2:275123/1‑149 (MQ=255)
                         cacaGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATCATTCGCGGCTTTGc                                                                                                                         >  5:119107/1‑149 (MQ=255)
                         cacaGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCg                                                                                                                        >  5:449686/1‑149 (MQ=255)
                          acaGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGAt                                                                                                                                        >  5:205494/1‑133 (MQ=255)
                          acaGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATCATTCGCGGCTTTGCg                                                                                                                        >  8:278988/1‑149 (MQ=255)
                          acaGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTCTTGGTACGCGGGGAGCATCCAACACCTATTGACCACCTGATCATTCGCGGCTTTGCg                                                                                                                        >  7:456087/1‑149 (MQ=255)
                          acaGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCGGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGAt                                                                                                                                        <  6:205494/133‑1 (MQ=255)
                           caGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATCATTCGCGGCTTTGCgg                                                                                                                       <  3:154575/149‑1 (MQ=255)
                             gTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCgg                                                                                                                       >  7:223867/1‑146 (MQ=255)
                              tCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGAt                                                                                                                                        <  6:199330/129‑1 (MQ=255)
                              tCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGAt                                                                                                                                        >  5:199330/1‑129 (MQ=255)
                               cTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTATAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTAGGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTTCGGAAg                                                                                                                    >  4:28874/1‑147 (MQ=255)
                                tATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAg                                                                                                                 >  5:152901/1‑149 (MQ=255)
                                  tttCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCCCGGCTTTGCGGCAGaa                                                                                                                  <  6:449686/146‑1 (MQ=255)
                                     cTGTTGTGACTGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATCATTCGCGGCTTTGCgg                                                                                                                       <  2:25600/139‑1 (MQ=255)
                                     cTGTTGTGACTGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATCATTCGCGGCTTTGCgg                                                                                                                       >  1:25600/1‑139 (MQ=255)
                                     cTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCa                                                                                                             >  4:144975/1‑148 (MQ=255)
                                        ttgtGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCaaa                                                                                                           <  3:264796/148‑1 (MQ=255)
                                        ttgtGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCa                                                                                                             >  4:264796/1‑146 (MQ=255)
                                        ttgtGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAg                                                                                                          <  5:154218/149‑1 (MQ=255)
                                        ttgtGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGc                                                                                                         >  3:191729/1‑149 (MQ=255)
                                           tgACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCa                                                                                                       <  6:154526/149‑1 (MQ=255)
                                            gACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATc                                                                                                                                       <  5:69622/116‑1 (MQ=255)
                                            gACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATc                                                                                                                                       >  6:69622/1‑116 (MQ=255)
                                               cgccgAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATCATTCGCGGCttt                                                                                                                           >  6:69016/1‑125 (MQ=255)
                                               cgccgAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATCATTCGCGGCttt                                                                                                                           <  5:69016/125‑1 (MQ=255)
                                                 ccgAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATaa                                                                                                 <  6:119107/149‑1 (MQ=255)
                                                   gAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACt                                                                                               >  4:393247/1‑148 (MQ=255)
                                                     agagATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTAt                                                                                             <  2:34878/149‑1 (MQ=255)
                                                     agagATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATc                                                                                            <  4:191729/149‑1 (MQ=255)
                                                      gagaTGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATc                                                                                            <  2:144327/149‑1 (MQ=255)
                                                      gagaTGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATc                                                                                            <  8:332749/149‑1 (MQ=255)
                                                         aTGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTg                                                                                          >  1:405052/1‑148 (MQ=255)
                                                          tGGCGCTGCAGCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGaa                                                                                        <  1:275123/149‑1 (MQ=255)
                                                           ggCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGCATg                                                                                      >  2:194431/1‑149 (MQ=255)
                                                           ggCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTAAACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAAt                                                                                       >  5:154783/1‑149 (MQ=255)
                                                             cgcTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGc                                                                                    <  5:70563/149‑1 (MQ=255)
                                                                  cATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTCCTACGCGTGGCGAATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGaa                                                                               <  3:28874/149‑1 (MQ=255)
                                                                   aTCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTg                                                                                 >  4:139817/1‑147 (MQ=255)
                                                                   aTCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGaa                                                                               <  2:405052/149‑1 (MQ=255)
                                                                      aGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGcc                                                                                   <  1:471530/142‑1 (MQ=255)
                                                                      aGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGcc                                                                                   >  2:471530/1‑142 (MQ=255)
                                                                      aGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGcc                                                                                                        >  4:190374/1‑120 (MQ=255)
                                                                      aGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGcc                                                                                                        <  3:190374/120‑1 (MQ=255)
                                                                        ccGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTg                                                                           >  1:103672/1‑147 (MQ=255)
                                                                         cGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCa                                                                         <  7:398209/149‑1 (MQ=255)
                                                                         cGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAg                                                                        <  8:287978/149‑1 (MQ=255)
                                                                           aaGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAgg                                                                       <  3:139817/149‑1 (MQ=255)
                                                                             gTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGc                                                                     <  7:278988/149‑1 (MQ=255)
                                                                              tCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGc                                                                     >  6:215023/1‑148 (MQ=255)
                                                                                ggTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTGTCTGAATGGCCTGAAGTTGCAGGGCAACt                                                                 >  7:322942/1‑149 (MQ=255)
                                                                                  tAAAAGAAGCGCTTAGCTACCCGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACtt                                                                >  1:272951/1‑149 (MQ=255)
                                                                                   aaaaGAAGCGCTTAGCTACACGCAGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACtt                                                                >  3:94964/1‑148 (MQ=255)
                                                                                       gaagCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATaa                                                          <  6:152901/149‑1 (MQ=255)
                                                                                       gaagCGCTTAGCTACACGCCGGGTGGCTCTGTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATCATTCGCGGCGTTGCGGCAGAAGGCCTAAGCCAGATTATCTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCtata                                                           <  8:456087/149‑1 (MQ=255)
                                                                                         agCGCTTAGCTACACGCCGTGTGTCTCTGTTGGTACGCGTGGCGCATCGAACACCTATTGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATTGCCTGAAGTTGCAGGGCATCTTCTATAAc                                                         <  4:94964/149‑1 (MQ=255)
                                                                                         agCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGGAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAAc                                                         <  8:223867/148‑1 (MQ=255)
                                                                                              ttAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGc                                                    >  3:173866/1‑149 (MQ=255)
                                                                                                   tACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCa                                               >  2:206231/1‑149 (MQ=255)
                                                                                                    acacGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCAt                                              >  5:304400/1‑149 (MQ=255)
                                                                                                      acGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTg                                            >  2:199562/1‑149 (MQ=255)
                                                                                                       cGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGa                                           <  6:304400/149‑1 (MQ=255)
                                                                                                        gCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGAc                                          >  7:360645/1‑149 (MQ=255)
                                                                                                          cGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGAccc                                        <  1:199562/149‑1 (MQ=255)
                                                                                                                  tctGTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGc                                >  1:360541/1‑149 (MQ=255)
                                                                                                                  tctGTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGc                                >  8:63675/1‑149 (MQ=255)
                                                                                                                  tctGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCt                               >  2:218085/1‑149 (MQ=255)
                                                                                                                  tctGTTGGTACGCGTGGCGCATCCAACACCTATGACCAACTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGCATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCt                               >  5:283769/1‑149 (MQ=255)
                                                                                                                  tctGTTGGTACGCGTGGAGCATCCAACACCTATTGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGc                                >  6:194876/1‑149 (MQ=255)
                                                                                                                     gTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTg                              >  2:401018/1‑148 (MQ=255)
                                                                                                                     gTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTg                              >  6:263126/1‑148 (MQ=255)
                                                                                                                      ttGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTgg                             >  7:49435/1‑148 (MQ=255)
                                                                                                                      ttGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTgg                             <  6:154783/148‑1 (MQ=255)
                                                                                                                      ttGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTg                              <  4:173866/147‑1 (MQ=255)
                                                                                                                        ggTACGCGTGGCGCATCCAACACCTATTGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGaa                           <  5:215023/148‑1 (MQ=255)
                                                                                                                             gcgTGGCGCATCCAACACCTATTGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGAc                                          >  7:391507/1‑128 (MQ=255)
                                                                                                                             gcgTGGCGCATCCAACACCTATTGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGAc                                          <  8:391507/128‑1 (MQ=255)
                                                                                                                                        ccAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGt         >  7:43693/1‑149 (MQ=255)
                                                                                                                                         cAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCATGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGt         <  3:393247/148‑1 (MQ=255)
                                                                                                                                          aaCACCTATTGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTg        >  6:275450/1‑149 (MQ=255)
                                                                                                                                          aaCACCTATTGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTg        >  8:314169/1‑149 (MQ=255)
                                                                                                                                          aaCACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAAtt               >  8:59195/1‑141 (MQ=255)
                                                                                                                                          aaCACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAAtt               <  7:59195/141‑1 (MQ=255)
                                                                                                                                               cTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCgg  >  5:251074/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
ACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGG  >  NC_000913/167697‑167988

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: