Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,256,061	0	C	T	34.5%	32.8 / 10.6	29	E136K (GAA→AAA)	plsB	glycerol‑3‑phosphate 1‑O‑acyltransferase
Reads supporting (aligned to +/- strand):  ref base C (6/13);  new base T (10/0);  total (16/13)
Fisher's exact test for biased strand distribution p-value = 4.14e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AAACAGATCCTGACGAGCAGGCAAACGTGGGCCTACGGCAGCCAGACGTTGACGGGCAAAGTGCATACGCGCCACGCGCGCCAGTTTCTGAGCGATAGTTTTATCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTC‑GCCTTTTTCAC‑GCCCCGGCGCGCG‑ACCAAA‑CATCACCGACACTGGCACCATCTGCACATCCAG‑A‑TTTGGGTTGCT‑ACGGTGCAAATCGAGATAGTCGTGGAA‑CAGCTTAATAGACTCTTCTTTCGGCGTGTAATAGGTGAACACACGCGGCCCGCCGTGAATGAACACATAGCGCGGCAGTAGCGTGCCGTCGATTTCCAGCGGCTCTAACGGGTCAGGCAAGTCATGTGCCAGACACTGGGCGCGCAACGTCAGC  >  NC_000913/4255816‑4256308                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                               aaccAGATCCTGACGAGCAGGCAAACGTGGGCCTACGGCAGCCAGACGTTGACGGGCAAAGTGCATACGCGCCACGCGCGCCAGTTTCTGAGCGATAGTTTTATCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGGTTACTTC‑TCCCTTTTTAC‑Gcc                                                                                                                                                                                                                                                          >  1:109119/4‑251 (MQ=255) aaaCAGATCCTGACGAGCAGGCAAACGTGGGCCTACGGCAGCCAGACGTTGACGGGCAAAGTGCATACGCGCCACGCGCGCCAGTTTCTGAGCGATAGTTTTATCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGGTCACTTT‑GCCCTTTTctcacc                                                                                                                                                                                                                                                           >  1:109121/1‑246 (MQ=255)         ccTGACGAGCAGGCAAACGTGGGCCTACGGCAGCCAGACGTTGACGGGCACAGTGCATACGCGCCACTCGCTCCAGTTTCTGAGCGATAGTTTTATCCGTCCCGAGTTCATCCGCCATACGGCGCAGCGAGACTGACGGCGAGAAACGCACAAAACTGTCACGACCGAGCCACAGTACAGCGAGAAATTTCTGTACGCCGTTAAGCATACTCAGCGGCTGATTTACTTG‑CCCTTTTTTAC‑Gcccccacgac                                                                                                                                                                                                                                                  >  1:82068/1‑245 (MQ=255)                 gCAGGCAAACGTGGGCCTACGGCAGCCAGACGTTGACGGGCAAAGTGCATACGCGCCACGCGCGCCAGTTTCTGAGCGATAGTTTTATCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCTGCGGGTTCACTTC‑GCCTTTTTCAC‑GCCCCGGCGCGCG‑ACCaaa                                                                                                                                                                                                                                         <  1:114992/251‑1 (MQ=255)                           gTGGGCCTACGGCAGCCAGACGTTGACGGGCAAAGTGCATACGCGCCACGCGCGCCAGTTTCTGAGCGATAGTTTTATCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACTGCCGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGATCTCATT‑GCCTTTTTCAc‑ccccccccccgcgccccaaccaaccacac                                                                                                                                                                                                                                >  2:74226/1‑227 (MQ=255)                                     ggCAGCCAGACGTTGACGGGCAAAGTGCATACGCGCCACGCGCGCCAGTTTCTGAGCGATAGTTTTATCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGGTCACTTC‑CCCTTTTTTCC‑GCCCCCGCGCGCGGCCCAAA‑CAtacccacccctggaccc                                                                                                                                                                                                                     >  1:355477/1‑235 (MQ=255)                                      gacgCCAGGCGTTGGCGGGCAAAGGTCATACGCCCCCCGCGCGCCAGTTTCTGAGCGATGGTTTTACCCGTGCCGGTTTCACCCGCCATACGGCGCAGCGAACCTGACGGCGAGAAACGCACAAACCTGCCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTC‑GCCTTTTGCAC‑GCCCCGGCGCGCG‑ACCAAA‑CACCACCGACACTGGCACCAt                                                                                                                                                                                                                   <  2:443598/248‑1 (MQ=255)                                       cggccagccttggACGGGAAAATTGCAAACGCGCCCCGCGCGCCAGTTTCTGAGCGATAGTTTTATCCGTCCCGTGGTCAGCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGCCCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTC‑GCCTTTTTCAC‑GCCCCGGCGCGCG‑ACCAAA‑CATCACCGACACTGGCACCATc                                                                                                                                                                                                                  <  1:97726/239‑1 (MQ=255)                                                gTTGACGGGCAAAGTGCATACGCGCCACGCGCGCCAGTTTCTGAGCGATAGTTTTATCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGGTTA‑TTCTGCCCTTTTTAC‑GCCCCGGCGCGCGGACAAAA‑CA‑CACCGaaccatgcaccaccttgacatcc                                                                                                                                                                                                         >  1:389013/1‑229 (MQ=255)                                                      gggCAAGGTGAATCCGCGCCACGCGCGCCAGTTTCTGGGCGATATTTTTATCCGTGCCGTGTTCTCCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACGGTCGCGCCCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTC‑GCCTTTTTCAC‑GCCCCGGCGCGCG‑ACCAAA‑CATCACCGACACTGGCACCATCTGCACATCCAG‑A‑ttt                                                                                                                                                                                                 <  1:15164/251‑1 (MQ=255)                                                           aaGTGCATCCGCGCCACGCGCGCCAGTTTCTGAGCGATAGTTTTACCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTC‑GCCTTTTTCAC‑GCCCCGGCGCGCG‑ACCAAA‑CATCACCGACACTGGCACCATCTGCACATCCAG‑A‑TTTGGGtt                                                                                                                                                                                            <  1:302683/251‑1 (MQ=255)                                                             gTGCATACGCGCCACGCGCGCCAGTTTCTGAGCGATAGTTTTATCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGGTCACTTC‑GCCTTTTTTTC‑CCCCCCGCGCGCG‑ACCAAA‑CATCACCCACACTGGCACCCCCTGCACATCCAGCT‑TTTGGGTTg                                                                                                                                                                                           >  1:107931/1‑251 (MQ=255)                                                                   ccgcgCCACGCGCGCCAGTTTCTGAGCGATAGTTTTATCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTC‑GCCTTTTTCAC‑GCCCCGGCGCGCG‑ACCAAA‑CATCACCGACACTGGCACCATCTGCACATCCAG‑A‑TTTGGGTTGCT‑ACGGt                                                                                                                                                                                   <  1:94216/250‑1 (MQ=255)                                                                     gcgcCAGGCGCGCAAGTTCCTGAGCGATAGTTTTATCCGTGCCGTGTCCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTC‑GCCTTTTTCAC‑GCCCCGGCGCGCG‑ACCAAA‑CATCACCGACACTGGCACCATCTGCACATCCAG‑A‑TTTGGGTTGCT‑ACGGTGc                                                                                                                                                                                 <  1:54145/251‑1 (MQ=255)                                                                      cgcCACGCGCGCCAGTTTCTGAGCGATAGTTTTATCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGGTCCCTTT‑CCCTTTTTTTC‑CCCCCGGCGCGGG‑GCCAAA‑CACAACAGCCACCGGCCCCACCTGCACATCCAG‑A‑TTTGGGTTGCT‑ACggggc                                                                                                                                                                                 >  1:303657/1‑247 (MQ=255)                                                                        cccggcgcgcCAGTTTCTGGGCGTTAGTTTTATCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGATCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTC‑GCCTTTTTCAC‑GCCCCGGCGCGCT‑ACCAAA‑CATCACCGACACTGGCACCATCTGCGCATCCAG‑A‑TTTGGGTTGCT‑ACGGTGCaaa                                                                                                                                                                              <  1:451068/247‑1 (MQ=255)                                                                              gcgcCAGTTTCTGAGCGATAGTTTTATCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTC‑GCCCTTTTTAC‑GCCCCCGCGCGCG‑CCCAAA‑CAACACCCACACCGGCACCAACTGCACCTCCAG‑AATTTGGGT‑GCT‑ACGGGGCaaaacg                                                                                                                                                                           >  1:224093/1‑245 (MQ=255)                                                                                      ttttGAGCGATAGTTTTATCCGTGCCTGGTTCATCCGCCATACGGCGCAGCGAAATTGACGGCGAGAAACGCACAAAACTGGCGTGACCGAGCCACAGTAAAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCATCGGCGGGTTCACTTC‑GCCTTTTTCAC‑GCCCCGGCGCGCG‑ACCAAA‑CATCACCGACACTGGCACCATCTGCACATCCAG‑A‑TTTGGGTTGCT‑ACGGTGCAAATCGAGATAGTCGt                                                                                                                                                                 <  1:515220/247‑1 (MQ=255)                                                                                                       aTCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGGTCACTTC‑GCCTTTTTCAC‑GCCCCGGCGCGCG‑CCCAAA‑AAACAACGACACCGGCACCCTCTCCCCATCCCG‑A‑TTTGGGTTGCT‑ACGGGGCAAATAGAGATAATCGTGGAA‑CAGCTTAATatacc                                                                                                                                              >  1:294989/1‑247 (MQ=255)                                                                                                           gtgccgtgTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGGTCACTTC‑GCCTTTTTTAC‑GCCCCCGCGCGCG‑GCCAAA‑CATCACCCACACTGGCACCCACTGCACCACCAG‑C‑TTTTGGGTGCT‑TCGGTGCAAAAAGGGATAAGCGTGGAA‑AAACTTAATAGACtcttt                                                                                                                                          >  1:137168/1‑250 (MQ=255)                                                                                                                         cGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGGTCACTTC‑GCCTTTTTTAC‑GCCCCCGCGCGCG‑GCCCAA‑AAACACCCACACCGGCCCCATCTTCCCATCCAC‑A‑GTTGTGTGGCT‑CCCGGGGAAATAGAGAAAGTCGTGGAA‑CAGCTTCATaaaaccttttttctgcggggaa                                                                                                                             >  1:123828/1‑229 (MQ=255)                                                                                                                               aCGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTC‑GCCTTTTTCAC‑GCCCCGGCGCGCG‑ACCAAA‑CATCACCGACACTGGCACCATCTGCACATCCAG‑A‑TTTGGGTTGCT‑ACGGTGCAAATCGAGATAGTCGTGGAA‑CAGCTTAATAGACTCTTCTTTCGGCGTGTAATAGGTg                                                                                                                       <  1:490876/250‑1 (MQ=255)                                                                                                                                 ggCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGGTCACTTT‑GCCTTTTTTTC‑CCCCCCGCGCGCG‑GCCCAA‑AAACACCGCCACCCGCACCATCTTCACATCCCG‑A‑TTTTGGGTGGT‑ACCGTGCAAAACGCGATAATCGTGGGA‑CAGCTTAATAGACTCTTCTTTCTGCGCGTGAAAGGTGaa                                                                                                                     >  1:109213/1‑250 (MQ=255)                                                                                                                                                                              gggcccaaagaaaagcaaaaatTTTCTGTCCGCCGTAAACAATCCGCAGGGGCGGTTTC‑CCTC‑CCCTTTTTCCCCGCCCGGGCGCGC‑‑ACCAAAACACCCCCGACCCTGGCCCCACCTGCACATCCAG‑T‑TTTGGGTTGCT‑ACGGTGCAAACCGAGATAGTCGTGGAA‑CAGCTTAATAGACTCTTCTTCCGGCGTGTAATAGGGGAACACACGCGGCCCGCCGTGAATGAACACATAGCGCGGCAGTAGCGTg                                                                       <  1:388948/234‑1 (MQ=255)                                                                                                                                                                                                                        tACGCAGCGGCGGGTTCACTTC‑GCCTTTTTCAC‑GCCCCGGCGCGCG‑ACCAAA‑CATCACCGACACTGGCACCATCTGCACATCCAG‑A‑TTTGGGTTGCT‑ACGGTGCAAATCGAGATAGTCGTGGAA‑CAGCTTAATAGACtctt                                                                                                                                           <  1:122992/141‑1 (MQ=255)                                                                                                                                                                                                                            cgggtgggggttcccgcc‑GCCTTTTCCAC‑GCCCCGCCGCGCG‑ACCAAG‑CGCCGCCGACTCTGGCACAAGCTCCGCATCAAG‑A‑TTTGG‑TTGCTTGCGGTGCAAAGCGAGAGAGTC‑TGAAAGCGGCTAAATAGACGCTTCTTTCGGCGTGTAATAGGTGAACGGACGCGGCCCGCCGTGAATGAACACATAGCGCGGCAGTAGCGTGGCGTCGATTTCCAGCGGCTCTAACGGGGCAGGCAAGTCATGTGCCa                         <  1:151978/234‑1 (MQ=255)                                                                                                                                                                                                                                     gTTCACTTC‑GCCTTTTTCAC‑GCCCCGGCGCGCG‑ACCAAA‑CATCACCGACACTGGCACCATCTGCACATCCAG‑A‑TTTGGGTTGCT‑ACGGTGCAAATCGAGATAGTCGTGGAA‑CAGCTTAATAGACTCTTCTTTCGGCGTGTAATAGGTGAACACACGCGGCCCGCCGTGAATGAACACATAGCGCGGCAGTAGCGTGCCGTCGATTTCCAGCGGCTCTAACGGGTCAGGCAAGTCATGTGCCAGACACTggg                >  1:128321/1‑251 (MQ=255)                                                                                                                                                                                                                                       tCACTTC‑GCCTTTTTCAC‑GCCCCGGCGCGCG‑ACCAAA‑CATCACCGACACTGGCACCATCTGCACATCCAG‑A‑TTTGGGTTGCT‑ACGGTGCAAATCGAGATAGTCGTGGAA‑CAGCTTAATAGACTCTTCTTTCGGCGTGTAATAGGTGAACACACGCGGCCCGCCGTGAATGAACACATAGCGCGGCAGTAGCGTGCCGTCGATTTCCAGCGGCTCTAACGGGTCAGGCAAGTCATGTGCCAGACACTGGgcg              >  1:396944/1‑251 (MQ=255)                                                                                                                                                                                                                                                ccTTTTTCAC‑GCCCCGGCGCGCG‑ACCAAA‑CATCACCGACACTGGCACCATCTGCACATCCAG‑A‑TTTGGGTTGCT‑ACGGTGCAAATCGAGATAGTCGTGGAA‑CAGCTTAATAGACTCTTCTTTCGGCGTGTAATAGGTGAACACACGCGGCCCGCCGTGAATGAACACATAGCGCGGCAGTAGCGTGCCGTCGATTTCCAGCGGCTCTAACGGGTCAGGCAAGTCATGTGCCAGACACTGGGCGCGCAACGt      >  1:145239/1‑251 (MQ=255)                                                                                                                                                                                                                                                     ttCAC‑GCCCCGGCGCGCG‑ACCAAA‑CATCACCGACACTGGCACCATCTGCACATCCAG‑A‑TTTGGGTTGCT‑ACGGTGCAAATCGAGATAGTCGTGGAA‑CAGCTTAATAGACTCTTCTTTCGGCGTGTAATAGGTGAACACACGCGGCCCGCCGTGAATGAACACATAGCGCGGCAGTAGCGTGCCGTCGATTTCCAGCGGCTCTAACGGGTCAGGCAAGTCATGTGCCAGACACTGGGCGCGCAACGTCAGc  >  1:357158/1‑250 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                               AAACAGATCCTGACGAGCAGGCAAACGTGGGCCTACGGCAGCCAGACGTTGACGGGCAAAGTGCATACGCGCCACGCGCGCCAGTTTCTGAGCGATAGTTTTATCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTC‑GCCTTTTTCAC‑GCCCCGGCGCGCG‑ACCAAA‑CATCACCGACACTGGCACCATCTGCACATCCAG‑A‑TTTGGGTTGCT‑ACGGTGCAAATCGAGATAGTCGTGGAA‑CAGCTTAATAGACTCTTCTTTCGGCGTGTAATAGGTGAACACACGCGGCCCGCCGTGAATGAACACATAGCGCGGCAGTAGCGTGCCGTCGATTTCCAGCGGCTCTAACGGGTCAGGCAAGTCATGTGCCAGACACTGGGCGCGCAACGTCAGC  >  NC_000913/4255816‑4256308

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