Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,732,139 G→T 6.6% intergenic (+16/‑62) xylF → / → xylG xylose ABC transporter periplasmic binding protein/xylose ABC transporter ATP binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,732,1390GT 6.6% 807.8 / 27.1 241intergenic (+16/‑62)xylF/xylGxylose ABC transporter periplasmic binding protein/xylose ABC transporter ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base G (122/102);  new base T (6/11);  total (128/113)
Fisher's exact test for biased strand distribution p-value = 1.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.10e-01

Only 100 of 249 total aligned reads displayed.

CAAAAGCAGATACCACACTGAATAATGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATAACGTCTGCTTGCGGTTGAATGCTG  >  NC_000913/3731990‑3732288
                                                                                                                                                     |                                                                                                                                                     
           aCCACACTGAATAATGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCCGCTCCGCGCtgggg                                                                                                                                            <  1:648543/150‑1 (MQ=255)
                acTGAATAATGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCCGCTCCGCGCTGGGGCGtaa                                                                                                                                       <  1:1541344/150‑3 (MQ=255)
                 cTGAATAATGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGAt                                                                                                                                      >  2:237301/1‑150 (MQ=255)
                  tGAATAATGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCCGCTCCGCGCTGGGGCGtaacg                                                                                                                                     <  2:1980596/150‑5 (MQ=255)
                  tGAATAATGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGAtt                                                                                                                                     <  2:1320955/150‑1 (MQ=255)
                     ataatGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCCGCTCCGCGCTGGGGCGtaacgctt                                                                                                                                  <  2:1867850/150‑8 (MQ=255)
                             cTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGc                                                                                                                          <  1:1490212/150‑1 (MQ=255)
                              tGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGAAACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGcc                                                                                                                         >  1:1430174/1‑150 (MQ=255)
                                 aaGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCcgcg                                                                                                                      >  2:1649681/1‑150 (MQ=255)
                                   gATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGtg                                                                                                                    <  1:16137/150‑1 (MQ=255)
                                       tCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGtgaatg                                                                                                                <  1:351140/150‑1 (MQ=255)
                                       tCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGtgaatg                                                                                                                >  2:1737333/1‑150 (MQ=255)
                                       tCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGtgaatg                                                                                                                >  2:357734/1‑150 (MQ=255)
                                        ccccTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGtgaatga                                                                                                               >  1:456704/1‑150 (MQ=255)
                                        ccccTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTATCCATGCCGCGtgaatga                                                                                                               >  2:846499/1‑150 (MQ=255)
                                           cTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAAtt                                                                                                            <  1:27144/150‑1 (MQ=255)
                                           cTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAAtt                                                                                                            <  1:419327/150‑1 (MQ=255)
                                             cccGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTgg                                                                                                          >  2:313112/1‑150 (MQ=255)
                                             cccGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTgg                                                                                                          >  2:1669996/1‑150 (MQ=255)
                                             cccGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTgg                                                                                                          >  2:1698451/1‑150 (MQ=255)
                                              ccGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGTCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGc                                                                                                         <  2:505292/150‑1 (MQ=255)
                                              ccGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATAAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGc                                                                                                         >  2:1852866/1‑150 (MQ=255)
                                                gCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCtt                                                                                                       >  1:1564161/1‑150 (MQ=255)
                                                         cacCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGt                                                                                              >  1:1739332/1‑150 (MQ=255)
                                                          acCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTc                                                                                             >  1:467130/1‑150 (MQ=255)
                                                           ccGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCg                                                                                            >  1:101339/1‑150 (MQ=255)
                                                           cagatcgatGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCg                                                                                            >  2:1496463/3‑150 (MQ=255)
                                                            cgatcgatGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGATACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGt                                                                                           >  1:1393096/1‑150 (MQ=255)
                                                              atcgatGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTa                                                                                         <  1:316436/150‑1 (MQ=255)
                                                                cgatGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATg                                                                                       <  1:480058/150‑1 (MQ=255)
                                                                 gatGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGc                                                                                      <  2:1364683/150‑1 (MQ=255)
                                                                  atGTCAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGACCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATCCCGCGTGAATGAATTCGCTTAGGTGGAGTCGTTATGca                                                                                     >  2:612951/1‑149 (MQ=255)
                                                                    gtgAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCtt                                                                                   <  2:101339/150‑1 (MQ=255)
                                                                    ctgAATAAAAACAAAATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCtt                                                                                   >  2:1394992/2‑150 (MQ=255)
                                                                     tgAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTa                                                                                  >  2:1787376/1‑150 (MQ=255)
                                                                      gAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGcgcg                                                                                                                                                     >  1:1585309/1‑82 (MQ=255)
                                                                      gAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGcgcg                                                                                                                                                     <  2:1585309/82‑1 (MQ=255)
                                                                      gAATAAAAACAACATCAAAGATACGATAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTAt                                                                                 >  2:1279407/1‑150 (MQ=255)
                                                                       aaTAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATc                                                                                >  2:647008/1‑150 (MQ=255)
                                                                       aaTAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATc                                                                                >  1:192060/1‑150 (MQ=255)
                                                                       aaTAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATc                                                                                >  1:80299/1‑150 (MQ=255)
                                                                        aTAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCt                                                                               >  2:1410880/1‑150 (MQ=255)
                                                                         tAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTa                                                                              >  2:883440/1‑150 (MQ=255)
                                                                         tAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTa                                                                              <  2:920094/150‑1 (MQ=255)
                                                                          aaaaaCAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTag                                                                             >  1:564147/1‑149 (MQ=255)
                                                                           aaaaCAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTACATGACGCGTGAATGAATTGGATTAGGTGGAGTCGTTATGCCTTATCTACt                                                                            >  2:385991/1‑150 (MQ=255)
                                                                             aacaacATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTg                                                                          >  2:1808382/1‑150 (MQ=255)
                                                                             aacaacATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTg                                                                          <  1:1121075/150‑1 (MQ=255)
                                                                               caacaTCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTg                                                                          <  2:228662/148‑1 (MQ=255)
                                                                               caacaTCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTg                                                                          >  1:228662/1‑148 (MQ=255)
                                                                                 acaTCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAAt                                                                      <  1:785799/150‑1 (MQ=255)
                                                                                 acaTCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAAt                                                                      <  1:700639/150‑1 (MQ=255)
                                                                                  caTCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATg                                                                     >  2:708526/1‑150 (MQ=255)
                                                                                    tCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCggggg                                                                                                                                            <  1:883349/77‑1 (MQ=255)
                                                                                     cAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGtgaatgaa                                                                                                              >  2:1338715/1‑106 (MQ=255)
                                                                                      aaaGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATgaaga                                                                 <  2:760193/150‑1 (MQ=255)
                                                                                      aaaGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATgaaga                                                                 >  2:1032518/1‑150 (MQ=255)
                                                                                      aaaGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATgaaga                                                                 >  1:1288598/1‑150 (MQ=255)
                                                                                      aaaGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATgaaga                                                                 >  1:1718382/1‑150 (MQ=255)
                                                                                      aaaGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATgaaga                                                                 >  1:1175436/1‑150 (MQ=255)
                                                                                       aaGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATgaagaa                                                                >  2:1604449/1‑150 (MQ=255)
                                                                                        aGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAAc                                                               <  1:1639888/150‑1 (MQ=255)
                                                                                         gATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACa                                                              <  2:104673/150‑1 (MQ=255)
                                                                                           tACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACAtt                                                            >  1:1804187/1‑150 (MQ=255)
                                                                                              ggTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTa                                                                                         >  1:581481/1‑118 (MQ=255)
                                                                                              ggTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTAcc                                                         >  2:590639/1‑150 (MQ=255)
                                                                                               gTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAGGAACATTACCa                                                        >  1:1995087/1‑150 (MQ=255)
                                                                                               gTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCa                                                        <  1:313112/150‑1 (MQ=255)
                                                                                                    tAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCTTGATTTCTCTCCATGCCGCGTGAATTAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAAcc                                                   <  1:385991/150‑1 (MQ=255)
                                                                                                    tAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCAATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAAcc                                                   <  2:974520/150‑1 (MQ=255)
                                                                                                     aaaGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCt                                                  >  2:152958/1‑150 (MQ=255)
                                                                                                      aaGACGGATTCCACAAAGAGAGTGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCtt                                                 >  2:1085938/1‑150 (MQ=255)
                                                                                                      aaGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTAAAATGAAGAACATTACCAAAACCtt                                                 >  1:738097/1‑150 (MQ=255)
                                                                                                          cGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGc                                             <  2:1635621/150‑1 (MQ=255)
                                                                                                          cGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGc                                             <  1:1649681/150‑1 (MQ=255)
                                                                                                          cGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGc                                             <  1:835105/150‑1 (MQ=255)
                                                                                                              ttCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAgtg                                         <  1:1032940/150‑1 (MQ=255)
                                                                                                                 cacaAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGa                                      >  2:454832/1‑150 (MQ=255)
                                                                                                                   caAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAg                                    >  2:1938791/1‑150 (MQ=255)
                                                                                                                         gagcgagcTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGAtt                              <  1:1308332/150‑1 (MQ=255)
                                                                                                                         gagcgagcTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGAtt                              <  1:1353251/150‑1 (MQ=255)
                                                                                                                          agcgagcTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTg                             <  2:80299/150‑1 (MQ=255)
                                                                                                                             gagcTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATa                          >  2:459119/1‑150 (MQ=255)
                                                                                                                                  gTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATAACGTc                     >  2:1440760/1‑150 (MQ=255)
                                                                                                                                  gTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATAACGTc                     >  1:288601/1‑150 (MQ=255)
                                                                                                                                  gTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATAACGTc                     >  2:333568/1‑150 (MQ=255)
                                                                                                                                   tAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATAACGTCt                    >  2:748927/1‑150 (MQ=255)
                                                                                                                                    aaGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATAACGTCTg                   >  2:492848/1‑150 (MQ=255)
                                                                                                                                     gaaatcaCGCCCCCGCTCCGCGCTGGGGCGTAATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATAACGTCTGc                  <  1:1937627/144‑1 (MQ=255)
                                                                                                                                        gTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATAACGTCTGCTTg               >  2:109952/1‑150 (MQ=255)
                                                                                                                                        atcaCGCCCCCGCTCCGCGCTGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACAATACCAAAACCTTCGGCAGTGTGAAGGCGATTGATAACGTCTGCTTg               >  2:1039674/4‑150 (MQ=255)
                                                                                                                                         ttACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATAACGTCTGCTTGc              <  1:1496762/150‑1 (MQ=255)
                                                                                                                                         tcaCGCCCCCGCTCCGCGCTGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATAACGTCTGCTTGc              >  1:1850364/3‑150 (MQ=255)
                                                                                                                                          tACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATAACGTCTGCTTGCg             <  1:1440760/150‑1 (MQ=255)
                                                                                                                                            cGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATg                                                                                       >  1:1218863/1‑74 (MQ=255)
                                                                                                                                            cGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATg                                                                                       <  2:1218863/74‑1 (MQ=255)
                                                                                                                                            cGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTg                             >  1:964463/1‑132 (MQ=255)
                                                                                                                                            cGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTg                             <  2:964463/132‑1 (MQ=255)
                                                                                                                                            cGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATAACGTCTGCTTGCGGt           >  2:112701/1‑150 (MQ=255)
                                                                                                                                               cccAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATAACGTCTGCTTGCGGTTGa        <  1:628314/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CAAAAGCAGATACCACACTGAATAATGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATAACGTCTGCTTGCGGTTGAATGCTG  >  NC_000913/3731990‑3732288

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