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

GCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCAC‑TGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGG‑CTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTATCGCGC  >  NC_000913/3744783‑3745245
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                   
gCCTCCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCGC‑TGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATgcgggcgg                                                                                                                                                                                                                                   <  1:818867/239‑1 (MQ=255)
gCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCAC‑TGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGAAATCATGAATgggggcgg                                                                                                                                                                                                                                   >  1:946984/1‑239 (MQ=255)
        aaGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCAC‑TGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGCGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGG‑CTTGTCaa                                                                                                                                                                                                                          <  1:768915/239‑1 (MQ=255)
               ttAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCAC‑TGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGG‑CTTGTCAAAACGCAt                                                                                                                                                                                                                   <  1:454309/239‑1 (MQ=255)
                      tGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCAC‑TGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATTCGGGCGG‑CTTGTCAAAACGCAttgttta                                                                                                                                                                                                             >  1:1228771/1‑236 (MQ=255)
                      tGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCAC‑TGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGGAAACATTGATTCGGGGGGCCTTGTCAAACCGcttg                                                                                                                                                                                                                  >  2:794543/1‑231 (MQ=255)
                      tGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCAC‑TGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGAAATCATGAATGCGGGGGG‑CTTGTCCAAACGCATTGTTGa                                                                                                                                                                                                             >  1:180324/1‑238 (MQ=255)
                         cTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCAC‑TGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGGGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGGGG‑CTTGTCAAAACGCATTGTTGACCt                                                                                                                                                                                                          <  1:794543/238‑1 (MQ=255)
                         cTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCAC‑TGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCGGGACATGTTTGATGCCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGG‑CTTGTCAAAACGCATTGTTGACCt                                                                                                                                                                                                          <  2:391720/238‑1 (MQ=255)
                           gtgCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCAC‑TGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGGAATCATGCATGCGGGGGG‑CCTGTTCAAACGCATTTTTTACCttc                                                                                                                                                                                                        >  1:871643/1‑236 (MQ=255)
                                        tgggtggtctgtTGGGTGGTATGGCTATCGGTTTGCTTTTTGCCGGGGCTC‑TGTGGTTGTGCGGGGGGGCGTTAATGTTCGGGCGGGACGTGTTTGATGTCCAGATCATGGCACAACCGCTGGTGAGCGGAGCCGATAGCTTCTCCCTGCGGGCGATTCCGTTCTTTGTGCTGGGGGGTGAAATCATGAATGCGGGCGG‑TTTGTAAAAACGCATTGTTGACCTGCCGATGAAACtggtg                                                                                                                                                                                          <  2:1338211/233‑1 (MQ=255)
                                               tctgtTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCAC‑TGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGAAATCATGAATTCGGGGGG‑GTTGTCAAAACGCCCTGTTGTCCTGCCGGCGAAACTGGGGGGGCCTa                                                                                                                                                                                   >  1:310541/1‑239 (MQ=255)
                                                 gtgggggTGGTGTTGCTTGTGGTTGTCCGTTTGCCGGGGCCT‑TGTTGTTGTGGGGGGCGGGGTTAATGTTCTGGCTGGACATGTTGGAGGTCCAATTCAGGGCAAAAACGTTGGTGAACGGGGCAGATAGCTTCTCCCGGCGGGCGATTCCGTTCTTTGTGCTGGCGGGGGAAATCATGAATGCGGGGGG‑CTTGTAAAAACGAATTGTTGACCTGCCGATGAAACtggt                                                                                                                                                                                           <  2:1281082/225‑1 (MQ=255)
                                                   ttGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCAC‑TGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGG‑Ct                                                                                                                                                                                                                                <  1:499809/190‑1 (MQ=255)
                                                      ggtggtATTGCTCTCGGTTTGCGTATTGCCTGGGCAC‑TGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCGTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAAGGCGGGCGG‑CTTGTCa                                                                                                                                                                                                                           <  1:1240614/192‑1 (MQ=255)
                                                           tatggatATCGGTTTGCCTTGGGCCTGGGCTCGTGTT‑TTGTGCGGGGCGGGGTTAGTGTTGGGGCGGGACATTTTTGATGTCCAAGTCAGGGCACAACCGTGGGGGAGCGGGGCCGATAGCTTCTCCCGGCGGGGGATCCGGTTTTTTGTGCTGGCGGGGAAAATCAGGAATGCGGGGGG‑TTTGTCAAAACGCATTGTGGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCggg                                                                                                                                                                        <  2:1234661/232‑1 (MQ=255)
                                                           tATTGCTATCGGTTTGCCTATTGCCTGGGCAC‑TGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCGGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGG‑CTTGTAAAAACGCATTGTTGGCCTGCCGATGAAACTGGTGGGGCATAAACCAGGCggg                                                                                                                                                                        <  1:557543/238‑1 (MQ=255)
                                                           tATTGCTATCGGTTTGCCTATTGCCTGGGCAC‑TGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACAAGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGAAATCCTGAATTCGGGCGG‑CTTGTCAAAACGCATTGTTTACCTGCCGGTTAAACATGTgggggctaca                                                                                                                                                                                 >  1:1281082/1‑223 (MQ=255)
                                                                    cGGTTTGCCTATTGCCTGGGCAC‑TGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGAAATTCTGAATTCGGGCGG‑CCTGTCCACACCCCTTGTTTACCCGCCGAGTGAACTGGTGGGGGATAAACACGGCgggccggggtcc                                                                                                                                                               >  2:1335160/1‑234 (MQ=255)
                                                                     ggTTTGCCTATTGCCTGGGCAC‑TGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGGAAACATGAATTCGGGGGG‑GTTTGCCAAAAGCACTTTTTTGCTTGCGCTGTAACTGGGGGGGGCGAAAACAGGCggggtgggggtcc                                                                                                                                                              >  2:963882/1‑233 (MQ=255)
                                                                                tGCCTGGGCAC‑TTTTGTGGTGCGGGGGGGCGTTATTGTTCTGGCTGGACATGTTTGTTGTCGAGATCGTGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGTTTCCGTTCTTTGTGCGGGCGGGTGAAATCATGAATGCGGGCGG‑TTTGTCAAAACGATTTGTTGCCCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTgg                                                                                                                                                  <  2:946984/239‑1 (MQ=255)
                                                                                 gCCTGGGCAC‑TGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCTCAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGGTTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGG‑CTTGTCa                                                                                                                                                                                                                           <  1:528530/165‑1 (MQ=255)
                                                                                    tGGGCAC‑TGTGTTTGTGGGGGGCGGCTTTAATTTTCTGGCTGGACATTTTTGATGTCCAGATCATGGCACAACCGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGG‑CTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGc                                                                                                                                              <  1:954289/239‑1 (MQ=255)
                                                                                            tgttgttggGGGGGGGGGCTTTAATGTTCGGGCGGGACATGTTTGAGGCCCAGATCGTGGCACAAACGCTGGTGAACGGGGCCGATAGCTTCTCCCTGCTGGCGTTTCCGTTTTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGG‑TTTGTCAAAACGCATTGTTGCCCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTgg                                                                                                                                                  <  2:776777/228‑1 (MQ=255)
                                                                                                    ggcggggcgggTTTATTGTTTGGGCTGGCCTTGTTTGTGGTCCAGATCATGCCACAAGGGCTGGGGAAGGGGGCTGCTAGTTTGTCCCTGCGGGCGATTCCGTTCTTTGGGCGGGCGGGGGAAATCGGGAATGGGGGGGG‑GTGGTAAAAACGCTTTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAg                                                                                                                               <  2:606829/238‑1 (MQ=255)
                                                                                                      cGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGAAATCATGAATTGGGGGGG‑CCTGTCAAAACGCATTGTTTACCTGCCGATGAAACTGGGGGGGGATAAAACAGGCGGGCTGGGCTACGCCGGCGTGCTGGGCGCTATGAt                                                                                                                                        >  1:776777/1‑228 (MQ=255)
                                                                                                        ggggtgggTTAATGTTCGGGCGTGACATGTTTGTTTTCCAGTTCATGGCCAAACCGCTGGTGACCGGGGCCGGTAGCTTCTCCCTGCTGGCGTTTCCGTTCTTTGTGCTGGCGGGTGAAATCATGATTGCGGGGGG‑TTTGTCAAAACGCTTTGTTGCCCTGCCGATGAAACtggt                                                                                                                                                                                           <  2:104008/168‑1 (MQ=255)
                                                                                                                  aaTGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCAAGAATGCGGGGGG‑CTTGGCAAAACGCATTGTTTAACTGCCGATGAAAATGGTGGGGGATaaa                                                                                                                                                                                 >  1:104008/1‑175 (MQ=255)
                                                                                                                           gctggACATGTTTGGTGTCCTGATCATGGCACAAACGCTGGTGACCGGAGCCGATAGCTTCTCCCTGCTGGCGTTTCCGTTCTTTGTGCTGGCGGGTGAAATCATAAATGCGGGCGG‑CTTGTCAAACCGCATTGTTGCCCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTa                                                                                                                                                                <  2:665362/183‑1 (MQ=255)
                                                                                                                                     tttGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGAAAACATGAATGCGGGCGG‑CTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCCGGGGTACGgcggcggg                                                                                                                                                      >  1:665362/1‑181 (MQ=255)
                                                                                                                                            tCCAGATAATGGCACAACCGTTGTTGAAGGGGGCCGATAGCTTCCCCTTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTAAAATCATGAATGCGGGGGG‑CTTGTCAAAACGCTTTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCt                                                                                                                                             <  2:837984/185‑1 (MQ=255)
                                                                                                                                              cAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGAAATTATGAATGCGGGCGG‑CTTGGCAAAAAGCACTTTTGTGCTCCCGGCGAAAATGGGGGGGGCTaaa                                                                                                                                                                                 >  1:33845/1‑147 (MQ=255)
                                                                                                                                                       gCACAACCGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGG‑CTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTgg                                                                             <  1:963882/238‑1 (MQ=255)
                                                                                                                                                         acaAGCGCTGGTGAACGGAGCCGATTGCTTCCCCCGGCTGGCGTTTCCGTTTTTTGTGCGGGGGGGTAAAATCATGAATGGGGGGGG‑CTTGTCAAAACGCTTTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGc                                                                          <  1:240723/239‑1 (MQ=255)
                                                                                                                                                           gataacgggTGAACGGGGCCGATAGCTTCTCCCTGCGGGCGATTCCGTTTTTGGTGCGGGCGGGTGAAATAATAAAGGCGGGCGG‑CTTGTAAAACCGCATTGTGGCCCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTg                                                                                                                                                      <  2:36200/154‑1 (MQ=255)
                                                                                                                                                                     gAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAAACATTAATGCGGGGGG‑CTTGTCAAAACGCATTGTTTACCTGCCGGTGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGGGCTGGGAGCTa                                                                                                                                            >  1:36200/1‑161 (MQ=255)
                                                                                                                                                                          gggCCGATAGCTTCCCCCTGCTGGCGATTCCGTTCTTTGTGCTGGGGGGTGAAACCATGAATGCGGGGGG‑CTTGTCAAAACGCATTTTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAgcgc                                                          <  1:993941/236‑1 (MQ=255)
                                                                                                                                                                           agcggttggCTTCTCCCTGCGGGCAATTCCGTTTTTTGTGTGGGGGGGTAAAATCAGGAATGCGGGGGG‑CTTGTAAAAACGCATTGTTGGCCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAgc                                                            <  1:462104/227‑1 (MQ=255)
                                                                                                                                                                             ccGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGG‑CTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATAccgccgc                                                                                              <  1:141728/199‑1 (MQ=255)
                                                                                                                                                                                   gtttCTCCCTGCGGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGG‑TTTGCAAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCg                                                                                                                                                           <  2:741863/130‑1 (MQ=255)
                                                                                                                                                                                    cTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCAGGAAGGCGGGCGG‑TTTGTAAAAACGCTTTGTGGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCgggc                                                                                                                                                                       <  2:394246/119‑1 (MQ=255)
                                                                                                                                                                                     ttCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGG‑CTTGTCAAAACGCATTGTTGACCTGc                                                                                                                                                                                                        <  1:984149/85‑1 (MQ=255)
                                                                                                                                                                                             gctgGCGATTCCGTTCTTTGTGCTGGCGGGGGGAATCATGAATTCGGGCGG‑CTTTTCAAAACACATTGTTGACCTGCCGGTGAAAATGGGGGGGGCTAAAACAGGCGGGGTGGGCTACGTTGGCGGGCttgg                                                                                                                                                 >  1:741863/1‑129 (MQ=255)
                                                                                                                                                                                              ctgGCGTTTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGGGG‑CTTGTCAAAACGCATTGTTGACCTGCCGAt                                                                                                                                                                                                    <  1:666470/80‑1 (MQ=255)
                                                                                                                                                                                              ctgGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCCTGAATGCGGGCGG‑CTTTTCCAAACGCCTTGTTGACCTGCCGGTGAAACTGGGGGGGGATAAACCAGGCGGGCTGGGGTACgg                                                                                                                                                             >  1:394246/1‑118 (MQ=255)
                                                                                                                                                                                               tgGCGATTCCGTTCTTTGTGCTGGCGGGGGGAAAAATGAATGGGGGGGG‑CTTTTCAAAACGCATTGTTGGCCTGCCGATGaaaat                                                                                                                                                                                              >  2:984149/1‑83 (MQ=255)
                                                                                                                                                                                                   gATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGG‑TTTGTCAAAACGCTTTTTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGGGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGgcggcg                                 <  1:833124/238‑1 (MQ=255)
                                                                                                                                                                                                                             aaaTCAGGAAGGCGGGGGG‑TTTGTAAAACCGCATTGTTACCCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTAt       <  1:1069062/238‑1 (MQ=255)
                                                                                                                                                                                                                                cggagagtGCGGGGGG‑TTTGTAAAAGCGGTTGTTTGCCCTGCGGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTATcgcg   <  2:180324/232‑1 (MQ=255)
                                                                                                                                                                                                                                 cATGAAGGGGGGCGG‑TTTGTCAAACCGCATTGTTGCCCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTATCgcgc  <  1:155206/239‑1 (MQ=255)
                                                                                                                                                                                                                                             cgg‑cTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAgg                                                                                                                                                                            >  2:683163/1‑57 (MQ=255)
                                                                                                                                                                                                                                             cgg‑cTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAgg                                                                                                                                                                            >  2:652678/1‑57 (MQ=255)
                                                                                                                                                                                                                                             cgg‑cTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAgg                                                                                                                                                                            <  1:683163/57‑1 (MQ=255)
                                                                                                                                                                                                                                             cgg‑cTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAgg                                                                                                                                                                            <  1:652678/57‑1 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                   
GCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCAC‑TGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGG‑CTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTATCGCGC  >  NC_000913/3744783‑3745245

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