Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,745,0030AG40.7% 82.1 / 17.4 59E67G (GAA→GGA) yiaN2,3‑diketo‑L‑gulonate TRAP transporter large permease protein
Reads supporting (aligned to +/- strand):  ref base A (8/27);  new base G (24/0);  total (32/27)
Fisher's exact test for biased strand distribution p-value = 5.12e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCCTCT  >  NC_000913/3744766‑3745229
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                         
aGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGa                                                                                                                                                                                                                                           <  1:308444/238‑1 (MQ=255)
       ctctAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGa                                                                                                                                                                                                                                  <  1:2252654/239‑1 (MQ=255)
            aCGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGGAATCATGGATtcg                                                                                                                                                                                                                            >  1:1139043/1‑239 (MQ=255)
               tCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGGAAACAttcatggtggg                                                                                                                                                                                                                          >  2:2031599/1‑229 (MQ=255)
                 gCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGGAATCATGAATtcggggg                                                                                                                                                                                                                        >  1:396176/1‑236 (MQ=255)
                  ccACCGCAAGGAGATTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGGAATCATGAATTCCgggggc                                                                                                                                                                                                                      >  1:787617/1‑239 (MQ=255)
                      cGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGGAAACATGAATgggggggg                                                                                                                                                                                                                       >  1:1031667/1‑234 (MQ=255)
                         aaGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCa                                                                                                                                                                                                              <  1:859539/238‑1 (MQ=255)
                                     cATGGCGGTGCTGATTTTTCTGGGCTGTTGGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCGGGCAACTGTTGTTGGGCGGGGCGGCGTAAATGTTCTGACGGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGACGGGAGCGGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATAAAGGCGGGGGGCTTGTAAAAACGCAttgttg                                                                                                                                                                                                <  2:1988278/239‑1 (MQ=255)
                                       tGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGGAATCATTAAAGCGGGCGGGTTTTCCAAAAGCCTTGTTGa                                                                                                                                                                                               >  1:571517/1‑238 (MQ=255)
                                       tGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATgcggg                                                                                                                                                                                                                          <  1:613957/214‑1 (MQ=255)
                                       tGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGAAATCATGAATgcggg                                                                                                                                                                                                                          >  2:613957/1‑214 (MQ=255)
                                        gtctttgctattttgtCGGGGTTTTCTGTGGGGTGTTATGGCTGTGGTTTGCCCTATGGCCTGGGCCTTTTTGTTGGGGGGGGCGGCTTTAATGTTCTGCCTGAATATGTTTGATGTCCAGTTCAGGGCAAAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGGGGGGGAAATCATGAAGGCGGGCGGCTTGTCAAACCGCTTGGTTGAc                                                                                                                                                                                              <  2:1537822/227‑1 (MQ=255)
                                         gCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGGAAATATGAATGCGGGCGGGCTGTCAAAACGCACTGTTGAccc                                                                                                                                                                                            >  1:1442205/1‑238 (MQ=255)
                                         gCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGGAAACATGAAATCGGGGGGCTTGTTAAAACGCCatgttta                                                                                                                                                                                                >  1:1150385/1‑234 (MQ=255)
                                             tgCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGGAATCATGACTGCTGGGGGCTGTGTCAAAAGCCttgtttgaccgcc                                                                                                                                                                                          >  2:373854/1‑230 (MQ=255)
                                             tgCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGGAAACATTAATGAGGGCGGGCGGTTCAAACGCAttgtttgcgctgcc                                                                                                                                                                                         >  2:27500/1‑230 (MQ=255)
                                             tgCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGGAAACATGAATGGGGGGGGCTTGTCAAAAACCATTTTTTACCTGCCg                                                                                                                                                                                        >  1:1272209/1‑239 (MQ=255)
                                             tgCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGAAAACATGAATGGGGGGGGGCTGTCTACACACCTTGTTGGCCGGcac                                                                                                                                                                                        >  2:1358660/1‑237 (MQ=255)
                                             tgCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCCATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGGAATTATGGCTGGCGGGGGCGTCTCCAAACCCAttgcgtt                                                                                                                                                                                                >  2:1699544/1‑228 (MQ=255)
                                                        cTGGGTTGTTTTTGGGGTGGTTTTGCTATGGTTGTGCTTATTGCCGGGGCAGTTTTGTTGGGCGGGGCGGCGTAAATGTTCTGGCGGGACTTTTTTGATGTCCAGATCATGGCACAAACGCTGGTGACCGGGGCGAATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGGGGGTGAAATCATGAAGGGGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACtggt                                                                                                                                                                            <  2:1812731/239‑1 (MQ=255)
                                                         tGGGTTGTCTGTGGGGTGGTTTTGTTATCGGTTGGCTTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGAATTGTTTGATGTCCAGATCATGGCACAAGCGCTGGTGAAGGGGGCCGATGGCTTCTCCCGGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCAGGAATGCGGGGGGTTTTTAAAACCGCATTGTTGACCTGCCGATGAAACtggtg                                                                                                                                                                           <  2:1218543/239‑1 (MQ=255)
                                                         tGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGGAATCATGAATGCGGGCGGCTTTTTAAAACGCATTGTTTACCCGGCGAGTAAACtggtg                                                                                                                                                                           >  1:236177/1‑239 (MQ=255)
                                                                       ggtggtATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGGAATCATGTAGACGGGCGGGCTGTCAAAAACGAATTGTGATGTGGCTGTGAAAACGGTGGGGGATAAAACaag                                                                                                                                                             >  2:913842/1‑237 (MQ=255)
                                                                          ggtATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGGAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCgg                                                                                                                                                          >  1:1304870/1‑239 (MQ=255)
                                                                          ggtATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGGAATCATGAATGCGGGCGGCCTGTCCAAAAGCATTGTTGACCTGCCGGTGAAACTGGTGGGGGGTAAAACAGGCgg                                                                                                                                                          >  1:1425572/1‑239 (MQ=255)
                                                                          ggtATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAAACATGAATGCGGGCGGCCTGTCAAAAAGCCTTGTTGACCTGCCGATGAAACTGGGGGGGCATAAACCAGGCgg                                                                                                                                                          >  1:2023029/1‑239 (MQ=255)
                                                                          ggtATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGGAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTTGACTGCCCATGAAACTGGTGGGGGATAAACCAGGCgg                                                                                                                                                          >  1:1486618/1‑239 (MQ=255)
                                                                          ggtATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGGAAACCTGGATGGGGGGGGGGTTGTCAAACTCCATTTTGGTCTGCCGGTGGAACATGTGGGGGGTATAACAGGCgc                                                                                                                                                          >  2:1159249/1‑238 (MQ=255)
                                                                          ggtATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGCGGAATCCTGAATGCGGGCGGCTTGTCAAAAAGCCTTGTTGTCCTGCCGATGAAAATGGTGGGGCATAAACCAGGCgg                                                                                                                                                          >  1:247038/1‑239 (MQ=255)
                                                                            tATTGCTATCGTTTGTCCTATTGCCGGGGCACTGTTGTTGGGGGGGGGGGCGTTAATGTTCTGGCTGGAAATGTTTGATGTCCAGATCAGGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTTTTTTTGCTGGCGGGTAAAATCATGAATGCGGGCGGTTTGTAAAACCGCTTGTTGGCCCTGCCGATGAAACTGGTGGGGCATAAACCAGGCggg                                                                                                                                                         <  2:189801/238‑1 (MQ=255)
                                                                              ttGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGGAAACATTAATTGGGGCGGCCTGTCGAAAACCATTTTTTGCCCTGCCATGAAACATGGGGGGGGAAAAAcagggggg                                                                                                                                                          >  2:1237199/1‑229 (MQ=255)
                                                                                   ctggttttGCTTTTTGCCGGGGCACTGTTGTTGTGCGGGGCGGCTTTAATGTTCGGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGTGATTCCGTTTTTTGTGTTGGCGGGTGAAATCATGAATGCGGGGGGTTTGTCAAAACGCTTTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTa                                                                                                                                                <  2:1150385/234‑1 (MQ=255)
                                                                                   aTCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGGAATCATGAATGCGGGGGGCTTGTTCAAAAGCAGTTTTGACCTGCCGATGAAACTGGGGGGGCATAAACCAAGCggggcggggg                                                                                                                                                  >  1:444237/1‑231 (MQ=255)
                                                                                   aTCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGGGGGTTGTTAAAAAGCATTTTTGACCTTCCGATGAAAATGGGGGGGGATAAACCAGGCgggccggggg                                                                                                                                                  >  1:189801/1‑236 (MQ=255)
                                                                                                   gCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGc                                                                                                                                 <  1:270864/239‑1 (MQ=255)
                                                                                                    ccTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACg                                                                                                                                                                                                         <  1:1220860/169‑1 (MQ=255)
                                                                                                     cTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACg                                                                                                                                                                                                         <  1:1984062/168‑1 (MQ=255)
                                                                                                       gggTACTGTGGTGGTGGGGGGCGGGGTTAATGTTCGGGGTGGACAGGTTTGATGTCAGGATCGGGGCACAAAGGCGGGTGAGGGGAGCAAATAGCTTCTCCCTGCTGGCACTTCCTTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGGGGTTTGACAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCt                                                                                                                             <  2:247038/239‑1 (MQ=255)
                                                                                                              gttgttgtGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGGAATTATGAATGCGGGCGGCCTGTCAAAACGCAttgttggc                                                                                                                                                                                              >  2:1220860/1‑167 (MQ=255)
                                                                                                              gttgttgtGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGGAATCATGAATGCGGGCGGCCTGTTAAAACGCCTTGTTTACCTGCCGATGAAACTGGGGGGGCCTAAACCAGGCGGGCCGGGGTGACGCCGCGTGCTGGCAGcctaggata                                                                                                                      >  2:699374/1‑231 (MQ=255)
                                                                                                               ttgttgtGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGAAATCATGAATGCGGGGGGCTTGTTATAAACCCTTGTTGaa                                                                                                                                                                                              >  2:1984062/1‑167 (MQ=255)
                                                                                                                     tgCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTTTTTGTGCTGGGGGGTGAAATCATGAAGGCGGGCGGCTTGTCAAAACGCATTGTTGCCCTGCCGATGAAACTGGTGGGGCATAAACCAGGCggg                                                                                                                                                         <  1:585637/198‑1 (MQ=255)
                                                                                                                       cGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCAt                                                                                                                                                                                                      >  1:2191071/1‑153 (MQ=255)
                                                                                                                       cGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCAt                                                                                                                                                                                                      <  2:2191071/153‑1 (MQ=255)
                                                                                                                                     gtttCTGGCTGGTCATGTTTGTGTTCAGTCTCTGGCCACAACCGCTGGTAACCGTGGCCGATAGCTTCTCCCTGCTGGCTATTCCGTTTTTTGTGCTGGGGGGGGAAACCATAAATGCGGGCGGTTTGTCAAAACGCATTGTGGCCCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTTCGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGcc                                                                                               <  2:485629/237‑1 (MQ=255)
                                                                                                                                            gctggACATGTTTGATGTCCAGATCATGGCACAAACGCGGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGa                                                                                                                                                                                    <  1:187489/149‑1 (MQ=255)
                                                                                                                                                             tCCAAATCATGGCACAAACGCTGTGGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAACCGCATTGTTGCCCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTc                                                                                                         <  1:576918/205‑1 (MQ=255)
                                                                                                                                                                        ggCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGTTTCCGTTTTTTGTGCTGGCGGGTAAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGCCCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATAccgccg                                                                               <  1:1497330/221‑1 (MQ=255)
                                                                                                                                                                           ccaagctggtgTGACCGGCCCCGATCGTTTCTCCCTGGTCGCACTTCTTTTTTTTGTGTGGGGGGTTAAACTCATGAAGGCGGACGCTTTTTAAACACGCATTGTTGCCCTGCCGATAAAACTGGTGGGGAATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGGGCCTTTCTGGTTCTGCCGTTGCGGATGCCGCCGCCGTTGCCGCTCTACTGGTGc                                                          <  2:1441573/229‑1 (MQ=255)
                                                                                                                                                                                          ggggcacGATAGTTCCCCCCGGGTGGCATTCCGTTTTTTTGTGGTGGGGGATGAAGTAAGGAAGGGGGGGGTTTTGTAAAAACGCTTTGTTGCCCTGCCGATGAACCTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAg                                             <  2:1442205/231‑1 (MQ=255)
                                                                                                                                                                                                      cTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCaaa                                                                                                                                                                                                            <  1:314071/69‑1 (MQ=255)
                                                                                                                                                                                                       ttCCCCGTGGTGGCTATCCCGTTTTTTGGGTGGGGGGTTAAAATCAGGATGGCGGGGGGTTTGTAAAAACGCATTGTTGCCCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGt                              <  1:1159249/239‑1 (MQ=255)
                                                                                                                                                                                                              tgctgGCGATTCCGTTCTTTGTGCTGGGGGGTAAAATCATGAAGGCGGGGGGTTTGTAAAAACGCATTTTTGACCTGCCGATGAAACtggt                                                                                                                                                                            <  2:2105686/91‑1 (MQ=255)
                                                                                                                                                                                                                ctgGCGATTCCGTTCTTTGTGCTGGCGGGGGGAAACATGAATGGCGGCGGCTTGTTAAAAACCAttgtt                                                                                                                                                                                                  >  2:314071/1‑69 (MQ=255)
                                                                                                                                                                                                                    cGATTCCGTTCTTTGTGCTGGCGGGGGAAATCATGa                                                                                                                                                                                                                                  >  1:1595451/1‑36 (MQ=255)
                                                                                                                                                                                                                       ttCCGTTTTTGGTGCGGGGGGGTAAAATCATGAATGCGGGGGGCTTGTCAAACCGCTTTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGt                                                                           <  1:609937/178‑1 (MQ=255)
                                                                                                                                                                                                                        tCCGTTCTTTGTGCTGGCGGGGGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGGGGGGCATaaa                                                                                                                                                                  >  1:2105686/1‑91 (MQ=255)
                                                                                                                                                                                                                           gTTTTTTGTGCGGGGGGGTAAAATCATGAGTGCGGGGGGCTTGTCAACACGCATTGTGGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGa          <  1:1358660/239‑1 (MQ=255)
                                                                                                                                                                                                                                     gCTGGGGGGAAAAATCATGAAGGGGGGGGGTTTGTCAAAACGCTTTGTTGCCCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATCCCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCctct  <  1:1705772/238‑1 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                         
AGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCCTCT  >  NC_000913/3744766‑3745229

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