Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,732,139 G→T 7.0% intergenic (+16/‑62) xylF → / → xylG D‑xylose transporter subunit/D‑xylose ABC transporter dual domain ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,732,1390GT 7.0% 411.8 / 13.2 136intergenic (+16/‑62)xylF/xylGD‑xylose transporter subunit/D‑xylose ABC transporter dual domain ATPase
Reads supporting (aligned to +/- strand):  ref base G (70/56);  new base T (4/6);  total (74/62)
Fisher's exact test for biased strand distribution p-value = 5.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 140 total aligned reads displayed.

CAAAAGCAGATACCACACTGAATAATGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATAACGTCTGCTTGCGGTTGAATGCT  >  NC_000913/3731990‑3732287
                                                                                                                                                     |                                                                                                                                                    
       aGATACCACACTGAATAATGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCg                                                                                                                                               >  2:293883/1‑150 (MQ=255)
        gATACCACACTGAATAATGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACTCCGATCGATGTGAATAAAAACAACATCAAAGATACGGTATTTTAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCgg                                                                                                                                              <  1:1812026/150‑1 (MQ=255)
          tACCACACTGAATAATGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGcggggcgggg                                                                                                                                            >  1:1895194/1‑150 (MQ=255)
           aCCACACTGAATAATGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCggggg                                                                                                                                           >  2:1116464/1‑150 (MQ=255)
            ccACACTGAATAATGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGc                                                                                                                                          >  2:1278493/1‑150 (MQ=255)
             cacaCTGAATAATGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCCGCTCCGCGCTGGGGCg                                                                                                                                         <  1:1730712/150‑1 (MQ=255)
               cacTGAATAATGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGCTGTGACTAAAACCAACATCAAAGATACGGTAATTAAAGACGGATGCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCCGCTCCGCGCTGGGGCGta                                                                                                                                       <  1:1054584/150‑2 (MQ=255)
               cacTGAATAATGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTg                                                                                                                                       >  2:1625556/1‑150 (MQ=255)
                acTGAATAATGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGa                                                                                                                                      >  2:744038/1‑150 (MQ=255)
                acTGAATAATGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGAACGATCTGAATAAAAACATCATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTATGAGTTACGCCCCAGCGCGGAGCGGGGGCGTGa                                                                                                                                      >  1:865524/1‑150 (MQ=255)
                 cTGAATAATGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGAt                                                                                                                                     >  2:961237/1‑150 (MQ=255)
                 cTGAATAATGGCCTGAAAGATGGCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGCGCGAGCTGTAAGCGTTACGCCCCCGCTCCGCGCTGGGGCGtaac                                                                                                                                     <  1:1637412/150‑4 (MQ=255)
                     ataatGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTtct                                                                                                                                 <  2:932610/150‑1 (MQ=255)
                      taatGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTtctc                                                                                                                                <  2:1383674/150‑1 (MQ=255)
                       aatGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCTATCTATGTGAATAAATACAACATCAAAGATACGGTTATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTtctct                                                                                                                               <  2:1911373/150‑1 (MQ=255)
                         tGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTcc                                                                                                                             <  1:1668728/150‑1 (MQ=255)
                          gggcTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACTGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCa                                                                                                                            <  2:1669606/147‑1 (MQ=255)
                          ggCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTCAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCa                                                                                                                            >  2:650161/1‑150 (MQ=255)
                          ggCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCCGCGCGGCGCGGGGGCGTGATTTGTCTCCa                                                                                                                            >  2:125768/1‑150 (MQ=255)
                          ggCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGACCTTTAAGCGTTACGCCCCAGCGCGCAGCGGGGGCGTGATTTCTCTCCa                                                                                                                            >  2:728935/1‑150 (MQ=255)
                           gCCTGAAAGATGTCCCCTCCCGCCTCCTGCCACCGATCGATGTGAATAAAAACAACATCAAAGCTACGGTAATTAAAGACGGATTCCCCAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCAt                                                                                                                           >  1:1100343/1‑150 (MQ=255)
                            ccTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATg                                                                                                                          >  2:65078/1‑150 (MQ=255)
                                   gATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGtg                                                                                                                   >  1:2011477/1‑150 (MQ=255)
                                     tGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGtgaa                                                                                                                 >  2:827544/1‑150 (MQ=255)
                                      gTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGtgaat                                                                                                                <  2:1746970/150‑1 (MQ=255)
                                       tCCCCTCCCGCCTCCTGACACCGCTCGATGTGCTTAAAAACATCATGAAAGATACGGTAATTAAATACGGATTCCACAATGATAGAGAGCGGTAAGCGGTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCTCGtgaatg                                                                                                               <  2:277898/150‑1 (MQ=255)
                                        ccccTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGtgaatga                                                                                                              >  2:178888/1‑150 (MQ=255)
                                        ccccTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGtgaatga                                                                                                              >  1:25614/1‑150 (MQ=255)
                                         cccTCCCGCCTCCTGCCACCGATCGATGTGAATAAATACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGtgaatgat                                                                                                             >  2:246643/1‑149 (MQ=255)
                                             cccGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTgg                                                                                                         <  2:1437608/150‑1 (MQ=255)
                                               cGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCt                                                                                                       >  2:807268/1‑150 (MQ=255)
                                                 cctcccGACACCGAGCGATTTGAAAAAAAACAACCTCAAAGATACGGTAATTAAAGACGGATTCCACAAAGCGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTa                                                                                                     <  2:910496/150‑1 (MQ=255)
                                                     ctGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTg                                                                                                 >  1:1712983/1‑150 (MQ=255)
                                                     ctGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTg                                                                                                 <  2:429515/150‑1 (MQ=255)
                                                     ctGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTg                                                                                                 >  1:1636614/1‑150 (MQ=255)
                                                     ctGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTg                                                                                                 <  1:1424356/150‑1 (MQ=255)
                                                     ctGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTg                                                                                                 <  1:671276/150‑1 (MQ=255)
                                                      tGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTgg                                                                                                >  2:250088/1‑150 (MQ=255)
                                                      tGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTgg                                                                                                >  1:1692424/1‑150 (MQ=255)
                                                         cacCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGt                                                                                             >  1:759194/1‑150 (MQ=255)
                                                            cgatcgatGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGt                                                                                          >  2:1259073/1‑150 (MQ=255)
                                                             gatcgatGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGtt                                                                                         <  1:1069879/150‑1 (MQ=255)
                                                             gatcgatGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGtt                                                                                         >  2:2018758/1‑150 (MQ=255)
                                                             gatcgatGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGtt                                                                                         >  1:166442/1‑150 (MQ=255)
                                                              atcgatGTGAATAAAAACAACATCAAAGATACTGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTa                                                                                        <  2:1681294/150‑1 (MQ=255)
                                                                 gatGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGc                                                                                     >  1:527872/1‑150 (MQ=255)
                                                                    gtgAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCtt                                                                                  >  1:1765820/1‑150 (MQ=255)
                                                                       aaTAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCCTGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTAta                                                                               >  2:302957/1‑149 (MQ=255)
                                                                       aaTAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATc                                                                               >  2:557737/1‑150 (MQ=255)
                                                                        aTAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCt                                                                              >  2:1967237/1‑150 (MQ=255)
                                                                         tAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGTCTTATCTa                                                                             <  1:1385988/150‑1 (MQ=255)
                                                                          aaaaaCAACATCAAAGCTACGGTAATTAAAGACGGCTTCCACAAAGAGAGCGAGCTGGAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTAc                                                                            <  1:1754013/150‑1 (MQ=255)
                                                                           aaaaCAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAATCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACt                                                                           >  2:307883/1‑150 (MQ=255)
                                                                           aaaaCAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACt                                                                           >  2:1564390/1‑150 (MQ=255)
                                                                           aaaaCAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACt                                                                           >  1:580371/1‑150 (MQ=255)
                                                                            aaaCAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACtt                                                                          <  1:780401/150‑1 (MQ=255)
                                                                             aacaacATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGCGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTg                                                                         <  2:1338831/150‑1 (MQ=255)
                                                                                aacaTCAAAGATACGGTAAATAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGaaa                                                                      <  2:1987375/150‑1 (MQ=255)
                                                                                 acaTCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAAt                                                                     <  2:1471304/150‑1 (MQ=255)
                                                                                   aTCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATga                                                                   >  1:782750/1‑150 (MQ=255)
                                                                                    tCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCAAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGAGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATgaa                                                                  >  1:1877093/1‑150 (MQ=255)
                                                                                     cAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATgaag                                                                 <  1:1283466/150‑1 (MQ=255)
                                                                                     cAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCATAGGTGGAGTCGTTATGCCTTATCTAATTGAAATgaag                                                                 >  1:1205165/1‑150 (MQ=255)
                                                                                      aaaGATACGGTAATTAAAGAGGGATTCCACAAAGAGAGCTAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATgaaga                                                                <  1:1982235/150‑1 (MQ=255)
                                                                                       aaGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATgaagaa                                                               <  1:1203979/150‑1 (MQ=255)
                                                                                         gATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCGACTTGAAATGAAGAACa                                                             >  1:742336/1‑150 (MQ=255)
                                                                                           tACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGGGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCGTTATCTACTTGAAATGAAGAACAtt                                                           >  1:1285504/1‑150 (MQ=255)
                                                                                             cGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTAc                                                         >  1:536616/1‑150 (MQ=255)
                                                                                             cGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTAc                                                         <  2:1660169/150‑1 (MQ=255)
                                                                                               gTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCa                                                       >  2:979921/1‑150 (MQ=255)
                                                                                                     aaaGACGGATTCCACAAAGATATCGAGCTGTAAGCGTTACGCCACAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCTTTATGCCTTATCTACTTGTAATGAATAACATTTCAAAAACCt                                                 >  2:1282743/1‑150 (MQ=255)
                                                                                                        gACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCg                                              >  1:498853/1‑150 (MQ=255)
                                                                                                            gATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAg                                          <  2:680829/150‑1 (MQ=255)
                                                                                                               tCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAgtgt                                       <  2:580371/150‑1 (MQ=255)
                                                                                                                  acaAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGaa                                    >  2:1188960/1‑150 (MQ=255)
                                                                                                                    aaaGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGGTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAgg                                  >  1:1575607/1‑150 (MQ=255)
                                                                                                                        tgagcgagcTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGAt                              <  1:822360/149‑1 (MQ=255)
                                                                                                                        agagCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGAt                              <  1:744038/150‑1 (MQ=255)
                                                                                                                           gcgagcTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCAATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGa                           >  2:1400031/1‑150 (MQ=255)
                                                                                                                            cgagcTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGACTGATTTGGCTTAGGTGGAGTCGGTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGAt                          >  1:170448/1‑150 (MQ=255)
                                                                                                                            cgagcTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTTTGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGAt                          >  1:1861041/1‑150 (MQ=255)
                                                                                                                              agcTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTTTGAAGGCGATTGATaa                        <  2:25614/150‑1 (MQ=255)
                                                                                                                                cTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTGTCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATAACg                      <  1:654059/150‑1 (MQ=255)
                                                                                                                                 tGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATAACgg                     >  1:244452/1‑149 (MQ=255)
                                                                                                                                  gTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGGCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATAACGTc                    <  2:1852486/150‑1 (MQ=255)
                                                                                                                                      gCGTTACGCCCCAGCGCTGAGAGGGGGGGTGATTTCTGTCCATGCCGCGTGTATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATAAAGTCTGCt                <  2:1257965/150‑1 (MQ=255)
                                                                                                                                       cGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGACCATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATAACGTCTGCTTg              >  1:916302/1‑150 (MQ=255)
                                                                                                                                        gTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATAACGTCTGCTTg              >  2:1376808/1‑150 (MQ=255)
                                                                                                                                        atcaCGCCCCCGCTCCGCGCTGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGGCTTATCTACTTGAAATGAAGAACCTTTCCAAACCCTTCGGCAGTGTGAAGGCGATTGATTACGTCTGCTTg              >  1:1786558/4‑150 (MQ=255)
                                                                                                                                        atcaCGCCCCCGCTCCGCGCTGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGCTAACGTCTGCTTg              >  2:398560/4‑150 (MQ=255)
                                                                                                                                        atcaCGCCCCCGCTCCGCGCTGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATAACGTCTGCTTg              >  2:1323497/4‑150 (MQ=255)
                                                                                                                                        atcaCGCCCCCGCTCCGCGCTGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATAACGTATGCTTg              >  1:1491281/4‑150 (MQ=255)
                                                                                                                                         ttACGCCCCAGCGCTGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAAGGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATAACGTCTGCTTGc             <  1:148070/150‑1 (MQ=255)
                                                                                                                                           aCGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAAACTTCGGCAGTGTGAAGGCGATTGATAACGTCGGCTTGCgg           >  2:530496/1‑150 (MQ=255)
                                                                                                                                             gCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAATTGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATAACGTCTGCTTGCGGtt         <  1:238187/150‑1 (MQ=255)
                                                                                                                                              ccccAGCGCGGGGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATAACGTCTGCTTGCGGTTg        <  2:1895194/150‑1 (MQ=255)
                                                                                                                                              ccccAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTTTGGCTTATCTACGTGAAATGAAGAACATTACCAAAACCTTGCGCAGTGTGAAGGCGAATGATAACCTCTGCTTGCGGttt        >  2:913075/1‑149 (MQ=255)
                                                                                                                                              ccccAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATAACGTCTGCTTGCGGTTg        >  1:389826/1‑150 (MQ=255)
                                                                                                                                                   gcgcGGAGCGGGGGCGTGATTCCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATAACGTCTGCTTGCGGTTGAATGc   >  1:1267956/1‑150 (MQ=255)
                                                                                                                                                    cgcgGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATAACGTCTGCTTGCGGTTGAATGCt  <  2:1487489/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CAAAAGCAGATACCACACTGAATAATGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATAACGTCTGCTTGCGGTTGAATGCT  >  NC_000913/3731990‑3732287

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