Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,745,0190CG35.7% 55.6 / 16.5 56G72G (GGC→GGG) yiaN2,3‑diketo‑L‑gulonate TRAP transporter large permease protein
Reads supporting (aligned to +/- strand):  ref base C (18/18);  new base G (9/11);  total (27/29)
Fisher's exact test for biased strand distribution p-value = 7.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.90e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTG‑CCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAG‑ATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCG‑TTCTTTGTGCTGGCGGGTGAAA‑TCATGAAT‑GCGGGCGGCT‑TGTC‑AAAACGCATTGTT‑GACCTGCCGATGAAAC‑TGGTGGGGCATAAACCAGGCGGGCTGGGCT‑ACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTATCGCGCCAATTATT  >  NC_000913/3744784‑3745253
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                               
ccACCGCAAGGAGATTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTG‑CCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAG‑ATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCG‑TTCTTTGTGCTGGCGGGTGGAA‑TCATGAAT‑TCCgggggc                                                                                                                                                                                                                                              >  1:787617/1‑239 (MQ=255)
    cGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTG‑CCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAG‑ATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCG‑TTCTTTGTGCTGGCGGGGGGAA‑ACATGAAT‑gggggggg                                                                                                                                                                                                                                               >  1:1031667/1‑234 (MQ=255)
       aaGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTG‑CCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAG‑ATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCG‑TTCTTTGTGCTGGCGGGTGAAA‑TCATGAAT‑GCGGGCGGCT‑TGTC‑a                                                                                                                                                                                                                                      <  1:859539/238‑1 (MQ=255)
                   cATGGCGGTGCTGATTTTTCTGGGCTGTTGGTTGGGTGGTATTGCTATCGGTTTG‑CCTATTGCCGGGCAACTGTTGTTGGGCGGGGCGGCGTAAATGTTCTGACGGGACATGTTTGATGTCCAG‑ATCATGGCACAAACGCTGGTGACGGGAGCGGATAGCTTCTCCCTGCTGGCGATTCCG‑TTCTTTGTGCTGGCGGGTGAAA‑TCATAAAG‑GCGGGGGGCT‑TGTA‑AAAACGCAttgtt‑g                                                                                                                                                                                                                        <  2:1988278/239‑1 (MQ=255)
                     tGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTG‑CCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAG‑ATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCG‑TTCTTTGTGCTGGCGGGTGGAA‑TCATTAAA‑GCGGGCGGGT‑TTTC‑CAAAAGCCTTGTT‑Ga                                                                                                                                                                                                                       >  1:571517/1‑238 (MQ=255)
                      gtctttgctattttgtCGGGGTTTTCTGTGGGGTGTTATGGCTGT‑GGTTTGCCCTATGGCCTGGGCCTTTTTGTTGGGGGGGGCGGCTTTAATGTTCTGCCTGAATATGTTTGATGTCCAG‑TTCAGGGCAAAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCG‑TTCTTTGTGCTGGGGGGGGAAA‑TCATGAAG‑GCGGGCGGCT‑TGTC‑AAACCGCTTGGTT‑GAc                                                                                                                                                                                                                      <  2:1537822/227‑1 (MQ=255)
                       gCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTG‑CCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAG‑ATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCG‑TTCTTTGTGCTGGCGGGTGGAA‑ATATGAAT‑GCGGGCGGGC‑TGTC‑AAAACGCACTGTT‑GAccc                                                                                                                                                                                                                    >  1:1442205/1‑238 (MQ=255)
                       gCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTG‑CCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAG‑ATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCG‑TTCTTTGTGCTGGCGGGGGGAA‑ACATGAAA‑TCGGGGGGCT‑TGTT‑AAAACGCCatgttta                                                                                                                                                                                                                        >  1:1150385/1‑234 (MQ=255)
                           tgCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTG‑CCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAG‑ATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCG‑TTCTTTGTGCTGGCGGGTGGAA‑TCATGACT‑GCTGGGGGCTGTGTC‑AAAA‑GCCttgtttgaccgcc                                                                                                                                                                                                                  >  2:373854/1‑230 (MQ=255)
                           tgCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTG‑CCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAG‑ATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCG‑TTCTTTGTGCTGGCGGGTGGAA‑ACATTAAT‑GAGGGCGGGC‑GGTT‑CAAACGCAttgtttgcgctgcc                                                                                                                                                                                                                 >  2:27500/1‑230 (MQ=255)
                           tgCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTG‑CCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAG‑ATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCG‑TTCTTTGTGCTGGCGGGGGGAA‑ACATGAAT‑GGGGGGGGCT‑TGTC‑AAAAACCATTTTT‑TACCTGCCg                                                                                                                                                                                                                >  1:1272209/1‑239 (MQ=255)
                           tgCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTG‑CCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAG‑ATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCG‑TTCTTTGTGCTGGCGGGGGAAA‑ACATGAAT‑GGGGGGGGGC‑TGTC‑TACACACCTTGTT‑GGCCGGcac                                                                                                                                                                                                                >  2:1358660/1‑237 (MQ=255)
                           tgCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCCATCGGTTTG‑CCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAG‑ATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCG‑TTCTTTGTGCTGGCGGGGGGAA‑TTATGGCT‑GGCGGGGGCG‑TCTC‑CAAACCCAttgcgtt                                                                                                                                                                                                                        >  2:1699544/1‑228 (MQ=255)
                                      cTGGGTTGTTTTTGGGGTGGTTTTGCTATGGTTGTG‑CTTATTGCCGGGGCAGTTTTGTTGGGCGGGGCGGCGTAAATGTTCTGGCGGGACTTTTTTGATGTCCAG‑ATCATGGCACAAACGCTGGTGACCGGGGCGAATAGCTTCTCCCTGCTGGCGATTCCG‑TTCTTTGTGCTGGGGGGTGAAA‑TCATGAAG‑GGGGGCGGCT‑TGTC‑AAAACGCATTGTT‑GACCTGCCGATGAAAC‑tggt                                                                                                                                                                                                    <  2:1812731/239‑1 (MQ=255)
                                       tGGGTTGTCTGTGGGGTGGTTTTGTTATCGGTTGG‑CTTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGAATTGTTTGATGTCCAG‑ATCATGGCACAAGCGCTGGTGAAGGGGGCCGATGGCTTCTCCCGGCTGGCGATTCCG‑TTCTTTGTGCTGGCGGGTGAAA‑TCAGGAAT‑GCGGGGGGTT‑TTTA‑AAACCGCATTGTT‑GACCTGCCGATGAAAC‑tggtg                                                                                                                                                                                                   <  2:1218543/239‑1 (MQ=255)
                                       tGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTG‑CCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAG‑ATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCG‑TTCTTTGTGCTGGCGGGTGGAA‑TCATGAAT‑GCGGGCGGCT‑TTTT‑AAAACGCATTGTT‑TACCCGGCGAGTAAAC‑tggtg                                                                                                                                                                                                   >  1:236177/1‑239 (MQ=255)
                                                     ggtggtATTGCTATCGGTTTG‑CCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAG‑ATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCG‑TTCTTTGTGCTGGCGGGGGGAA‑TCATGTAG‑ACGGGCGGGC‑TGTCAAAAACGAATTG‑T‑GATGTGGCTGTGAAAA‑CGGTGGGGGATAAAACaag                                                                                                                                                                                     >  2:913842/1‑237 (MQ=255)
                                                        ggtATTGCTATCGGTTTG‑CCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAG‑ATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCG‑TTCTTTGTGCTGGCGGGTGGAA‑TCATGAAT‑GCGGGCGGCT‑TGTC‑AAAACGCATTGTT‑GACCTGCCGATGAAAC‑TGGTGGGGCATAAACCAGGCgg                                                                                                                                                                                  >  1:1304870/1‑239 (MQ=255)
                                                        ggtATTGCTATCGGTTTG‑CCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAG‑ATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCG‑TTCTTTGTGCTGGCGGGTGGAA‑TCATGAAT‑GCGGGCGGCC‑TGTC‑CAAAAGCATTGTT‑GACCTGCCGGTGAAAC‑TGGTGGGGGGTAAAACAGGCgg                                                                                                                                                                                  >  1:1425572/1‑239 (MQ=255)
                                                        ggtATTGCTATCGGTTTG‑CCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAG‑ATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCG‑TTCTTTGTGCTGGCGGGTGAAA‑ACATGAAT‑GCGGGCGGCC‑TGTC‑AAAAAGCCTTGTT‑GACCTGCCGATGAAAC‑TGGGGGGGCATAAACCAGGCgg                                                                                                                                                                                  >  1:2023029/1‑239 (MQ=255)
                                                        ggtATTGCTATCGGTTTG‑CCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAG‑ATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCG‑TTCTTTGTGCTGGCGGGGGGAA‑TCATGAAT‑GCGGGCGGCT‑TGTC‑AAAACGCATTGTT‑TGACTGCCCATGAAAC‑TGGTGGGGGATAAACCAGGCgg                                                                                                                                                                                  >  1:1486618/1‑239 (MQ=255)
                                                        ggtATTGCTATCGGTTTG‑CCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAG‑ATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCG‑TTCTTTGTGCTGGCGGGGGGAA‑ACCTGGATGGGGGGGGGGT‑TGTC‑AAACTCCATT‑TT‑GGTCTGCCGGTGGAAC‑ATGTGGGGGGTATAACAGGCgc                                                                                                                                                                                  >  2:1159249/1‑238 (MQ=255)
                                                        ggtATTGCTATCGGTTTG‑CCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAG‑ATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCG‑TTCTTTGTGCTGGCGGGCGGAA‑TCCTGAAT‑GCGGGCGGCT‑TGTC‑AAAAAGCCTTGTT‑GTCCTGCCGATGAAAA‑TGGTGGGGCATAAACCAGGCgg                                                                                                                                                                                  >  1:247038/1‑239 (MQ=255)
                                                          tATTGCTATCGTTTGT‑CCTATTGCCGGGGCACTGTTGTTGGGGGGGGGGGCGTTAATGTTCTGGCTGGAAATGTTTGATGTCCAG‑ATCAGGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCG‑TTTTTTTTGCTGGCGGGTAAAA‑TCATGAAT‑GCGGGCGGTT‑TGTA‑AAACCGC‑TTGTTGGCCCTGCCGATGAAAC‑TGGTGGGGCATAAACCAGGCggg                                                                                                                                                                                 <  2:189801/238‑1 (MQ=255)
                                                            ttGCTATCGGTTTG‑CCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAG‑ATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCG‑TTCTTTGTGCTGGCGGGGGGAA‑ACATTAAT‑TGGGGCGGCC‑TGTCGAAAACCATTTTTT‑GCCCTGCC‑ATGAAACATGGGGGGGGAAAAAcagggggg                                                                                                                                                                                  >  2:1237199/1‑229 (MQ=255)
                                                                 ctggttttG‑CTTTTTGCCGGGGCACTGTTGTTGTGCGGGGCGGCTTTAATGTTCGGGCTGGACATGTTTGATGTCCAG‑ATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGTGATTCCG‑TTTTTTGTGTTGGCGGGTGAAA‑TCATGAAT‑GCGGGGGGTT‑TGTC‑AAAACGCTTTGTT‑GACCTGCCGATGAAAC‑TGGTGGGGCATAAACCAGGCGGGCTGGGCT‑a                                                                                                                                                                        <  2:1150385/234‑1 (MQ=255)
                                                                 aTCGGTTTG‑CCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAG‑ATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCG‑TTCTTTGTGCTGGCGGGTGGAA‑TCATGAAT‑GCGGGGGGCT‑TGTT‑CAAAAGCAGTTTT‑GACCTGCCGATGAAAC‑TGGGGGGGCATAAACCAAGCggggcggggg                                                                                                                                                                          >  1:444237/1‑231 (MQ=255)
                                                                 aTCGGTTTG‑CCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAG‑ATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCG‑TTCTTTGTGCTGGCGGGTGAAA‑TCATGAAT‑GCGGGGGGGT‑TGTT‑AAAAAGCATTTTT‑GACCTTCCGATGAAAA‑TGGGGGGGGATAAACCAGGCgggccggggg                                                                                                                                                                          >  1:189801/1‑236 (MQ=255)
                                                                                 gCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAG‑ATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCG‑TTCTTTGTGCTGGCGGGTGAAA‑TCATGAAT‑GCGGGCGGCT‑TGTC‑AAAACGCATTGTT‑GACCTGCCGATGAAAC‑TGGTGGGGCATAAACCAGGCGGGCTGGGCT‑ACGTCGGCGTGCTGGc                                                                                                                                                         <  1:270864/239‑1 (MQ=255)
                                                                                  ccTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAG‑ATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCG‑TTCTTTGTGCTGGCGGGTGAAA‑TCATGAAT‑GCGGGCGGCT‑TGTC‑AAAACg                                                                                                                                                                                                                                 <  1:1220860/169‑1 (MQ=255)
                                                                                   cTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAG‑ATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCG‑TTCTTTGTGCTGGCGGGTGAAA‑TCATGAAT‑GCGGGCGGCT‑TGTC‑AAAACg                                                                                                                                                                                                                                 <  1:1984062/168‑1 (MQ=255)
                                                                                     gggTACTGTGGTGGTGGGGGGCGGGGTTAATGTTCGGGGTGGACAGGTTTGATGTCAGG‑ATCGGGGCACAAAGGCGGGTGAGGGGAGCAAATAGCTTCTCCCTGCTGGCACTTCCT‑TTCTTTGTGCTGGCGGGTGAAA‑TCATGAAT‑GCGGGGGGTT‑TGAC‑AAAACGCATTGTT‑GACCTGCCGATGAAAC‑TGGTGGGGCATAAACCAGGCGGGCTGGGCT‑ACGTCGGCGTGCTGGCAGCt                                                                                                                                                     <  2:247038/239‑1 (MQ=255)
                                                                                            gttgttgtGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAG‑ATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCG‑TTCTTTGTGCTGGCGGGGGGAA‑TTATGAAT‑GCGGGCGGCC‑TGTC‑AAAACGCAttgtt‑ggc                                                                                                                                                                                                                      >  2:1220860/1‑167 (MQ=255)
                                                                                            gttgttgtGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAG‑ATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCG‑TTCTTTGTGCTGGCGGGGGGAA‑TCATGAAT‑GCGGGCGGCC‑TGTT‑AAAACGCCTTGTT‑TACCTGCCGATGAAAC‑TGGGGGGGCCTAAACCAGGCGGGCCGGGGTGACG‑CCGCGTGCTGGCAGcctaggata                                                                                                                                              >  2:699374/1‑231 (MQ=255)
                                                                                             ttgttgtGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAG‑ATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCG‑TTCTTTGTGCTGGCGGGGGAAA‑TCATGAAT‑GCGGGGGGCT‑TGTT‑ATAAACCCTTGTT‑Gaa                                                                                                                                                                                                                      >  2:1984062/1‑167 (MQ=255)
                                                                                                   tgCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAG‑ATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCG‑TTTTTTGTGCTGGGGGGTGAAA‑TCATGAAG‑GCGGGCGGCT‑TGTC‑AAAACGCATTGTT‑GCCCTGCCGATGAAAC‑TGGTGGGGCATAAACCAGGCggg                                                                                                                                                                                 <  1:585637/198‑1 (MQ=255)
                                                                                                     cGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAG‑ATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCG‑TTCTTTGTGCTGGCGGGTGAAA‑TCATGAAT‑GCGGGCGGCT‑TGTC‑AAAACGCAt                                                                                                                                                                                                                              <  2:2191071/153‑1 (MQ=255)
                                                                                                     cGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAG‑ATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCG‑TTCTTTGTGCTGGCGGGTGAAA‑TCATGAAT‑GCGGGCGGCT‑TGTC‑AAAACGCAt                                                                                                                                                                                                                              >  1:2191071/1‑153 (MQ=255)
                                                                                                                   gtttCTGGCTGGTCATGTTTG‑TGTTCAGTCTCTGGCCACAACCGCTGGTAACCGTGGCCGATAGCTTCTCCCTGCTGGCTATTCCG‑TTTTTTGTGCTGGGGGGGGAAA‑CCATAAAT‑GCGGGCGGTT‑TGTC‑AAAACGCATTGTG‑GCCCTGCCGATGAAAC‑TGGTGGGGCATAAACCAGGCGGGCTGGGCT‑TCGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGcc                                                                                                                       <  2:485629/237‑1 (MQ=255)
                                                                                                                          gctggACATGTTTGATGTCCAG‑ATCATGGCACAAACGCGGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCG‑TTCTTTGTGCTGGCGGGTGAAA‑TCATGAAT‑GCGGGCGGCT‑TGTC‑AAAACGCATTGTT‑GACCTGCCGATGa                                                                                                                                                                                                            <  1:187489/149‑1 (MQ=255)
                                                                                                                                           tCCAA‑ATCATGGCACAAACGCTGTGGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCG‑TTCTTTGTGCTGGCGGGTGAAA‑TCATGAAT‑GCGGGCGGCT‑TGTC‑AAACCGCATTGTT‑GCCCTGCCGATGAAAC‑TGGTGGGGCATAAACCAGGCGGGCTGGGCT‑ACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTc                                                                                                                                 <  1:576918/205‑1 (MQ=255)
                                                                                                                                                      ggCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGTTTCCG‑TTTTTTGTGCTGGCGGGTAAAA‑TCATGAAT‑GCGGGCGGCT‑TGTC‑AAAACGCATTGTT‑GCCCTGCCGATGAAAC‑TGGTGGGGCATAAACCAGGCGGGCTGGGCT‑ACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATAccgccg                                                                                                       <  1:1497330/221‑1 (MQ=255)
                                                                                                                                                         ccaagctggtgTGACCGGCCCCGATCGTTTCTCCCTGGTCGCACTTCTT‑TTTTTTGTG‑TGGGGGGTTAAACTCATGAAG‑GCGGACGCTT‑TTTA‑AACACGCATTGTT‑GCCCTGCCGATAAAAC‑TGGTGGGGAATAAACCAGGCGGGCTGGGCT‑ACGTCGGCGTGCTGGCAGCTATGATTATGGCGGGCCTTTCTGGTTCTGCCGTTGCGGATGCCGCCGCCGTTGCCGCTCTACTGGTGc                                                                                  <  2:1441573/229‑1 (MQ=255)
                                                                                                                                                                        ggggcacGATAGTTCCCCCCGGGTGGC‑ATTCCGTTTTTTTGTGGTGGGGGATGAAG‑TAAGGAAG‑GGGGGGGTTT‑TGTA‑AAAACGCTTTGTT‑GCCCTGCCGATGAACC‑TGGTGGGGCATAAACCAGGCGGGCTGGGCT‑ACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAg                                                                     <  2:1442205/231‑1 (MQ=255)
                                                                                                                                                                                    cTTCTCCCTGCTGGCGATTCCG‑TTCTTTGTGCTGGCGGGTGAAA‑TCATGAAT‑GCGGGCGGCT‑TGTC‑aaa                                                                                                                                                                                                                                    <  1:314071/69‑1 (MQ=255)
                                                                                                                                                                                     ttCCCCGTGGTGGCTATCCCG‑TTTTTTGGGTGGGGGGTTAAAA‑TCAGGATG‑GCGGGGGGTT‑TGTA‑AAAACGCATTGTT‑GCCCTGCCGATGAAAC‑TGGTGGGGCATAAACCAGGCGGGCTGGGCT‑ACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGt                                                      <  1:1159249/239‑1 (MQ=255)
                                                                                                                                                                                            tgctgGCGATTCCG‑TTCTTTGTGCTGGGGGGTAAAA‑TCATGAAG‑GCGGGGGGTT‑TGTA‑AAAACGCATTTTT‑GACCTGCCGATGAAAC‑tggt                                                                                                                                                                                                    <  2:2105686/91‑1 (MQ=255)
                                                                                                                                                                                              ctgGCGATTCCG‑TTCTTTGTGCTGGCGGGGGGAA‑ACATGAAT‑GGCGGCGGCT‑TGTT‑AAAAACCAttgtt                                                                                                                                                                                                                          >  2:314071/1‑69 (MQ=255)
                                                                                                                                                                                                     ttCCG‑TTTTTGGTGCGGGGGGGTAAAA‑TCATGAAT‑GCGGGGGGCT‑TGTC‑AAACCGCTTTGTT‑GACCTGCCGATGAAAC‑TGGTGGGGCATAAACCAGGCGGGCTGGGCT‑ACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGt                                                                                                   <  1:609937/178‑1 (MQ=255)
                                                                                                                                                                                                      tCCG‑TTCTTTGTGCTGGCGGGGGAAA‑TCATGAAT‑GCGGGCGGCT‑TGTC‑AAAACGCATTGTT‑GACCTGCCGATGAAAC‑TGGGGGGGCATaaa                                                                                                                                                                                          >  1:2105686/1‑91 (MQ=255)
                                                                                                                                                                                                         g‑TTTTTTGTGCGGGGGGGTAAAA‑TCATGAGT‑GCGGGGGGCT‑TGTC‑AACACGCATTGTG‑GACCTGCCGATGAAAC‑TGGTGGGGCATAAACCAGGCGGGCTGGGCT‑ACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGa                                  <  1:1358660/239‑1 (MQ=255)
                                                                                                                                                                                                                   gCTGGGGGGAAAAA‑TCATGAAG‑GGGGGGGGTT‑TGTC‑AAAACGCTTTGTT‑GCCCTGCCGATGAAAC‑TGGTGGGGCATAAACCAGGCGGGCTGGGCT‑ACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATCCCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCctct                          <  1:1705772/238‑1 (MQ=255)
                                                                                                                                                                                                                                    aTGAAG‑GCGGGGGGCT‑TGTC‑AAAACGCTTTTTT‑GCCCTGCCGATGAAAC‑TGGTGGGGCATAAACCAGGCGGGCTGGGCT‑ACGTCGGCGTGCTGGCAGCTATGATTATGGCGAg                                                                                                                                       <  2:12385/113‑1 (MQ=255)
                                                                                                                                                                                                                                      gAAT‑GCGGGCGGCT‑TGTC‑AAAACGCATTGTT‑GACCTGCCGATGAAAC‑TGGTGGGGCATAAACCAGGCGGGCTGGGCT‑ACGTCGGCGTGCt                                                                                                                                                            >  2:1284013/1‑90 (MQ=255)
                                                                                                                                                                                                                                      gAAT‑GCGGGCGGCT‑TGTC‑AAAACGCATTGTT‑GACCTGCCGATGAAAC‑TGGTGGGGCATAAACCAGGCGGGCTGGGCT‑ACGTCGGCGTGCt                                                                                                                                                            <  1:1284013/90‑1 (MQ=255)
                                                                                                                                                                                                                                       aaT‑GCGGGCGGCT‑TGTC‑AAAACGCATTGTT‑GACCTGCCGATGAAAC‑TGGTGGGGCATAAACCAGGCGGGCTGGGCT‑ACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTg                                                                                      >  1:1287004/1‑159 (MQ=255)
                                                                                                                                                                                                                                       aaT‑GCGGCCGGCT‑TGTA‑AAACCGCATTGTT‑GCCCTGCCGATGAAAC‑TGGTGGGGCATAAACCAGGCGGGCTGGGCT‑ACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTATCGCGcc         <  1:625925/236‑1 (MQ=255)
                                                                                                                                                                                                                                         t‑GCGGGCGGCT‑TGTA‑AAAACGCTTTGTG‑GCCCTGCCGATGAAAC‑TGGTGGGGCATAAACCAGGCGGGCTGGGCT‑ACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTATCGCGCCAatta    <  1:1956349/239‑1 (MQ=255)
                                                                                                                                                                                                                                            gggggggTT‑TGTC‑AAACCTCATTGTT‑GCCCTGCCGATGAAAC‑TGGTGGGGCATAAACCAGGCGGGCTGGGCT‑ACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTATCGCGcc         <  1:373854/231‑1 (MQ=255)
                                                                                                                                                                                                                                            ggggcggCT‑TTTA‑AAACCGCATTGTT‑GACCTGCCGATGAAAC‑TGGTGGGGCATAAACCAGGCGGGCTGGGCT‑ACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTATCGCGCCAattat   <  1:1651926/237‑1 (MQ=255)
                                                                                                                                                                                                                                             gggcggTT‑TGTA‑AAACCGCTTTGTT‑GCCCTGCCGATGAAAC‑TGGTGGGGCATAAACCAGGCGGGCTGGGCT‑ACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTATCGCGCCAattatt  <  1:27500/238‑1 (MQ=255)
                                                                                                                                                                                                                                               gcggcT‑TGTC‑AAAACGCATTGTT‑GACCTGCCGATGAAAC‑TGGTGGGGCATAAACCAGGCGGGCTGGGCT‑ACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGt                                                                                                                             >  1:12385/1‑113 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                               
CCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTG‑CCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAG‑ATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCG‑TTCTTTGTGCTGGCGGGTGAAA‑TCATGAAT‑GCGGGCGGCT‑TGTC‑AAAACGCATTGTT‑GACCTGCCGATGAAAC‑TGGTGGGGCATAAACCAGGCGGGCTGGGCT‑ACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTATCGCGCCAATTATT  >  NC_000913/3744784‑3745253

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