Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,745,0220CT32.6% 34.4 / 20.1 55G73G (GGC→GGTyiaN2,3‑diketo‑L‑gulonate TRAP transporter large permease protein
Reads supporting (aligned to +/- strand):  ref base C (16/15);  new base T (0/15);  total (24/31)
Fisher's exact test for biased strand distribution p-value = 5.52e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.62e-02
Rejected as consensus: Frequency below/above cutoff threshold.

AAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTATCGCGCCA  >  NC_000913/3744791‑3745247
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                     
aaGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCa                                                                                                                                                                                                                                 <  2:1086898/238‑1 (MQ=255)
     gCTTAGTCATGGCTGGTGTGATTTTTCTGGGCGGTCTGTTGGGTGGTTTTGCTATCGGTTTGCCTATTGCCTGGGCTATGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGAGTTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGGTTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAAGGCGGGCGGCTTGGCAAAACg                                                                                                                                                                                                                           <  2:2144000/238‑1 (MQ=255)
     gCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACg                                                                                                                                                                                                                           <  2:968452/238‑1 (MQ=255)
        tAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGAAAACATGAATGCGGGCGGCTTGTCAAAAAGCaat                                                                                                                                                                                                                       >  1:1894285/1‑237 (MQ=255)
                    tgCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGATGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGGAATTCTGAAAGCGGGGGGGTTGTCAAAACGCATTGTTTGCCTGCCg                                                                                                                                                                                                          >  1:953653/1‑239 (MQ=255)
                     tttgttttgtCTGGGTTGTTTGTTGGGTGTAATTGGTAGCGTTTTCCGTTTGGCTTGGGCTTTGTTTTTGTGCGGGGCGGTGTTATTGTTCTGGCTGAACTTGTTTGATGTCCAGATCAGGGCACAAACGCTGGTGAACGGAGCCGTTAGCTTCCCCCTGCTGGCGTTTCCGTTTTTTGTGCGGGCGGGTAAAAACATAAATGGGGGGGGTTTGTCAAACCCATTTGTTGTCCTGCCGa                                                                                                                                                                                                         <  2:2356693/234‑1 (MQ=255)
                                tGGGTTGTTTGTGGGGTGGTTTTGGTATGGGTTTGCCTGGTGCTGGGGCCTTGTTGTGGTGGGGGAGGGGGAAAGTGTTGTGACTGGAAATGTTTGATGTCGAAAGAATGGCAAAACCGGTGGTGACGGGAGCGGATAGCTTCCCCCGGCTGGCAATTCCGTTTTTTGTGTTGGCGGAAGAAATAATAAATGGGGGGGGTTTAAAAAAACGAATTGTTGACCTGCCGATGAAACtggtg                                                                                                                                                                                             <  2:277850/239‑1 (MQ=255)
                                gggtctgtctgtGGGGTGGTATTGCTATCGGTTTGCCTATTGCGGGGGCACTGTTGTTGTGGGGGGCGGCGTTAATGTTCTGGGTGGACAGGTTTGATGTCAAGATCATGGCACAAACCCTGGTGAACGGAGCCGATAGTTTCTCCCTGTGGGCGATTCCGTTCTTTGTGCTGGCGGGTAAAATCATAAAGGCGGGGGGTTTGTAAAAACGCATTGTGGACCTGCCGATGAAACtggtg                                                                                                                                                                                             <  2:2186473/235‑1 (MQ=255)
                                ggggCTGTTTTGTGGGGGGTTTGGCTATCGGTTTGCCTATTGCCGGGGCACTGTTGTTGTGCGGGGCGGCGTAAATGTTCGGGCGGAACATGTTTGATGTCCAGATCATGGCACAACGGCGGGTGAAGGGGGCCGATAGTTTCTCCCTGCTGGCGATTCCGTTTTTTGTGCTGGCGGGTGAAATCATGAATGGGGGGGGTTTGTAAAACCGCATTGTTGACCTGCCGATGAAACtggtg                                                                                                                                                                                             <  2:1952341/238‑1 (MQ=255)
                                       tctgtTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGAAATCATGAATGCGGGGGGGTTGTCAAAAAGCATTGTTGACCTGCCGGTGGAACTGGTggggggta                                                                                                                                                                                      >  1:2186473/1‑235 (MQ=255)
                                                tggtATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACg                                                                                                                                                                                                                           <  1:161036/195‑1 (MQ=255)
                                                 ggtATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCCGCCGATGAAACTGGTGGGGCATAAACCAGGCgg                                                                                                                                                                           >  1:572611/1‑239 (MQ=255)
                                                  gtATTGCTTTCGGTTGTCCTATGGCCTGGGCTTTGTTGTTGGGGGGGGCGGGGTTAATGTTCTGGCTGGACATGTTTGATGTCAAGATCATGGCACAACCGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGGGGTTTGTCAAAACGCATTGTTGCCCTGCCGATGAAACTGGTGGGGCATAAACCAgg                                                                                                                                                                              <  2:226875/235‑1 (MQ=255)
                                                        cTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGGAAACATGAAAGCGGGGGGCTTGTCCAAAAGCCTTTTTTACCTGCCGGTGAAAATGGTGGGGCATAAAACAGGCgggctgggg                                                                                                                                                                    >  1:27048/1‑238 (MQ=255)
                                                          aTCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATTAATGCGGGGGGCTTGTCAAAAAGCATTGTTTACCTGCCGCTGTAACTGGTGGGGGCTAAACCAGGCgggcgggggt                                                                                                                                                                   >  1:255321/1‑236 (MQ=255)
                                                          aTCGGTCTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGGAAACACGGAATGGGGGGGCTTTCTCAAACCAATAGTTTGACCGCCGGTGAAAACAGGGGGGGCAAAAACAGGGCgggttgggg                                                                                                                                                                    >  2:2000359/1‑237 (MQ=255)
                                                             ggTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGGAATAATGAATGCGGGGGGGTTGTGCAAAAGCCTTTTTGATCTGCCGGTGAAAAAGGTGGGGGGTAAAACAAGGGgggcgggggacg                                                                                                                                                                 >  2:2340885/1‑234 (MQ=255)
                                                             ggTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGGAATTATGAATTCGGGGGGCCTGTCAAAAAGCATTTTTGACCTGCCGGTGAAACTGGTGGGGGATAAAACAGGCGGGCTGGGCTACg                                                                                                                                                                >  1:1172428/1‑238 (MQ=255)
                                                             ggTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGAAATCATGAATGCGGGCGGCTTGTCAAAAAGCATTTTTGACCTGCCGGTGAAAATGGGGGGGGGTAAAACAGGCggggggggggcttt                                                                                                                                                               >  1:226875/1‑228 (MQ=255)
                                                             ggTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGAAATCATGAAAGCGGGGGGCTTGTCAAAAAGCATTGTTGGACTGCCGATGAAACTGGGGGGGGATTAAACAGGCggggcgggggtcg                                                                                                                                                                >  1:827339/1‑228 (MQ=255)
                                                                       ttGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGGAATCATGAATGTGGGGGGGTTGTCAAAACACCATTTTGTCCTGCCCGTGTAAAAGGGGGGGGGTAAACCCAGCGGGGTGGGGTGCGTCgggggggtg                                                                                                                                                     >  2:2350586/1‑232 (MQ=255)
                                                                         gCCTGGGCACTGTTGTTGTGGGGGGCGGGGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGTTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGc                                                                                                                                                   <  1:2350586/239‑1 (MQ=255)
                                                                                     tttttgtGGGGGGGGGTGTTATTTTTCGGGCGGAACATGTTGGATGCCAAGATAAGGGCACAAACGTTGGGGAAGGGGCCCGATAGTTTCTCCCTGCTGGCATTCCCGTTCTTTGTGCTGGGGGGTGAAATCATGAAGGGGGGGGGTTTGACAAAACGCTTTGTGGCCCTGCCGATGAAACTGGTGGGGCa                                                                                                                                                                                        <  2:2088257/188‑1 (MQ=255)
                                                                                           tgCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCg                                                                                                                                                                                                          <  1:627867/168‑1 (MQ=255)
                                                                                               gggCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCa                                                                                                                                                                                                                         <  2:895385/150‑1 (MQ=255)
                                                                                               gggCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCa                                                                                                                                                                                                                         >  1:895385/1‑150 (MQ=255)
                                                                                               gggCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATAATGAAAGCGGGCGGGTTGTCAAAACGCATTGTTGACCCGCCGGTGAAACTGGTGGGGGATAAACCAGGc                                                                                                                                                                             >  1:2088257/1‑191 (MQ=255)
                                                                                                     cgTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGGAATCATGAAAGCGGGGGGGTTTTCCAAAAGCATTGTTGTCCTGCCGATGAAAAtgg                                                                                                                                                                                               >  2:627867/1‑168 (MQ=255)
                                                                                                        tAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCt                                                                                                                                                                                                                                      <  1:55948/129‑1 (MQ=255)
                                                                                                                  gctggACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGAAATCATGAATGCGGGGGGGTTGTCaaaaag                                                                                                                                                                                                                            >  2:55948/1‑127 (MQ=255)
                                                                                                                      gACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGCCCt                                                                                                                                                                                                              <  1:595109/137‑1 (MQ=255)
                                                                                                                       aCATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCAAGAATGCGGGGGGCTTGTCAAAACGCCTTGTTGTCCTTCCGATGAAACAGGTGGGGGATTAAACAGGCgggctgggg                                                                                                                                                                    >  1:733846/1‑175 (MQ=255)
                                                                                                                                aTGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGAAATTATGAATGCGGGCGGGCTGTTTAAACGCATTGTTGGATTGCCGGTGaaa                                                                                                                                                                                                   >  2:595109/1‑137 (MQ=255)
                                                                                                                                 ttcccAGTTCGTGCCAAAAGCGCTGGGAACGGGGCCGATTAGCTCCCCCGGGCGGGCGATCCGGTTTTTTGTGTGGGGGGTGAAAATAAAAATGGGGGGGGGTTGTCAAAAGCGCATGTTTGCCCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTg                                                                                           <  2:255321/236‑1 (MQ=255)
                                                                                                                                     cAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACg                                                                                                                                                                                                                           <  2:1473829/110‑1 (MQ=255)
                                                                                                                                     cAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACg                                                                                                                                                                                                                           >  1:1473829/1‑110 (MQ=255)
                                                                                                                                     cAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAAACATGAATGCGGGGGGCTTTTCAAAACGCATTGTTGACCTGCCGATGAAACtggtg                                                                                                                                                                                             >  1:548598/1‑138 (MQ=255)
                                                                                                                                     cAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTAAAATCATGAATGCGGGGGGCTTTTCAAAACGCATTGTTGACCTGCCGATGAAACtggtg                                                                                                                                                                                             <  2:548598/138‑1 (MQ=255)
                                                                                                                                         tCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGGGGGTGAAATCATGAAGGCGGGCGGCTTGTAAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGc                                                                                       <  1:1701398/235‑1 (MQ=255)
                                                                                                                                             ggCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGGAATCATGAATGCGGGCGGCTTGTCAAAAAGCATTGTTTACCTGCCGATGAAAATGGGGGGGGATAAAACAGGCgggctggggcacgg                                                                                                                                                               >  1:1012585/1‑153 (MQ=255)
                                                                                                                                                        tGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGTTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCt                                                                                                                                               <  1:1061303/164‑1 (MQ=255)
                                                                                                                                                        tGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCa                                                                                                                                                                                                                                 <  1:1681060/86‑1 (MQ=255)
                                                                                                                                                                  aGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGGGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCgggc                                                                                                                                                                         <  1:1397987/128‑1 (MQ=255)
                                                                                                                                                                  aGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGGAATCCTGAATGCGGGGGGCTTGTCAAAAAGCCttg                                                                                                                                                                                                                      >  2:1681060/1‑86 (MQ=255)
                                                                                                                                                                           gCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCa                                                                                                                                                                                                                         <  1:33838/75‑1 (MQ=255)
                                                                                                                                                                                  ccTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGGGGGGGGTTTGTCAAACCGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATgcg                                                                  <  1:1069451/216‑1 (MQ=255)
                                                                                                                                                                                      ctgGCGATCCCTTTTTTTGGGGTGGGGGGTAAAATCAGGAAGGGGGGGGGTTTGTAAAAACGCATTTTTGCCCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCggg                                        <  1:835685/238‑1 (MQ=255)
                                                                                                                                                                                             ttCCGTTCTTTGTGCTGGCGGGTGAAATCATGAAGGGGGGGGGTTTGTCAAAACGCATTGTTGACCTGCCGATGaaa                                                                                                                                                                                                   <  1:893155/77‑1 (MQ=255)
                                                                                                                                                                                                          ggtGGCGGGGAAAATCGGGAAGGGGGGGGGTTGGGCAAAACGCGTTGTAGCCCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCctct                    <  1:2273074/236‑1 (MQ=255)
                                                                                                                                                                                                                       aTCATGAAGGGGGGCGTTTTGAAAAACCGAGTTTTTGCCCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTATcgc      <  2:827339/239‑1 (MQ=255)
                                                                                                                                                                                                                            tgatggcggggggTTTGTAAAAACGTTTTTTAAACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGTTACGGCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTATCGCGCCa  <  1:2000359/234‑1 (MQ=255)
                                                                                                                                                                                                                              agggggggcggTTTGAAAAAGCGCTTTGTGGGCCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTATCgcgc    <  1:1137495/230‑1 (MQ=255)
                                                                                                                                                                                                                              aaTGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCg                                                                                                                                                             >  1:400216/1‑82 (MQ=255)
                                                                                                                                                                                                                                  cgggcggCTTGTCAAAACGCATTGTTGACCTGCCGTTGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGc                                                                                                                                                <  1:1845803/91‑1 (MQ=255)
                                                                                                                                                                                                                                  cgggcggCTTGTCAAAACGCATTGTTGACCTGCCGTTGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGc                                                                                                                                                >  2:1845803/1‑91 (MQ=255)
                                                                                                                                                                                                                                    ggcggcTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAg                                                                                                                                 <  1:1447000/104‑1 (MQ=255)
                                                                                                                                                                                                                                    ggcggcTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAg                                                                                                                                 >  2:1447000/1‑104 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                     
AAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTATCGCGCCA  >  NC_000913/3744791‑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: