Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,268,970	A→G	T186A (ACA→GCA)	yhaC →	uncharacterized protein YhaC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,268,970	0	A	G	94.7%	63.1 / NA	19	T186A (ACA→GCA)	yhaC	uncharacterized protein YhaC
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (6/12);  minor base T (0/1);  total (6/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

AAAAACCAGGAAATATCCACAGGTAATATCAATATTAATTGCCCATTAAAAGATATTTACAATGAAATCAGGAGGTTAAACGTAATTTTTAGTTGTGGGACTGGAGATATCGTTGATCTATCCTCTCTGGACTTACGTAATGTCGATTTAGATTATTATGATTTCACAGATAAACATATGGCTAATACTATTTTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTAGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCATTAGCGATATGTATTCTGAAGAAAAAATCAAGTTTGTCGAAAGTTGTTTAAATGAGCCTGGAGACATTACCGAAGAACAATTAAAAATTATGAGATTTGCAATTATAAAATCTATACCAAGGGCAACTCTTACAGA  >  NC_000913/3268754‑3269214                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                     aaaaaCCAGGAAATATCCACAGGTAATATCAATATTAATTGACCATTAAAAGATATTTACAATGAAATCAGGAGGTTAAACGTAATTTTTAGTTGTGGGACTGGAGATATCGTTGATCTATCCTCTCTGGACTTACGTAATGTCGATTTAGATTATTATGATTTCACAGATAAACATATGGCTAATACTATTTTAAATCCTTTTAAATTGAATTCAGCAAATTTTACTAATGCCAACATGTTTCAGGTTaa                                                                                                                                                                                                                    <  2:135329/251‑1 (MQ=255)                            atCAATATTAATTGCCCATTAAAAGATATTTACAATGAAATCAGGAGGTTAAACGTAATTTTTAGTTGTGGGACTGGAGATATCGTTGATCTATCCTCTCTGGACTTACGTAATGTCGATTTAGATTATTATGATTTCACAGATAAATATATGGCTAATACTATTTTAAATCCTTTTAAATTGAATTCAGCAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTAGTTCAACACAAAACGCcac                                                                                                                                                                                         <  1:322631/251‑1 (MQ=255)                                      aTTGCCCATTAAAAGATATTTACAATGAAATCAGGAGGTTAAACGTAATTTTTAGTTGTGGGACTGGAGATATCGTTGATCTATCCTCTCTGGACTTACGTAATGTCGATTTAGATTATTATGATTTCACAGATAAACATATGGCTAATACTATTTTAAATCCTTTTAAATTGAATTCAGCAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTAGTTCAACACAAAACGCCACAATCTCCTgg                                                                                                                                                                               <  1:196276/251‑1 (MQ=255)                                                       atTTCCAATGAAATCAGGAGGTTAATCGTAATTTTTAGTTGTGGGACTGGAGATATCGTTGATCTATCCTCTCTGGACTTACGTAATGTCGATTTAGATTATTATGATTTCACAGATAAACATATGGCTAATACTATTTTAAATCCTTTTAAATTGAATTCAGCAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTAGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAAt                                                                                                                                                              >  2:2952/1‑251 (MQ=255)                                                                  aaTCAGGAGGTTAAACGTAATTTTTAGTTGTGGGACTGGAGATATCGTTGATCTATCCTCTCTGGACTTACGTAATGTCGATTTAGATTATTATGATTTCACAGATAAACATATGGCTAATACTATTTTAAATCCTTTTAAATTGAATTCAGCAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTAGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGttt                                                                                                                                                   >  1:248910/1‑251 (MQ=255)                                                                                 cGTAATTTTTAGTTGTGGGACTGGAGATATCGTTGATCTATCCTCTCTGGACTTACGTAATGTCGATTTAGATTATTATGATTTCACAGATAAACATATGGCTAATACTATTTTAAATCCTTTTAAATTGAATTCAGCAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTAGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCATTAGCg                                                                                                                                    <  1:283157/251‑1 (MQ=255)                                                                                               gtgGGACTGGAGATATCGTTTATCTATCCTCTCTGGACTTTCGTAATGTCGATTTAGATTATTATGATTTCACAGATAAACATATGGCTAATACTAGTTTAAATCCTTTTAAATTGAATTCATCAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGGTAGTTCAACACAAACCGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCATTAGCGATATGTATTCTga                                                                                                                       <  1:267815/250‑1 (MQ=255)                                                                                               gtgGGACTGGAGATATCGTTGATCTATCCTCTCTGGACTTACGTAATGTCGATTTAGATTATTATGATTTCACAGATAAACATATGGCTAATACTATTTTAAATCCTTTTAAATTGAATTCAGCAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTAGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCATTAGCGATATGTATTCTga                                                                                                                       <  1:267806/250‑1 (MQ=255)                                                                                                             atCGTTGATCTATCCTCTCTGGACTTACGTAATGTCGATTTAGATTATTATGATTTCACAGATAAACATATGGCTAATACTATTTTAAATCCTTTTAAATTGAATTCAGCAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTAGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCATTAGCGATATGTATTCTGAAGAAAAAATCAAGtt                                                                                                        >  2:67015/1‑251 (MQ=255)                                                                                                             atCGTTGATCTATCCTCTCTGGACTTACGTAATGTCGATTTAGATTATTATGATTTCACAGATAAACATATGGATAATACTATTTTAAATCCTTTTAAATTGAATTCAGCAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTAGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCATTAGCGATATGTATTCTGAAGAAAAAATCAAGtt                                                                                                        >  2:22949/1‑251 (MQ=255)                                                                                                                   gATCTTTCCTCTCTGGAATTACGTAATGTCGATTTAGATTATTATGATTTCACAGATAAACATATGGCTAATACTAATTTAAATCCTTTTAAATTGAATTCAGCAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTAGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCATTAGCGATATGTATTCTGAAGAAAAAATCAAGTTTGTCg                                                                                                   <  2:282790/250‑1 (MQ=255)                                                                                                                              tctGGACTTACGTAATGTCGATTTAGATTATTATGATTTCACAGATAAACATATGGCTAATACTATTTTAAATCCTTTTAAATTGAATTCAGCAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTAGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCATTAGCGATATGTATTCTGAAGAAAAAATCAAGTTTGTCGAAAGTTGTTTa                                                                                        >  1:73183/1‑250 (MQ=255)                                                                                                                                  gACTTACGTAATGTCGATTTAGATTATTATGATTTCACAGATAAACATATGGCTAATACTATTTTAAATCCTTTTAAATTGAATTCAGCAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTAGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCATTAGCGATATGTAGTCAGAAGAAAAAATCAAGTTTGTCGAAAGTTGTTTAAATGa                                                                                   <  1:368743/251‑1 (MQ=255)                                                                                                                                                 gATTTAGATTATTATGATTTCACAGATAAACATATGGCTAATACTATTTTAAATCCTTTTAAATTGAATTCAGCAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTAGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCATTAGCGATATGTATTCTGAAGAAAAAATCAAGTTTGTCGAAAGTTGTTTAAATGAGCCTGGAGACATTAc                                                                    <  1:195754/251‑1 (MQ=255)                                                                                                                                                           attatGATTTCACAGATAAACATATGGCTAATACTATTTTAAATCCTTTTAAATTGAATTCAGCAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTAGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCATTAGCGATATGTATTCTGAAGAAAAAATCAAGTTTGTCGAAAGTTGTTTAAATGAGCCTGGAGACATTACCGAAGAACa                                                           <  1:353881/250‑1 (MQ=255)                                                                                                                                                                              aCATATGGCTAATACTATTTTAAATCCTTTTAAATTGAATTCAGCAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTAGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCATTAGCGATATGTATTCTGAAGAAAAAATCAAGTTTGTCGAAAGTTGTTTAAATGAGCCTGGAGACATTACCGAAGAACAATTAAAAATTATGAGATTTg                                       >  1:211525/1‑251 (MQ=255)                                                                                                                                                                                           aCTATTTTAAATCCTTTTAAATTGAATTCAGCAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTAGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCATTAGCGATATGTATTCTGAAGAAAAAATCAAGTTTGTCGAAAGTTGTTTAAATGAGCCTGGAGACATTACCGAAGAACAATTAAAAATTATGAGATTTGCAATTATAAAAt                           <  1:67015/250‑1 (MQ=255)                                                                                                                                                                                                               aTTGAATTCAGCAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTAGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCATTAGCGATATGTATTCTGAAGAAAAAATCAAGTTTGTCGAAAGTTGTTTAAATGAGCCTGGAGACATTACCGAAGAACAATTAAAAATTATGAGATTTGCAATTATAAAATCTATACCAAGGGCAACTCTTa      <  1:58647/251‑1 (MQ=255)                                                                                                                                                                                                                   aaTTCAGCAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTAGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCATTAGCGATATGTATTCTGAAGAAAAAATCAAGTTTGTCGAAAGTTGTTTAAATGAGCCTGGAGACATTACCGAAGAACAATTAAAAATTATGAGATTTGCAATTATAAAATCTATACCAAGGGCAACTCTTActga  <  2:364070/251‑4 (MQ=255)                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                     AAAAACCAGGAAATATCCACAGGTAATATCAATATTAATTGCCCATTAAAAGATATTTACAATGAAATCAGGAGGTTAAACGTAATTTTTAGTTGTGGGACTGGAGATATCGTTGATCTATCCTCTCTGGACTTACGTAATGTCGATTTAGATTATTATGATTTCACAGATAAACATATGGCTAATACTATTTTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTAGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCATTAGCGATATGTATTCTGAAGAAAAAATCAAGTTTGTCGAAAGTTGTTTAAATGAGCCTGGAGACATTACCGAAGAACAATTAAAAATTATGAGATTTGCAATTATAAAATCTATACCAAGGGCAACTCTTACAGA  >  NC_000913/3268754‑3269214

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