Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,634,3700TC54.5% ‑2.1 / 40.2 33V172V (GTA→GTGtfaQputative tail fiber assembly protein TfaQ
Reads supporting (aligned to +/- strand):  ref base T (12/3);  new base C (14/4);  total (26/7)
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.82e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGTCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAAGAGCCAAAGCAAAACTTCAGACACCTGACATATAAAAAATAATCTCGGTGTGAGATGCTTTACGTCTTCCAAGCCCCCTTCCTTGCCGTAAATGGAAAGATACATCTAATTATAGAATTTATATGTTTTACCCTACGGCAGTGCTGGCCATTCAATATCCTGTGCAGTTGACGTATCAACACGGTTCAGCAATACCCGATACTTTTTCCATGCTTCCAGCAACGAGATTTCTTCCTCCGTTGCAATTTCCAGATCTGCAGCATCCTGAAGCGGCGCAATATGCTCACTGGCTACCTGCATCAGGTTGTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCC  >  NC_000913/1634072‑1634670
                                                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                                            
gtgtTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGTCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAAGAGCCAAAGCAAAACTTCAGACACCTGACATATAAAAAATAATCTCGGTGTGAGATGCTTTACGTCTTCCAAGCCCCCTTCCTTGCCGTAAATGGAAAGATACATCTAATTATAGAATTTATATGTTTTACCCTACGGCAGTGCTGGCCATTCAATATCCTGTGCAGTTGACGTATCAACACGGTTCAGCAATa                                                                                                                                                                                                                                                                                                             >  1:864695/1‑300 (MQ=31)
          aagaaCAGAAACAGTCGATATTCGAACGAATTAAGTCAGGGGTAAGGATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGACAAACCATTTTAAGAGCCAAAGCAAAACTTCAGACACCTGACATATAAAAAATAATCTCGGTGTGAGATGCTTTACGTCTTCCAAGCCCCCTTCCTTGCCGTAAATGGAAAGATACATCTAATTATAGAATTTATATGTTTTACCCTACGGCAGTGCTGGCCATTCAATATCCTGTGCAGTTGACGTATCAACACGGTTCAGCAATACCCGagac                                                                                                                                                                                                                                                                                                     <  2:176197/298‑4 (MQ=9)
            gaaCAGAACCAGGCGGTATTCGAACGAATTAAGTCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAAGAGCCAAAGCAAAACTTCAGACACCTGACATATACAAAATAATCTCGGTGTGAGATGCTTTACGTCTTCCAAGCCCCCTTCCTTGCCGTAAATGGAAAGATACATCTAATTATAGAATTTATATGTTTTACCCTCCGGCAGTGCTGGCCATTCAATATCCTGTGCAGTTGACGTATCAACAAGGTTCAGCAATACCCGCTACtttt                                                                                                                                                                                                                                                                                                 >  1:376584/1‑300 (MQ=9)
            gaaCAGAAACAGGCGGTATTCGAACGAATTAAGTCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAAGAGCCAAAGCAAAACTTCAGACACCTGACATATAAAAAATAATCTCGGTGTGAGATGCTTTACGTCTTCCAAGCCCCCTTCCTTGCCGTAAATGGAAAGATACATCTAATTATAGAATTTATATGTTTTACCCTACGGCAGTGCTGGCCATTCAATATCCTGTGCAGTTGACGTATCAACACGGTTCAGCAATACCCGATACtttt                                                                                                                                                                                                                                                                                                 >  1:577759/1‑300 (MQ=31)
                 gAAACAGGCGGTATTCGAACGAATTAAGTCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAAGAGCCAAAGCAAAACTTCAGACACCTGACATATAAAAAATAATCTCGGTGTGAGATGCTTTACGTCTTCCAAGCCCCCTTCCTTGCCGTAAATGGAAAGATACATCTAATTATAGAATTTATATGTTTTACCNTACGGCAGTGCTGGCCATTCAATATCCTGTGCAGTTGACGTATCAACACGGTTCAGCAATACCCGATACTTTTTCCAt                                                                                                                                                                                                                                                                                            >  1:187471/1‑300 (MQ=14)
                                                                                                                                                                                                                                                                             gCAGTTGATGTATCAACACGGTTCAGCAACACCCGATACTTCTTCCAGGCTTCCAGCAACGAGGTTTCTTCCTCCGTTGCAATTTCCAGATCTGCAGCATCCTGAAGCGGCGCAATATGCTCACTGGCTACCTGCATCAGGCTTTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGATATTCCCATACCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTc                                 >  2:402219/1‑299 (MQ=11)
                                                                                                                                                                                                                                                                              cAGTTGATGTATCAACACGGTTCAGCAACACCCGATACTTCTTCCAGGCTTCCAGCAACGAGGTTTCTTCCTCCGTTGCAATTTCCAGATCTGCAGCATCCTGAAGCGGCGCAATATGCTCACTGGCTACCTGCATCAGGCTTTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGc                               >  1:801820/1‑300 (MQ=11)
                                                                                                                                                                                                                                                                                gTTGATGTATCAACACGGTTCAGCAACACCCGATACTTCTTCCAGGCTTCCAGCAACGAGGTTTCTTCCTCCGTTGCAATTTCCAGATCTGCAGCATCCTGAAGCGGCGCAATATGCTCACTGGCTACCTGCATCAGGCTTTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGcc                              >  1:158435/1‑299 (MQ=11)
                                                                                                                                                                                                                                                                                       tATCAACACGGTTCAGCAATACCCGATACTTTTTCCATGCTTCCAGCAACTAGATTTCTTCCTCCGTTGCAATTTCCAGATCTGCAGCATCCTTAAGCGGCACAATATGCTCACTGGCTACCTGCATCAGTTTGTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGc                        <  2:190987/298‑1 (MQ=9)
                                                                                                                                                                                                                                                                                       tATCAACACGGTTCAGCAACACCCGATACTTCTTCCAGGCTTCCAGCAACGAGGTTTCTTCCTCCGTTGCAATTTCCAGATCTGCAGCATCCTGAAGCGGCGCAATATGCTCACTGGCTACCTGCATCAGGCTTTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCa                      >  1:547180/1‑300 (MQ=11)
                                                                                                                                                                                                                                                                                       tATCAACACGGTTCAGCAACACCCGATACTTCTTCCAGGCTTCCAGCAACGAGGTTTCTTCCTCCGTTGCAATTTCCAGATCTGCAGCATCCTGAAGCGGCGCAATATGCTCACTGGCTACCTGCATCAGGCTTTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCAGTTCCACTTCTGATATTCCACTACCGGCGATAACCAGGTAAAATTTTCAGGTAACGGAACGAGTTCAGAAATAAATAACGCGTCGAAGGAAGcaa                      >  2:480605/1‑298 (MQ=11)
                                                                                                                                                                                                                                                                                       tATCAACACGGTTCAGCAACACCCGATACTTCTTCCAGGCTTCCAGCAACGAGGTTTCTTCCTCCGCTGCAATTTCCAGATCTGCAGCATCCTGAAGCGGCGCAATATGCTCACTGGCTACCTGCATCAGGCTTTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTACCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGACGGAAGcc                       >  1:836347/1‑299 (MQ=11)
                                                                                                                                                                                                                                                                                            acacGGTTCAGCAATACCCGATACTTTTTCCATGCTTCCAGCAACGAGATTTCTTCCTCCGTTGCAATTTCCAGATCTGCAGCATCCTGAAGCGGCGCAATATGCTCACTGGCTACCTGCATCAGGTTGTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCa                 >  1:162129/1‑300 (MQ=31)
                                                                                                                                                                                                                                                                                            acacGGTTCAGCAACACCCGATACTTCTTCCAGGCTTCCAGCAACGAGGTTTCTTCCTCCGTTGCAATTTCCAGATCTGCAGCATCCTGAAGCGGCGCAATATGCTCACTGGCTACCTGCATCAGGCTTTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCa                 >  1:888456/1‑300 (MQ=11)
                                                                                                                                                                                                                                                                                               cGGTTCAGCAACACCCGATACTTCTTCCAGGCTTCCAGCAACGAGGTTTCTTCCTCCGTTGCAATTTCCAGATCTGCAGCATCCTGAAGCGGCGCAATATGCTCACTGGCTACCTGCATCAGGCTTTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATa               >  2:900745/1‑299 (MQ=11)
                                                                                                                                                                                                                                                                                               cGGTTCAGCAACACCCGATACTTCTTCCAGGCTTCCAGCAACGAGGTTTCTTCCTCCGTTGCAATTTCCAGATCTGCAGCATCCTGAAGCGGCGCAATATGCTCACTGGCTACCTGCATCAGGCTTTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAAAGCGTCGCCGGAAGCCACGTCATa               >  2:455703/1‑299 (MQ=11)
                                                                                                                                                                                                                                                                                                 gTTCAGCAATACCCGATACTTTTTCCATGCTTCCAGCAACGAGATTTCTTCCTCCGTTGCAATTTCCAGATCTGCAGCATCCTGAAGCGGCGCAATATGCTCNCTGGCTACCTGCATCAGGTTGTTTTTTGTTTCNTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCAAGTCATAAAc            >  2:711227/1‑300 (MQ=14)
                                                                                                                                                                                                                                                                                                 gTTCAGCAATACCCGATACTTTTTCCATGCTTCCAGCAACGAGATTTCTTCCTCCGTTGCAATTTCCAGATCTGCAGCATCCTGAAGCGGCGCAATATGCTCACTGGCTACCTGCATCAGGTTGTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGACACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAAc            >  1:115682/1‑300 (MQ=14)
                                                                                                                                                                                                                                                                                                 gTTCAGCAATACCCGATACTTTTTCCATGCTTCCAGCAACGAGATTTCTTCCTCCGTTGCAATTTCCAGATCTGCAGCATCCTGAAGCGGCGCAATATGCTCACTGGCTACCTGCATCAGGTTGTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAAc            >  2:171706/1‑300 (MQ=31)
                                                                                                                                                                                                                                                                                                 gTTCAGCAATACCCGATACTTTTTCCATGCTTCCAGCAACGAGATTTCTTCCTCCGTTGCAATTTCCAGATCTGCAGCATCCTGAAGCGGCGCAATATGCTCACTGGCTACCTGCATCAGGTTGTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAAc            >  1:12785/1‑300 (MQ=31)
                                                                                                                                                                                                                                                                                                 gTTCAGCAATACCCGATACTTTTTCCATGCTTCCAGCAACGAGATTTCTTCCTCCGTTGCAATTTCCAGATCTGCAGCATCCTGAAGCGGCGCAATATGCTCACTGGCTACCTGCATCAGGTTGTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCAGGAAGCCACGTCATAAAc            >  1:542028/1‑300 (MQ=14)
                                                                                                                                                                                                                                                                                                 gTTCAGCAATACCCGATACTTTTTCCATGCTTCCAGCAACGAGATTTCTTCCTCCGTTGCAATTTCCAGATCTGCAGCATCCTGAAGCGGCGCAATATGCTCACTGGCTACCTGCATCAGGTTGTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTATGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAAc            >  1:443958/1‑300 (MQ=14)
                                                                                                                                                                                                                                                                                                 gTTCAGCAATACCCGATACTTTTTCCATGCTTCCAGCAACGAGATTTCTTACTCCGTTGCAATTTCCAGATCTGCAGCATCCTGAAGCGGCGCAATATGCTCACTGGCTAACTGCATCAGGTTGTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCAGTTCCACTTATGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAATGGACCGAGTTCAGAAATAAATAACGCGTCGACGGAAGCAACGTAATAAAc            >  2:214966/1‑300 (MQ=11)
                                                                                                                                                                                                                                                                                                 gTTCAGCAACACCCGATACTTCTTCCATGCTTCCAGCAACGAGGTTTCTTCCTCCGTTGCAATTTCCAGATCTGCAGCATCCTGAAGCGGCGCAATATGCTCACTGGCTACCTGCATCATGCTTTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCTGCGATAACCAGGTAAAATTTTCCGGTAACGGATCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATagac            <  1:880219/300‑4 (MQ=11)
                                                                                                                                                                                                                                                                                                 gTTCAGCAACACCCGATACTTCTTCCAGGCTTCCAGCAACTAGGTTTCTTCCTCCTTTGCAATTTCCAGATCTGCAACATCCTGAAGCGGCGCAATATGCTCACTGGCTACCTGCATCAGGCTTTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATagac            <  2:487244/300‑4 (MQ=11)
                                                                                                                                                                                                                                                                                                 gTTCAGCAACACCCGATACTTCTTCCAGGCTTCCAGCAACGAGGTTTCTTCCTCCGTTGCAATTTCCAGATCTGCAGCATCCTGAAGCGGCGCAATATGCTCNCTGGCTACCTGCATCAGGCTTTTTTTTGTTTCNTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATagac            >  2:713252/1‑297 (MQ=11)
                                                                                                                                                                                                                                                                                                 gTTCAGCAACACCCGATACTTCTTCCAGGCTTCCAGCAACGAGGTTTCTTCCTCCGTTGCAATTTCCAGATCTGCAGCATCCTGAAGCGGCGCAATATGCTCACTGGCTACCTGCATCAGGCTTTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATagac            >  1:141630/1‑297 (MQ=11)
                                                                                                                                                                                                                                                                                                 gTTCAGCAACACCCGATACTTCTTCCAGGCTTCCAGCAACGAGGTTTCTTCCTCCGTTGCAATTTCCAGATCTGCAGCATCCTGAAGCGGCGCAATATGCTCACTGGCTACCTGCATCAGGCTTTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGACTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCAAGTCAtcgac            >  2:409832/1‑296 (MQ=11)
                                                                                                                                                                                                                                                                                                 gTTCAGCAACACCCGATACTTCTTCCAGGCTTCCAGCAACGAGGTTTCTTCCTACGTTGCAATTTCCAGATCTGCAGCATCCTGAAGCGGCGCAATATGCTCACTGGCTACCTGCATCAGGCTTTTTTTTGTTTCTTCCGCCTCACGGATCCGGAACCGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCAGGAAGCCACGTCATaga             >  2:251341/1‑297 (MQ=11)
                                                                                                                                                                                                                                                                                                 gTTCAGCAACACCCGATACTTCTTCCAGGCTTCCAGCAACGAGGTTTCTTCATCCGTTGCAATTTCCAGATCTGCAGCATCCTGAAGCGGCGCAATATGCTCACTGGCTACCTGCATCAGGCTTTTTTTTGTTTCTTCCGCCTCCCGGATCCGGCACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCAACGTCATagac            >  2:408375/1‑297 (MQ=11)
                                                                                                                                                                                                                                                                                                  ttCAGCAATACCCGATACTTTTTCCATGCTTCCAGCAACTAGATTTCTTCCTCCGTTGCAATTTCCAGATCTGCAGCATCCTGAAGCGGCGCAATATGCTCACTGGCTACCTGCATCAGTTTGTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAAc            <  2:550804/299‑1 (MQ=14)
                                                                                                                                                                                                                                                                                                     aGCAACACCCGATACTTCTTCTAGGCTTCCAGCAACGAGGTTTCTTCCTCCGTTGCAATTTCCAGATCTGCAGCATCCTGAAGCGGCGCAATATACTCACTGGCTACCTGCATCAGGCTTTTTTTTGTTCCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATagacggt         <  2:513407/299‑7 (MQ=11)
                                                                                                                                                                                                                                                                                                     aGCAACACCCGATACTTCTTCCAGGCTTCCAGCAACGAGGTTTCTTCCTCCGTTGCAATTTCCAGATCTGCAGCATCCTGAAGCGGCGCAATATGCTCACTGGCTACCTGCATCAGGCTTTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATagac            <  2:487537/296‑4 (MQ=11)
                                                                                                                                                                                                                                                                                                          tACCCGATACTTTTTCCATGCTTCCAGCAACTAGATTTCTTCCTCCGTTGCAATTTCCAGATCTGCAGCATCCTGAAGCGGCGCAATATGCTCACTGGCTACCTGCATCAGGTTGTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGTTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTAccc  <  1:131557/300‑1 (MQ=14)
                                                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                                            
GTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGTCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAAGAGCCAAAGCAAAACTTCAGACACCTGACATATAAAAAATAATCTCGGTGTGAGATGCTTTACGTCTTCCAAGCCCCCTTCCTTGCCGTAAATGGAAAGATACATCTAATTATAGAATTTATATGTTTTACCCTACGGCAGTGCTGGCCATTCAATATCCTGTGCAGTTGACGTATCAACACGGTTCAGCAATACCCGATACTTTTTCCATGCTTCCAGCAACGAGATTTCTTCCTCCGTTGCAATTTCCAGATCTGCAGCATCCTGAAGCGGCGCAATATGCTCACTGGCTACCTGCATCAGGTTGTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCC  >  NC_000913/1634072‑1634670

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