Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,745,0190CG43.6% 29.7 / 15.9 39G72G (GGC→GGG) yiaN2,3‑diketo‑L‑gulonate TRAP transporter large permease protein
Reads supporting (aligned to +/- strand):  ref base C (12/10);  new base G (8/9);  total (20/19)
Fisher's exact test for biased strand distribution p-value = 7.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATG‑TTCTGGC‑TGGACATGTTTGATGTCCAG‑ATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAA‑TCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGC‑TGGGCTACGTCGGCGTGCTGGCAGCTA‑TGATTATGGCGAG‑CCTTTCTGGTTCTGCCG‑TTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGA‑TGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTATCGCGCCAATTATTCC  >  NC_000913/3744796‑3745255
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                 
gCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATG‑TTCTGGC‑TGGACATGTTTGATGTCCAG‑ATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAA‑TCATGAATGCGGGCGGCTTGTCAAAACg                                                                                                                                                                                                                                     <  1:1125288/238‑1 (MQ=255)
               tgCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATG‑TTCTGGC‑TGGACATGTTTGATGTCCAG‑ATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGGAA‑ACATTAATGTGGGGGGCTTTTCCAAAAGCCTTGTTTAACTGCCg                                                                                                                                                                                                                     >  1:466351/1‑239 (MQ=255)
                 cTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATG‑TTCTGGC‑TGGACATGTTTGATGTCCAG‑ATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGGAA‑TTATGAAATCGGGGGGGCTTTCAAAAAGCATTGTTTACCCGcc                                                                                                                                                                                                                      >  1:1513335/1‑236 (MQ=255)
                                tttttgtTGGGTGGTTTTGCTTTCGGTTTGCCTATTGCCGGGGCATTGTTGTTGTGGGGGGCGGCGTTAAGG‑TTCGGGC‑TGGACATTTTTGATGTCCAG‑ATCATGGCACAACCGCTGGTGAACGGGGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGGGGGTAAAA‑TAATGAAGGCGGGCGGTTTGTAAAAACGCATTGTTGCCCTGCCGATGAAACTGGTGGGGc                                                                                                                                                                                                     <  2:703244/234‑1 (MQ=255)
                                 tgtttgtgggtgggATTGCTAGCGTTTTGCCTATTGCCTGGGCGCTGTTGTTGTGCGGGGGGGCTGTAATG‑TTCTGGC‑TGGACATGTTTGATG‑CCAGAATCAAGCCACAACCGTGGGAGAGGGGAGCCATTAGCTTCCCCCGGCTGGCTATTCCGTTCTTTGTGCTGGAGGGTAAAA‑TAAAGAAGGGGGGGGGTTTGTAAAAACGCTTTGTTGCCCTGCCGATGAAACTGGTGGGGCa                                                                                                                                                                                                    <  2:471016/232‑1 (MQ=255)
                                            ggtATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATG‑TTCTGGC‑TGGACATGTTTGATGTCCAG‑ATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAA‑TCATGAATTCGGGCGGCTTGTCCAAACGCCTTGTTGACCTGCCGCTGAAACTGGTGGGGCATTAACCAAGCgg                                                                                                                                                                                        >  1:1446171/1‑239 (MQ=255)
                                            ggtATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATG‑TTCTGGC‑TGGACATGTTTGATGTCCAG‑ATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAA‑TCATGAAAGCGGGCGGCTTGTCAAAAAGCCTTTTTGACCTGCCGAGGAAACTGGTGGGGGCTAAAACAGGCgg                                                                                                                                                                                        >  1:1256752/1‑239 (MQ=255)
                                            ggtATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATG‑TTCTGGC‑TGGACATGTTTGATGTCCAG‑ATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGAAA‑ACATGAAATCGGGGGGGTTTTCAAAAAGCATTTTTTGCCTGCCGGTGGAACTGGGGGGGGATAAAACAAGCgg                                                                                                                                                                                        >  1:194245/1‑239 (MQ=255)
                                             gtATTGTTATGTGTTTCACTATTGCCTGGGCTCTGTTGTTGTGCGGGGCGGCTTTAATT‑TTCTGGC‑TGGACAGTTTTGATTGCAAG‑ATTATGCCACAAGCGCTGGTGAGAGGAGCAGATATCTTCTCCCTGCTGGATATTGCGTTTTTTGTGCTGGGGTGGGAAA‑CTATGCTCGGGGGCGGTTGCGCAAAACGCGTTTTGGACCTGCCGATGAATCTGGTGGGGCATAAACCAGGCggg                                                                                                                                                                                       <  2:1016262/239‑1 (MQ=255)
                                              tATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATG‑TTCTGGC‑TGGACATGTTTGATGTCCAG‑ATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGAAA‑ACATGAATGCGGGCGGCTTGTCAAAACGCATTGTTTGACTGCCGGGGGAACTGGGTGGGGATTAACCCggggggg                                                                                                                                                                                      >  1:717720/1‑238 (MQ=255)
                                              gtttGCTTTTGTTTTGCCTTTTCCGTGGGCACTGTTGTTGTGGGGGGCGGGTTTAATG‑TTCTGGC‑TGGACATGTTTGATGTCCAG‑ATCATGGCACAAACGCTGGTGACCGGAGCCGATAGCTTCTCCCTGCTGGCGTTTCCGTTCTTTGTGCTGGGGGGTGAAA‑TCATGAATGCGGGGGGTTTGTCAAACCGCTTTTTTGCCCTGCCGATGAAACTGGTGGGGCATAAACCAGGCgggc                                                                                                                                                                                      <  2:566770/237‑1 (MQ=255)
                                                ttGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATG‑TTCTGGC‑TGGACATGTTTGATGTCCAG‑ATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGAAA‑TCATGAATGCGGGGGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCgggg‑tg                                                                                                                                                                                   >  1:334331/1‑239 (MQ=255)
                                                             gCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATG‑TTCTGGC‑TGGACATGTTTGATGTCCAG‑ATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGGAA‑ACATGAATGCGGGGGGCTTGTCAAAAAGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGG‑TGGGCGACGTCGGCg                                                                                                                                                                      >  1:1023419/1‑239 (MQ=255)
                                                                       ggggCATTGTTGTTGGGGGGGGCGGCTTAAATG‑TTCTGGC‑TGGACATTTTTGATGTCCAG‑ATCATGCCACAAACGTTGGTAAACGGGGCCGATAGCTTCTCCCTGCGGGCGTTTCCTTTTTTTTTGTTGGCGGGTAAAA‑TCATGAATGCGGGCGGCTTGTAAAAACGCATTGTTGCCCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGC‑TGGGCTACGTCGGCGTGCTGGCAGc                                                                                                                                                            <  2:194245/238‑1 (MQ=255)
                                                                        gggCACTTTTGTTGGGGGGGGCGTCGTAAATG‑TT‑TGGCGTGTACATTGTTTATGCCCAG‑ATCATGGCACAACCGCTGGTGAACGGGGCCGATGGCTTCCCCCGGCGGGCGTTTCCGTTTTTTGTGCGGGGGGGTGAAA‑TCATGAAGGGGGGGGGTTTGTCAAACCGCATTGTGGCCCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGC‑TGGGCTACGTCGGCGTGCTGGCAGc                                                                                                                                                            <  2:178354/238‑1 (MQ=255)
                                                                           ccattgtgttttGGGGGGGGGGGTTTATGTTTCTGGT‑GGGACTTTTTTGATGCCCAG‑TTGTTGGCACAAACGTGGGGGAACGGAGCCGATAGTTTCTCCCTGGGGGCGATTCCGTTTTTTGTGGTGGGGGGTAAAA‑TCAGGAAGGCGGGGGGTTTGTCAAAACGAATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGC‑TGGGCTACGTCGGCGTGCTGGCAGCTa                                                                                                                                                          <  2:1513335/232‑1 (MQ=255)
                                                                              tgttgttgtGCGGGGCGGCGTTAATG‑TTCTGGC‑TGGACATGTTTGATGTCCAG‑ATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAA‑TCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGGGGGGCATAAACCAGGCGGGC‑TGGGGTACGTCCGCGTGCCGGCAGGTa‑ttatt                                                                                                                                                    >  1:1023704/1‑239 (MQ=255)
                                                                                ttgttgtGCGGGGCGGCGTTAATG‑TTCTGGC‑TGGACTTGTTTGATGCCCAG‑ATCATGGCCCAACCGCTGGTGAAGGGAGCCGATAGCTTCCCCCTGCTGGCGATTCCGTTTTTTGTGCTGGCGGGTAAAA‑TCATAAATGCGGGGGGTTTGTCAAAACGCATTGTTGCCCTGCCGATGAAACTGGTgggg                                                                                                                                                                                                      <  2:1154066/189‑1 (MQ=255)
                                                                                 tgttgtGCGGGGCGGCGTTAATG‑TTCTGGC‑TGGACATGTTTGATGTCCAG‑ATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAA‑TCATGAATGCGGGCGGCTTGTCAAAACg                                                                                                                                                                                                                                     <  1:1503251/157‑1 (MQ=255)
                                                                                          gggCGGCGTTAATG‑TTCTGGC‑TGGACATGTTTGATGTCCAG‑ATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAA‑TCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGAGGGAAATGGTGGGGGATAAAACAg                                                                                                                                                                                            >  1:1154066/1‑189 (MQ=255)
                                                                                           ggcggcGTTAATG‑TTCTGGC‑TGGACATGTTTGATGTCCAG‑ATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAA‑TCATGAATGCGGGGGGCCTGTCAAAAAGCCtttttgtt                                                                                                                                                                                                                           >  2:1503251/1‑155 (MQ=255)
                                                                                                         ttCTGGC‑TGGACATGTTTGATGTCCAG‑ATCATGGCACAACCGCTGGTGAACGGAACCGAGAGCTTCTCCCTGCTGGCGATCCCGTTCTTTGTGCTGGCGGGTGAAA‑TCATGAATGCGGGCGGCTTGTCAAAGCGAATTGTTAACCTGCCGATGAAACtggtg                                                                                                                                                                                                         <  1:1113463/162‑1 (MQ=255)
                                                                                                              gc‑tggACATGTTTGATGTCCAG‑ATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAA‑TCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGc                                                                                                                                                                                                                       <  1:1011733/143‑1 (MQ=255)
                                                                                                                   gACATGTTTGATGTCCAG‑ATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAA‑TCATGAATGCGGGCGGCTTGTCAAAACg                                                                                                                                                                                                                                     >  2:992893/1‑125 (MQ=255)
                                                                                                                   gACATGTTTGATGTCCAG‑ATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAA‑TCATGAATGCGGGCGGCTTGTCAAAACg                                                                                                                                                                                                                                     <  1:992893/125‑1 (MQ=255)
                                                                                                                         tttGATGTCCAG‑ATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGAAA‑TCATGAATGCGGGGGGCTTGTCAAAACGCATTTTTGTCCTGCCGGTTAAACt                                                                                                                                                                                                             >  2:1011733/1‑143 (MQ=255)
                                                                                                                            atggcccAG‑ATACTGCCAAAAACGTGGGTGAACGGGCCCGATGTTTTCTCCTTGGTGGCG‑TTCCTGTTTTTGGGCGGGGGGGGGAAATTAATGAGGGCGGGGGGTTTGAAAAACCGTTGTGTTGCCCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGC‑TGGGCTACGTCGGCGTGCTGGCAGCTA‑TGATTATGGCGAG‑CCTTTCTGGTTCTGCCG‑TTGCGGATACCGCCGCCGt                                                                                                      <  2:334331/234‑1 (MQ=255)
                                                                                                                                tCCAG‑ATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAA‑TCATGAATGCGGGCGGCTTGTCa                                                                                                                                                                                                                                          <  1:1125208/107‑1 (MQ=255)
                                                                                                                                  cAG‑ATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGAAA‑TCATGAATGCGGGGGGGCTGTCAAAACCCCTTGTTGAGCTGCCGATGAAAATGGTGGGGGATAAACCAGGCGGGG‑CGGGGTCAGTCGGGGTGCGGGGGGCTATTGGTTATGGGGGG‑CCTTTCTGGTTTTCCCT‑TGGGGGAGAccggcgccgggcgcc                                                                                                 >  2:55487/1‑233 (MQ=255)
                                                                                                                                    g‑ATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGAAA‑TCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTTACCTGCCGATGAAACTGGTGGGGGATAAACCAGGCGGGC‑TGGGCCACCTCGGCGGGCGGGCAGCTA‑TGATTATTGCGGG‑GCCTTCCTGTTCTGCCG‑GTGCGGGTACCGCCGCCGCTggccgc                                                                                               >  2:1306000/1‑233 (MQ=255)
                                                                                                                                        cATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGGAA‑TCATGAATGCGGGCGGCTTGTCAAAACCCATTGTTGACCTGCCGATGAAACTGGGGGGGCATAAAACAAGCGGGC‑GTGGGTACGCCGGCGCGCTGGCGGCTA‑TGATTTTGTCGAGCCCTTTTTGGTT‑TGCCC‑TGGTGCGGACCGCCGCCGCTGCCGcg                                                                                               >  2:198171/1‑234 (MQ=255)
                                                                                                                                           ggCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAA‑TCCTGCATGCGGGCGGCGTGTTCAAACGCCTTGTTGAACTGCCCGTGAAACTGGCGGGGCCTAAACCCGGCGGGCGTGGGC‑CCCTCCCGGGGCCGGGAGGCG‑TGATTATTGCGGG‑GCCTTCCGGTTTTGCCGTTTGCGGATCCCCCCACCGTTGCCGCCCTAcgt                                                                                         >  2:613943/1‑237 (MQ=255)
                                                                                                                                             cacaAGCGCGGGTGACCGAAGCCGATAACTTCCCCCGGCTGGCGTTTCCTTTTTTGGTGCTGGCGGGTAAAA‑TCATAAATGGGGGGGGCTTGTAAAAACGCATTGTTGACCTGCCGATGAAACtggtg                                                                                                                                                                                                         <  2:1225710/128‑1 (MQ=255)
                                                                                                                                                       ggTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGAAA‑ACATGAATGCGGGCGGGTTGTCAAAACGCCTTGTTGTACTCCCGATGAAACTGGTGGGGCATAAAc                                                                                                                                                                                               >  1:1225710/1‑128 (MQ=255)
                                                                                                                                                                                 tgctgGCGATTCCGTTCTTTGTGCTGGCGGGTGAAA‑TCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGCCCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGC‑TGGGCTACGTCGGCGTGCTgg                                                                                                                                                                <  1:462266/132‑1 (MQ=255)
                                                                                                                                                                                          tcccGTTTTTTGTGCGGGCGGGTGAAA‑TCATGAAGGGGGGGGGTTTGTCAAACCGCTTTGTTGCCCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGC‑TGGGCTACGTCGGCGTGCTGGCAGCTA‑TGATTATGGCGAG‑CCTTTCTGGTTCTGCCG‑TTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGA‑TGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCgggg                                        <  1:198171/237‑1 (MQ=255)
                                                                                                                                                                                                                        tagatggggggcgtCTTGCAAAACCGCTTTTTTACCCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGC‑TGGGCTACGTCGGCGTGCTGGCAGCTA‑TGATTATGGCGAG‑CCTTTCTGGTTCTGCCG‑TTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGA‑TGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTATCGCGcc           <  1:1283601/231‑1 (MQ=255)
                                                                                                                                                                                                                         tGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGC‑TGGGCTACGTCGGCGTGCTGGCAGCTA‑TGATTATGGCGAG‑CCTTTCTGGTTCTGCCG‑TTGCGGATACCGCCGCCGTTGCCGCTCTACCGGGGCCGAGGGATG‑GCCGCGCCAACTATCCGGTCAACCGGGGGGGGGGGCTGGATGCCTCTGGGGGCCTTTTCGCGCCa          >  2:1544722/1‑239 (MQ=255)
                                                                                                                                                                                                                          gAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGC‑TGGGCTACGTCGGCGTGCTGGCAGCTA‑TGATTATGGCGAG‑CCTTTCTGGTTCTGCCG‑GTGCGGATACCGCCGCCGGTGCCGCTCTACCGGGGCCGA‑TGATGCGCCGCGCCCACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCCTCCGGCGGCCTTATCGCGccca         >  2:1437773/1‑237 (MQ=255)
                                                                                                                                                                                                                             tGCGGGGGGCTGGAAAAACCGCTTTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGC‑TGGGCTACGTCGGCGTGCTGGCAGCTA‑TGATTATGGCGAG‑CCTTTCTGGTTCGGCCG‑TTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGA‑TGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTATCGCGCCAatta      <  1:613943/239‑1 (MQ=255)
                                                                                                                                                                                                                                  gcggcTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGC‑TGGGCTACGTCGGCGTGCTGGCAGCTA‑TGATTATGGCGAG‑CCTTTCTGGGTCTGCCG‑TTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGA‑TGATGCGCAGCGCCAACTATCCGGGCAACCGGGCGGCGGGGCTGGTTGCCTCTGGCGGGCTTTTCGCGCCCATTATTcc  >  2:853003/1‑238 (MQ=255)
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                 
GCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATG‑TTCTGGC‑TGGACATGTTTGATGTCCAG‑ATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAA‑TCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGC‑TGGGCTACGTCGGCGTGCTGGCAGCTA‑TGATTATGGCGAG‑CCTTTCTGGTTCTGCCG‑TTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGA‑TGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTATCGCGCCAATTATTCC  >  NC_000913/3744796‑3745255

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