Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,107,8660TG51.4% 32.4 / 54.6 35G105G (GGT→GGG) pfkA6‑phosphofructokinase I
Reads supporting (aligned to +/- strand):  ref base T (3/14);  new base G (18/0);  total (21/14)
Fisher's exact test for biased strand distribution p-value = 2.93e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.22e-06

GACATGATCAACCGTGGCGGTACGTTCCTCGGTTCTGCGCGTTTCCCGGAATTCCGCGACGAGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTA‑CATGGGTGCAATGCGTCT‑GACC‑G‑AAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATCGACCGTCTGCGT  >  NC_000913/4107729‑4108016
                                                                                                                                         |                                                                                                                                                          
gACATGATCAACCGTGGCGGTACGTTCCTCGGTTCTGCGCGTTTCCCGGAATTCCGCGACGAGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTA‑CATGGGt                                                                                                                                               <  1:755369/150‑1 (MQ=255)
        cAACCGTGGCGGTACGTTCCTCGGTTCTGCGCGTTTCCCGGAATTCCGCGACGAGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTA‑CATGGGTGCAATGCg                                                                                                                                       <  1:889712/150‑1 (MQ=255)
         aaCCGTGGCGGTACGTTCCTCGGTTCTGCGCGTTTCCCGGAATTCCGCGACGAGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGAAGGGGCCTA‑CAAGGGGGCaaaggg                                                                                                                                       >  2:40768/1‑145 (MQ=255)
                 cGGTACGTTCCTCGGTTCTGCGCGTTTCCCGGAATTCCGCGACGAGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGAGGAGGGTCCTT‑CCTGGGGGGAAAGCGTCT‑GGcc                                                                                                                               >  2:555393/1‑148 (MQ=255)
                      cGTTCCTCGGTTCTGCGCGTTTCCCGGAATTCCGCGACGAGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGGTTCTC‑CATGGGGGCAATGCGTCT‑GGCC‑G‑AAAAggg                                                                                                                     >  2:638452/1‑151 (MQ=255)
                       gTTCCTCGGTTCTGCGCGTTTCCCGGAATTCCGCGACGAGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGGTGCTT‑CATGAGTGCAAAAAGGCG‑GACC‑G‑aaaaaggg                                                                                                                    >  1:759649/1‑146 (MQ=255)
                                            cccGGAATTCCGCGACGAGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTA‑CATGGGTGCAATGCGTCT‑GACC‑G‑AAATGGGCTTCCCGTGCATCGGTCTGcc                                                                                                <  2:163145/150‑1 (MQ=255)
                                            cccGGAATTCCGCGACGAGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTA‑CATGGGTGCAATGCGTCT‑GACC‑G‑AAATGGGCTTCCCGTGCATCGGTCTGCCg                                                                                               <  1:461860/151‑1 (MQ=255)
                                              cGGAATTCCGCGACGAGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGGCGGTTCCTC‑CATCGGGGCAAAGCGGCT‑GACC‑G‑AAAAGGGGGTCCCCTGCCCCGGGCTGCCggg                                                                                             >  1:767835/1‑151 (MQ=255)
                                               ggAATTCCGCGACGAGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGGCGGGTTCTA‑CATGGGGGGAAAGGGGCT‑GGCC‑G‑AAAAAGGGGTCCCCCGCCTCCGTCTGcc                                                                                                >  1:511578/1‑147 (MQ=255)
                                                gAATTCCGCGACGAGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTA‑CATGGGTGCAATGCGTCT‑GACC‑G‑AAATGGGCTTCCCGTGCATCGGTCTGCCGGGc                                                                                            <  1:258693/150‑1 (MQ=255)
                                                     ccGCGACGAGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGGCGGGTGCTT‑CCTGGGGGGAAAGCGGCT‑TAAC‑GCAAATGGG‑TTCCCCTGCATCTGGCGGCCggggccttt                                                                                       >  2:965540/1‑144 (MQ=255)
                                                      cgcgACGAGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTA‑CATGGGTGCAATGCGTCT‑GACC‑G‑AAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCg                                                                                     <  1:405548/151‑1 (MQ=255)
                                                        cgacgaGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTTACGGGTCCTACCATGGG‑GGAAAGAGGCTGGACG‑G‑AAAAGGG‑TTCCCCTGCCTCGGGCTGCCGGGCGCTATcccc                                                                                   >  1:298840/1‑148 (MQ=255)
                                                               aaCATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTA‑CATGGGTGCAATGCGTCT‑GACC‑G‑AAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACa                                                                            <  2:493822/151‑1 (MQ=255)
                                                                  aTCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGGGACGGGTCCTA‑CACCGGGGGAATGGGGCT‑GACC‑C‑GAATAGGGGTCCCCCGCATCGGTTTGGCCGGGGCTATCTCCAAAGACaac                                                                          >  2:352175/1‑148 (MQ=255)
                                                                    ccGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTA‑CATGGGTGCAATGCGTCT‑GACC‑G‑AAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCaa                                                                        <  2:964291/150‑1 (MQ=255)
                                                                     cgcgCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGGTTCTC‑CATGGGTGCAATGCGTCT‑GACC‑G‑AAAAGGGGTTCCCCGGCCTCCGTCTGCCGGGGGCCATCGACAACGACAACaaa                                                                       >  2:66623/1‑150 (MQ=255)
                                                                        gcCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTA‑CATGGGTGCAATGCGTCT‑GACC‑G‑AAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGc                                                                    <  2:40609/150‑1 (MQ=255)
                                                                          cGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTA‑CATGGGTGCAATGCGTCT‑GACC‑G‑AAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACt                                                                 <  2:32242/151‑1 (MQ=255)
                                                                               cTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGGTCCTA‑CATTGGGGGAAAGCGGCT‑GACC‑G‑AAAAGGGCTTCCCCTGCATCGGGCTGCCGGGGACTATCGACAACGAAATCAAAAGCACTg                                                                >  2:583360/1‑147 (MQ=255)
                                                                                         ccTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTA‑CATGGGTGCAATGCGTCT‑GACC‑G‑AAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGt                                                  <  1:813546/151‑1 (MQ=255)
                                                                                           tGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGAAGGGTTCTA‑CATTGGGGGAATTCGGCT‑GAAC‑C‑AAAAGGGGTTTCCGGGCATCCGGTTGGCGGGGCCTATCGACAACCAAAACAAAAGCACCTGCTACACCATTGGgtt                                                >  2:416879/1‑151 (MQ=255)
                                                                                             aaaaaaCGTGGTATCGACGCGCTGGTGGTTATCGGCGGGGAAGGGTCCCA‑CAAAGGGGCAATGCGGCT‑GGCCGG‑AAAGGGGCTCCCCG‑CCCTCGGTCTGCCCGGGACTAACGACAACGAAACCAAAAGAACAGACTACACAATccgggtgt                                              >  1:217283/1‑144 (MQ=255)
                                                                                                aaaCGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGAAGGGTTCCA‑CAAGGGGGGAATGCGGCT‑GTCC‑G‑GAAAGGGGCTCCCGTGCCCCGGTCTGCCGGGCCCTATCGACAACCACACAAAAGACACCGGCTACACTATAGGGttctt                                             >  2:122042/1‑149 (MQ=255)
                                                                                                   cGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTA‑CATGGGTGCAATGCGTCT‑GACC‑G‑AAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTg                                        <  1:510189/151‑1 (MQ=255)
                                                                                                     tGGTATCGACGCGCTGGTGGTTATCGGCGGTGGCGGGTCCTA‑AAATGGGGCAAATCGTCT‑GACCGG‑AAATGGG‑TTCCCCTGCATCGGGCTGCCGGGGGCCATCGACAACAACAACAAAAAAACTGACTACAACATCCGGTTTTTCAc                                          >  2:818679/1‑147 (MQ=255)
                                                                                                     tGGTATCGACGCGCTGGTGGTTATCGGCGGAGACGGGTACTA‑TATAGGGGGAAAGGGTCT‑GACC‑G‑AAAAAGGCGTTCCCTGCATCGGTGTGGCCGGGGCTATAGACCAAGACACCAAAAGGACCGACTACACTATCCGTttcttt                                            >  2:68654/1‑144 (MQ=255)
                                                                                                      ggTATCGACGCGCTGGTGGTTATCGGCGGTGGCGGTTCCTT‑CATTGGGGGGATGCGTCT‑TCCC‑G‑CAAAGGGGTTCCCCTGCCGCGGTGTGCCGGGGACCATCGGCAACCACATAAAAAGAACAGACTACACTATTGGGTTTTTTTCTgc                                       >  2:702881/1‑149 (MQ=255)
                                                                                                           cGACGCGCTGGTGGTTATCGGCGGTGGCGGGTCCTA‑CCTGGGGGCGATGCGGCC‑GACC‑G‑AAAAGGGGGTTCCCCGCATTGGTCTGCCGGGCACTATCGACAACGACATCAAAAGAACTGACTACACCAACGGGTTTTTTCCTGCGCTGAg                                 >  2:782264/1‑150 (MQ=255)
                                                                                                               gcgcTGGTGGTTATCGGCGGTGACGGGTTCTA‑CACTGGGGCAATGCGTCT‑TACC‑G‑GAAAGGGGGTCCCCGGCAACGGTCTGCCCGGGACCATCTACAAAAACATCAAAAGCACTTCCTACACCATCGGGttttttccctcgctggg                                 >  2:152883/1‑134 (MQ=255)
                                                                                                                 gcTGGTGGTTATCGGCGGTGACGGGTTCTA‑CATTGGGGCGAAGCGGCT‑GACC‑G‑CAAAGGGGGTCCCGTGCCTCGGGCTGCCGGGGGCTATCGACAAAGACATCAAAAGCACTGACTACACTATCTGGTTTTTTACTtcgcgggggacc                             >  2:747869/1‑140 (MQ=255)
                                                                                                                               gcggTGACGGTTCCTA‑CATGGGTGCAATGCGTCT‑GACC‑G‑AAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATCg             >  2:878156/1‑150 (MQ=255)
                                                                                                                               gcggTGACGGTTCCTA‑CATGGGTGCAATGCGTCT‑GACC‑G‑AAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATCGa            <  1:303598/151‑1 (MQ=255)
                                                                                                                                     aCGGTTCCTA‑CATGGGTGCAATGCGTCT‑GACC‑G‑AAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATCGACCGTCt      <  1:429313/151‑1 (MQ=255)
                                                                                                                                         ttCCTA‑CATGGGTGCAATGCGTCT‑GACC‑G‑CAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATCGACCGTCTGCGt  >  2:520062/1‑151 (MQ=255)
                                                                                                                                         |                                                                                                                                                          
GACATGATCAACCGTGGCGGTACGTTCCTCGGTTCTGCGCGTTTCCCGGAATTCCGCGACGAGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTA‑CATGGGTGCAATGCGTCT‑GACC‑G‑AAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATCGACCGTCTGCGT  >  NC_000913/4107729‑4108016

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