Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,547,958	G→T	M186I (ATG→ATT)	fdnG →	formate dehydrogenase N subunit alpha

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,547,958	0	G	T	94.7%	57.3 / NA	19	M186I (ATG→ATT)	fdnG	formate dehydrogenase N subunit alpha
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (10/8);  total (10/9)
Fisher's exact test for biased strand distribution p-value = 4.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCGCATTAGCTGGGAAGAAGCATTCTCCCGTATTGCGAAGCTGATGAAAGCTGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGCGTAACGGTAAACCGTTGGCTTTCTACCGGTATGCTGTGTGCCTCCGGTGCCAGCAACGAAACCGGGATGCTGACCCAGAAATTTGCCCGCTCCCTCGGGATGCTGGCGGTAGACAACCAGGCGCGCGTCTGACACGGACCAACGGTAGCAAGTCTTGCTCCAACATTTGGTCGCGGTGCGATGACCAACCACTGGGTGGATATCAAAAACGCTAACGTCGTGATGGTGATGGGCGGTAACGCTGCTGAAGCGCATCCCGTCGGTTTCCGCTGGGCGATGGAAGCGAAAAACAACAACGACGCAACCTTGATCGTTGTCGATCCCCGTTTTACGCGTACC  >  NC_000913/1547765‑1548195                                                                                                                                                                                                  |                                                                                                                                                                                                                                              aGCGCATTAGCTGGGAAGAAGCATTCTCCCGTATTGCGAAGCTGCTGAAAGCTGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGCGTAACGGTAAACCGTTGGCTTTCTACCGGTATGCTGTGTGCCTCCGGTGCCAGCAACGAAACCGGGATGCTGACCCAGAAATTTGCCCGCTCCCTCGGGATTCTGGCGGTAGACAACCAGGCGCGCGTCTGACACGGACCAACGGTAGCAAGTCTTGCt                                                                                                                                                                                      >  2:16517/1‑251 (MQ=255)                                                            gACGCTAACTTTATTGAAAAGAACGAGCAGGGCGTAACGGTAAACCGTTGGCTTTCTACCGGTATGCTGTGTGCCTCCGGTGCCAGCAACGAAACCGGGATGCTGACCCAGAAATTTGCCCGCTCCCTCGGGATTCTGGCGGTAGACAACCAGGCGCGCGTCTGACACGGACCAACGGTAGCAAGTCTTGCTCCAACATTTGGTCGCGGTGCGATGACCAACCACTGGGTGGATATCAAAAACGCTAAcgt                                                                                                                           >  1:246342/1‑251 (MQ=255)                                                             aCGCTAACTTTATTGAAAAGAACGAGCAGGGCGTAACGGTAAACCGTTGGCTTTCTACCGGTATGCTGTGTGCCTCCGGTGCCAGCAACGAAACCGGGATGCTGACCCAGAAATTTGCCCGCTCCCTCGGGATTCTGGCGGTAGACAACCAGGCGCGCGTCTGACACGGACCAACGGTAGCAAGTCTTGCTCCAACATTTGGTCGCGGTGCGATGACCAACCACTGGGTGGATATCAAAAACGCTAAcgtc                                                                                                                          <  1:157387/251‑1 (MQ=255)                                                                            aaaaGAACGAGCAGGGCGTAACGGTAAACCGTTGGCTTTCTACCGGTATGCTGTGTGCCTCCGGTGCCAGCAACGAAACCGGGATGCTGACCCAGAAATTTGCCCGCTCCCTCGGGATTCTGGCGGTAGACAACCAGGCGCGCGTCTGACACGGACCAACGGTAGCAAGTCTTGCTCCAACATTTGGTCGCGGTGCGATGACCAACCACTGGGTGGATATCAAAAACGCTAACGTCGTGATGGTGATGGGc                                                                                                           <  1:16517/251‑1 (MQ=255)                                                                             aaaGAACGAGCAGGGCGTAACGGTAAACCGTTGGCTTTCTACCGGTATGCTGTGTGCCTCCGGTGCCAGCAACGAAACCGGGATGCTGACCCAGAAATTTGCCCGCTCCCTCGGGATTCTGGCGGTAGACAACCAGGCGCGCGTCTGACACGGACCAACGGTAGCAAGTCTTGCTCCAACATTTGGTCGCGGTGCGATGACCAACCACTGGGTGGATATCAAAAACGCTAACGTCGTGATGGTGATGGGCg                                                                                                          >  1:26856/1‑251 (MQ=255)                                                                              aaGAACGAGCAGGGCGTAACGGTAAACCGTTGGCTTTCTACCGGTATGCTGTGTGCCTCCGGTGCCAGCAACGAAACCGGGATGCTGACCCAGAAATTTGCCCGCTCCCTCGGGATTCTGGCGGTAGACAACCAGGCGCGCGTCTGACACGGACCAACGGTAGCAAGTCTTGCTCCAACATTTGGTCGCGGTGCGATGACCAACCACTGGGTGGATATCAAAAACGCTAACGTCGTGATGGTGATGGGCgg                                                                                                         >  1:162045/1‑251 (MQ=255)                                                                                     aGCAGGGCGTAACGGTAAACCGCTGGCTTTCTACCGGCATGCTGTGAGCCAACGGTGCCAGCAACGAAACCGGGACGCTGACCCAGCAATTTGCCCGCTCCCTCGGGATTCTGGCGGTAGACAACCAGGCGCGCGTCTGCCACGGAACAACGGGAGCAAGTCCTGCTCCAACATTTGGTCGCGGTGCGATGACCCACCACTAGGTGGATCTCAAACACGCTAACGCCGTGATGGTGATGGGCGGTCACgct                                                                                                  >  2:134005/1‑251 (MQ=37)                                                                                           gCGTAACGGTAAACCGTTGGCTTTCTACCGGTATGCTGTGTGCCTCCGGTGCCAGCAACGAAACCGGGATGCTGACCCAGAAATTTGCCCGCTCCCTCGGGATTCGGGCGGGAGACAACCAGGCGCGCGTCGGACACGGACCAACGGTAGCAAGTCTTGCTCCAACATTTGGTCGCGGTGCGATGACCAACCACTGGGTGGATATCAAAAACGCTAACGTCGTGATGGTGATGGGCGGTAACGCTGCTGaa                                                                                            <  2:199829/251‑1 (MQ=255)                                                                                                  ggTAAACCGTTGGCTTTCTACCGGTATGCTGTGTGCCTCCGGTGCCAGCAACGAAACCGGGATGCTGACCCAGAAATTTGCCCGCTCCCTCGGGATTCTGGCGGTAGACAACCAGGCGCGCGTCTGACACGGACCAACGGTAGCAAGTCTTGCTCCAAAATTTGGTCGCGGTGCGATGAACCACCACTGGGTGGATATCAAAAACGCTAACGTCGTGATGGTGATGGGCGGTAACGCTGCTGAAGCGCATc                                                                                     >  2:245298/1‑251 (MQ=255)                                                                                                                                    gCCTCCGGTGCCAGCAACGAAACCGGGATGCTGACCCAGAAATTTGCCCGCTCCCTCGGGATTCTGGCGGTAGACAACCAGGCGCGCGTCTGACACGGACCAACGGTAGCAAGTCTTGCTCCAACATTTGGTCGCGGTGCGATGACCAACCACTGGGTGGATATCAAAAACGCTAACGTCGTGATGGTGATGGGCGGTAACGCTGCTGAAGCGCATCCCGTCGGTTTCCGCTGGGCGATGGAAGCGaaaa                                                    <  1:9504/250‑1 (MQ=255)                                                                                                                                      cTCCGGTGCCAGCAACGAAACCGGGATGCTGACCCAGAAATTTGCCCGCTCCCTCGGGATTCTGGCGGTAGACAACCAGGCGCGCGTCTGACACGGACCAACGGTAGCAAGTCTTGCTCCAACATTTGGTCGCGGTGCGATGACCAACCACTGGGTGGATATCAAAAACGCTAACGTCGTGATGGTGATGGGCGGTAACGCTGCTGAAGCGCATCCCGTCGGTTTCCGCTGGGCGATGGAAGCGAAAaac                                                  >  1:28079/1‑250 (MQ=255)                                                                                                                                              ccAGCAACGAAACCGGGATGCTGACCCAGAAATTTGCCCGCTCCCTCGGGATTCTGGCGGTAGACAACCAGGCGCGCGTCTGACACGGACCAACGGTAGCAAGTCTTGCTCCAACATTTGGTCGCGGTGCGATGACCAACCACTGGGTGGATATCAAAAACGCTAACGTCGTGATGGTGATGGGCGGTAACGCTGCTGAAGCGCATCCCGTCGGTTTCCGCTGGGCGATGGAAGCGAAAAACAACAacgac                                         >  1:167227/1‑251 (MQ=255)                                                                                                                                                              gATGCTGACCCAGAAATTTGCCCGCTCCCTCGGGATTCTGGCGGTAGACAACCAGGCGCGCGTCTGACACGGACCAACGGTAGCAAGTCTTGCTCCAACATTTGGTCGCGGTGCGATGACCAACCACTGGGTGGATATCAAAAACGCTAACGTCGTGATGGTGATGGGCGGTAACGCTGCTGAAGCGCATCCCGTCGGTTTCCGCTGGGCGATGGAAGCGAAAAACAACAACGACGCAACCTTGATCGTTg                         <  1:266676/251‑1 (MQ=255)                                                                                                                                                                       ccAGAAATTTGCCCGCTCCCTCGGGATTCTGGCGGTAGACAACCAGGCGCGCGTCTGACACGGACCAACGGTAGCAAGTCTTGCTCCAACATTTGGTCGCGGTGCGATGACCAACCACTGGGTGGATATCAAAAACGCTAACGTCGTGATGGTGATGGGCGGTAACGCTGCTGAAGCGCATCCCGGCGGTTTCCGCTGGGCGCTGGAAGCGAAAAAAAACAACGCAGCCACCTTGATCGTTGTCGATcccc                >  2:1550/1‑251 (MQ=25)                                                                                                                                                                        cagtaatgttcccGCTACCTCGGGGTTCTGGGGGTAGACAACCAGGCTCGCGGCTGACACGGACCAACGGTAGCAAGTCTTGCGCGAACATTTGGGCGCGGTGCGATGACCAACCACTGGGTGGATATCAAAATCGCTAACGTCGTGATGGTGATGGGCGGTAACGCTGCTGAAGCGCATCCCGTCGGTTTCCGCTGGTCGATGGAAGCGAAAAACAACAACGACGCAACCTTGTTCGTTGTCGGTCCCCg               <  2:134756/241‑1 (MQ=21)                                                                                                                                                                                     gCTCCCTCGGGATTCTGGCGGTAGACAACCAGGCGCGCGTCTGACACGGACCAACGGTAGCAAGTCTTGCTCCAACATTTGGTCGCGGTGCGATGACCAACCACTGGGTGGATATCAAAAACGCTAACGTCGTGATGGTGATGGGCGGTAACGCTGCTGAAGCGCATCCCGTCGGt                                                                             >  1:295508/1‑176 (MQ=255)                                                                                                                                                                                     gCTCCCTCGGGATTCTGGCGGTAGACAACCAGGCGCGCGTCTGACACGGACCAACGGTAGCAAGTCTTGCTCCAACATTTGGTCGCGGTGCGATGACCAACCACTGGGTGGATATCAAAAACGCTAACGTCGTGATGGTGATGGGCGGTAACGCTGCTGAAGCGCATCCCGTCGGt                                                                             <  2:295508/176‑1 (MQ=255)                                                                                                                                                                                     gCTCCCTCGGGATTCTGGCGGTAGACAACCAGGCGCGCGTCTGACACGGACCAACGGTAGCAAGTCTTGCTCCAACATTTGGTCGCGGTGCGATGACCAACCACTGGGTGGATATCAAAAACGCTAACGTCGTGATGGTGATGGGCGGTAACGCTGCTGAAGCGCATCCCGTCGGTTTCCGCTGGGCGATGGAAGCGAAAAACAACAACGACGCAACCTTGATCGTTGTCGATCCCCGTTTTACGCGTAcc  <  2:246342/251‑1 (MQ=255)                                                                                                                                                                                     gCTCCCTCGGGATGCTGGCGGTAGACAACCAGGCGCGCGTCTGACCCGGACCAACGGTAGCAAGTCGTGCTCCAACATGTGGTCGCGTTGCTATGACCAACCACTGTGTGGATATCAAAAACGCTAACGACGTTATGGTGATGGGCGGTAACGCTGCTGAAGCGCCCCCCGTCGGTTTCCGCTGGGCTATTGCTGCGAAAAACAACAACGACGCCACGTTGATCGGTGTCGATCCCCGCTTTACGCGttcc  <  2:51706/251‑4 (MQ=255)                                                                                                                                                                                                  |                                                                                                                                                                                                                                              AGCGCATTAGCTGGGAAGAAGCATTCTCCCGTATTGCGAAGCTGATGAAAGCTGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGCGTAACGGTAAACCGTTGGCTTTCTACCGGTATGCTGTGTGCCTCCGGTGCCAGCAACGAAACCGGGATGCTGACCCAGAAATTTGCCCGCTCCCTCGGGATGCTGGCGGTAGACAACCAGGCGCGCGTCTGACACGGACCAACGGTAGCAAGTCTTGCTCCAACATTTGGTCGCGGTGCGATGACCAACCACTGGGTGGATATCAAAAACGCTAACGTCGTGATGGTGATGGGCGGTAACGCTGCTGAAGCGCATCCCGTCGGTTTCCGCTGGGCGATGGAAGCGAAAAACAACAACGACGCAACCTTGATCGTTGTCGATCCCCGTTTTACGCGTACC  >  NC_000913/1547765‑1548195

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