Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913933,5120GT25.9% 79.0 / 11.5 55A97S (GCC→TCC) ftsKcell division DNA translocase FtsK
Reads supporting (aligned to +/- strand):  ref base G (10/30);  new base T (14/0);  total (25/30)
Fisher's exact test for biased strand distribution p-value = 6.04e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGCCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTACACCATTCCCGTCATTATTGTCGGCGGTTGTT‑GGTTTGCCTGGCGTCATCAGTCCAGCGACGAAT‑ACATTGATTATTTTGCCGTTT‑CGCTACGCATCATT‑GGCGTTTTGGCGCTCATCCTTACCTCCTGTGGT‑CTGGCGGCAAT‑CAACGCTGACGATATCTGGTATTTTGCCTCCGGTGGCGTCATTGGCAGCTTACTAAGCACTACGCTAC‑AACCACTGCTACACAGTAGCGGGGGAACTATTGCGCTGCTCTGCGTTTGGGCAGCGGGCCTGACGTTGTTCACCG  >  NC_000913/933263‑933761
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                               
aCAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTGCGGACCCCAGCTGGTCGCAAACGGCCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTTGCGTGATGGCTTACTCCATTCCCGGCATTATTGTCGGCGGTTGTT‑GGTTTGc                                                                                                                                                                                                                                                                <  1:51686/251‑1 (MQ=255)
         aGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGCCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTACACCATTCCCGTCATTATTGTCGGCGGTTGTT‑GGTTTTCCTTGCggc                                                                                                                                                                                                                                                        >  1:302440/1‑248 (MQ=255)
         aGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGAACCCAGCTGGTCGCAAACGGCCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGAGTGATGGCTTACACCATTCCCGTCATTATTGTCGGCGGTATTT‑GGTTTTCCttgggtct                                                                                                                                                                                                                                                       >  1:48843/1‑244 (MQ=255)
              cGGCCGCCGCCTTCTGGATGCTTTGCTGATCCTTATTTTCCTGTTTGCCGTCTGGGTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGCCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTACACCATTCCCGTCATTATTGTCGGCGGTTGTT‑GGTTTGCCTGGGGTCATCAGt                                                                                                                                                                                                                                                  <  2:48843/251‑1 (MQ=255)
                  cgccgccTTCTGGAGGCGTTGCTGCTCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGCCTGGCATGAACCTATCCATAATTTAGGTGGGCTGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTACACCATTCCCGTCATTCTTTTCGGCGGTTGTT‑GGTTTgtcttgcggattcagttccg                                                                                                                                                                                                                                              >  2:313748/1‑232 (MQ=255)
                  cgccgccTTCTGGAGGCGTTGCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGCCTGGCATGAACCTATCCATACTTTAGGTGGGATGCACGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTACACCATTCCCGTCATTATTGTCGGCGGTTTTT‑GGGTTGccgtgcggtctcagtcc                                                                                                                                                                                                                                                >  2:211374/1‑234 (MQ=255)
                  cgccgccTTCTGGAGGCGTTGCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGCCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTACACCATTCCCGTCATTATTGTCGGCGGTTGTT‑GGTTTTCCTGGCGGCATCAGTcc                                                                                                                                                                                                                                                >  1:474868/1‑249 (MQ=255)
                   gccgccTTCTGGAGGCGTTGCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGCCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTACACCATTCCCGTCATTATTGTCGGCGGTTGTT‑GGTTTGCCTGGCGTCATCAGTCCAGc                                                                                                                                                                                                                                             <  1:211374/251‑1 (MQ=255)
                      gccTTCTGGAGGCGTTGCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGCCCCCAGCTGGTCGCAAACGGCCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTACACCATTCCCGTCATTATTGTCGGCGGTTGGT‑GGTTTTCCTGTGGTCTTCAGTctcgtggc                                                                                                                                                                                                                                          >  2:434459/1‑244 (MQ=255)
                        cgtCTGGAGGCGTTGCTGTTCCGTATTGTCCTGGTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGCCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTACACCATTCCCGTCATTATTGTCGGCGGTTGTT‑GGTTTGCCTGGCGTCATCAGTCCAGcgacga                                                                                                                                                                                                                                        <  2:413429/249‑1 (MQ=255)
                              gAGGCGTTGCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGCCCCCAGCTGGTCGCACACGGCCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTACACCATTCCCGTCATTATTGTCGGCGGCTGTT‑GGTTTGCCTGGCGTCATCAGTCCAGAGACGAAT‑ACAt                                                                                                                                                                                                                                 >  1:524817/1‑251 (MQ=255)
                               aGGCGTTGCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTCAGCTTTAACCCTTCGGACCCCAGCTGGTCGCCAACGGCCTGGCATGAACCTATCCATCATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTACACCATTCCCGTCATTATTGTCGGCGGTTGTT‑GTTTTTCCTTGCGTCATCACTcctccggcgaaaccatt                                                                                                                                                                                                                                 >  2:50410/1‑236 (MQ=255)
                                    ttGCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCCGCTGGTCGCCACCGGCCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTACACCATTCCCGTCATTATTGTCGGCGGTTGTT‑GGTTTTCCTGGCGTCATCAGTCCAacgacggaaaaatttattt                                                                                                                                                                                                                            >  1:262563/1‑238 (MQ=255)
                                      gCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGCCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTACACCATTCCCGTCATTATTGTCGGCGGTTGTT‑GGTTTGCCTGGCGTCATCAGTCCAGCGACGAAT‑ACATTGattat                                                                                                                                                                                                                          <  1:366577/250‑1 (MQ=255)
                                                   gtgcTGTTTGCCGTCTGGTTGAGGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGCCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGGGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTACACCATTCCCGTCATTATTGTCGGCGGTTGTT‑TGTTTGCCTGGCGTCATCAGTCCAGCGACGAAT‑ACATTGATTATTTTGCCGTTT‑CGCt                                                                                                                                                                                                           <  2:384259/248‑1 (MQ=255)
                                                     ccTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGCCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTACACCATTCCCGTCATTATTGTCGGCGGTTGTT‑GGTTTGCCTGGCGTCATCAGTCCAGCGACGAAT‑ACATTGATTATTTTGCCGTTT‑CGCTAc                                                                                                                                                                                                         <  1:400398/251‑1 (MQ=255)
                                                             ccGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGCCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTACACCATTCCCGTCATTATTGTCGGCGGTTGTT‑GGTTTGCCTGGCGTCATCAGTCCAGCGACGAAT‑ACATTGATTATTTTGCCGTTT‑CGCTACGCATCAtt                                                                                                                                                                                                 <  1:55668/251‑1 (MQ=255)
                                                                      tGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGCCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTACACCATTCCCGTCATTATTGTCGGCGGTTGTT‑GGTTTTCCTGGCGGCATCCGTCCCGCGCCGAAA‑ACCTTTATTTTTTTTCCGTTT‑CGCTTCGCCTCATT‑GGCGGTTTg                                                                                                                                                                                       >  2:384599/1‑251 (MQ=255)
                                                                              gCCTTCCTAAGCTTTACCCCTTCGGACCCCAGCTGGTCGCAAACGGCCGGGCATGAACCTAGCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGGTGTTCTGTATTGTTGGCGTGATGGCTTACACCATTCCCGTCATTATTGTCGGCTGTTGTT‑GGTTTGCCTGGCGTCATCATTCCAGCGACGAAT‑ACGTTGATTGTTTTGCCGTTG‑TTCTGCGCATCAGG‑GGCGTTTTGGCGCTc                                                                                                                                                                                 <  2:329975/249‑1 (MQ=255)
                                                                              gCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGCCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTACACCATTCCCGTCATTATTGTCGGCGGTTGTT‑GGGTTTCCTGGCGGCATCAGTCCCGCGACGGAT‑ACATTGGTTATTTTGCCGGTT‑CGCTTCGCATCATT‑GGCGTTTTTGGGCTCAt                                                                                                                                                                               >  2:105622/1‑251 (MQ=255)
                                                                                  tACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGCCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTACACCATTCCCGTCATTATTGTCGGCGGTTGTT‑GGTTTGCCTGGCGTCATCAGTCCAGCGACGAAT‑ACATTGATTATTTTGCCGTTT‑CGCTACGCATCATT‑GGCGTTTTGGCGCTCATcc                                                                                                                                                                             <  1:447549/249‑1 (MQ=255)
                                                                                                ccTTCGGGCCCCAGCTGGTCGCAAACGGCCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTACACCATTCCCGTCATTATTGTCGGCGGTTGTT‑GGTTTGCCTGGCGTCATCAGTCCAGCGACGAAT‑ACATTGATTATTTTGCCGTTT‑CGCTACGCATCATT‑GGCGTTTTGGCGCTCATCCTTCCCTCCTGTGGT‑Ct                                                                                                                                                            <  2:410767/251‑1 (MQ=255)
                                                                                                ccTTCGGACCCCAGCTTGACGCAAACGGCCTGGCATGAACCTATCCATAGTGTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTCATTTTTGGCGTGATGGCTTACACCATTCCCGTCATTATTGTCGGCGGTTGTG‑GGTTTGCCTGGCGTCATCAGTCCAGCGACGAAT‑ACATTGATTATTTTGCCGGTT‑CTCTACGCATCATT‑GGCGTTTTGGCGCTCATCCTTACCTCCTGTGGT‑Ct                                                                                                                                                            <  1:230706/251‑1 (MQ=255)
                                                                                                          ccAGCTGGTCGCAAACGGCCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTACACCATTCCCGTCATTATTGTCGGCGGTTGTT‑GGTTTGCCTGGCGTCATCAGTCCAGCGACGAAT‑ACATTGATTATTTTGCCGTTT‑CGCTACGCATCATT‑GGCGTTTTGGCGCTCATCCTTACCTCCTGTGGT‑CTGGCGGCAAT‑c                                                                                                                                                 <  2:406975/251‑1 (MQ=255)
                                                                                                                gggcGCAAACGGCCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTACACCATTCCCGTCATTATTGTCGGCGGTTGTT‑GGTTTGCCTGGCGTCATCAGTCCAGCGACGAAT‑ACATTGATTATTTTGCCGTTT‑CGCTACGCATCATT‑GGCGTTTTGGCGCTCATCCTTACCTCCTGTGGT‑CTGGCGGCAAT‑CAACGCt                                                                                                                                           <  2:147291/248‑1 (MQ=255)
                                                                                                                        aCGGCCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTACACCATTCCCGTCATTATTGTCGGCGGTTTTT‑GGTTTTCCTGGCGTCTTCAGTTCCGCGACGAAT‑AAATTTATTTTTTTGCCGTTT‑CGCTCCGCCTCATT‑GGCGTTTTTGCGCTCATCCTTTCCTTCCGTGTT‑CTGTCGGCAAttcacgcttgccgtat                                                                                                                                    >  2:402018/1‑236 (MQ=255)
                                                                                                                                  aTGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTACACCATTCCCGTCATTATTGTCGGCGGATGTT‑GGTTTTCCTGGCGGAATCAGTCCAGCGGAGACA‑ACATTGATTATTTTGCCGTTT‑TGCTACGCACCATT‑GGCGGTTTGGCGCTCATACTTACCTCCTGTGGT‑CTTGGGGCAAT‑CCAAGCTGGCGAGATATGGTCtttt                                                                                                                         >  2:444718/1‑251 (MQ=255)
                                                                                                                                    gAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGATGTTCTTTATTTTTGGCGTGATGGCTTACAACATTCCCGTCATTATTGTCGGCGGTTGTT‑GGGTTGCCTGGCGTCATTAGTCCAGCGAGCAAT‑ACATTGATTATTTTGCCGTTTTCG‑TACGCACCATTTGGCGTTTTG‑CGGTCATCCTTTCCTCCTGTGGT‑CTGGGGGCAAT‑CAAGACTGACGATATCTGGGATTTTGc                                                                                                                       >  1:343024/1‑251 (MQ=255)
                                                                                                                                     aaCCTATCGATCATTTAGGTGGGATGCCAGGTGCGTCGTTGGCCGATATGCTGTTCTTTATTTTTGGCGTGATGGCTTACACCGTTAAGGTCATTATTGACGGCGGTTGTT‑GGTTTGCCTGGCGATATCTTTACAGCCACGAATGCCATT‑ATTATTTTAACGTTT‑CCATACGCAACCAT‑GGCGTTTTGTCGCTCCCCCTTACCTCCTGTTGT‑CTGGCGGCCAA‑CAAAGCGGACGatgattggtgtttttcc                                                                                                                      >  1:356655/1‑236 (MQ=255)
                                                                                                                                       ccTATCCATAATTTAGGTGGGATGCCCGGTTCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTACACCATTCCCGTCATTATTGTCGGCGGTTGTT‑GGTTTGCCTGGCGTCATCAGTCCAGCGACGAAT‑ACATTGATTATTTTGCCGTTT‑CGCTACGCATCATT‑GGCGTTTTGGCGCTCATCCTTACCTCCTGTGGT‑CTGGCGGCAAT‑CAACGCTGACGATATCTGGTATTTTGCCTc                                                                                                                    <  1:458894/251‑1 (MQ=255)
                                                                                                                                       ccTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTACACCATTCCCGTCATTATTGTCGGCGGTTGTT‑GGTTTGCCTGGCGTCATCAGTCCAGCGACGAAT‑ACATTGATTATTTTTCCGTTT‑CGCTACGCCTCATT‑GGCGTTTTTGCGCTCATCCTTACCTCCTGTGGT‑CTGGCGGCAAT‑CAACGCTGACGATATCTGGTATTTTGCCTc                                                                                                                    >  1:525699/1‑251 (MQ=255)
                                                                                                                                            ccATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTACACCATTCCCGTCATTATTGTCGGCGGTTGTT‑GGTTTGCCTGGCGTCATCAGTCCAGCGACGAAT‑ACATTGATTATTTTGCCGTTT‑CGCTACGCATCATT‑GGCGTTTTGGCGCTCATCCTTACCTCCTGTGGT‑CTGGCGGCAAT‑CAACGCTGACGATATCTGGTATTTTGCCTCCGGt                                                                                                                <  1:203890/250‑1 (MQ=255)
                                                                                                                                                      ggTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTACACCATTCCCGTCATTATTGTCGGCGGTTGTT‑GGTTTGCCTGGCGTCATCAGTCCAGCGACGAAT‑ACATTGATTATTTTGCCGTTT‑CGCTACGCATCATT‑GGCGTTTTGGCGCTCATCCTTACCTCCTGTGGT‑CTGGCGGCAAT‑CAACGCTGACGATATCTGGTATTTTGCCTCCGGTGGCGTCATTgg                                                                                                     <  1:8312/251‑1 (MQ=255)
                                                                                                                                                         gggATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTACACCATTCCCGTCATTATTGTCGGCGGTTGTT‑GGTTTGCCTGGCGTCATCAGTCCAGCGACGAAT‑ACATTGATTATTTTGCCGTTT‑CGCTACGCATCATT‑GGCGTTTTGGCGCTCATCCTTACCTCCTGTGGT‑CTGGCGGCAAT‑CAACGCTGACGATATCTGGTATTTTGCCTCCGGTGGCGTCATTGGCa                                                                                                   <  1:5916/250‑1 (MQ=255)
                                                                                                                                                           gATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGCTGGCTTACACCATTCCCGTCATTATTGTCGGCGGTTGTT‑GGGTTGCCTGGCGTCATCAGTCCAGCGCCGGAA‑AAATTGATTTTTTTGCCGTTT‑CGCTACGCCTCATT‑GGCGGTTTTGCGGTCCTCCTTACCTCCTGTGGT‑CTGGCGGCAAT‑CAACGCCGACGATATCTGGTAtttttccagatcggaagagcggcgcgg                                                                                                  >  2:47840/1‑226 (MQ=255)
                                                                                                                                                           gATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTACACCATTCCCGTCATTATTGTCGGCGGTTGTT‑GGTTTGCCTGGCGTCATCAGTCCAGCGACGAAT‑ACATTGATTATTTTGCCGTTT‑CGCTACGCATCAGT‑GGCGTTTTGGCGCTCATCCGTACCGCCTGTGGT‑CTGGCGGCAAT‑CAACGCTGACGATATCTGGTATTTTGcc                                                                                                                      <  1:47840/229‑1 (MQ=255)
                                                                                                                                                                  ggTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTACACCATTCCCGTCATTATTGTCGGCGGTTGTT‑GGTTTGCCTGGCGTCATCCGTCCAGCGACGAAT‑ACATTGTTTATTTTTCCGTTT‑TGCTACGCATCATT‑TGCGTTTTTGCGCTCATCCTTACCTCCTGTGGG‑CTGGCGGCCAT‑CCACGCTGGCGATATCTGGGTTTTTTCCTCCGGGTGCGTCCTTGGCAGCTTACTa                                                                                           >  2:110821/1‑249 (MQ=255)
                                                                                                                                                                          gTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTACACCATTCCCGTCATTATTGTCGGCGGTTGTTTGGTTT‑CCCGGCGTCCTCAGTTCCGCGACGAAT‑ACATTGATTTTTTTTCCGTTT‑TGCTACCCCTCATT‑GGCGTTTTTGGGCTCATACTTACCTCCAGTGGTACTG‑CGGCCAT‑CAACGCTGCCGCTTTCTTGTATTTTTCCTCCGGTGGCGTCCTTGGACGCTTACTaaaccctcagc                                                                                 >  2:451149/1‑242 (MQ=255)
                                                                                                                                                                             ggCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTACACCATTCCCGTCATTATTGTCGGCGGTTGTT‑GGTTTTTCTGGCGTCATTCGTCCCGCGAACAAT‑ACATTTATTTTTTTGCCGTTT‑CGCTACCCCTTATT‑GGCGTTTTTGCGCCCATCCTTTCCTCCTGTGGT‑CTGGCGGCAAT‑CAACGCTGCCGACATCTGGTCTTTTGCCTCCCGGGGCGTCATTGGCAGCTTACTTAGCCCTACGccac                                                                              >  2:519878/1‑248 (MQ=255)
                                                                                                                                                                               cAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTACACCATTCCCGTCATTATTGTCGGCGGTTGTT‑GGGTTTACTGGCGTCATCAGTCCCGCGACGAAT‑ACATTGATTATTTTGCCGTTT‑TGCTACGCATCCTT‑GGCGTTTTTGCGCTCATCCTTACCTCCTGTGGG‑CTGGCGGCAAT‑CAACGCTGACGATATCTGGGATTTTTCCTCCGGGGGCGTCATTGGCAGCTTACTAAGCACTACGCTAC‑a                                                                            >  1:206600/1‑250 (MQ=255)
                                                                                                                                                                                         gTTCTTTATTTTTGGCGTGATGGCTTACACCATTCCCGTCATTATTGTCGGCGGTTGTT‑GGGTTGCCTGGCGTCATCAGTCCCGCGACGAAT‑ACATTGATTATTTTGCCGTTT‑TGCTACGCACCCTT‑GGCGTTTTGGCCCCCATCCTTACCTCCTGTGGG‑CTGGCGGCAAT‑CAACCCTGACGATATCTGGTATTTTGCCTCCCGTGGCGTCATTGGCCGCCTACTAAGCACCCAGCTAC‑AAccccccctac                                                                 >  2:140861/1‑243 (MQ=255)
                                                                                                                                                                                            cgttATTTTGGGCGTGATGGCTTACACCATTCCCGTCATTATTGTCGGCGGTTGTT‑GGTTTGCCTGGCGTCATCAGTCCAGCGACGAAT‑ACATTGATTATTTTGCCGTTT‑CGCTACGCATCATT‑GGCGTTTTGGCGCTCATCCTTACCTCCTGTGGT‑CTGGCGGCAAT‑CAACGCTGACGATATCTGGTATTTTGCCTCCGGTGGCGTCATTGGCAGCTTACTAAGCACTACGCTAC‑AACCACTGCTacac                                                               <  2:55518/248‑1 (MQ=255)
                                                                                                                                                                                                  ttttGGCGTGATGGCTTACACCATTCCCGTCATTATTGTCGGCGGTTGTT‑GGTTTGCCTGGCGTCATCAGTCCAGCGACGAAT‑ACATTGATTATTTTGCCGTTT‑CGCTACGCATCATT‑GGCGTTTTGGCGCTCATCCTTACCTCCTGTGGT‑CTGGCGGCAAT‑CAACGCTGACGATATCTGGTATTTTGCCTCCGGTGGCGTCATTGGCAGCTTACTAAGCACTACGCTAC‑AACCACTGCTACACAGTAGc                                                         <  1:444718/250‑1 (MQ=255)
                                                                                                                                                                                                  ttttGGCGTGATGGCTTACACCATGCCCGTCATTATTGTCGGCGGTTGTT‑GGTTTGCCTGGCGTCATCAGTCCAGCGACGAAT‑ACATTGATTATTTTGCCGTTT‑CGCTACGCATCATT‑GGCGTTTTGGCGCTCATCCTTACCTCCTGTGGT‑CTGGCGGCAAT‑CAACGCTGACGATATCTGGTATTTTGCCTCCGGTGGCGTCATTGGCAGCTTACTAAGCACTACGCTAC‑AACCACTGCTACACAGTAGc                                                         <  2:302440/250‑1 (MQ=255)
                                                                                                                                                                                                    ttGGCGTGATGGCTTACACCATTCCCGTCATTATTGTCGGCGGTTGTT‑GGTTTGCCTGGCGTCATCAGTCCAGCGACGAAT‑ACATTGATTATTTTGCCGTTT‑CGCTACGCATCATT‑GGCGTTTTGGCGCTCATCCTTACCTCCTGTGGT‑CTGGCGGCAAT‑CAACGCTGACGATATCTGGTATTTTGCCTCCGGTGGCGTCATTGGCAGCTTACTAAGCACTACGCTAC‑AACCACTGCTACACAGTAGCgg                                                       <  1:110821/250‑1 (MQ=255)
                                                                                                                                                                                                              ggCTTACACCATTCCCGTCATTATTGTCGGCGGTTGTT‑GGTTTGCCTGGCGTCATCAGTCCAGCGACGAAT‑ACATTGATTATTTTGCCGTTT‑CGCTACGCATCATT‑GGCGTTTTGGCGCTCATCCTTACCTCCTGTGGT‑CTGGCGGCAAT‑CAACGCTGACGATATCTGGTATTTTGCCTCCGGTGGCGTCATTGGCAGCTTACTAAGCACTACGCTAC‑AACCACTGCTACACAGTAGCGGGGGAACTATTg                                            <  1:402018/251‑1 (MQ=255)
                                                                                                                                                                                                              ggCTTACACCATTCCCGTCATTATTGTCGGCGGATGTT‑GGGTTTCCTGGCGGTATCAGTCCAGCGACGAAT‑ACATTGTTTATTTTTGCGTTT‑TGCTCCGCATCATT‑GGCGGTTTGGGGCTCATCCTTACCTCCTGTGGT‑CTGGCGGCAAT‑CAACCCTGACGGTATCTGGGTTTTTGCCTCCGGTGGCGTCAATTGCAGCTTTCTAAGCACTACGCTACCAACCAC‑GCTAACCAGTAGGGGGGGGACtcatt                                            >  2:264249/1‑247 (MQ=255)
                                                                                                                                                                                                              ggCTTACAACATTCCCGTCATTATTGTCGGCGGTTGTT‑GGTTTGCCTTGCGTCATCAGTCCAGCGACGAAT‑ACATTGATTATTTTGCCGTTT‑CGCTACGCATCATT‑GGCGTTTTGGCGCTCATCCTTACCTCCTGTGGT‑CTGGCGGCAAT‑CAACTCTGACGATATCTGGTATTTTGCCTCCGGTGGCGTCATTGGCAGCTTACTAAGCACTACGCTAC‑AACCACTGCTACACAGTAGCGGGGGAACTATTg                                            <  1:317407/251‑1 (MQ=255)
                                                                                                                                                                                                               gCTTACACCATTCCCGTCATTATTGTCGGCGGTTGTT‑GGTTTGCCTGGCGTCATCAGTCCAGCGACGAAT‑ACATTGATTATTTTGCCGTTT‑CGCTACGCATCATT‑GGCGTTTTGGCGCTCATCCTTACCTCCTGTGGT‑CTGGCGGCAAT‑CAACGCTGACGATATCTGGTATTTTGCCTCCGGTGGCGTCATTGGCAGCTTACTAAGCACTACGCTAC‑AACCACTGCTACACAGTAGCGGGGGAACTATTgc                                           <  1:237622/251‑1 (MQ=255)
                                                                                                                                                                                                                      ccATTCCCGTCATTATTGTCGGCGGTTGTT‑GGTTTGCCTGGCGTCATCAGTCCAGCGACGAAT‑ACATTGATTATTTTGCCGTTT‑CGCTACGCATCATT‑GGCGTTTTGGCGCTCATCCTTACCTCCTGTGGT‑CTGGCGGCAAT‑CAACGCTGACGATATCTGGTATTTTGCCTCCGGTGGCGTCATTGGCAGCTTACTAAGCACTACGCTAC‑AACCACTGCTACACAGTAGCGGGGGAACTATTGCGCTGctc                                    <  1:105622/251‑1 (MQ=255)
                                                                                                                                                                                                                         ttCCCGTCATTATTGTCGGCGGTTGTT‑GGTTTTCCTGGCGTCCTCACTCCCGCGACGGAT‑TCATTGATTATTTTGCAGTTT‑CGCTACGCATCATT‑GGGGCTTTGGCGCTCATCCTTACCTCCTGTGGT‑CTGGCGGCAATTCAA‑GCTGGCGATATCTGGGATTTTGCCTCCGGTGGGGTCAATGGCAGCTTACTAAGCACAACGATAC‑AACCACTGCTACCCAGTAGCGGGGGAACAATTGCGCTGctcctc                                 >  2:465061/1‑248 (MQ=255)
                                                                                                                                                                                                                             cGTCATTATTGTCGGCGGTTGTT‑GGTTTGCCTGGCGTCATCAGTCCAGCGACGAAT‑GCATTGATTATTTTGCCGTTT‑CGCTACGCATCATT‑GGCGTTTTGGCGCTCATCCTTGCCTCCTGTGGT‑CTGGCGGCAAT‑CAACGCTGACGATATCTGGTATTTTGCCTCCGGTGGCGTCATTGGCAGCTTACTAAGCACTGCGCTAC‑AACCACTGCTACACAGTAGCGGGGGAACTATTGCGCTGCTCTGCGttt                             <  2:262563/251‑1 (MQ=255)
                                                                                                                                                                                                                             cGTCATTATTGTCGGCGGTTGTT‑GGTTTGCCTGGCGTCAGCAGTCCAGCGACGAAT‑ACATTGATGATGTTGCCGTTT‑CGCTACGCATCATT‑GGCGTTTTTGCGCTCATCCTTACCTCCTGTGGT‑CTGGCGGCAAT‑CAACGCTGACGATATCTTGTATTTTGCCTCCGGTGGCTTCATTTGCATGTTACTAAGCACTACGCTAC‑AACTACTTCTACACAGTAGCGGGGGAACTATTGCGCTGCTCTGCGttt                             <  1:342153/251‑1 (MQ=255)
                                                                                                                                                                                                                                 attattGTCGGCGGTTGTT‑GGTTTGCCTGGCGTCATCACGCCCGCGAAGAAT‑ACATTGATTATTTTGCCGTTT‑CGCTACGCATCATT‑GGCGCTTTTGCGATCATCCTTTCCTCCTCTGGT‑CTGGCGGCCAT‑AAACGCTGACGCAATCTGGCCTTTTGCCTCCGGAGGCGGCATTGGCAGCTTACTAAGCACTACGCTAC‑AACCACTGCTTCACAGTAGCGGGGGGACTATTTCGCCGCTTTGCGGTTGGGc                         >  1:351309/1‑251 (MQ=255)
                                                                                                                                                                                                                                       gTCGGCGGTTGTT‑GGTTTGCCTGGCGTCATCAGTCCAGCGACGAAT‑ACATTGATTATTTTGCCGTTT‑CGCTACGCATCATT‑GGCGTTTTGGCGCTCATCCTTACCTCCTGTGGT‑CTGGCGGCAAT‑CAACGCTGACGATATCTGGTATTTTGCCTCCGGTGGCGTCATTGGCAGCTTACTAAGCACTACGCTAC‑AACCACTGCTACACAGTAGCGGGGGAACTATTGCGCTGCTCTGCGTTTGGGCAGCgg                    <  1:478017/250‑1 (MQ=255)
                                                                                                                                                                                                                                                         tGCCTGGCGTCATCAGTCCAGCGACGAAT‑ACATTGATTATTTTGCCGTTT‑CGCTACGCATCATT‑GGCGTTTTGGCGCTCATCCTTACCTCCTGTGGT‑CTGGCGGCAAT‑CAACGCTGACGATATCTGGTATTTTGCCTCCGGTGGCGTCATTGGCAGCTTACTAAGCACTACGCTAC‑AACCACTGCTACACAGTAGCGGGGGAACTATTGCGCTGCTCTGCGTTTGGGCAGCGGGCCTGACGTTGTTCACCg  >  1:114564/1‑251 (MQ=255)
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                               
ACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGCCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTACACCATTCCCGTCATTATTGTCGGCGGTTGTT‑GGTTTGCCTGGCGTCATCAGTCCAGCGACGAAT‑ACATTGATTATTTTGCCGTTT‑CGCTACGCATCATT‑GGCGTTTTGGCGCTCATCCTTACCTCCTGTGGT‑CTGGCGGCAAT‑CAACGCTGACGATATCTGGTATTTTGCCTCCGGTGGCGTCATTGGCAGCTTACTAAGCACTACGCTAC‑AACCACTGCTACACAGTAGCGGGGGAACTATTGCGCTGCTCTGCGTTTGGGCAGCGGGCCTGACGTTGTTCACCG  >  NC_000913/933263‑933761

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