Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,745,0190CG52.9% 24.0 / 66.1 51G72G (GGC→GGG) yiaN2,3‑diketo‑L‑gulonate TRAP transporter large permease protein
Reads supporting (aligned to +/- strand):  ref base C (9/15);  new base G (14/13);  total (23/28)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.00e-03
Rejected as consensus: Frequency below/above cutoff threshold.

TCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGA‑TAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCA‑TGAATGCGGGCGGCTTGTCAAAA‑CGCATTG‑TTGACCTGCCGATG‑AAACTGGTGGGGCAT‑AAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTATCGCGCCA  >  NC_000913/3744781‑3745247
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                        
tCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGA‑TAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAAACC‑TGAAAgcgggc                                                                                                                                                                                                                                          >  1:277850/1‑239 (MQ=255)
          aaGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGA‑TAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCA‑TGAATGCGGGCGGCTTGTCa                                                                                                                                                                                                                                 <  2:1086898/238‑1 (MQ=255)
               gCTTAGTCATGGCTGGTGTGATTTTTCTGGGCGGTCTGTTGGGTGGTTTTGCTATCGGTTTGCCTATTGCCTGGGCTATGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGAGTTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGA‑TAGCTTCTCCCTGCTGGCGGTTCCGTTCTTTGTGCTGGCGGGTGAAATCA‑TGAAGGCGGGCGGCTTGGCAAAA‑Cg                                                                                                                                                                                                                           <  2:2144000/238‑1 (MQ=255)
               gCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGA‑TAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCA‑TGAATGCGGGCGGCTTGTCAAAA‑Cg                                                                                                                                                                                                                           <  2:968452/238‑1 (MQ=255)
                  tAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGA‑TAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGAAAACA‑TGAATGCGGGCGGCTTGTCAAAA‑AGCaat                                                                                                                                                                                                                       >  1:1894285/1‑237 (MQ=255)
                              tgCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGATGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGA‑TAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGGAATTC‑TGAAAGCGGGGGGGTTGTCAAAA‑CGCATTG‑TTTGCCTGCCg                                                                                                                                                                                                          >  1:953653/1‑239 (MQ=255)
                               tttgttttgtCTGGGTTGTTTGTTGGGTGTAATTGGTAGCGTTTTCCGTTTGGCTTGGGCTTTGTTTTTGTGCGGGGCGGTGTTATTGTTCTGGCTGAACTTGTTTGATGTCCAGATCAGGGCACAAACGCTGGTGAACGGAGCCGT‑TAGCTTCCCCCTGCTGGCGTTTCCGTTTTTTGTGCGGGCGGGTAAAAACA‑TAAATGGGGGGGGTTTGTCAAAC‑CCATTTG‑TTGTCCTGCCGa                                                                                                                                                                                                         <  2:2356693/234‑1 (MQ=255)
                                          tGGGTTGTTTGTGGGGTGGTTTTGGTATGGGTTTGCCTGGTGCTGGGGCCTTGTTGTGGTGGGGGAGGGGGAAAGTGTTGTGACTGGAAATGTTTGATGTCGAAAGAATGGCAAAACCGGTGGTGACGGGAGCGGA‑TAGCTTCCCCCGGCTGGCAATTCCGTTTTTTGTGTTGGCGGAAGAAATAA‑TAAATGGGGGGGGTTTAAAAAAA‑CGAATTG‑TTGACCTGCCGATG‑AAACtggtg                                                                                                                                                                                             <  2:277850/239‑1 (MQ=255)
                                          gggtctgtctgtGGGGTGGTATTGCTATCGGTTTGCCTATTGCGGGGGCACTGTTGTTGTGGGGGGCGGCGTTAATGTTCTGGGTGGACAGGTTTGATGTCAAGATCATGGCACAAACCCTGGTGAACGGAGCCGA‑TAGTTTCTCCCTGTGGGCGATTCCGTTCTTTGTGCTGGCGGGTAAAATCA‑TAAAGGCGGGGGGTTTGTAAAAA‑CGCATTG‑TGGACCTGCCGATG‑AAACtggtg                                                                                                                                                                                             <  2:2186473/235‑1 (MQ=255)
                                          ggggCTGTTTTGTGGGGGGTTTGGCTATCGGTTTGCCTATTGCCGGGGCACTGTTGTTGTGCGGGGCGGCGTAAATGTTCGGGCGGAACATGTTTGATGTCCAGATCATGGCACAACGGCGGGTGAAGGGGGCCGA‑TAGTTTCTCCCTGCTGGCGATTCCGTTTTTTGTGCTGGCGGGTGAAATCA‑TGAATGGGGGGGGTTTGTAAAAC‑CGCATTG‑TTGACCTGCCGATG‑AAACtggtg                                                                                                                                                                                             <  2:1952341/238‑1 (MQ=255)
                                                 tctgtTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGA‑TAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGAAATCA‑TGAATGCGGGGGGGTTGTCAAAA‑AGCATTG‑TTGACCTGCCGGTG‑GAACTGGTggggggta                                                                                                                                                                                      >  1:2186473/1‑235 (MQ=255)
                                                          tggtATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGA‑TAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCA‑TGAATGCGGGCGGCTTGTCAAAA‑Cg                                                                                                                                                                                                                           <  1:161036/195‑1 (MQ=255)
                                                           ggtATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGA‑TAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGAAATCA‑TGAATGCGGGCGGCTTGTCAAAA‑CGCATTG‑TTGACCCGCCGATG‑AAACTGGTGGGGCAT‑AAACCAGGCgg                                                                                                                                                                           >  1:572611/1‑239 (MQ=255)
                                                            gtATTGCTTTCGGTTGTCCTATGGCCTGGGCTTTGTTGTTGGGGGGGGCGGGGTTAATGTTCTGGCTGGACATGTTTGATGTCAAGATCATGGCACAACCGCTGGTGAACGGAGCCGA‑TAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCA‑TGAATGCGGGGGGTTTGTCAAAA‑CGCATTG‑TTGCCCTGCCGATG‑AAACTGGTGGGGCAT‑AAACCAgg                                                                                                                                                                              <  2:226875/235‑1 (MQ=255)
                                                                  cTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGA‑TAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGGAAACA‑TGAAAGCGGGGGGCTTGTCCAAA‑AGCCTTT‑TTTACCTGCCGGTG‑AAAATGGTGGGGCAT‑AAAACAGGCgggctgggg                                                                                                                                                                    >  1:27048/1‑238 (MQ=255)
                                                                    aTCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGA‑TAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCA‑TTAATGCGGGGGGCTTGTCAAAA‑AGCATTG‑TTTACCTGCCGCTG‑TAACTGGTGGGGGCT‑AAACCAGGCgggcgggggt                                                                                                                                                                   >  1:255321/1‑236 (MQ=255)
                                                                    aTCGGTCTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGA‑TAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGGAAACACGGAATG‑GGGGGGCTTTCTCAAA‑CCAATAGTTTGACC‑GCCGGTGAAAACAGGGGGGGCAA‑AAACAGGGCgggttgggg                                                                                                                                                                    >  2:2000359/1‑237 (MQ=255)
                                                                       ggTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGA‑TAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGGAATAA‑TGAATGCGGGGGGGTTGTGCAAA‑AGCCTTT‑TTGATCTGCCGGTG‑AAAAAGGTGGGGGGTAAAACAAGGGgggcgggggacg                                                                                                                                                                 >  2:2340885/1‑234 (MQ=255)
                                                                       ggTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGA‑TAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGGAATTA‑TGAATTCGGGGGGCCTGTCAAAA‑AGCATTT‑TTGACCTGCCGGTG‑AAACTGGTGGGGGAT‑AAAACAGGCGGGCTGGGCTACg                                                                                                                                                                >  1:1172428/1‑238 (MQ=255)
                                                                       ggTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGA‑TAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGAAATCA‑TGAATGCGGGCGGCTTGTCAAAA‑AGCATTT‑TTGACCTGCCGGTG‑AAAATGGGGGGGGGT‑AAAACAGGCggggggggggcttt                                                                                                                                                               >  1:226875/1‑228 (MQ=255)
                                                                       ggTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGA‑TAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGAAATCA‑TGAAAGCGGGGGGCTTGTCAAAA‑AGCATTG‑TTGGACTGCCGATG‑AAACTGGGGGGGGAT‑TAAACAGGCggggcgggggtcg                                                                                                                                                                >  1:827339/1‑228 (MQ=255)
                                                                                 ttGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGA‑TAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGGAATCA‑TGAATGTGGGGGGGTTGTCAAAA‑CACCATT‑TTGTCCTGCCCGTG‑TAAAAGGGGGGGGGT‑AAACCCAGCGGGGTGGGGTGCGTCgggggggtg                                                                                                                                                     >  2:2350586/1‑232 (MQ=255)
                                                                                   gCCTGGGCACTGTTGTTGTGGGGGGCGGGGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGA‑TAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCA‑TGAATGCGGGCGGTTTGTCAAAA‑CGCATTG‑TTGACCTGCCGATG‑AAACTGGTGGGGCAT‑AAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGc                                                                                                                                                   <  1:2350586/239‑1 (MQ=255)
                                                                                               tttttgtGGGGGGGGGTGTTATTTTTCGGGCGGAACATGTTGGATGCCAAGATAAGGGCACAAACGTTGGGGAAGGGGCCCGA‑TAGTTTCTCCCTGCTGGCATTCCCGTTCTTTGTGCTGGGGGGTGAAATCA‑TGAAGGGGGGGGGTTTGACAAAA‑CGCTTTG‑TGGCCCTGCCGATG‑AAACTGGTGGGGCa                                                                                                                                                                                        <  2:2088257/188‑1 (MQ=255)
                                                                                                     tgCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGA‑TAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCA‑TGAATGCGGGCGGCTTGTCAAAA‑CGCATTG‑TTGACCTGCCg                                                                                                                                                                                                          <  1:627867/168‑1 (MQ=255)
                                                                                                         gggCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGA‑TAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCA‑TGAATGCGGGCGGCTTGTCAAAA‑CGCa                                                                                                                                                                                                                         <  2:895385/150‑1 (MQ=255)
                                                                                                         gggCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGA‑TAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCA‑TGAATGCGGGCGGCTTGTCAAAA‑CGCa                                                                                                                                                                                                                         >  1:895385/1‑150 (MQ=255)
                                                                                                         gggCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGA‑TAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATAA‑TGAAAGCGGGCGGGTTGTCAAAA‑CGCATTG‑TTGACCCGCCGGTG‑AAACTGGTGGGGGAT‑AAACCAGGc                                                                                                                                                                             >  1:2088257/1‑191 (MQ=255)
                                                                                                               cgTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGA‑TAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGGAATCA‑TGAAAGCGGGGGGGTTTTCCAAA‑AGCATTG‑TTGTCCTGCCGATG‑AAAAtgg                                                                                                                                                                                               >  2:627867/1‑168 (MQ=255)
                                                                                                                  tAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGA‑TAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCA‑TGAATGCGGGCGGCt                                                                                                                                                                                                                                      <  1:55948/129‑1 (MQ=255)
                                                                                                                            gctggACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGA‑TAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGAAATCA‑TGAATGCGGGGGGGTTGTCaaaaag                                                                                                                                                                                                                            >  2:55948/1‑127 (MQ=255)
                                                                                                                                gACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGA‑TAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCA‑TGAATGCGGGCGGCTTGTCAAAA‑CGCATTG‑TTGCCCt                                                                                                                                                                                                              <  1:595109/137‑1 (MQ=255)
                                                                                                                                 aCATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGA‑TAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCA‑AGAATGCGGGGGGCTTGTCAAAA‑CGCCTTG‑TTGTCCTTCCGATG‑AAACAGGTGGGGGAT‑TAAACAGGCgggctgggg                                                                                                                                                                    >  1:733846/1‑175 (MQ=255)
                                                                                                                                          aTGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGA‑TAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGAAATTA‑TGAATGCGGGCGGGCTGTTTAAA‑CGCATTG‑TTGGATTGCCGGTG‑aaa                                                                                                                                                                                                   >  2:595109/1‑137 (MQ=255)
                                                                                                                                           ttcccAGTTCGTGCCAAAAGCGCTGG‑GAACGGGGCCGATTAGCTCCCCCGGGCGGGCGATCCGGTTTTTTGTG‑TGGGGGGTGAAAATA‑AAAATGGGGGGGGGTTGTCAAAAGCGCATGT‑TTGCCCTGCCGATG‑AAACTGGTGGGGCAT‑AAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTg                                                                                           <  2:255321/236‑1 (MQ=255)
                                                                                                                                               cAGATCATGGCACAAACGCTGGTGAACGGAGCCGA‑TAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCA‑TGAATGCGGGCGGCTTGTCAAAA‑Cg                                                                                                                                                                                                                           <  2:1473829/110‑1 (MQ=255)
                                                                                                                                               cAGATCATGGCACAAACGCTGGTGAACGGAGCCGA‑TAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCA‑TGAATGCGGGCGGCTTGTCAAAA‑Cg                                                                                                                                                                                                                           >  1:1473829/1‑110 (MQ=255)
                                                                                                                                               cAGATCATGGCACAAACGCTGGTGAACGGAGCCGA‑TAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAAACA‑TGAATGCGGGGGGCTTTTCAAAA‑CGCATTG‑TTGACCTGCCGATG‑AAACtggtg                                                                                                                                                                                             >  1:548598/1‑138 (MQ=255)
                                                                                                                                               cAGATCATGGCACAAACGCTGGTGAACGGAGCCGA‑TAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTAAAATCA‑TGAATGCGGGGGGCTTTTCAAAA‑CGCATTG‑TTGACCTGCCGATG‑AAACtggtg                                                                                                                                                                                             <  2:548598/138‑1 (MQ=255)
                                                                                                                                                   tCATGGCACAAACGCTGGTGAACGGAGCCGA‑TAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGGGGGTGAAATCA‑TGAAGGCGGGCGGCTTGTAAAAA‑CGCATTG‑TTGACCTGCCGATG‑AAACTGGTGGGGCAT‑AAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGc                                                                                       <  1:1701398/235‑1 (MQ=255)
                                                                                                                                                       ggCACAAACGCTGGTGAACGGAGCCGA‑TAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGGAATCA‑TGAATGCGGGCGGCTTGTCAAAA‑AGCATTG‑TTTACCTGCCGATG‑AAAATGGGGGGGGAT‑AAAACAGGCgggctggggcacgg                                                                                                                                                               >  1:1012585/1‑153 (MQ=255)
                                                                                                                                                                  tGGTGAACGGAGCCGA‑TAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCA‑TGAATGCGGGCGGTTTGTCAAAA‑CGCATTG‑TTGACCTGCCGATG‑AAACTGGTGGGGCAT‑AAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCt                                                                                                                                               <  1:1061303/164‑1 (MQ=255)
                                                                                                                                                                  tGGTGAACGGAGCCGA‑TAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCA‑TGAATGCGGGCGGCTTGTCa                                                                                                                                                                                                                                 <  1:1681060/86‑1 (MQ=255)
                                                                                                                                                                            aGCCGA‑TAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCA‑TGAATGCGGGGGGCTTGTCAAAA‑CGCATTG‑TTGACCTGCCGATG‑AAACTGGTGGGGCAT‑AAACCAGGCgggc                                                                                                                                                                         <  1:1397987/128‑1 (MQ=255)
                                                                                                                                                                            aGCCGA‑TAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGGAATCC‑TGAATGCGGGGGGCTTGTCAAAA‑AGCCttg                                                                                                                                                                                                                      >  2:1681060/1‑86 (MQ=255)
                                                                                                                                                                                     gCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCA‑TGAATGCGGGCGGCTTGTCAAAA‑CGCa                                                                                                                                                                                                                         <  1:33838/75‑1 (MQ=255)
                                                                                                                                                                                            ccTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCA‑TGAATGGGGGGGGTTTGTCAAAC‑CGCATTG‑TTGACCTGCCGATG‑AAACTGGTGGGGCAT‑AAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATgcg                                                                  <  1:1069451/216‑1 (MQ=255)
                                                                                                                                                                                                ctgGCGATCCCTTTTTTTGGGGTGGGGGGTAAAATCA‑GGAAGGGGGGGGGTTTGTAAAAA‑CGCATTT‑TTGCCCTGCCGATG‑AAACTGGTGGGGCAT‑AAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCggg                                        <  1:835685/238‑1 (MQ=255)
                                                                                                                                                                                                       ttCCGTTCTTTGTGCTGGCGGGTGAAATCA‑TGAAGGGGGGGGGTTTGTCAAAA‑CGCATTG‑TTGACCTGCCGATG‑aaa                                                                                                                                                                                                   <  1:893155/77‑1 (MQ=255)
                                                                                                                                                                                                                    ggtGGCGGGGAAAATCG‑GGAAGGGGGGGGGTTGGGCAAAA‑CGCGTTG‑TAGCCCTGCCGATG‑AAACTGGTGGGGCAT‑AAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCctct                    <  1:2273074/236‑1 (MQ=255)
                                                                                                                                                                                                                                 aTCA‑TGAAGGGGGGCGTTTTGAAAAAC‑CGAGTTT‑TTGCCCTGCCGATG‑AAACTGGTGGGGCAT‑AAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTATcgc      <  2:827339/239‑1 (MQ=255)
                                                                                                                                                                                                                                      tgatggcggggggTTTGTAAAAA‑CGTTTTT‑TAAACCTGCCGATG‑AAACTGGTGGGGCAT‑AAACCAGGCGGGCTGGGTTACGGCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTATCGCGCCa  <  1:2000359/234‑1 (MQ=255)
                                                                                                                                                                                                                                        agggggggcggTTTGAAAAAG‑CGCTTTG‑TGGGCCTGCCGATG‑AAACTGGTGGGGCAT‑AAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTATCgcgc    <  1:1137495/230‑1 (MQ=255)
                                                                                                                                                                                                                                        aaTGCGGGCGGCTTGTCAAAA‑CGCATTG‑TTGACCTGCCGATG‑AAACTGGTGGGGCAT‑AAACCAGGCGGGCTGGGCTACGTCg                                                                                                                                                             >  1:400216/1‑82 (MQ=255)
                                                                                                                                                                                                                                            cgggcggCTTGTCAAAA‑CGCATTG‑TTGACCTGCCGTTG‑AAACTGGTGGGGCAT‑AAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGc                                                                                                                                                <  1:1845803/91‑1 (MQ=255)
                                                                                                                                                                                                                                            cgggcggCTTGTCAAAA‑CGCATTG‑TTGACCTGCCGTTG‑AAACTGGTGGGGCAT‑AAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGc                                                                                                                                                >  2:1845803/1‑91 (MQ=255)
                                                                                                                                                                                                                                              ggcggcTTGTCAAAA‑CGCATTG‑TTGACCTGCCGATG‑AAACTGGTGGGGCAT‑AAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAg                                                                                                                                 <  1:1447000/104‑1 (MQ=255)
                                                                                                                                                                                                                                              ggcggcTTGTCAAAA‑CGCATTG‑TTGACCTGCCGATG‑AAACTGGTGGGGCAT‑AAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAg                                                                                                                                 >  2:1447000/1‑104 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                        
TCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGA‑TAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCA‑TGAATGCGGGCGGCTTGTCAAAA‑CGCATTG‑TTGACCTGCCGATG‑AAACTGGTGGGGCAT‑AAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTATCGCGCCA  >  NC_000913/3744781‑3745247

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