Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	71,569	T→C	I73I (ATT→ATC)	yabI →	DedA family protein YabI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	71,569	0	T	C	96.2%	81.6 / ‑6.1	27	I73I (ATT→ATC)	yabI	DedA family protein YabI
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (7/18);  total (9/18)
Fisher's exact test for biased strand distribution p-value = 3.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCCACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAG‑TTTCTGGCACGCCTGGCTGGCAGGGATTATTGGCTGCTTGATGGGCGACT‑GGATTTCTTTCTGGCTGGGTTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCACTACTTGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTA  >  NC_000913/71330‑71814                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                      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)                                                                                                                                                                                                                 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)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                      CGTCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCCACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAG‑TTTCTGGCACGCCTGGCTGGCAGGGATTATTGGCTGCTTGATGGGCGACT‑GGATTTCTTTCTGGCTGGGTTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCACTACTTGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTA  >  NC_000913/71330‑71814

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: ‑