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

CACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTC‑TCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGC‑GGGC‑GGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGT‑GGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGC‑TATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTATCGCGCC  >  NC_000913/3744785‑3745246
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                      
cACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTC‑TCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGC‑GGGC‑GGCt                                                                                                                                                                                                                                   <  2:294740/239‑1 (MQ=255)
          aGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTC‑TCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGC‑GGGC‑GGCTTGTCAAAAc                                                                                                                                                                                                                          <  2:1668585/238‑1 (MQ=255)
           gCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTC‑TCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGC‑GGGC‑GGCTTGTCAAAACg                                                                                                                                                                                                                         <  1:347839/238‑1 (MQ=255)
            cTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTC‑TCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGAAATCATGAAAGC‑GGGC‑GGCCTTTCAAAAAGCa                                                                                                                                                                                                                       >  1:375476/1‑239 (MQ=255)
                     ggCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTC‑TCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGC‑GGGC‑GGCTTGTCAAAAc                                                                                                                                                                                                                          <  1:957869/227‑1 (MQ=255)
                        ttttttgttttttCTGGGTTGTTTGTTGGGTGGTATTGCTTTGGGTTTGCCTTTTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGTGGGACATGTTTGATTTCCAGATCATGGCACAAGCGGTGTTGAAGGGGGCCGATATCTTC‑TCCCTGCTGGCGTTTCCTTTCTTTTTGCTGGGGGGTGAAATCATGAATGC‑GGGG‑GTTTTGTCAAAACGCATTTTGGCCCTg                                                                                                                                                                                                            <  2:353766/230‑1 (MQ=255)
                          tgCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTC‑TCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGGAAACATGGATGC‑GGGG‑GGCTTGTTCAAAAGCATTGTTGACCTGCCg                                                                                                                                                                                                         >  1:1392344/1‑239 (MQ=255)
                          tgCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTC‑TCCCTGCAGGCGATTCCGTTCTTTGTGCTGGCGGGGGGAAACATGAATGC‑GGGCGGGCTTGTC‑AAACGCATTGGTGTCCTGCCg                                                                                                                                                                                                         >  2:751848/1‑239 (MQ=255)
                            tgtaTTTTTCTGGGCTGTTGGTTGGGGTGTTTGGCTTTCGTTTTGCCTTTTGCCGGGGCAGTGTTGTTGGGCGGGGCGGCGTTAATGTTCTGGCTGGACGTTTTTGATGTCCAGTTCATGGCACAAACGCGGGTGAACGGAGCCGATAGCTTC‑CCCCTGCTGGCGATTCCTTTCTTTGTGCTGGCGGGTGAAATCATGAAGGC‑GGGG‑GGCTTGTCAAAACGCATTGTTGACCTGCCGa                                                                                                                                                                                                        <  2:479658/235‑1 (MQ=255)
                               aTTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTC‑TCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGGAAACATGAAATG‑GGGG‑GGGTTGTCTACACGCAttgttgt                                                                                                                                                                                                                >  2:957869/1‑226 (MQ=255)
                                    tgtgggttttttgtTGGGTGGTATTGCTTTCGTTTTGCCTATGGCCGGGGCACTGTTGTTGTGCGGGGCGGAGTTATTGTTCGGGCTGGACATGTTTAATGTCCAGATACGGGCACAACCCCTGGTGAGCGG‑GCCTGTTGCCTCATCCCTGCTGGCGGTCCCGTTTTTTGTGCTGGGGGGTAAAGTCGAAAATGC‑GGGC‑GGTTTAGCAAAACACTTTGTTGACCTGCCGATGAAACtgg                                                                                                                                                                                               <  2:813942/228‑1 (MQ=255)
                                    tCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTC‑TCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAAAAAAGAATGCGGGGG‑GGGCTGCCAAAAC‑CATTTCTGACCTGCGGAGGACACtgg                                                                                                                                                                                               >  2:222374/1‑239 (MQ=255)
                                     gtGGGTTTTCTGTTGGGGGGTTTTTTGTTGGGTTTGCCTTTGGCCGGGGCTGTGTTGTTGTGGGGGGGGGCTTTAATTTTCTGGCTGGACATGTTTGATGTCCAGATCAGGGCACAACCGGTGGTGAAGGGGGCCGATAGCTTC‑CCCCTGCTGGCGTTTCCGTTTTTTGGGTTGGGGGGGAAAATCATGAAGGC‑GGGG‑GGTTTGTCAAAACGCTTTGTTGACCTGCCGATGAAACtggt                                                                                                                                                                                              <  2:782973/238‑1 (MQ=255)
                                      ggggCTGTCTTTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTC‑TCCCTGCTGGCAATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGC‑GGGC‑GGTTTGTCAAAACGCATTGTTGACCTGCCGATGAAACtggt‑g                                                                                                                                                                                            <  1:751848/238‑1 (MQ=255)
                                             tctgtTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTC‑TCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGC‑GGGC‑GGCt                                                                                                                                                                                                                                   <  1:925774/194‑1 (MQ=255)
                                                     ttggtTTTGCTGTCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCTTTAATGTTTTGGCTGGACATGTTAGATGTCAAGATCATGGCAGAAGCGCTGGTGAACGGAGCAGATAGCTTC‑TCCCTGCTGGCGATTCCGTTCTTTGGGCTGGCGGGTGAAATCATGAATGC‑GGGG‑GGTTTGTAAAAACGCATTGTTGACCTGCCGATGAAACTGGT‑GGGGCATAAACCAgg                                                                                                                                                                              <  1:222374/237‑1 (MQ=255)
                                                       ggtATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTC‑TCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGCTGAAATACTGAATGC‑GGGC‑GGGGTGTCTAAACg                                                                                                                                                                                                                         >  2:925774/1‑194 (MQ=255)
                                                         tATTGCTATCGGTTGGCCTATTGCCGGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCGGGCTGGACATGTTTGATGTCCAGATCATGGCACAAGCGCTGGTGAACGGAGCCGATAGCTTC‑TCCCTGCTGGCGATTCCATTCTTTGTGCTGGGGGGGAAAATCAGGAATGC‑GGGC‑GGCTTGTCAAAACGCATTGTTGCCCTGCCGATGAAACTGGT‑GGGGCATAAACCAGGCggg                                                                                                                                                                          <  1:1288470/238‑1 (MQ=255)
                                                                   ggTTTGCCTATCGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTC‑TCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGC‑GGGG‑GGCTTGTCAAAACGCATTGTTTACCTGCCGATGAAACTGGG‑GGGGCATAAACCCGGCgggctgg                                                                                                                                                                      >  1:1734944/1‑232 (MQ=255)
                                                                      ttGCCTATTGCTTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTC‑TCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGC‑GGGC‑GGGTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGG‑GGGGCATAAACCAGGCGGGCCGGGCTACGTCg                                                                                                                                                             >  1:1744173/1‑238 (MQ=255)
                                                                           tATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTC‑TCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGGAAACATGAATGC‑GGGG‑GGCTTGTCAAAAAGCAATTTTGACCTGCCGGTGAAACTGGT‑GGGGCATTAACCAGGCGGGGCGGGCTACggtcgtct                                                                                                                                                         >  1:1021845/1‑230 (MQ=255)
                                                                              tGCCTGGGCACTGTTGTTGGGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGGGCCGATAGCTTC‑TCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGC‑GGGC‑GGCTTGTCAAAACGCATTGTTGCCCTGCCGATGAAACTGGT‑GGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTgg                                                                                                                                                    <  1:155762/239‑1 (MQ=255)
                                                                               gCCTGGGCATTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTC‑TCCCTGCTGGCGATTCCGTTCTTGGTGCTGGCGGGTGAAATCATGAATGC‑GGGC‑GGCTTGTCAAAACGCATTGTTGACCTGCCAATGAAACtggt                                                                                                                                                                                              <  1:1761921/197‑1 (MQ=255)
                                                                                           ttgtggtGCGGGGCGGCGTTAATGTTCTGGCGGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAAGGGAGCCGATAGCTTC‑TCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGC‑GGGC‑GGCTTGTCAAACCGCATTGTTGACCTGCCGATGAAACTGGT‑GGGGCATAAACCAGGCggg                                                                                                                                                                          <  1:1031377/204‑1 (MQ=255)
                                                                                                     gggCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTC‑TCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGGAATCATGAATGC‑GGGG‑GGCCTGTCAAAACGCCTTGTTGACCTGCCGAAGAAAATGGG‑GGGGCCTAAAACACGGggggtggggcctg                                                                                                                                                                >  2:1031377/1‑199 (MQ=255)
                                                                                                     gggCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTC‑TCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGAAATCATGAATGC‑GGGC‑GGCTTGTCAAAACGCATTGTTGACCTGCCGATGaaaa                                                                                                                                                                                                  >  1:1536534/1‑170 (MQ=255)
                                                                                                       gcgggtttAATGTTCTGGCTGGACATTTTTGATGTCCAGATCAGGGCACAAACGCTGGGAAGCGGAGCCGATAGCTTC‑TCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGC‑GGGC‑GGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGT‑GGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGC‑TATGATTATGGCGAGCCttt                                                                                                                           <  1:205382/232‑1 (MQ=255)
                                                                                                                 tGTTCTGGCGGGACATGTTTGATGTCCAGATCATGGCACAAACGGTGGTGAACGGAGCCGATAGCTTC‑TCCCTGCGGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGC‑GGGC‑GGCTTGTCAAAACGCATGGTTGCCCTGCCGATGAAACTGGT‑GGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGC‑TATGATTATGGCGAGCCTTTCTGGTTCTGcc                                                                                                                <  1:438347/239‑1 (MQ=255)
                                                                                                                            gACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTC‑TCCCTGCTGGCGATTCCGTTTTTTGTGCGGGCGGGTGAAATCATGAATGC‑GGGC‑GGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGT‑GGGGCATaa                                                                                                                                                                                    <  1:1402100/161‑1 (MQ=255)
                                                                                                                                      aTGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTC‑TCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGAAAATATGAATTC‑TGGG‑GGCTCGTCAAAACGCATTGGTGGCCTGTCGCTGGAAATGGTGGGGGGATAAA‑CAGGCggg                                                                                                                                                                          >  2:1402100/1‑161 (MQ=255)
                                                                                                                                       tGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTC‑TCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGC‑GGGC‑GGCTTGTCAAAACGCAt                                                                                                                                                                                                                      <  1:1636847/117‑1 (MQ=255)
                                                                                                                                              aTCATGACAAAAACCCTGGTGAAGGGAGCAAATAGCTTC‑ACCCTGCTGACATTTCCGTTTTTTGGGCTGGCGGGTAAATACAAGAAGGC‑GGGG‑GGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACtggt                                                                                                                                                                                              <  1:1407706/134‑1 (MQ=255)
                                                                                                                                               tCATGGCAAAAACGCTGGTGAACGGAGCCGATAGCTTC‑TCCCTGCTGGCGATTCCTTTCTTTGTGCTGGCGGGTAAAATCATGAAGGC‑GGGC‑GGTTTGTCAAAACGCATTTTTGACCTGCCGATGAAACTGGT‑GGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGC‑TATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGct                                                                                     <  1:1598555/236‑1 (MQ=255)
                                                                                                                                               tCAGGGCACAAACGCTGGTGACCGGAGCCGATAGCTTC‑TCCCTGCTGGCGTTTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGC‑GGGC‑GGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGT‑GGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGC‑TATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCATTGCCGCTCTa                                                                                  <  1:1777298/239‑1 (MQ=255)
                                                                                                                                                cATGGCACAACGGGTGTGGAGGGGAGCCGATAGCTTC‑TCCCTGCTGGCGATTCCGTTTTTTGTGTTGGGGGGTAAAATCATGAATGC‑GGGC‑GGCTTGTCAAACCGCATTGTTGACCTGCCGATGAAACTGGT‑GGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGC‑TATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTAc                                                                                 <  1:165803/239‑1 (MQ=255)
                                                                                                                                                   ggCACAAACGCTGGTGAACGGAGCCGATAGCTTC‑TCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGAAAAAATGAATGC‑GGGG‑GGCCTGGCAAAAAGCATTGTTGACCTGCCGATGGAACTGGG‑GGGGCAAAAACCAGGCGGGGTGGGGCACGTCGGCGGGCGGGGAGCTTATGTTAGTGGGGGGCCTTGTTGGTTTTGCCGTGGCGGAGCCCGCGCCCGGTGCCggcctatgg                                                                               >  2:1703637/1‑231 (MQ=255)
                                                                                                                                                        aaaCGCTGGTGAACGGAGCCGATAGCTTC‑TCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGACATCCTTAATTC‑GGGC‑GGGTTGTCCAAAAGCCATGTTGGCCTGCCGATGAAACTGGT‑GGGGCCTaaa                                                                                                                                                                                   >  2:1407706/1‑134 (MQ=255)
                                                                                                                                                                          cggATAGCTTC‑CCCCGGGTGGCTTTTGTGTTCTTTGTGGGGGGGGGTCAAATGATGAGGGG‑GGGG‑GGGTTGACAAAGCGGTTTTTTGCCCTGCCGATGAAACTGGT‑GGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGC‑TATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGAATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCaa                                                        <  2:388955/236‑1 (MQ=255)
                                                                                                                                                                                       cccTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGC‑GGGC‑GGCTTGTCAAAACGCATTGTTGACCCGCCGATGAAACTGGT‑GGGGCATAAACCAGGCGGGCTTGGCTACGGCGGCGGGCCGGCAGG‑TATGAt                                                                                                                                         >  1:54608/1‑145 (MQ=255)
                                                                                                                                                                                              gggttttCCTTTCTTTGTGCTGGGGGGTAAAATTTATATTGG‑GGGC‑GGTTTGTAAAACCGCTTGTTTGCCCTGCCGATGAAACTGGT‑GGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGC‑TATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGggcgg                                    <  1:708149/234‑1 (MQ=255)
                                                                                                                                                                                                       gTTCTTTGTGCTGGCGGGTGAAATCATGAATGC‑GGGC‑GGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGT‑GGGGCATAAACCAGGCgggc                                                                                                                                                                         <  1:494705/98‑1 (MQ=255)
                                                                                                                                                                                                                                aTGAATGC‑GGGC‑GGCTTGCCAAAACGCATGTTTGACCTGCCGATGAAACTGGT‑GGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGC‑TATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTATCGCGcc  <  1:968237/239‑1 (MQ=255)
                                                                                                                                                                                                                                         gggc‑ggCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGT‑GGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGC‑TATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGGGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCCTCTGGCGGCAt           >  1:1508791/1‑222 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                      
CACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTC‑TCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGC‑GGGC‑GGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGT‑GGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGC‑TATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTATCGCGCC  >  NC_000913/3744785‑3745246

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: ‑