Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	71,548	C→T	G66G (GGC→GGT)	yabI →	DedA family protein YabI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	71,548	0	C	T	91.3%	63.8 / ‑4.0	23	G66G (GGC→GGT)	yabI	DedA family protein YabI
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (5/16);  minor base G (2/0);  total (7/16)
Fisher's exact test for biased strand distribution p-value = 8.30e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

AGAATCCCTGCTTCGTCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCCACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAG‑TTTCTGGCACGCCTGGCTGGCAGGGATTATTGGCTGCTTGATGGGCGACT‑GGATTTCTTTCTGGCTGGGTTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCACTACTTGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAAT  >  NC_000913/71317‑71768                                                                                                                                                                                                                                         |                                                                                                                                                                                                                              aGAATCCCTGCTTCGTCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAG‑TTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACT‑GGa                                                                                                                                                                                                           <  1:94349/251‑1 (MQ=255)              cGTCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCTACCGTGTATTCATTGATTGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAG‑TTTCTGGCACGCCTTGCTGGCAGGGATTGTCGGTTGTCTGCTGGTCGACT‑GGATCTCCTTctggct                                                                                                                                                                                              <  2:263874/251‑1 (MQ=255)                           ggCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTAATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAG‑TTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACT‑GGATCTCCTTCTGGCTGGGCTGGCGtttt                                                                                                                                                                                 <  1:204832/251‑1 (MQ=255)                                  tATGCAAGCATTGCTGGAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAG‑TTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACT‑GGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAg                                                                                                                                                                           <  1:20470/250‑1 (MQ=255)                                  tATGCAAGCATTGCTGGAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCGGGGAGCGCTGATTTGCAGCGGGGAGGTAAGGTTTC‑GGCACGCCTGGCTGGCAGGGAGTGTCGGGTGTCTGATGGGCGAAT‑GGAGATCCCTTTGGGTGGGGTGGGGTTTTAAAAAGc                                                                                                                                                                          >  1:234329/1‑251 (MQ=255)                                        agcgttGCTGGAACACTTTATTACCCAATCTCCCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGTTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAG‑TTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACT‑GGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGc                                                                                                                                                                    <  2:248672/247‑1 (MQ=255)                                        aGCATTGCTGGAACACTTTATTACCCAATCTACCGTGTATTCATTGATGTCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAG‑TTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACT‑GGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGc                                                                                                                                                                    <  2:248679/251‑1 (MQ=255)                                             tGCTGGAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGGTAAG‑TTTCTTGCACGCCTGGCTGGGCGGGGTTTTCGGGTGTCTGCTGGGCGACTTGGTCTTCTT‑CTGGCTGGGGTGGGGTTTTTAAAAACCGTTGCAttgc                                                                                                                                                               >  1:232006/1‑248 (MQ=255)                                                     cacTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGGCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAG‑TTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACT‑GGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCAtt                                                                                                                                                       <  2:259914/251‑1 (MQ=255)                                                             taCCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAG‑TTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACT‑GGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTgaag                                                                                                                                                <  1:259353/250‑1 (MQ=255)                                                                                            gtggtgTTGGGGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAG‑TTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACT‑GGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACa                                                                                                                <  2:160785/251‑1 (MQ=255)                                                                                                    ggtggCCGTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAG‑TTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGTGACT‑GGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTg                                                                                                         <  1:134846/250‑1 (MQ=255)                                                                                                                     cgctggcgctggTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAG‑TTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACT‑GGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGtt                                                                                        <  2:98859/250‑1 (MQ=255)                                                                                                                              tggTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAG‑TTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACT‑GGATCTCCTTCTGGCTGGGCCGGCGGTTTAAAAAGACGTTGCATCGCTGGGCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAAAATGCGTTGCCTCAACACAGCATGTTCACCATTCtt                                                                              >  1:309092/1‑250 (MQ=255)                                                                                                                               ggTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAG‑TTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACT‑GGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCt                                                                               <  1:198597/249‑1 (MQ=255)                                                                                                                                     tttGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAG‑TTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGGGGCT‑GGATCTCCTTCTGGCTGGGCTGGCGGTTTAAAAAGCCGGTGCATTGCTGGTCATTTCTGAAGAAAAACAAATCATTACTCGGTAAAAATGAACCTGGGGTGCGTCAACACAGCATGTTCACCATTCTggtcggtt                                                                       >  1:169595/1‑250 (MQ=255)                                                                                                                                                   gTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAG‑TTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACT‑GGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCg                                                         <  1:301356/251‑1 (MQ=255)                                                                                                                                                     aCGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAG‑TTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACT‑GGATCTCCTTCTGGCTGGGCTGGCGTATTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGAc                                                       <  2:221656/251‑1 (MQ=255)                                                                                                                                                                           gAGCGCTGATTGGCAGCGGCGAGTTAAG‑TTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACT‑GGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAAtggtg                                 >  1:301877/1‑251 (MQ=255)                                                                                                                                                                             gcgcTGATTGGCAGCGGCGAGTTAAG‑TTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACT‑GGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGc                               >  1:197601/1‑251 (MQ=255)                                                                                                                                                                                 tGATTGGCAGCGGCGAGTTAAG‑TTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACT‑GGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCgggg                           >  2:152674/1‑250 (MQ=255)                                                                                                                                                                                                    aaG‑TTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACT‑GGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTg         <  2:309092/250‑1 (MQ=255)                                                                                                                                                                                                             ggCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACT‑GGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAAt  <  1:236100/249‑1 (MQ=255)                                                                                                                                                                                                                                         |                                                                                                                                                                                                                              AGAATCCCTGCTTCGTCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCCACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAG‑TTTCTGGCACGCCTGGCTGGCAGGGATTATTGGCTGCTTGATGGGCGACT‑GGATTTCTTTCTGGCTGGGTTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCACTACTTGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAAT  >  NC_000913/71317‑71768

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