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→GGGpfkA6‑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

GACATGATCAACCGTGGCGGTACGTTCCTCGGTTCTGCGCGTTTCCCGGAATTCCGCGACGAGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCGAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATCGACCGTCTGCGT  >  NC_000913/4107729‑4108016
                                                                                                                                         |                                                                                                                                                          
gACATGATCAACCGTGGCGGTACGTTCCTCGGTTCTGCGCGTTTCCCGGAATTCCGCGACGAGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGt                                                                                                                                               <  1:755369/150‑1 (MQ=255)
        cAACCGTGGCGGTACGTTCCTCGGTTCTGCGCGTTTCCCGGAATTCCGCGACGAGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCg                                                                                                                                       <  1:889712/150‑1 (MQ=255)
         aaCCGTGGCGGTACGTTCCTCGGTTCTGCGCGTTTCCCGGAATTCCGCGACGAGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGAAGGGGCCTACAAGGGGGCaaaggg                                                                                                                                       >  2:40768/1‑145 (MQ=255)
                 cGGTACGTTCCTCGGTTCTGCGCGTTTCCCGGAATTCCGCGACGAGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGAGGAGGGTCCTTCCTGGGGGGAAAGCGTCTGGcc                                                                                                                               >  2:555393/1‑148 (MQ=255)
                      cGTTCCTCGGTTCTGCGCGTTTCCCGGAATTCCGCGACGAGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGGTTCTCCATGGGGGCAATGCGTCTGGCCGAAAAggg                                                                                                                     >  2:638452/1‑151 (MQ=255)
                       gTTCCTCGGTTCTGCGCGTTTCCCGGAATTCCGCGACGAGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGGTGCTTCATGAGTGCAAAAAGGCGGACCGaaaaaggg                                                                                                                    >  1:759649/1‑146 (MQ=255)
                                            cccGGAATTCCGCGACGAGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCGAAATGGGCTTCCCGTGCATCGGTCTGcc                                                                                                <  2:163145/150‑1 (MQ=255)
                                            cccGGAATTCCGCGACGAGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCGAAATGGGCTTCCCGTGCATCGGTCTGCCg                                                                                               <  1:461860/151‑1 (MQ=255)
                                              cGGAATTCCGCGACGAGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGGCGGTTCCTCCATCGGGGCAAAGCGGCTGACCGAAAAGGGGGTCCCCTGCCCCGGGCTGCCggg                                                                                             >  1:767835/1‑151 (MQ=255)
                                               ggAATTCCGCGACGAGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGGCGGGTTCTACATGGGGGGAAAGGGGCTGGCCGAAAAAGGGGTCCCCCGCCTCCGTCTGcc                                                                                                >  1:511578/1‑147 (MQ=255)
                                                gAATTCCGCGACGAGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCGAAATGGGCTTCCCGTGCATCGGTCTGCCGGGc                                                                                            <  1:258693/150‑1 (MQ=255)
                                                     ccGCGACGAGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGGCGGGTGCTTCCTGGGGGGAAAGCGGCTTAACGCAAATGGGTTCCCCTGCATCTGGCGGCCggggccttt                                                                                       >  2:965540/1‑144 (MQ=255)
                                                      cgcgACGAGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCGAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCg                                                                                     <  1:405548/151‑1 (MQ=255)
                                                        cgacgaGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTTACGGGTCCTACCATGGGGGAAAGAGGCTGGACGGAAAAGGGTTCCCCTGCCTCGGGCTGCCGGGCGCTATcccc                                                                                   >  1:298840/1‑148 (MQ=255)
                                                               aaCATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCGAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACa                                                                            <  2:493822/151‑1 (MQ=255)
                                                                  aTCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGGGACGGGTCCTACACCGGGGGAATGGGGCTGACCCGAATAGGGGTCCCCCGCATCGGTTTGGCCGGGGCTATCTCCAAAGACaac                                                                          >  2:352175/1‑148 (MQ=255)
                                                                    ccGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCGAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCaa                                                                        <  2:964291/150‑1 (MQ=255)
                                                                     cgcgCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGGTTCTCCATGGGTGCAATGCGTCTGACCGAAAAGGGGTTCCCCGGCCTCCGTCTGCCGGGGGCCATCGACAACGACAACaaa                                                                       >  2:66623/1‑150 (MQ=255)
                                                                        gcCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCGAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGc                                                                    <  2:40609/150‑1 (MQ=255)
                                                                          cGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCGAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACt                                                                 <  2:32242/151‑1 (MQ=255)
                                                                               cTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGGTCCTACATTGGGGGAAAGCGGCTGACCGAAAAGGGCTTCCCCTGCATCGGGCTGCCGGGGACTATCGACAACGAAATCAAAAGCACTg                                                                >  2:583360/1‑147 (MQ=255)
                                                                                         ccTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCGAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGt                                                  <  1:813546/151‑1 (MQ=255)
                                                                                           tGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGAAGGGTTCTACATTGGGGGAATTCGGCTGAACCAAAAGGGGTTTCCGGGCATCCGGTTGGCGGGGCCTATCGACAACCAAAACAAAAGCACCTGCTACACCATTGGgtt                                                >  2:416879/1‑151 (MQ=255)
                                                                                             aaaaaaCGTGGTATCGACGCGCTGGTGGTTATCGGCGGGGAAGGGTCCCACAAAGGGGCAATGCGGCTGGCCGGAAAGGGGCTCCCCGCCCTCGGTCTGCCCGGGACTAACGACAACGAAACCAAAAGAACAGACTACACAATccgggtgt                                              >  1:217283/1‑144 (MQ=255)
                                                                                                aaaCGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGAAGGGTTCCACAAGGGGGGAATGCGGCTGTCCGGAAAGGGGCTCCCGTGCCCCGGTCTGCCGGGCCCTATCGACAACCACACAAAAGACACCGGCTACACTATAGGGttctt                                             >  2:122042/1‑149 (MQ=255)
                                                                                                   cGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCGAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTg                                        <  1:510189/151‑1 (MQ=255)
                                                                                                     tGGTATCGACGCGCTGGTGGTTATCGGCGGTGGCGGGTCCTAAAATGGGGCAAATCGTCTGACCGGAAATGGGTTCCCCTGCATCGGGCTGCCGGGGGCCATCGACAACAACAACAAAAAAACTGACTACAACATCCGGTTTTTCAc                                          >  2:818679/1‑147 (MQ=255)
                                                                                                     tGGTATCGACGCGCTGGTGGTTATCGGCGGAGACGGGTACTATATAGGGGGAAAGGGTCTGACCGAAAAAGGCGTTCCCTGCATCGGTGTGGCCGGGGCTATAGACCAAGACACCAAAAGGACCGACTACACTATCCGTttcttt                                            >  2:68654/1‑144 (MQ=255)
                                                                                                      ggTATCGACGCGCTGGTGGTTATCGGCGGTGGCGGTTCCTTCATTGGGGGGATGCGTCTTCCCGCAAAGGGGTTCCCCTGCCGCGGTGTGCCGGGGACCATCGGCAACCACATAAAAAGAACAGACTACACTATTGGGTTTTTTTCTgc                                       >  2:702881/1‑149 (MQ=255)
                                                                                                           cGACGCGCTGGTGGTTATCGGCGGTGGCGGGTCCTACCTGGGGGCGATGCGGCCGACCGAAAAGGGGGTTCCCCGCATTGGTCTGCCGGGCACTATCGACAACGACATCAAAAGAACTGACTACACCAACGGGTTTTTTCCTGCGCTGAg                                 >  2:782264/1‑150 (MQ=255)
                                                                                                               gcgcTGGTGGTTATCGGCGGTGACGGGTTCTACACTGGGGCAATGCGTCTTACCGGAAAGGGGGTCCCCGGCAACGGTCTGCCCGGGACCATCTACAAAAACATCAAAAGCACTTCCTACACCATCGGGttttttccctcgctggg                                 >  2:152883/1‑134 (MQ=255)
                                                                                                                 gcTGGTGGTTATCGGCGGTGACGGGTTCTACATTGGGGCGAAGCGGCTGACCGCAAAGGGGGTCCCGTGCCTCGGGCTGCCGGGGGCTATCGACAAAGACATCAAAAGCACTGACTACACTATCTGGTTTTTTACTtcgcgggggacc                             >  2:747869/1‑140 (MQ=255)
                                                                                                                               gcggTGACGGTTCCTACATGGGTGCAATGCGTCTGACCGAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATCg             >  2:878156/1‑150 (MQ=255)
                                                                                                                               gcggTGACGGTTCCTACATGGGTGCAATGCGTCTGACCGAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATCGa            <  1:303598/151‑1 (MQ=255)
                                                                                                                                     aCGGTTCCTACATGGGTGCAATGCGTCTGACCGAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATCGACCGTCt      <  1:429313/151‑1 (MQ=255)
                                                                                                                                         ttCCTACATGGGTGCAATGCGTCTGACCGCAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATCGACCGTCTGCGt  >  2:520062/1‑151 (MQ=255)
                                                                                                                                         |                                                                                                                                                          
GACATGATCAACCGTGGCGGTACGTTCCTCGGTTCTGCGCGTTTCCCGGAATTCCGCGACGAGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCGAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATCGACCGTCTGCGT  >  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: