Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,126,2100TG35.3% 100.6 / 12.9 68T201P (ACC→CCC) priAprimosomal protein N'
Reads supporting (aligned to +/- strand):  ref base T (19/25);  new base G (24/0);  total (43/25)
Fisher's exact test for biased strand distribution p-value = 6.50e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GTCGTTCAGGCCGGAATGCAGAACTTCCACGGGGGCATTAAAACGTTCACGAAAACGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAGTTCGTTCGCCAGTCGCTAAACTCTGGTGTTTCACTTGCTAAATCACACAGACCTTTTTTGCGTAGCGCCTGCAACGCGGCATCATTAAATTCGAGCGTGGCGACCTGGTCGCGCCAGATTTTGCCTTGCCGTAACGCCGCCAGCGCCTGTTGTTGCTTGGGGGAGCGTTTCAGGCTGTTCAGATCCACCGCCTGGCCTTGTTCAGTGGCAAACC  >  NC_000913/4125965‑4126458
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                              
gTCGTTCAGGCCGGAATGCAGAACTTCCACGGGGGCATTAAAACGTTCACGAAAACGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGGCTGTTTGGTATTcc                                                                                                                                                                                                                                                           >  1:218485/1‑250 (MQ=255)
                  aagaccTTCCACGGGGGCATTAAAACGTTCCGGAAAACGGGGGATAGTTTGCGGTGCCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGGTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAg                                                                                                                                                                                                                                        <  2:373016/246‑1 (MQ=255)
                        ttCCACGGGGGCATTAAAACGTTCACGAAAACGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGACTTGCGCCAACGGCGGCGGACTGTTCGGTATTTCATTGCAAACGCTCACCAGCAACgg                                                                                                                                                                                                                                 >  1:178034/1‑251 (MQ=255)
                         tCCACGGGGGCATTAAAACGTCCACGAAAACGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGGGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACg                                                                                                                                                                                                                                  <  1:439/249‑1 (MQ=255)
                            aCGGGGGCATTAACACGTTCACGAAAACGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCCCACCCGAAACGGCAt                                                                                                                                                                                                                             >  2:336050/1‑250 (MQ=255)
                              gggggCATTAAAACGTTCACGAAAACGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAgt                                                                                                                                                                                                                          <  2:150362/251‑1 (MQ=255)
                                  gCATTAAAACGTTCACGAAAACGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAgttcgt                                                                                                                                                                                                                      <  1:319522/251‑1 (MQ=255)
                                  gCATTAAAACGTTCACGAAAACGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGGATTCCACCGCAACCGCTCCCCCGAAACAGCAttggtgcg                                                                                                                                                                                                                      >  2:17063/1‑244 (MQ=255)
                                              tCACGAAAACGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCCGGTCGGGATTCAATCGCAACCGCTCACCAGAAACGGCATAggtcgtatgcccgtcgcc                                                                                                                                                                                                          >  1:231154/1‑239 (MQ=255)
                                                aCGAAAACGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGGATTCAATCGCAACCCCCCCCCAGAAAACGCATAGTTCGTTCCCCAGTCGcga                                                                                                                                                                                                        >  1:487950/1‑248 (MQ=255)
                                                    aaaCGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGACAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGGCTGTTTGGGATTCAATATCAACCGCCCCCCACAAACGGGCAAATTCCTTCCCCAGTCCCTAAActc                                                                                                                                                                                                   >  2:103451/1‑251 (MQ=255)
                                                       cGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGGATTCAATCTCAACCGCTCACCAGAAACGGCATAGTTCGTTCCCCCATCGCTAAACTCTgg                                                                                                                                                                                                >  1:317920/1‑251 (MQ=255)
                                                           cGATAGTTGGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAGTTCGTTCGCCAGTCGCTAAACTCTGGTGtt                                                                                                                                                                                            <  1:112300/251‑1 (MQ=255)
                                                            gATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAGTTCGTTCGCCAGTCGCTAAACTCTGGTGttt                                                                                                                                                                                           <  1:458910/251‑1 (MQ=255)
                                                                  ttGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAACCCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAGTTCGTTCGCCAGTCGCTAAACTCTGGTGTTTCACtt                                                                                                                                                                                      <  1:232647/250‑1 (MQ=255)
                                                                   tGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGGCCGTTCGGTATTCAATCGCAACCGCCCCCCACAAACGGCATAGGTCGTTCCCCCGTCGCTAACCTCTGGCGTTTCCCTTGc                                                                                                                                                                                    >  2:53615/1‑251 (MQ=255)
                                                                      ggTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGGGGCCTGTTCGGTATTCAATCGCAACCCCTCCCCCGAAACCGCATAATTCGTTCGCCAGTCGCTAAAACCTGGGGTTTCACTTTCTa                                                                                                                                                                                  >  2:412200/1‑250 (MQ=255)
                                                                       gtgtCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAGTTCGTTCGCCAGTCGCTAAACTCTGGTGTTTCACTTGCTaaa                                                                                                                                                                                <  1:218275/251‑1 (MQ=255)
                                                                                cGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGGCCGTTTGGTGTTCTATCCCAAACGCGCACCCGAAACTGGATAGTTCGTTCGCCCGTCGCTCACCTCTGGTGTTTTTCTTGCTTAATcacacac                                                                                                                                                                        >  2:185780/1‑249 (MQ=255)
                                                                                    ttCCGGCACAATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTCCCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAGTTCGTTCGCCAGTCGCTAAACTCTGGTGTTTCACTTGCTAAATCACACAGACCt                                                                                                                                                                   <  2:436542/250‑1 (MQ=255)
                                                                                     tCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGGATTCAATCGCAACCCCTCACCAGAAACGGCATAGATCGTTCGCCAGTCGCTAAACTCTGGGGTTTCACTTTCTAAATCACACAGACCttt                                                                                                                                                                 >  1:490474/1‑251 (MQ=255)
                                                                                      ccGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGGATTCAATCGCAACCGCTCACCAGAAACGGCATAGTTCGTTCCCCAGTCGCTAAACTCTGGTGTTTCCCTTTCTAAATCACACAGACCtttt                                                                                                                                                                >  1:126365/1‑251 (MQ=255)
                                                                                        ggCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGGATTCAATCGCAACCCCTCACCCGAAAAGGCATAGTTTTTTTCCCAGTCGCCAAACTCTGGTGGTTCCCTTGCTAAATCCCACAAACCtttttt                                                                                                                                                              >  1:132621/1‑251 (MQ=255)
                                                                                          cACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGATGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGGGGCCTGGTTGGGATTCAATCGCAAAAGCTCAACAGGAAACGGCTAGTTCGTTTTCCACTCGCTAAAATCTGGTGGTTCAACTGCTAAATCACACAGACCtttttttt                                                                                                                                                            >  1:368681/1‑249 (MQ=255)
                                                                                              aTCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCCGTTCGGGCTTTAATCGCAACCGCTCACCAGCAACGGCATAGTTCGTTCGCCACTCGCTAAACTCTTGTGTTTCACTTTCTAAAACAAACAGACCCTTTTTGCGt                                                                                                                                                          >  2:107829/1‑249 (MQ=255)
                                                                                              aTCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAACTGCGCCAACGGCGGTGGCCCGTTCGGGATTCAAACGCAAACCCTCACCAGCAAAGGCATAGTTCGTTCGCCACGCGCTAAACCCTGGTGGTTCACTTTCTCAATCACACAAAACCttttttcgc                                                                                                                                                           >  2:397679/1‑246 (MQ=255)
                                                                                                cACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCGGAGATAAACCTCCGTTTTACCGGAGCCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACTGCATAGTTCGTTCGCCAGTCGCTAAACTCTGGTGTTTCACTTGCTAAATCACACAGACCTTTTTTGCGTAgcg                                                                                                                                                      <  1:74724/251‑1 (MQ=255)
                                                                                                  ccAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCCGGTCGGTATTCAATCGCAACCGCTCACCAGGAACGGCATAGTTTGTTCGCCAATACCCTAACACTGGTGGTTCACTTGCTAAAACACAAAGAACTTTTTTGCGGAGCGcc                                                                                                                                                    >  2:486356/1‑251 (MQ=255)
                                                                                                              tGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAGTTCGTTCGCCAGTCGCTAAACTCTGGTGTTTCACTTGCTAAATCACACAGACCTTTTTTGCGTAGCGCCTGCAACGCGGc                                                                                                                                         <  1:296342/250‑1 (MQ=255)
                                                                                                                 ccTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCCCCAGAAACGGCATAGTTCGTTCGCCAGGCGCTAAACTCTGGTGTTTTACTTGCTAAATCACACAGACCTTTTTTGCGTAGCGCCTGCAACGCGGcctc                                                                                                                                      >  1:491881/1‑247 (MQ=255)
                                                                                                                  cTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCCTAGTTCGTTCGCCAGGCGCTAAACTCTGGTGTTTCACTTTCTAAATCACACAGACCTTTTTTTCGTAGCGCCTGCAACGCGGcatcat                                                                                                                                    >  1:185209/1‑251 (MQ=255)
                                                                                                                   tGAGCGAGCAAGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGGATTCAATCCCAACCCCCCACCACAAACGGCATAGTTCGTTCGCCAGTCGCTAAACTCTGGTGGTTCACTTGCTAAATCAAACAGACCTTTTTTTCGTAGCGCCTGCAACGCGGCATCAtt                                                                                                                                   >  1:178753/1‑251 (MQ=255)
                                                                                                                    gagcgagcACGTTTTCCAGTACGCTGAGATAACCCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAGTTCGTTCGCCAGTCGCTAAACTCTGGTGTTTCACTTGCTAAATCACACAGACCTTTTTTGCGTAGCGCCTGCAACGCGGcatcat                                                                                                                                    <  2:255433/249‑1 (MQ=255)
                                                                                                                          gcACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAGTTCGTTCGCCAGTCGCTAAACTCTGGTGTTTCACTTGCTAAATCACACAGACCTTTTTTGCGTAGCGCCTGCAACGCGGCATCATTAAATTCg                                                                                                                            <  1:170244/251‑1 (MQ=255)
                                                                                                                            aCGTTTTCAAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAGTTCGTTCGCCAGTCGCTAAACTCTGGTGTTTCACTTGCTAAATCACACAGACCTTTTTTGCGTAGCGCCTGCAACGCGGCATCATTAAATTCGAg                                                                                                                          <  2:437205/251‑1 (MQ=255)
                                                                                                                             cGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTTGGGATTCAATCGCAACCGCTCACCAGAAACGGCATAGTTCGTTCCCCAGTCGCTAAAATCTGGTGTTTCACTTGCTAAATCACACAGACCTTTTTTTCGTAGCGCCTGCAACGCGGCATCATTAAATTCGAGc                                                                                                                         >  1:79622/1‑251 (MQ=255)
                                                                                                                                 ttCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGGATTCCATCGCAAACGCTCACCAGAAACGGCATAGTTCGTTCCCCCGTCGCTAAACTCTTGTGGTTCACTTGCTAAATCACACAGACCTTTTTTTCGTAGCGCCTGCAACGCGGCATCATTAAATTCGATCGTgg                                                                                                                     >  2:198072/1‑251 (MQ=255)
                                                                                                                                  tCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAGTTCGTTCGCCAGTCGCTAAACTCTGGTGTTTCACTTGCTAAATCACACAGACCTTTTTTGCGTAGCGCCTGCAACGCGGCATCATTAAATTCGAGCGTgg                                                                                                                     <  1:121127/250‑1 (MQ=255)
                                                                                                                                   ccAGTACGCTGAGATAACCCTCCGTTTTCCCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAGTTCGTTCGCCAGTCGCTAAACTCTGGTGTTTCACTTGCTAAATCACACAGACCTTTTTTGCGTAGCGCCTGCAACGCGGCATCATTAAATTCGAGCGTGGCg                                                                                                                   <  2:341665/251‑1 (MQ=255)
                                                                                                                                                aTAAACCTACGATTTACCGGAACCGGTAACGCCAGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTACGCCAACGGCGGTTGGCCGGTCGGTATTCAATCGCAACCCCTCACCAGAAACGCCATAGATCGTTCGCCCTTCGCTAAACTCTTGTGTTTCACTTGCTACATACCACACACCTTTTTTGCTGAGTCCCTGCAACTCGGCCACATTAAATTAGAGCCTGGCgcaccgggcgcgcc                                                                                                      >  2:294269/1‑238 (MQ=255)
                                                                                                                                                   aaCCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTTGGGATTCAATCGCAACCGCTCACCCGAAAAGGCATAGTTCGTTCGCCCGTCTCTAAACTCTGGGGTTTTACCTGCTAAATCACAAAGACCTTTTTTGCGTAGCGCCTTCAACGGTGCATCATTAAATTCGAGCGTGGCGAACTGGTCGCGcccg                                                                                                    >  2:583400/1‑248 (MQ=255)
                                                                                                                                                    aCCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCCGTTCGGGATTCAAACGCCAACGCCCACCAGAAACGGCAAAGTTCGTTCGCCCGGCCCTAAAATCTTGGGTTTCACTTGCTTAATTACCCAGACCTTTTTTTCGTAGGGGCCGCCAAACCGCGACATTAAATTCGAGCGTGGCGACCTGGTCGCGcccgct                                                                                                  >  2:47444/1‑247 (MQ=255)
                                                                                                                                                     ccTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGGGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCCTAGTTCGTTCCCCAGTCGCTCAACCCTGGTGTTTCACTTGCTAAATCCCACAGACCTTTTTTGCGTAACGCCCGCAAAGCGGCATCCTTAAATTTGAGCGTGGGGACCTGGGAGCGCCAGAtt                                                                                                 >  1:281317/1‑251 (MQ=255)
                                                                                                                                                       cccGTTTTACCGGACCCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAGTTCGTTCGCCAGTCGCTAAACTCTGGTGTTTCACTTGCTAAATCACACAGACCTTTTTTGCGTAGCGCCTGCAACGCGGCATCATTAAATTCGAGCGTGGCGACCTGGTCGCGCCAGAtttt                                                                                               <  2:470086/250‑1 (MQ=255)
                                                                                                                                                                ccGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGGATTTAATCGCAACCGCTCCCCAGAAACGGCCTAGTTCTTTCCCCCGTCGCTAAACTCTGGTGTTTTACTTGCTAAATCACACCAACCTTTTTTTCGTCGCGCCTGCAACGCGGCACCATTAAATTCGCGCGGGGCGGCCTGGTCGCGCCAGATTtttcctttccc                                                                                     >  2:155947/1‑247 (MQ=255)
                                                                                                                                                                 cGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCCGTTCGGTTTTCAATCGCAACCGCCCACCAGAAACGGCATAGTTCGTTCTCCAGTCGCTTAACTCTGGGGTTTTACTTGCTAAATCACACAGACCTTTTTTGCGTAGCGCCTGCAACGCGGCATCATTAAATTCGAGCGTGGCGACCTGGTCGCGCCCGATTTTGCCTTGCCgg                                                                                    >  1:300565/1‑250 (MQ=255)
                                                                                                                                                                    aaCCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAGTTCGTTCGCCAGTCGCTAAACTCTGGTGTTTCACTTGCTAAATCACACAGACCTTTTTTTCGTAGCGCCTGCAACGCGGCATCATTAAATTCGAGCGTGGCGACCTGGTCGCGCCCGATTTTGCCTTGCCGTAAc                                                                                 >  1:293909/1‑251 (MQ=255)
                                                                                                                                                                       cGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGACTGTTCGGGATTCAATCGCAACCGCTCACCAGAAACGGCATAGTTCGTTCGCCAGTCGCTAAACTCTGGGGTTTCCCTTGCTAAAACACACAAACCTTTTTTGCGTAGCGCCTGCAACTCGGCATCATTAAATTCGAGCGTGTGGACCTGGGGGCGCAAGAGTTTTCCTTTCAGTAAcggc                                                                              >  1:170642/1‑249 (MQ=255)
                                                                                                                                                                           aaCGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGGTCGGTATTTAATCGCAAACGCTCACCAGAAACGGCATCGGTCGTTCGCCAGTCGCTAAACTCTGGGGTTTCACTTGCTTAATCTCACACACCCTTTTTGCGGAGCGCCCGCAACGCGGCATCATTACATTCGAGCGTGGCGACCTGTTCGCTCCAGATTTTGCCTTGCCCTAACgccgcc                                                                           >  2:455604/1‑250 (MQ=255)
                                                                                                                                                                              gcccgccAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACCGCGGCGGACAGATCGGGCTTTAATCGCAACCGCTCCCCCGGCAACGCCTAGTTCTTTTTCCCGCCGCCAAACTCTGCTGTTTTCCTTTCTAAATCACACACACCTTTTTTTCGTTGCGCCTTCCCCTCCGCATCATTTAATTTGAGCGTTGCGACCTCGCCCCGCCCGCTTTTGCCTTTCCTTTACgccgcccgcc                                                                       >  1:56982/1‑247 (MQ=255)
                                                                                                                                                                                ccgccAGCAGCCAGGCAGAACAAGGATCTGCCGCGCTATGAATTGCGCCAACGGCGGGTGCCTTTTCGGGATTCACTCGCCACCGCTCACCAGAAACGCCGTAGTTTGTTTGCCCGCCGCCTCACCCTTGTGTTTCACCTCCTAAATCACACAGAACTTTTTTGCGTAGCGCCCGCAACGCGGCACCATTAAATTCGCGCGTGGCGACCTTGTCGCGCCCGATTTTTCCTTGCCGTAACACCCCCAGCGcc                                                                     >  2:450742/1‑251 (MQ=255)
                                                                                                                                                                                  gccAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAGTTCGTTCGCCAGTCGCTAAACTCTGGTGTTTCACTTGCTAAATCACACAGACCTTTTTTGCGTAGCGCCTGCAACGCGGCATCATTAAATTCGAGCGTGGCGACCTGGTCGCGCCAGATTTTGCCTTGCCGTAACGCCGCCAGCGCCtg                                                                   <  2:126365/251‑1 (MQ=255)
                                                                                                                                                                                       cagcCAGGCAGAAAAAGTATCTGCCGCGCTATGCATTGCGCCAACGGCGGTGGCCCGTTCGGTATTCCATCGCAACCGCCCCCCAGAAACGGCATAGTTCGTTCGCCCGGCGCTAAACTCCGGTGGTTCACTTGCTAAATCAACCACACCTTTTTTTCGTAGCGCCTGCCACGCGGCATCATTAAATTCGAGCGTGGCGACCTGGTCGCGCCAGATTTTGCCTTGCCGTAACGCCGCCACCGCCTgttgtt                                                              >  2:526717/1‑251 (MQ=255)
                                                                                                                                                                                                 gAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGGATTCAATCGCAACCGCCCACCAGAAACGGCATAGTTCGTTCGCCAGTCGCTAAACTCTGGTGTTTCACTTGCTAAATCACACAGACCTTTTTTGCGTAGCGCCTGCAACGCGGCATCATTAAATTCGAGCGGGGCGACCTGGTCGCGCCAGATTTTGCCTTGCCGTAACGCCGCCAGCGCCCGTTGTTGCTTtggggg                                                    >  1:151661/1‑250 (MQ=255)
                                                                                                                                                                                                  aaaaaGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGGCCGTTCGGTATTCCATCCCCAACGCTCACCCGCAAAGGCAGAGTTCGTTCGCCAGGCGCTAAACTCTGGTGTTTCACCTTCTAAATCACACCGACCTTTTTTGCGTCGCGCCCGCAACGCGGCATCCTTAAATTCGAGGGGGGCGACCTGGTCCCGCCCGATTTTGGCCTGGCCTTAAGACGCCCCCGCCCGGTGTTGCTTggggg                                                     >  2:135144/1‑249 (MQ=255)
                                                                                                                                                                                                  aaaaaGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAGTTCGTTCGCCAGTCGCTAAACTCTGGTGTTTCACTTGCTAAATCACACAGACCTTTTTTGCGTAGCGCCTGCAACGCGGCATCATTAAATTCGAGCGTGGCGACCTGGTCGCGCCAGAtt                                                                                                 <  2:378828/206‑1 (MQ=255)
                                                                                                                                                                                                  aaaaaGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAGTTCGTTCGCCAGTCGCTAAACTCTGGTGTTTCACTTGCTAAATCACACAGACCTTTTTTGCGTAGCGCCTGCAACGCGGCATCATTAAATTCGAGCGTGGCGACCTGGTCGCGCCAGATTTTGCCTTGCCGTAACGCCGCCAGCGCCTGTTGTTGCTTggggg                                                     >  1:114219/1‑249 (MQ=255)
                                                                                                                                                                                                  aaaaaGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAAAGGCATAGTTCGTTCCCCCGTCGCTAAACTCTGGTGTTTCACTTGCTAAATCACACAAACCTTTTTTGCGTAGCGCCTGCAACGCGGCACCATTAAATTCGCGCGTGGCGACCTGGTCGCGCCAGAtt                                                                                                 >  1:378828/1‑206 (MQ=255)
                                                                                                                                                                                                       gTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCCGTTCGGGATTCAAACGCAACCGCTCACCCGAAACGGCATAGTTTGTTCGCCAGTCGCTAAACTCTGGTGTTTCACTTGCTAAATCACCCAGACCTTTTTTTCGTTGCGCCTGCAACGCGGCATCATTAAATTCCAGCGGGGCGACCTGGTCGCGCCCGATTTTGTCTTGCCGTAACGCCGCCCGCGCCTGTTGGTTCCTGGGGGAGCGttt                                              >  2:474914/1‑251 (MQ=255)
                                                                                                                                                                                                                   cTATGAATTGCGCCAACGGCGGTGGCCCGTTCGGGATTTAATCGCAACCGCTCACCAGAAACGGCATAGTTCGGTCGCCCGTCGCTAAACTCTGGTGTTTTCCTTGCTAAATCACACAGACCTTTTTTGCGTAGCGCCCGCAACGCGGCATCATTAAATTCGAGCGTGGCGACCTGGTCGCGCCAGATTTTGCCTTTCCGTCACGCCGCCAGCGCCCGTTTTTGCTTGGGGGGGCGGTTCAGGCTGTTc                                    >  2:269785/1‑249 (MQ=255)
                                                                                                                                                                                                                       gAATTGCGCCAACGGCGGTGGCCTGTTCGGGATTCAATCGCAACCGCTCACCAGCAACGGGATAGTTCGTTAGCCAGTCGCTAAACTCTGGTGTTTCACTTGCTAAATCACACAGACCTTTTTTGCGTAGCGCCTGCAACGCGGCATCATTAAATTCGAGCGTGGCGACCTGGTCGCGCCAGATTTTTCCTTTCCGTAACGCCGCCCGCGCCTGGTGTTGCTTGGGGGAGCGGTTCAGGCTGTTCAGATcc                              >  1:362109/1‑251 (MQ=255)
                                                                                                                                                                                                                           ggcgcCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAGTTCGTTCGCCAGTCGCTAAACTCTGGTGTTTCACTTGCTAAATCACACAGACCTTTTTTGCGTAGCGCCTGCAACGCGGCATCATTAAATTCGAGCGTGGCGACCTGGTCGCGCCAGATTTTGCCTTGCCGTAACGCCGCCAGCGCCTGTTGTTGCTTGGGGGAGCGTTTCAGGCTGTTCAGATCCACCg                          <  2:182879/250‑1 (MQ=255)
                                                                                                                                                                                                                            gcgcCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAGTTCGTTCGCCAGTCGCTAAACTCTGGTGTTTCACTTGCTAAATCACACAGACCTTTTTTGCGTAGCGCCTGCAACGCGGCATCATTAAATTCGAGCGTGGCGACCTGGTCGCGCCAGATTTTGCCTTGCCGTAACGCCGCCAGCGCCTGTTGTTGCTTGGGGGAGCGTTTCAGGCTGTTCAGATCCACCGc                         <  1:412200/251‑1 (MQ=255)
                                                                                                                                                                                                                            gcgcCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAGTTCGTTCGCCAGTCGCTAAACTCTGGTGTTTCACTTGCTAAATCACACAGACCTTTTTTGCGTAGCGCCTGCAACGCGGCATCATTAAATTCGAGCGTGGCGACCTGGTCGCGCCAGATTTTGCCTTGCCGTAACGCCGCCAGCGCCTGTTGTTGCTTGGGGGAGCGTTTCAGGCTGTTCAGATCCACCGc                         <  1:572590/251‑1 (MQ=255)
                                                                                                                                                                                                                                    ggcggTGGCCTGTTCGTTATTCAATCGCAACCGCTCACCAGAAACGGCATAGTTCGTTCGCCAGTCGCTAAACTCTGGTGTTTCACTTGCTAAATCACACAGACCTTTTTTGCGTAGCGCCTGCAACGCGGCATCATTAAATTCGAGCGTGGCGACCTGGTCGCGCCAGATTTTGCCTTGCCGTAACGCCGCCAGCGCCTGTTGTTGCTTGGGGGAGCGTTTCAGGCTGTTCAGATCCACCGCCTGGCCtt                 <  2:300160/251‑1 (MQ=255)
                                                                                                                                                                                                                                    ggcggTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAGTTCGTTCGCCAGTCGCTAAACTCTGGTGTTTCACTTGCTAAATCACACAGACCTTTTTTGCGTAGCGCCTGCAACGCGGCATCATTAAATTCGAGCGTGGCGACCTGGTCGCGCCAGATTTTGCCTTGCCGTAACGCCGCCAGCGCCTGTTGTTGCTTGGGGGAGCGTTTCAGGCTGTTCAGATCCACCGCCTGGCCtt                 <  2:406349/251‑1 (MQ=255)
                                                                                                                                                                                                                                    ggaggTGGCCTGTTCGGTATGCAATCGCACCCGCGCACCAGAAACGGCATAGTTCGTTCGCCAGTCGCTAAACTCTGGTGTTTCACTTGCTAAATCACACAGACCTTTTTTGCGTAGCGCCTGCAACGCGGCATCATTAAATTCGAGCGTGGCGACCTGGTCGCGCCAGATTTTGCCTTGCCGTAACGCCGCCAGCGCCTGTTGTTGCTTGGGGGAGCGTTTCAGGCTGTTCAGATCCACCGCCTGGCCtt                 <  2:500617/248‑1 (MQ=255)
                                                                                                                                                                                                                                                    ggTATTCAATCGCAACCGCTCACCAGAAACGGCATAGTTCGTTCGCCAGTCGCTAAACTCTGGTGTTTCACTTGCTAAATCACACAGACCTTTTTTGCGTAGCGCCTGCAACGCGGCATCATTAAATTCGAGCGTGGCGACCTGGTCGCGCCAGATTTTGCCTTGCCGTAACGCCGCCAGCGCCTGTTGTTGCTTGGGGGAGCGTTTCAGGCTGTTCAGATCCACCGCCTGGCCTTGTTCAGTGGCAAAcc  >  2:574614/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                              
GTCGTTCAGGCCGGAATGCAGAACTTCCACGGGGGCATTAAAACGTTCACGAAAACGGGCGATAGTTTGCGGTGTCAGGCCGATTTCCGGCACCATCACCAGCGCCTGTTTGCCCTGAGCGAGCACGTTTTCCAGTACGCTGAGATAAACCTCCGTTTTACCGGAACCGGTAACGCCCGCCAGCAGCCAGGCAGAAAAAGTATCTGCCGCGCTATGAATTGCGCCAACGGCGGTGGCCTGTTCGGTATTCAATCGCAACCGCTCACCAGAAACGGCATAGTTCGTTCGCCAGTCGCTAAACTCTGGTGTTTCACTTGCTAAATCACACAGACCTTTTTTGCGTAGCGCCTGCAACGCGGCATCATTAAATTCGAGCGTGGCGACCTGGTCGCGCCAGATTTTGCCTTGCCGTAACGCCGCCAGCGCCTGTTGTTGCTTGGGGGAGCGTTTCAGGCTGTTCAGATCCACCGCCTGGCCTTGTTCAGTGGCAAACC  >  NC_000913/4125965‑4126458

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