Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,725,576	C→A	M19I (ATG→ATT)	ydhF ←	putative oxidoreductase YdhF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,725,576	0	C	A	92.6%	82.4 / ‑5.7	27	M19I (ATG→ATT)	ydhF	putative oxidoreductase YdhF
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base A (10/15);  total (12/15)
Fisher's exact test for biased strand distribution p-value = 1.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGTAATGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATACCGCATTTACTGACGATTTCCATCCGTTCACGCAGGTGAGGTGCCAGTTTCAGTGCCTCGCCAAACGCCGCTTCGCACTGATAGCCACCATAAATATCAGCATGGTCCACGGTGGTCACGCCGAGATCCAGATGCTCTTCAATAAAACTGACCAGCTGGCGGGCGGACATATTCCAGTCCATCAATCGCCAGTAGCCCA‑TCACAAAACGGGAAAACTCCGGGCCTTGCGGCGCAATAGTAATACGCTGAACCATAATCGCTTCCTCTTATCAGATATGAGAGGAGTATACGCAAGATTAGGTTCAAAAGAGTGATGGTTGCTCCGGTTCGTCTGATGACGCTGGCTTATTTGCGCGTAATTTGCGCATTAATCGCTGCCGACAAAGGCGCAGCACCTCTTGTTTTTCGCCATCGCTCATTTTATTCCAGTTAAAACG  >  NC_000913/1725334‑1725815                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                 aTGTAATGACCAATGCCGTTTTCTTCACGCGCGGTCGTCGCGATACCGCATTTACTGACGATTTCCATCCGTTCACGCAGGTGAGGTGCCAGTTTCAGTGCCTCGCCAAACGCCGCTTCGCACTGATAGCCACCATAAATATCAGCATGGTCCACGGTGGTCACGCCGAGATCCAGATGCTCTTCAATAAAACTGACCAGCTGGCGGGCGGACATATTCCAGTACATCAATCACCAGTAGCCAA‑ACACaa                                                                                                                                                                                                                                          >  2:85241/1‑250 (MQ=255)     aaTGACCAATGACGTATTCTGCACGCGCGGTCGTCGCGATACCGCATTTACTGACGATTTCCATCCGTTCACGCAGGTGAGGTGCCAGTTTCAGTGCCTCGCCAAACGCCGCTTCGCACTGATAGCCACCATAAATATCAGCATGGTCCACGGTGGGCACGCCGAGATCCAGATGCTCTTCAATAAAACTGACCAGCTGGCGGGCGGACATATTCCAGTCCATCAATCGCCAGTAGCCAA‑TCACAAAACgg                                                                                                                                                                                                                                     <  2:157111/251‑1 (MQ=255)           cAATGACGTTTTCTTCACGCGCGGTCGTCGCGATACCGCATTTACTGACGATTTCCATCCGTTCACGCAGGTGAGGTGCCAGTTTCAGTGCCTCGCCAAACGCCGCTTCGCACTGATAGCCTCCATAAATATCAGCATGGTCCACGGTGGTCACGCCGAGATCCAGATGCTCTTCAATAAAACTGACCAGCTGGCGGGCGGACATATTCCAGTCCATCAATCGCCAGTAGCCAA‑TCACAATACGGGAAAAc                                                                                                                                                                                                                               <  2:95582/251‑1 (MQ=255)                                cgGTCGTCGCGATACCGCATTTACTGGCGATTTCCATCCGTTCACGCAGGTGAGGTGCCAGTTTCAGTGCCTCGCCAAACGCCGCTTCGCACTGATAGCCACCAGAAATATCAGCATGGTCCACGGTGGTCACGCCGAGATCCAGATGCTCTTCAATAAAACTGACCAGCTGGCGGGCGGACATATTCCAGTCCATCAATCGCCAGTAGCCAA‑TCACAAAACGGGAAAACTCCGGGCCTTGCGGCGCAATa                                                                                                                                                                                                          <  2:154328/251‑1 (MQ=255)                                    cgtcgCGATACCGCATTTACTGACGCTTTCCATCCGTTCACGCAGGGGAGGTGCCAGTTTCAGTGCCTCGCCAAACTCCGCTTCGCACTGATAGCCACCATAAATATCAGCCTGGTCCACGGTGGTCACGCCGAGATCCAGATGCTCTTCAATACAACTGACCAGCTGGCGGTCGGCCATTTTCCAGTCCATCAATCGCCCGTA‑CCCACTCCCAAAACGGGACACCTCCGGTCCTTgccccgcactacta                                                                                                                                                                                                       >  2:101875/1‑238 (MQ=255)                                      tcgCGATACCGCATTTACTGACGATTTCCATCCGTTCACGCAGGTGAGGTGCCAGTTTCAGTGCCTCGCCAAACGCCGCTTCGCACTGATAGCCACCATAAATATCAGCATGGTCCACGGTGGTTACGCCGAGATCCAGATGCTCTTCAATAAAACTGACCAGCTGGCGGGCGGACATATTCCAGTCCATCAATCGCCAGTAGCCAA‑TCACAAAACGGGAAAACTCCGGGCCTTGCGGCGCAATAGTAAt                                                                                                                                                                                                     >  2:168806/1‑250 (MQ=255)                                        gcgATACCGCATTTACTGACGATTTCCATCCGTTCACGCAGGTGAGGTGCCAGTTTCAGTGCCTCGCCAAACGCCGCTTCGCACTGATAGCCACCATAAATATCAGCATGGTCCACGGTGGTCACGCCGAGATCCAGATGCTCTTCAATAAAACTGACCAGCTGGCGGGCGGACATATTCCAGTCCATCAATCGCCAGTAGCCAA‑TCACAAAACGGGAAAACTCCGGGCCTTGCGGCGCAATAGTAATACg                                                                                                                                                                                                  <  1:126743/251‑1 (MQ=255)                                                    ttACTGACGATTTCCATCCGTTCACGCAGGTGAGGTGCCAGTTTCAGTGCCTCGCCAAACGCCGCTTCGCACTGATAGCCACCATAAATATCAGCATGGTCCACGGTGGTCACGCCGAGATCCAGATGCTCTTCAATAAAACTGACCAGCTGGCGGGCGGACATATTCCAGTCCATCAATCGCCAGTAGCCAA‑TCACAAAACGGGAAAACTCCGGGCCTTGCGGCGCAATAGTAATACGCTGAACCATAAt                                                                                                                                                                                      >  1:190645/1‑251 (MQ=255)                                                      aCTGACGATTTCCATCCGTTCACGCAGGTGAGGTGCCAGTTTCAGTGCCTCGCCAAACGCCGCTTCGCACTGATAGCCACCATAAATATCAGCATGGTCCACGGTGGTCACGCCGAGATCCAGATGCTCTTCAATAAAACTGACCAGCTGGCGGGCGGACATATTCCAGTCCATCAATCGCCAGTAGCCAA‑TCACAAAACGGGAAAACTCCGGGCCTTGCGGCGCAATAGTAATACGCTGAACCATAATCg                                                                                                                                                                                    >  1:74342/1‑251 (MQ=255)                                                                            cGCAGGTTAGGTGCCAGTTTTAGTGACTCGCCAAACACCGCTTCGCACTGATAGCCACCGTAAATATCAGCATGGTCCACGGGGGTCACGCCGAGATCCAGATGCTCTTCAATAGAACTGACCAGCTGGCGGGCGGACATATTCCAGGCCATCAATCGCCAGTAGCCCA‑TCAGCAACCGGGAAAACTCCAGGCCGTGAGGCGACAGCGTTATACGCTGGAGCATCATCGCTTTCTCTTATCCGAGATgagg                                                                                                                                                              >  2:72340/1‑250 (MQ=255)                                                                            cGCAGGTGAGGTGCCAGTTTCAGTGCCTCGCCAAACGCCGCTTCGCACTGATAGCCACCATAAATATCAGCATGGTCCACGGTGGTCACGCCGAGATCCAGATGCGCTTCAATAAAACTGACCAGCTGGCGGGCGGACATATTCCAGTCCATCAATCGCCAGTAGCCAA‑TCACAAAACGGGAAAACTCCGGGCCTTGCGGCGCAATAGTAATACGCGGAACCATAATCGCTTCCTCTTATCAGATATgaga                                                                                                                                                              <  2:71402/251‑1 (MQ=255)                                                                                          cAGTTTCAGTGCCTCGCCAAACGCCGCTTCGCACTGATAGCCACCATAAATATCAGCATGGTCCACGGTGGTCACGCCGAGATCCAGATGCTCTTCAATAAAACTGACCAGCTGGCGGGCGGACATATTCCAGTCCATCAATCGCCAGTAGCCAA‑TCACAAAACGGGAAAACTCCGGGCCTTGCGGCGCAATAGTAATACGCTGAACCATAATCGCTTCCTCTTATCAGATATGAGAGGAGTATACGCAAg                                                                                                                                                <  1:95935/251‑1 (MQ=255)                                                                                              ttCAGTGCCTCGCCAAACGCCGCTTCGCACTGATAGCCACCATAAATATCAGCATGGTCCACGGTGGTCACGCCGAGATCCAGATGCTCTTCAATAAAACTGACCAGCTGGCGGGCGGACATATTCCAGTCCATCAATCGCCAGTAGCCAA‑TCACAAAAAGGGAAAACTACGGGCCTTGCGGCGCAATAGTAATAAGCTGAACCATAATCGCTTCCTCTTATCAGATATGAGAGGCGTATACGCAAGAtt                                                                                                                                             >  2:205256/1‑250 (MQ=255)                                                                                                  gggCCTGGCCAAACGCCTCTTCGCACTTATAGCCACCATAAATATCAGCATGGTCCACGGTGGTCACGCCGAGATCCAGATGCGCTTCAATAAAACTGACCAGCTGGCGGGCGGACATAGTCCAGTCCATCAATCGCCAGGAGCCAA‑TCACAAAACGGGAAAACTCCGGGCCTTGCGGCGCAAGATTAATACGCTGAACCATAATCGCTTCCTCTTATCAGAGATGAGAGGAGTGTACGCAAGATTAGGtt                                                                                                                                        <  2:114189/249‑1 (MQ=255)                                                                                                              ccgccgcGTCGCACTGATAGCCACCATAAATATCAGCATGGTCCACGGTGGTCACGCCGAGATCCAGATGCTCTTCAATAAAACTGACCAACTGGCGGGCGGACATATTCCAGTCCATCAATCGCCAGTAGCCAA‑TCACAAAACGGGAAAACTCCGGGCCTTGCGGCGCAATAGTAATACGCTGAACCATAATCGCTTCCTCTTATCAGATATGAGAGGAGTATACGCAAGATTAGGTTCAAAAGAGTGAt                                                                                                                            <  1:78305/250‑1 (MQ=255)                                                                                                                   gcTTCGCACTGATAGCCACCATAAATATCAGCATGGTCCACGGTGGTCACGCCGAGATCCAGATGCTCTTCAATAAAACTGACCAGCTGGCGGGCGGACATATTCCAGTCCATCAATCGCCAGTAGCCAA‑TCACAAAACGGGAAAACGCCGGGCCTTGCGGCGCAATAGTAATACGCTGAACCATAATCGCTTCCTCTTATCAGATATGAGAGGAGTATACTCAAGATTATGTTCCAAAGAGTGATGGTTg                                                                                                                       <  2:79350/251‑1 (MQ=255)                                                                                                                        gCACTGATAGCCACCATAAATATCAGCATGGTCCACGGTAGTCACGCCGAGATCCAGATGATCTTCAATAAAACTGACCAGCTGGCGGGCGGACATATTCAAGTCCATCAATCGCCAGTAGCCAA‑TCACAAAACGGGAAAACTCCGGGCCTTGCGGCGCAATAGTAATACGCTGAACCATAATAGCTTCCTCTTATCAGATATGAGAGGAGTATACGCAAGATTAGGTTCAAAAGAGTGATGGTTGCTCCg                                                                                                                  >  1:137197/1‑251 (MQ=255)                                                                                                                         cACTGATAGCCACCATAAATATCAGCATGGTCCACGGTGGTCACGCCGAGATCCAGATGCTCTTCAATAAAACTGACCAGCTGGCGGGCGGACATATTCCAGTCCATCAATCGCCAGTAGCCAA‑TCACAAAACGGGAAAACTCCGGGCCTTGCGGCGCAATAGTAATACGCTGAACCATAATCGCTTCCTCTTATCAGATATGAGAGGAGTATACGCAAGATTAGGTTCAAAAGAGTGATGGTTGCTCCgg                                                                                                                 >  1:142797/1‑251 (MQ=255)                                                                                                                                                     ggTCCACGGTGGTCACGCCGAGATCCAGATGCTCTTCAATAAAACTGACCAGCTGGCGGGCGGACATATTCCAGTCCATCAATCGCCAGTAGCCAA‑TCACAAAACGGGAAAACTCCGGGCCTTGCGGCGCAATAGTAATACGCTGAACCATAATCGCTTCCTCTTATCAGATATGAGAGGAGTATACGCAAGATTAGGTTCAAAAGAGTGATGGTTGCTCCGGTTCGTCTGATGACGCTGGCTTATTTgcg                                                                                     <  2:190645/251‑1 (MQ=255)                                                                                                                                                     ggTCCACGGTGGTCACGCCGAGATCCAGATGCTCTTCAATAAAACTGACCAGCTGGCGGGCGGACATATTCCAGTCCATCAATCGCCAGTAGCCAA‑TCACAAAACGGGAAAACTCCGGGCCTTGCGGCGCAATAGTAATACGCTGAACCATAATCGCTTCCTCTTATCAGATATGAGAGGAGTATACGCAAGATTAGGTTCAAAAGAGTGATGGTTGCTCCGGTTCGTCTGATGACGCTGGCTTATTTgcg                                                                                     <  2:215445/251‑1 (MQ=255)                                                                                                                                                                             ccAGATGCTCTTCAATAAAACTGACCAGCTGGCGGGCGGACATATTCCAGTCCATCAATCGCCAGTAGCCAA‑TCACAAAACGGGAAAACTCCGGGCCTTGCGGCGCAATAGTAATACGCTGAACCATAATCGCTTCCTCTTATCAGATATGAGAGGAGTATACGCAAGATTAGGTTCAAAAGAGTGATGGTTGCTCCGGTTCGTCTGATGACGCTGGCTTATTTGCGCGTAATTTGCGCATTAATCGCTg                                                              <  2:72761/250‑1 (MQ=255)                                                                                                                                                                                      ctTCAATAAAACTGACCAGCTGGCGGGCGGACATATTCCAGTCCATCAATCGCCAGTAGCCAA‑TCACAAAACGGGAAAACTCCGGGCCTTGCGGCGCAATAGTAATACGCTGAACCATAATCGCTTCCTCTTATCAGATATGAGAGGAGTATACGCAAGATTAGGTTCAAAAGAGTGATGGTTGCTCCGGTTCGTCTGATGACGCTGGCTTATTTGCGCGTAATTTGCGCATTAATCGCTGCCGACaaa                                                      <  1:47338/249‑1 (MQ=255)                                                                                                                                                                                              aaaCTGACCAGCTGGCGGGCGGACATATTCCAGTCCATCAATCGCCAGTAGCCAA‑TCACAAAACGGGAAAACTCCGGGCCTTGCGGCGCAATAGTAATACGCTGAACCATAATCGCTTCCTCTTATCAGATATGAGAGGAGTATACGCAAGATTAGGTTCAAAAGAGTGATGGTTGCTCCGGTTCGTCTGATGACGCTGGCTTATTTGCGCGTAATTTGCGCATTAATCGCTGCCGACAAAGGCGCAGCAc                                            >  1:135502/1‑251 (MQ=255)                                                                                                                                                                                                                          tCCAGTCCATCAATCGCCAGTAGCCAA‑TCACAAAACGGGAAAACTCCGGGCCTTGCGGCGCAATAGTAATACGCTGAACCATAATCGCTTCCTCTTATCAGATATGAGAGGAGTATACGCAAGATTAGGTTCAAAAGAGTGATGGTTGCTCCGGTTCGTCTGATGACGCTGGCTTATTTGCGCGTAATTTGCGCATTAATCGCTGCCGACAAAGGCGCAGCACCTCTTGTTTTTCGCCATCGCTCAtttt                 >  1:166482/1‑250 (MQ=255)                                                                                                                                                                                                                                    cAATCGCCAGTAGCCAA‑TCACAAAACGGGAAAACTCCGGGCCTTGCGGCGCAATAGTAATGCGCTGAACCATAATCGCTTCCTCTTATCAGATATGAGAGGAGTATACGCAAGATTAGGTTCAAAAGAGTGATGGTTGCTCCGGTTCGTCTGATGACGCTGGCTTATTTGCGCGTAATTTGCTCATTAATCGCTTCTGACAAAGGCGCAGCACCTCTTTTTTGTCGTCATCGGTCATTTTATTCCAGTTa       <  2:156645/250‑1 (MQ=255)                                                                                                                                                                                                                                      aTCGCCAGTAGCCAA‑TCACAAAACGGGAAAACTCCGGGCCTTGCGGCGCAATAGTAATACGCTGAACCATAATCGCTTACTCTTATCAGATATGAGAGGAGTATACGCAAGATTAGGTTCAAAAGAGTGATGGTTGCTCCGGTTCGTCTGATGACGCTGGCTTATTTGCGCGTAATTTGCGCATTAATCGCTGCCGACAAAGGCGCAGCACCTCTTGTTTTTCGCCATCGCTCATTTTATTCCAGTTaaa     >  1:86266/1‑250 (MQ=255)                                                                                                                                                                                                                                        cGCCAGTAGCCAA‑TCACAAAACGGGAAAACGCCGGGCCTTGCGGCGCAATAGTAATACGCTGAACCATAATCGCTTCCTCTTATCAGAGATGAGAGGAGTATACGCAAGATTAGGTTCAAAAGAGTGATGGTTGCTCCGGTTCGTCTGATGACGCTGGCTTATTTGCGCGTAATTTGCGCATTAATCGCTGCCGACAAAGGCGCAGCACCCCTTGTTTTTCGCCATCGCTCATTTTATTCCAGTTAAAACg  <  1:63471/251‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                 ATGTAATGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATACCGCATTTACTGACGATTTCCATCCGTTCACGCAGGTGAGGTGCCAGTTTCAGTGCCTCGCCAAACGCCGCTTCGCACTGATAGCCACCATAAATATCAGCATGGTCCACGGTGGTCACGCCGAGATCCAGATGCTCTTCAATAAAACTGACCAGCTGGCGGGCGGACATATTCCAGTCCATCAATCGCCAGTAGCCCA‑TCACAAAACGGGAAAACTCCGGGCCTTGCGGCGCAATAGTAATACGCTGAACCATAATCGCTTCCTCTTATCAGATATGAGAGGAGTATACGCAAGATTAGGTTCAAAAGAGTGATGGTTGCTCCGGTTCGTCTGATGACGCTGGCTTATTTGCGCGTAATTTGCGCATTAATCGCTGCCGACAAAGGCGCAGCACCTCTTGTTTTTCGCCATCGCTCATTTTATTCCAGTTAAAACG  >  NC_000913/1725334‑1725815

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