Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,718,6410GA37.5% 10.7 / 10.6 17A321V (GCC→GTC) anmKanhydro‑N‑acetylmuramic acid kinase
Reads supporting (aligned to +/- strand):  ref base G (6/4);  new base A (2/4);  total (9/8)
Fisher's exact test for biased strand distribution p-value = 6.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.75e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CCCCCAGTACCGTCTCCTGGCTTGCGCCAGTGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGGCATCGGTGGTGGTGACTTCTGTGCCTGGCAGTAATGCCGCCAGACGCGCCATGAGTAGCGGATTACGACTACCTCCACCACATACCATCAATCGTTCGCAGCCACCGCTCAACAAAACTTGTTCAGAAATGGTCA  >  NC_000913/1718516‑1718776
                                                                                                                             |                                                                                                                                       
cccccAGTACCGTCTCCTGGCTTGCGCCAGTGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGGCATCggtggtggt                                                                                                                            <  1:217334/139‑1 (MQ=255)
                 tGGCTTGCGCCAGTGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGACATCGGTGGTGGTGACTTCTGTGCCTGGCa                                                                                                           >  2:523233/1‑139 (MQ=255)
                                             ttCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGGCATCGGTGGTGGTGACTTCTGTGCCTGGCAGTAATGCCGCCAGACGCGCCATGAGTAg                                                                               <  2:292132/139‑1 (MQ=255)
                                                         ccGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGGCATCGGTGGTGGTGACTTCTGTGCCTGGCAGTAATGCCGCCAGACGCGCCATGAGTAGCGGATTACGACt                                                                   <  1:605843/139‑1 (MQ=255)
                                                         ccGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGACATCGGTGGTGGTGACTTCTGTGCCTGGCAGTAATGCCGCCAGACGCGCCATGAGTAGCGGATTACGACt                                                                   <  1:124336/139‑1 (MQ=255)
                                                                  tCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGGCATCGGTGGTGGTGACTTCTGTGCCTGGCAGTAATGCCGCCAGACGCGCCATGAGTAGCGGATTACGACTACCTccacc                                                          >  2:240550/1‑139 (MQ=255)
                                                                  tCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGCCATCGGTGGTGGTGACTTCTGTGCCTGGCAGTAATGCCGCCAGACGCGCCATGAGTAGCGGATTACGACt                                                                   >  2:347828/1‑130 (MQ=255)
                                                                  tCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGACATCGGTGGTGGTGACTTCTGTGCCTGGCAGTAATGCCGCCAGACGCGCCATGAGTAGCGGATTACGACt                                                                   <  1:347828/130‑1 (MQ=255)
                                                                   ccGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGGCATCGGTGGTGGTGACTTCTGTGCCTGGCAGTAATGCCGCCAGACGCGCCATGAGTAGCGGATTACGACTACCTccacca                                                         >  2:186707/1‑139 (MQ=255)
                                                                    cGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGACATCGGTGGTGGTGACTTCTGTGCCTGGCAGTAATGCCGCCAGACGCGCCATGAGTAGCGGATTACGACTACCTCcaccac                                                        <  1:326246/139‑1 (MQ=255)
                                                                             aGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGGCATCGGTGGTGGTGACTTCTGTGCCTGGCAGTAATGCCGCCAGACGCGCCATGAGTAGCGGATTACGACTACCTCCACCACATACCATCa                                               >  2:147990/1‑139 (MQ=255)
                                                                                  ggCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGGCATCGGTGGTGGTGACTTCTGTGCCTGGCAGTAATGCCGCCAGACGCGCCATGAGTAGCGGATTACGACTACCTCCACCACATACCATCAATCGt                                          >  1:544727/1‑139 (MQ=255)
                                                                                   gCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGGCATCGGTGGTGGTGACTTCTGTGCCTGGCAGTAATGCCGCCAGACGCGCCATGAGTAGCGGATTACGACTACCTCCACCACATACCATCAATCGtt                                         <  1:147990/139‑1 (MQ=255)
                                                                                                     aTGTCATCGCCACTAATGCCAACGACATCGGTGGTGGTGACTTCTGTGCCTGGCAGTAATGCCGCCAGACGCGCCATGAGTAGCGGATTACGACTACCTCCACCACATACCATCAATCGTTCGCAGCCACCGCTcaaca                       >  2:358591/1‑139 (MQ=255)
                                                                                                          aTCGCCACTAATGCCAACGGCATCGGTGGTGGTGACTTCTGTGCCTGGCAGTAATGCCGCCAGACGCACCATGAGTAGCGGATTACGACTACCTCCACCACATACCATCAATCGTTCGCAGCCACCGCTCAACAAAACt                  >  1:152071/1‑139 (MQ=255)
                                                                                                                         aaCGGCATCGGTGGTGGTGACTTCTGTGCCTGGCAGTAATGCCGCCAGACGCGCCATGAGTAGCGGATTACGACTACCTCCACCACATACCATCAATCGTTCGCAGCCACCGCTCAACAAAACTTGTTCAGAAATGGTc   >  1:191672/1‑139 (MQ=255)
                                                                                                                          aCGACATCGGTGGTGGTGACTTCTGTGCCTGGCAGTAATGCCGCCAGACGCGCCATGAGTAGCGGATTAAGACTACCTCCACCACATACCATCAATCGTTCGCAGCCACCGCTCAACAAAACTTGTTCAGAAATTGTCa  <  1:358591/139‑1 (MQ=255)
                                                                                                                             |                                                                                                                                       
CCCCCAGTACCGTCTCCTGGCTTGCGCCAGTGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGGCATCGGTGGTGGTGACTTCTGTGCCTGGCAGTAATGCCGCCAGACGCGCCATGAGTAGCGGATTACGACTACCTCCACCACATACCATCAATCGTTCGCAGCCACCGCTCAACAAAACTTGTTCAGAAATGGTCA  >  NC_000913/1718516‑1718776

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑