Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,208,1320GA68.2% 42.3 / 58.5 63intergenic (‑600/‑15)purH/rrsEbifunctional AICAR transformylase/IMP cyclohydrolase/16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base G (10/10);  new base A (21/22);  total (31/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CACACAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTT  >  NC_000913/4207889‑4208365
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                         
cacacaAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCa                                                                                                                                                                                                                                    <  2:150345/251‑1 (MQ=38)
       aCAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACtttt                                                                                                                                                                                                                             <  1:335210/251‑1 (MQ=38)
       aCAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACtttt                                                                                                                                                                                                                             <  1:218994/251‑1 (MQ=38)
           ccGGTTCGGTTGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAAt                                                                                                                                                                                                                         >  2:344700/1‑251 (MQ=38)
              gttcggttGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTg                                                                                                                                                                                                                       >  1:371542/1‑250 (MQ=38)
                              aaaTCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGc                                                                                                                                                                                                      >  2:25922/1‑251 (MQ=25)
                                     gAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGAt                                                                                                                                                                                                <  1:264410/250‑1 (MQ=25)
                                        aTTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGaa                                                                                                                                                                                            <  1:47750/251‑1 (MQ=25)
                                                  gACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTggcggc                                                                                                                                                                                  >  1:446974/1‑251 (MQ=25)
                                                           agagGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAac                                                                                                                                                                         <  2:240783/251‑1 (MQ=25)
                                                                        atatACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTggc                                                                                                                                                                                     >  1:357496/1‑226 (MQ=37)
                                                                        atatACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTggc                                                                                                                                                                                     <  2:357496/226‑1 (MQ=37)
                                                                             cGCCACCTCGCGACCGTTCGCTAAAGCGCGTCGCTATTGCTCTTTAACAATTTCTCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGCCGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGGTTAATTCTTTAAACATCAATCTTTTAAATAGAAGAGGTTGAGCATGACTCAGATTTAACGCTGAAGGCAGGCCTAAAACAGGCAAGACGAACggg                                                                                                                                                       >  2:350635/1‑250 (MQ=11)
                                                                               ccACCTCGTGACAGTGCGCTAAAGCGCGTCGCAACTGCTCGTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAATACGTAATTCATTACGAAGTTTAATTCTTTAAGGATCAAAGGTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGGTGGCGGCAGTCCTTACTCAGGCAAGTCGAACGGTaa                                                                                                                                                     <  1:173055/251‑1 (MQ=11)
                                                                                cACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAAc                                                                                                                                                    <  1:355454/251‑1 (MQ=25)
                                                                                cACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTGTTAAATTGAAGAGTTTGAACATGGCTCAGATTGAACGCTGGCGGCAGGCCTGACGCATTGAAGTCGAACGGTAAc                                                                                                                                                    <  2:464202/251‑1 (MQ=11)
                                                                                           aGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAgct                                                                                                                                         <  1:212330/251‑1 (MQ=37)
                                                                                            gTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAgctt                                                                                                                                        >  2:466040/1‑251 (MQ=25)
                                                                                                gcTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACGTTTAAATTGAATAGTGTGATCATTGCTCAGAGTGAACGCTGGCGGCAGGCCTAGCACATGCAAGTCGAACGTTAACATGTGGAAgcttgctt                                                                                                                                    <  2:53580/251‑1 (MQ=11)
                                                                                                             gCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACg                                                                                                                         >  2:4496/1‑249 (MQ=14)
                                                                                                             gCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACg                                                                                                                         <  1:4496/249‑1 (MQ=14)
                                                                                                               caCTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTg                                                                                                                     <  2:309531/250‑1 (MQ=14)
                                                                                                                   gCTCTTTAACAATTTATCAGACAATCTGTGTGGGCCCTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCgg                                                                                                                 >  1:259101/1‑251 (MQ=14)
                                                                                                                      ctTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACg                                                                                                              >  1:417904/1‑251 (MQ=31)
                                                                                                                          aaCAATTTATCAGACAATCTGTGTGGGCACTCGAAGATTCGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGt                                                                                                           <  2:371542/250‑1 (MQ=14)
                                                                                                                          aaCAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTg                                                                                                          >  2:465760/1‑251 (MQ=31)
                                                                                                                                             tgtggggtcACTCGAAGATACGGATTCTTACCGTCGCAAGGCGAAAAATGACTACCAAGTCTCAAGAGTTAACACGTAATTGATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTCAATTGAAGAGTTTGATCATGACTCAGATTCAACGCTGGCGGCAGGCCTATCACATGCAACTCAAACGGTTACAGGCAGAAGCTTGCTTCTTTGCTGACGAGTGACGGAAGGGTGTGTaagtgccgggaaccgt                                                                                       >  2:94939/9‑237 (MQ=11)
                                                                                                                                              gtgtgGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGaaa                                                                                          <  1:436958/247‑1 (MQ=14)
                                                                                                                                                gtgGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCt                                                                                    >  1:297356/1‑251 (MQ=14)
                                                                                                                                                    gCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAACTGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGCGAGTAATGTCTGGGAAACTGCCTGATg                                                                                >  1:30702/1‑251 (MQ=14)
                                                                                                                                                          gAAGATACGGATTCTTAACGTCGCAAGACGAAAAAAGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAATTTTAATTCTTTAAGCATCAAACTTGTAAAGTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCGTAAGACATGCAAGTCGAGCGGTAACAGGAAGAAGCTTGTTTGTTTGCTGATGAGTGGCGGAAGGGTGAGTAATGTCTGTGaaagtgtgtgaagcagggg                                                                          <  2:261730/251‑17 (MQ=11)
                                                                                                                                                              aTACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATa                                                                      <  1:343508/251‑1 (MQ=14)
                                                                                                                                                               tACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGGCTGGGAAACTGGCTGATGGAGGGGGATa                                                                      >  2:32034/1‑250 (MQ=14)
                                                                                                                                                                   gATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAactac                                                                 <  1:351799/251‑1 (MQ=31)
                                                                                                                                                                         tAACGTCGCAAGACGAAAAATGAATACCAAGTCGCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGTGATAACTACTGGaa                                                            <  2:220391/250‑1 (MQ=14)
                                                                                                                                                                           aCGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACg                                                         >  1:168963/1‑251 (MQ=14)
                                                                                                                                                                               cGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAg                                                     <  1:450635/251‑1 (MQ=31)
                                                                                                                                                                               cGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAg                                                     <  1:338278/251‑1 (MQ=31)
                                                                                                                                                                               cGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAg                                                     <  2:65020/251‑1 (MQ=14)
                                                                                                                                                                                 cAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCt                                                   <  1:259012/251‑1 (MQ=14)
                                                                                                                                                                                  aaGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCt                                                   >  1:287099/1‑250 (MQ=14)
                                                                                                                                                                                    gACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAAt                                                <  1:344700/251‑1 (MQ=14)
                                                                                                                                                                                            aTGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGAAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGGGTAATGTTTGGGAAACTGCCTGCTGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCAt                                         >  1:313642/1‑250 (MQ=9)
                                                                                                                                                                                              gAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAAc                                      <  1:85929/251‑1 (MQ=31)
                                                                                                                                                                                              gAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAAc                                      >  2:341224/1‑251 (MQ=14)
                                                                                                                                                                                                   ccAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGc                                 <  1:25922/251‑1 (MQ=14)
                                                                                                                                                                                                    cAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCa                                <  2:302942/251‑1 (MQ=31)
                                                                                                                                                                                                       gTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGa                             <  2:290595/251‑1 (MQ=31)
                                                                                                                                                                                                           cAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGCTGGAGGGGGGTAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGGCAGACCa                          >  2:334046/1‑250 (MQ=9)
                                                                                                                                                                                                            aaGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCaa                         >  1:412426/1‑250 (MQ=31)
                                                                                                                                                                                                             agagTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCaa                         >  1:237756/1‑249 (MQ=31)
                                                                                                                                                                                                                      acGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGAcc              >  2:137420/1‑251 (MQ=31)
                                                                                                                                                                                                                        gTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTTGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCtt            <  1:57593/251‑1 (MQ=14)
                                                                                                                                                                                                                           aTTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGTCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAAAGGAAGCTAATACCGCATAACGTCGCCAGACCAAAGCGGGGGACCCTCgg         >  1:104058/1‑251 (MQ=11)
                                                                                                                                                                                                                           aTTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCgg         >  2:161889/1‑251 (MQ=14)
                                                                                                                                                                                                                           aTTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCgg         <  1:398226/251‑1 (MQ=14)
                                                                                                                                                                                                                            ttCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGGGAGTAATGTCTGGGAAACTGACTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGACGCAAGACCAAAGAGGGGGACCTTCggg        >  2:415239/1‑251 (MQ=9)
                                                                                                                                                                                                                              cATTACGAAGTTTAATTCTTTAAGCATCCAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGcc      >  2:463683/1‑251 (MQ=14)
                                                                                                                                                                                                                               aTTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGGAAAGTCTGGGAAACTGCCTGATGGAGGGGGATAAATACTGGAAAAGGTAGCTAATAACGCATAACGTCGCAAGACCAAAGAGGGGGACCCTCgggact     >  2:347561/1‑248 (MQ=11)
                                                                                                                                                                                                                               aTTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGGGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCct     >  1:251464/1‑251 (MQ=14)
                                                                                                                                                                                                                               aTTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCTGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGCAATGTCTGGGAAACTGCCTGATGGAGGGGGCTAACTACTGGAAACGGTAGCTAATACCGCATTACGTCGCAAGACCCACGAGGGGGACCTTCgcgacc     >  2:9395/1‑246 (MQ=11)
                                                                                                                                                                                                                               aTTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCct     >  1:86850/1‑251 (MQ=14)
                                                                                                                                                                                                                                  aCGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCtt  <  1:197834/251‑1 (MQ=31)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                         
CACACAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTT  >  NC_000913/4207889‑4208365

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