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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	71,584	0	T	C	92.6%	77.7 / ‑5.6	27	G78G (GGT→GGC)	yabI	DedA family protein YabI
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (8/17);  minor base G (2/0);  total (10/17)
Fisher's exact test for biased strand distribution p-value = 1.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCCACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAG‑TTTCTGGCACGCCTGGCTGGCAGGGATTATTGGCTGCTTGATGGGCGACT‑GGATTTCTTTCTGGCTGGGTTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCACTACTTGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCC  >  NC_000913/71343‑71819                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                            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)                                                                                                                                                                                                    gggattgtcggttgtctgctGGGAGACT‑GGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCAGCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGATGGGAATTCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAAtattat                                  <  2:24695/232‑1 (MQ=255)                                                                                                                                                                                                     ggattgtcggttgtctgctGGGCGACT‑GGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATc                                 <  2:143371/233‑1 (MQ=255)                                                                                                                                                                                                         tgtcggttgtctgctGGGCGACT‑GGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGc                              <  1:244317/236‑1 (MQ=255)                                                                                                                                                                                                                 gtctgctGGGCGACT‑GGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAACCAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTTGATTTGTCGGTGGCAAAACTCCTTACTCCGAATCTTATTGGctgcccgcct                     >  2:237554/7‑246 (MQ=255)                                                                                                                                                                                                                               t‑GGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGCAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGGGTGGTGGATCTGCCGGTGGCAAAATTCATTACCCCGGATATTATCGTCTGCCTGCTTTGGCCACCGgtttt       >  2:230810/1‑250 (MQ=255)                                                                                                                                                                                                                                     ctccttctGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCCAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGCGGCCACCGTGTTACTTcc  >  2:258067/5‑251 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                            GGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCCACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAG‑TTTCTGGCACGCCTGGCTGGCAGGGATTATTGGCTGCTTGATGGGCGACT‑GGATTTCTTTCTGGCTGGGTTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCACTACTTGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCC  >  NC_000913/71343‑71819

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