Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,014,2220CA27.4% 90.9 / 18.5 62Q332K (CAA→AAA) fliFflagellar basal‑body MS‑ring and collar protein
Reads supporting (aligned to +/- strand):  ref base C (40/5);  new base A (0/17);  total (40/22)
Fisher's exact test for biased strand distribution p-value = 3.56e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ATCCTGTCGCCTATTGTTGGTAACGGTAATATTCACGCCCAGGTTACGGCGCAGCTGGACTTCGCCAGTAAAGAACAAACGGAAGAACAGTATCGCCCTAACGGTGATGAATCTCATGCGGCGCTTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGCGGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCGATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCAACGTCTGTCAGTCGCGGTCGTGGTGAATTACAAAACCTTGCCAGATGGCAAACCGTTGCCTCTCAGCAACGAACAGATGAAGCAAATTGAAGATCTGACCCGCGAGGCGATGGG  >  NC_000913/2013973‑2014466
                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                    
aTCCTGTCGCCTATTGTTGGTAACGGTAATATTCACGCCCAGGTTACGGCGCAGCTGGACTTCGCCAGTAAAGAACAAACGGAAGAACAGTATCGCCCTAACGGTGATGAATCTCATGCGGCGCTTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGCGGCGTCCCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTACGGCAAATc                                                                                                                                                                                                                                                      >  1:338643/1‑250 (MQ=255)
   cTGTCGCCTATTGTTGGTAACGGTAATATTCACGCCCAGGTTACGGCGCAGCTGGACTTCGCCAGTAAAGAACAAACGGAAGAACAGTATCGCCCTAACGGTGATGAATCTCATGCGGCGCTTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGCGGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCACATCaaa                                                                                                                                                                                                                                                   >  2:128687/1‑250 (MQ=255)
      tCGCCTATTGTTGGTAACGGTAATATTCACGCCCAGGTTACGGCGCAGCCGGACTTCGCCAGTAAAGAACAAACGGAAGAACAGTATCGCCCTAACGGTGATGAATCTCATGCGGCGCTTCGTTCACGCCAGCTTAATGAGAGAGAGCAAAGCGGTTCCGGTTATCCGGGCGGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATa                                                                                                                                                                                                                                                >  1:579596/1‑250 (MQ=255)
         ccTATTGTTGGTAACGGTAATATTCACGCCCAGGTTACGGCGCAGCTGGACTTCGCCAGTAAAGAACAAACGGAAGAACAGTATCGCCCTAACGGTGATGAATCTCATGCGGCGCTTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGGTATCCGGGCGGCGTACCGGGGGCGTTGTCGACTCAACCGGCCCCTGCGAATAACGCGCCAATCCGCACGCCTCCGGCAAATCAAAATAACCg                                                                                                                                                                                                                                            >  2:227529/1‑251 (MQ=255)
                         gTAATATTCACGCCCAGGTTACGGCGCAGCTGGACTTCGCCAGTAAAGAACAAACGGAAGAACAGTATCGCCCTAACGGTGATGAATCTCATGCGGCGCTTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGCGGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGc                                                                                                                                                                                                                            >  1:258752/1‑251 (MQ=255)
                            atatTCACGCCCAGGTTACGGCGCAGCTGGACTTCGCCAGTAAAGAACAAACGGAAGAACAGTATCGCCCTAACGGGGATGAATCTCATGCGGCGCTTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGCGGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCCGCAGCAGACGAGcacc                                                                                                                                                                                                                         >  2:370996/1‑251 (MQ=255)
                                    cccaaGGTTCCGGCGCAGCTGAACTTCCCCAGAAAAAAACAAACGGAAGAACAGTATCGCCCTACGGGTGATGAATCTCAGGCGGCGTTTCGTTCACGCCAGTTAAATTGTGGCGGGCAAAGCGGTCCGGGTTACCCGTGCGGCTTCCCGGGGGCTTTGCCGAACAACCCGGCACCTGCGAAAAACTCCCCAATCAGCCCGCCCCCGGCAAACAAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCaa                                                                                                                                                                                                                 <  1:544443/247‑1 (MQ=255)
                                       cAGGTTACGGCGCAGTTGGACTTCGCCAGTAAAGAAAAAACGGAAGAACAGTATCGCCCAAACGGTGATGAATCTCATGCGGCGCTTGGTCCACGCCAGCTTAATGAGAGCGAGCAAAGGGGTTCGGGTTATCCGGGCGGCGTACCGGGGGCGTTGCAGAATCAACCGGCCCCTGCAAAAAACGCGCAAACCAGCACGCCTCCGGCAAATCAAAAAAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAg                                                                                                                                                                                                              <  1:466206/251‑1 (MQ=255)
                                        aGGTTACGGCGCAGCTGGACTTCGCCAGTAAAGAACAAACGGAAGAACAGTATCGCCCTACCGGTGATGAATCTCATGCGGCGCTTCGTTCACGCAAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGCGGCGTACCGGGGGCGTTGTCGAATCAACCGGCACATGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCCATAGt                                                                                                                                                                                                             >  1:30854/1‑251 (MQ=255)
                                        aGGTTACGGCGCAGCTGGACTTCGCCAGTAAAGAACAAACGGAAGAACAGTATCGCCCTAACGGTGATGAATCTCATGCGGCGCTTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGCGGCGTACCGGGGGCGTTGTGGAATCACCCGGCACCTGCGAGTAACGCGCCACTCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGCGCACCACCAGCAATAGt                                                                                                                                                                                                             >  2:61669/1‑251 (MQ=255)
                                        aGGTTACGGCGCAGCTGGACTTCGCCAGTAAAGAACAAACGGAAGAACAGTATCGCCCTAACGGTGATGAATCTCATGCGGCGCTTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGCGGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCCCCCGCAATAgg                                                                                                                                                                                                             >  2:93698/1‑250 (MQ=255)
                                        aGGTTACGGCGCAGCTGGACTTCGCCAGTAAAGAACAAACGGAAGAACAGTATCGCCCTAACGGTGATGAATCTCATGCGGCGCTTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGCGGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGt                                                                                                                                                                                                             >  1:123185/1‑251 (MQ=255)
                                             aCGGCGCAGCTGGACTTCGCCAGTAAAGAACAAACGGAAGAACAGTATCGCCCTAACGGTGATGAATCTCATGCGGCGCTTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGCGGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGcc                                                                                                                                                                                                        >  2:376021/1‑251 (MQ=255)
                                               ggCGCAGCTGGACTTCGCCAGTAAAGAACAAACGGAAGAACAGTATCGCCCTAACGGTGATGAATCTCATGCGGCGCTTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGCGGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCcgc                                                                                                                                                                                                      >  1:333574/1‑251 (MQ=255)
                                                  gcAGCTGGACTTCGCCAGTAAAGAACAAACGGAAGAACAGTATCGCCCTAACGGTGATGAATCTCATGCGGCGCTTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGCGGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGt                                                                                                                                                                                                    >  2:107026/1‑250 (MQ=255)
                                                               gCCAGTAAAGAACAAACGGAAGAACAGTATCGCCCTAACGGTGATGAATCTCATGCGGCGCTTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGCGGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCCCCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGa                                                                                                                                                                                       >  2:544430/1‑250 (MQ=255)
                                                                 cagaaaaGAACAAACGAAAGAACATTATCGCCCAAACGGTGTTAAATCTCATGCGGCGTTCCGTTCCCGCCAGCTTAATGAGAGGGGGCAAAGCGGTTCCGGTTACCGGGGCGGCGTCCGGGGGGCGTTGTCGAATCACCCGGCACTTGCGAATAACGCGCCAACCAGCACGCCTCCGGCAAATAAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATg                                                                                                                                                                                    <  1:533467/247‑1 (MQ=255)
                                                                 cAGTAAAGAACAAACGGAAGAACAGTATCGCCCTAACGGTGATGAATCTCATGCGGCGCTTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGCGGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATg                                                                                                                                                                                    >  1:476310/1‑251 (MQ=255)
                                                                           aaaaCGGAAGAACAGTACCCCCCTAACGGTGATAAATCTCGTGCGGCGCTTCGTTCCCCCCAGCTAAATGAGAGCGAGCAAAGCGGTTCCGGTTACCCGGGGGGCGTCCGGGGGGCGTTGTCGAACCACCCGGCCCCTGCGAATACCGCGCCAATCAGCACGCCTCCGGAAAACAAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTaa                                                                                                                                                                          <  1:503516/250‑1 (MQ=255)
                                                                                   agaaCAGTATCGCCCTAACGGTGATGAATCTCATGCGGCGCTTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGCGGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGa                                                                                                                                                                    >  1:221725/1‑249 (MQ=255)
                                                                                       cAGTATCGCCCTAACGGTGATGAATCTCATGCGGCGCTTCGTTCACGCCCGCTTAATGAGAGCGAGCCAAGCGGTTCCGGTTATCCGGGCGGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGACTGAAACCAGTAACTACGAAGTCgc                                                                                                                                                              >  1:504584/1‑250 (MQ=255)
                                                                                          tATCGCCCTAACGGTGATGAATCTCATGCGGCGCTTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGCGGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCGATCg                                                                                                                                                           >  2:250660/1‑251 (MQ=255)
                                                                                            tCGCCCTAACGGTGATGAATCTCATGCGGCGCTTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGCGGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCCGTAACTACGAAGGCGATCGcc                                                                                                                                                         >  2:462274/1‑250 (MQ=255)
                                                                                                  tAACGGTGATGAATCTCATGCGGCGCTTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGAGGTTCCGGTTATCCGGGCGGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCGATCGCACCATTc                                                                                                                                                   >  1:310067/1‑251 (MQ=255)
                                                                                                   accGGTGATGAATCCCATGCGGCGCTTCGTTCCCGCCAGCTTAATGAGAGCGAGCAAGGCGGTTCCGGTAACCGGGGCGGCGTACGGGGGGCGTTGTCGAATCAACCGGCACCTGCAAATACCGCGCCAATCAGCACGCCTCCGGCAAATCAAAAAAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCGATCGCACCAtt                                                                                                                                                    <  1:494992/247‑1 (MQ=255)
                                                                                                   aacgggtatgaATCTCAGGCGGCCCTTCTTTCCCGCCAGTTTAATGAGAGCGAGAAAAGCGGTTCCGTTACTCCGGGCGGCGACCCGGGGGGGTTGCCGAACCAACCGGCCCCTGCGAATACCGCCCAAATCAGCACGCCCCCGGAAAATAAAAATACCCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCGATCGCACCATTc                                                                                                                                                   <  2:393939/243‑1 (MQ=255)
                                                                                                   aaCGGTGATGAATCTCATGCGGCGCTTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTACCCGGGGGGCGTACCGGGGGCGTTGCCGAATCAACCGCCCCCTGCGAATAACGCGCAAATCAGCACGCCCCCGGCAAATAAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCGATCGCACCATTCg                                                                                                                                                  <  1:376021/251‑1 (MQ=255)
                                                                                                           tgaATCTCATGCGGCGCTTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGCGGCGTCCCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCACTCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCGATCGCACCATTCGTCATAcc                                                                                                                                           >  2:240320/1‑250 (MQ=255)
                                                                                                                   cggcggggCTCCGTTCACGCCACCTTAATGAGAGGGAGCAAAGCGTTTCCGGTTACCCGGGGGTCGTACCGGGGGCTTTGTCGAATCAACCGGCCCCTGCGAATAACGCGCAAATCAGCACGCCTCCGGCAAAAAAAAAAAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCGATCGCACCATTCGTCATACCAAAATGAAc                                                                                                                                  <  1:28512/249‑1 (MQ=255)
                                                                                                                                 tCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGCGGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCGATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGc                                                                                                                     >  1:117961/1‑250 (MQ=255)
                                                                                                                                  cACGCCAGATTAAGGAGGGCGAGCAAAGCGGTCCCGGTTACCCGGGGGGCTTCCCGGGGGCGTTGTCGAACCAACCGGCACCTGCAAATAACGCGCAAATCAGCACGCCTCCGGCAAATAAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCGATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCaa                                                                                                                   <  1:227529/251‑1 (MQ=255)
                                                                                                                                            tAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGCGGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCGATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCAACGTCTggc                                                                                                           >  1:340809/1‑247 (MQ=255)
                                                                                                                                                     gcgagcAAAGCGGTTCCGGTAACCCGGGCGGCGTACCGGGGGCGTTGTCGAATCACCCGGCACCTGCGAAAAACCCCCCAATCAGCACCCCCCCGGCAAACAAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCGATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCAACGTCTGTCAGTCGCGGTc                                                                                                 <  1:544430/251‑1 (MQ=255)
                                                                                                                                                     gcgagcAAAGCGGGCCCGGTTATCCGGGCGGCGCACCGGGGGCGTTGTCCAATCAACCGGCACCTGCGAATAAAGCGCCAATCAGCACGCCTCCGGCACCTCAAAATACCCGCCAGCAGCAGGCGAGCACCACCAGAAATCGTGGGCCGCGTAGCACACAGCGGAATGAAAACAGTAACTACGAAGTCGATCGCACCATTCGTCATACCAAAATGAACGTGGGCGAGGTTGACCGTCCGTCAGTCGCggcc                                                                                                 >  2:263864/1‑249 (MQ=255)
                                                                                                                                                      ccgagtaaaccgggcccggtaacccGGGCGGCGTCCGGGGGGCGTGGTCGAACCAGCCGGCCCCTGCAAATAACGCGCCAATCAGCCCGCCTCCGACAAACAAAAAGAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCGATCGCACCATTCGTCATACCAAAATTAACGTGGGCGATGTGCAACGTCTGTCAGTCGCGGTc                                                                                                 <  2:30854/228‑1 (MQ=255)
                                                                                                                                                         gcAAAGCGGTTCCGTTTACCCGGGCGGCGTACGGGGGGCGTTGTCGAACCAACCGGCCCCTGCGAATACCGCGCCAATCAGCACCCCCCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCGATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCAACGTCTGTCAGTCGCGGTCgtgg                                                                                             <  1:250660/251‑1 (MQ=255)
                                                                                                                                                               cGGTGCCGGTTATCCGGGCGGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGGGAATAACGCGACAAGCAGCACGCCTCCGGCACATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGGGGGCCGCGGAGCACACAGCGGAATGAGACCAGTAACTACGAAGTCGATCGCACCATTCGTCATACCAAAATGAACGAGGGCGATGTGCAACGTCTGTCAGTCGCGGGCgtgggacatt                                                                                       >  2:80139/1‑245 (MQ=255)
                                                                                                                                                                   tCCGGTTACCCGGGCGGCGTCCCGGGGGCGTTGTCGAACCACCCGGCCCCGGCGAAAAACGCCCCAATCACCACGCCCCCGGCAAATAAAAAAAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCGATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCAACGTCTGTCAGTCGCGGTCGTGGTGAATTACa                                                                                    <  1:93698/250‑1 (MQ=255)
                                                                                                                                                                      ggTTATCCGGGGGGCGTGCGGGGGGCGTTGTCAAATCAACCGGCACCTGCAAATAACCGGCAAACCACCCCGCCCCCGGCAAATAAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCGATCGCACCATTCGGCAGACCAAAATGAACGTGGGCGATGTGCAACGTCTGTCAGTCGCGGTCGTGGTGAATTACAAAAc                                                                                <  1:552577/251‑1 (MQ=255)
                                                                                                                                                                        ttATCCGGGCGGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCGATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCAACGTCTGTCAGTCGCGGTCGTGGTGACTTACAAAACCt                                                                              >  2:574830/1‑251 (MQ=255)
                                                                                                                                                                          aTCCGGGCGGGGTACGGGGGGCGTTGGCAAATCACCCGCCCCCGGCGATAACCGCCCCAATCAGCCCCCCTCCGGCAAAAAAAAATATCCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCGATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCAACGTCTGTCAGTCGCGGTCGTGGTGAATTACAAAACCTTg                                                                            <  1:462274/251‑1 (MQ=255)
                                                                                                                                                                          aTCCGGGCGGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAAAGCGCCAATCAGCACGCCTCCGGCCAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCGATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCAACGTCTGTCAGTCGCGGTCGTGGTGAATTACAAAACCTTg                                                                            >  1:62485/1‑251 (MQ=255)
                                                                                                                                                                             ccGGGGGGCGACCGGGGGGGGTGTACGAAACAACCGGCACCTGCAAAAAACCCCCCAACCAGCACGCCCCCGGAAAATAAAAAAAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCGATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCAACGTCTGTCAGTCGCGGTCGTGGTGAATTACAAAACCTTGcc                                                                          <  1:370996/251‑1 (MQ=255)
                                                                                                                                                                               gggCGGCGTCCCGGGGGCGTTGCCGAATCAACCGCCCCCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCGATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCAACGTCTGTCAGTCGCGGTCGTGGTGAATTACAAAACCTTGCCAg                                                                        <  1:128687/251‑1 (MQ=255)
                                                                                                                                                                                ggcggcGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAACTCAAAATACCCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCGATCGCACCAt                                                                                                                                                     >  1:235539/1‑173 (MQ=255)
                                                                                                                                                                                ggcggcGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCGATCGCACCATTCGTCATACCAAAATGAACGTGGGc                                                                                                                            >  1:377949/1‑198 (MQ=255)
                                                                                                                                                                                 gcggcgTACGGGGGGCGTTGTCGAATCAACCGGCCCCTGCGAAAAACGCGCCAATCAGCGCGCCTCCGCAAAATAAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCGATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCAACGTCTGTCAGTCGCGGTCGTGGTGAATTACAAAACCTTGCCAGAt                                                                      <  1:574830/251‑1 (MQ=255)
                                                                                                                                                                                            ggggCGTTGTCGAATCACCCGGCCCCTGCGAATACCGCGCCAACCAGCACGCCTCCGGCAAATAAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCGATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCAACGTCTGTCAGTCGCGGTCGTGGTGAATTACAAAACCTTGCCAGATGGCAAACCGtt                                                           <  1:345480/251‑1 (MQ=255)
                                                                                                                                                                                                 gTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGAACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGAGTAACACACAGCGGAATGAACCCAGTAACTACGAAGTCGATCGAACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCAACGTCTGTCAGTCGCGGTCGAGGTGAATTACAAAACCTTGCCAGATGGCAAACCGTTGCctc                                                      >  1:535331/1‑251 (MQ=255)
                                                                                                                                                                                                  ttGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCCCACAGCGGAATGAAACCAGTAACTACGAAGTCGATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCAACGTCTGTCAGTCGCGGTCGTGGTGAATTACAAAACCTTGCCAGATGGCAAACCGTTGCctct                                                     >  2:85537/1‑251 (MQ=255)
                                                                                                                                                                                                        ggaaccaacggcccccTGCGAATACCCCGCCAAGCAGCTCGCCTCGGGCAAGTCAAAATAGCCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCGATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCAACGTCTGTCAGTCGCGGTCGGGGTGAATTACAAAACCTTGCCAGATGGCAAACCGTTGCCTCTCAg                                                  <  1:251090/236‑1 (MQ=255)
                                                                                                                                                                                                             aaCCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCGATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCAACGTCTGTCAGTCGCGGTCGTGGTGAATTACAAAACCTTGCCAGATGGCAAACCGTTGCCTCTCAGCAACGa                                            >  1:278947/1‑250 (MQ=255)
                                                                                                                                                                                                                cGGCACCTGCGAATAACGCGCCAATCAGCAAGCCTCCGGCACATCAAACTAACCGCCAGCGGCACGCGATCACCACCAGGGATAGTGGGCCCCGGAACACACACCGGCATGAAACCAGTAACTACGAAAACGATCGCACCATTCGCCATACCAGAATGAACGTGTGCGATGTGCCACGTCTGTCAGTCGCGGACGTGGGCGATTACAAAACCTTGCCAGATGGCACACCGTGGCCTCTCAGCAACGGACac                                        >  1:339659/1‑250 (MQ=255)
                                                                                                                                                                                                                  ggCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCGATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCAACGTCTGTCAGTCGCGGTCGTGGTGAATTACAAAACCTTGCCAGATGGCAAACCGTTGCCTCTCAGCAACGAAc                                          >  1:251772/1‑248 (MQ=255)
                                                                                                                                                                                                                    cACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAAACGCAAGCAGCAGGCGCGCCCCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCGATCGCACCATTCGTCATCCCAAAATGAACGTGGGCGATGTGCAACGTCTGTCAGTCGCGGTCGTGGTGAATTACAAAACCTTGCCAGATGGCAAACCGTTGACTCTCAGCAACGAACAGATg                                     >  1:188700/1‑251 (MQ=255)
                                                                                                                                                                                                                                   cccaaccccagcccccccccGGCAAATAAAAAAAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCGATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCAACGTCTGTCAGTCGCGGTCGTGGTGAATTACAAAACCTTGCCAGATGGCAAACCGTTGCCTCTCAGCAACGAACAGATGAAGCAAATTGaa                         <  2:449966/231‑1 (MQ=255)
                                                                                                                                                                                                                                         aTCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCGATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCAACGTCTGTCAGTCGCGGTCGTGGTGAATTACAAAACCTTGCCAGATGGCAAACCGTTGCCTCTCAGCAACGAACAGATGAAGCAAATTGAAGATCTGAc                 >  1:217312/1‑251 (MQ=255)
                                                                                                                                                                                                                                           cAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCGATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCAACGTCTGTCAGTCGCGGTCGTGGTGAATTACAAAACCTTGCCAGATGGCAAACCGTTGCCTCTCAGCAACGAACAGATGAAGCAAATTGaa                         >  1:449966/1‑241 (MQ=255)
                                                                                                                                                                                                                                                  ccTCCGACAAAAAAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCGATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCAACGTCTGTCAGTCGCGGTCGTGGTGAATTACAAAACCTTGCCAGATGGCAAACCGTTGCCTCTCAGCAACGAACAGATGAAGCAAATTGAAGATCTGACCCGCGAGGc        <  1:240320/251‑1 (MQ=255)
                                                                                                                                                                                                                                                   cTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCCCCCCCAGCAATAGTGGGCCGCGTAGAACACAGCGGAATGAAACCAGTAACTACGAAGTCGATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCAACGTCTGTCAGTCGCGGTCGTGGTGAATTACAAAACCTTGCCAGATGGCAAACCGTTGCCTCTCAGCAACGAACAGATGAAGCAAATTGAAGATCTGACCCGCGAGGCg       >  2:280923/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   cTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCGATCGCACCATTCGTCATACCAAAATGCACGTGGGCGATGTGCACCGTCTGTCAGTCGCGGTCGTGGTGAATTACAAAACCTTGCCAGATGGCAAACCGTTGCCTCTCAGCAACGAACAGATGAAGCAAATTGAAGATCTGACCCGCGAGGCg       >  1:571201/1‑251 (MQ=255)
                                                                                                                                                                                                                                                    tCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCGATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCAACGTCTGTCAGTCGCGGTCGTGGTGAATTACAAAACCTTGCCAGATGGCAAACCGTTGCCTCTCAGCAACGAACAGATGAAGCAAATTGAAGATCTGACCCGCGAGGCg       >  1:250087/1‑250 (MQ=255)
                                                                                                                                                                                                                                                        gCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCGATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCAACGTCTGTCAGTCGCGGTCGTGGTGAATTACAAAACCTTGCCAGATGGCAAACCGTTGCCTCTCAGCAACGAACAGATGAAGCAAATTGAAGATCTGACCCGCGAGGCGATggg  >  2:46801/1‑251 (MQ=255)
                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                    
ATCCTGTCGCCTATTGTTGGTAACGGTAATATTCACGCCCAGGTTACGGCGCAGCTGGACTTCGCCAGTAAAGAACAAACGGAAGAACAGTATCGCCCTAACGGTGATGAATCTCATGCGGCGCTTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGCGGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCGATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCAACGTCTGTCAGTCGCGGTCGTGGTGAATTACAAAACCTTGCCAGATGGCAAACCGTTGCCTCTCAGCAACGAACAGATGAAGCAAATTGAAGATCTGACCCGCGAGGCGATGGG  >  NC_000913/2013973‑2014466

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