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	95.8%	76.6 / ‑6.4	24	T186A (ACA→GCA)	yhaC	uncharacterized protein YhaC
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (13/10);  total (13/11)
Fisher's exact test for biased strand distribution p-value = 4.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

GTGAACTTACCATTAATGAACCGTTAAAAAACCAGGAAATATCCACAGGTAATATCAATATTAATTGCCCATTAAAAGATATTTACAATGAAATCAGGAGGTTAAACGTAATTTTTAGTTGTGGGACTGGAGATATCGTTGATCTATCCTCTCTGGACTTACGTAATGTCGATTTAGATTATTATGATTTCACAGATAAACATATGGCTAATACTATTTTAAATCCTTTTAAATTGAATTCAACAAA‑TTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTAGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCATTAGCGATATGTATTCTGAAGAAAAAATCAAGTTTGTCGAAAGTTGTTTAAATGAGCCTGGAGACATTACCGAAGAACAATTAAAAATTATGAGATTTGCAATTATAAAATCTATACCAAGGGC  >  NC_000913/3268728‑3269202                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                          gtgAACTTACCATTAATGAACCGTTAAAAAACCAGGAAATATCCACAGGTAATATCAATATTAATTGCCCATTAAAAGATATTTACAATGAAATCAGGAGGTTAAACGTAATTTTTAGTTGTGGGACTGGAGATATCGTTGATCTATCCTCTCTGGACTTACGTAATGTCGATTTAGATTATTATGATTTCACATATAAACATATGGCTAATACTATTTTAAATCCTTTTAAAGTGAATTCAGCAAA‑ttt                                                                                                                                                                                                                                   <  2:182783/250‑1 (MQ=255) gtgAACTTACCATTAATGAACAGTTAAAAAACCAGGAAATATCCACAGGTAATAACAATATTAATTGCCCATTAAAAGATATTTACAATGAAATTAGGAGGTTAAACGTAAGTTTTAGTTGTGGGACTGGAGATATCGTTGATCTATCCGCTCTGGACTTACGTAATGTCGATGTAGATTATTATGATTTCACAGATAAACATATTGCTAATACTGTTTTAAATCCTTTTAAAGTTAATTCAACAAA‑tttt                                                                                                                                                                                                                                  <  2:213917/251‑1 (MQ=255)        tACCATTAATGAACCGTTAAAAAACCAGGAAATATCCACAGGTAATATCAATATTAATTGCCCATTAAAAGATATTTACAATGAAATCAGGAGGTTAAACGTAATTTTTAGTTGTGGGACTGGAGATATCGTTTATCTATCCTCTCTGGACTTACGTAATGTCGATTTAGATTATTATGATTTCACAGATACACATATGGCTAATACTATTTTAAATCCTTTTCAATTTATTTCAGCAAA‑TTTTTCTaa                                                                                                                                                                                                                             <  2:203765/249‑1 (MQ=255)                            aaaaCCAGGAAATATCCACAGGTAATATCAATATTAATTGCCCATTAAAAGATATTTACAATGAAATCAGGATGTTAAACGTAATTTTTAGTTGTGGGACTGGAGATATCGTTGATCTATCCTCTCTGGACTTACGTAATGTCGATTTAGATTATTATGATTTCACAGATAAACATATGGCTAATACTATTTTAAATCCTTTTAAATTGAATTCAGCAAA‑TTTTACTAATGCCAACATGTTTCAGGTTAAt                                                                                                                                                                                                       <  2:372761/251‑1 (MQ=255)                              aaCCAGGAAATATCCACAGGTAATATCAATATTAATTGCCCATTAAAAGATATTTACAATGAAATCAGGAGGTTAAACGTAATTTTTAGTTGTGGGACTGGAGATATCGTTGATCTATCCTCTCTGGACTTACGTAATGTCGATTTAGATTATTATGATTTCACAGATAAACATATGGCTAATACTATTTTAAATCCTTTTAAATTGAATTCAGCAAA‑TTTTACTAATGCCAACATGTTTCAGGTTAAttt                                                                                                                                                                                                     >  1:108757/1‑251 (MQ=255)                                                        cAATATTAATTGCCCATTAAAAGATATTTACAATGAAATCAGGAGGTTAAACGTAATTTTTAGTTGTGGGACTGGAGATATCGTTGATCTATCCTCTCTGGACTTACGTAATGTCGATTTAGATTATTATGATTTCACAGATAAACATATGGCTAATACTATTTTAAATC‑TTTTAAATTTAATTTAGCAAAATTTTA‑TAATTCCAACATGTTTTAGGTTTATTTTggttgtttaaccacaaacgcccaaa                                                                                                                                                                           >  1:20900/1‑226 (MQ=255)                                                              tAATTGCCCATTAAAAGATATTTACAATGAAATCAGGAGGTTAAACGTAATTTTTAGTTGTGGGACTGGAGATATCGTTGATCTATCCTCTCTGGACTTACGTAATGTCGATTTAGATTATTATGATTTCACAGATAAACATATGGCTAATACTATTTTAAATCCTTTTAAATTGAATTCAGCAAA‑TTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTAGTTCAACACAAAACGCCACAAtctcat                                                                                                                                                                     >  1:58001/1‑249 (MQ=255)                                                               aaTTGCCCATTAAAAGATATTTACAATGAAATCAGGAGGTTAAACGTAATTTTTAGTTGTGGGACTGGAGATATCGTTGATCTATCCTCTCTGGACTTACGTAATGTCGATTTAGATTATTATGATTTCACAGATAAACATATGGCTAATACTATTTTAAATCCTTTTAAATTGAATTCAGCAAA‑TTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTAGTTCAACACAAAACGCCACAATCTCCTg                                                                                                                                                                    >  2:135597/1‑251 (MQ=255)                                                                  tGCCCATTAAACGATATTTACAATGAAATCAGGAGGTTAAACGTAATTTTTAGTTGTGGGACTGGAGATATCGTTGATCTATCCTCTCTGGACTTACGTAATGTCGATTTAGATTATTATGATTTCACAGATAAACATATGGCTAATACTATTTTAAATCCTTTTAAATTGAATTCAGCAAA‑TTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTAGTTCAACACAAAACGCCACAATCTCCTggg                                                                                                                                                                  >  1:31307/1‑250 (MQ=255)                                                                             aGATATTTACAATGAAATCAGGAGGTTAAACGTAATTTTTAGTTGTGGGACTGGAGATATCGTTGATCTATCCTCTCTGGACTTACGTAATGTCGATTTAGATTATTATGATTTCACAGATAAACATATGGCTAATACTATTTTAAATCCTTTTAAATTGAATTCAGCAAA‑TTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTAGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTa                                                                                                                                                       >  1:322906/1‑250 (MQ=255)                                                                                           aaaTCATGATGTTAAACGTAATTTTTAGTTGTGGGACTGGAGATATCGTTGATCTATCCTCTCTGGACTTACGTAATGTCGATTTAGATTATTATGATTTCACAGATAAACATATGGCTAATACTATTTTAAATCCTTTTAAATTGACTTCAGCAAA‑TTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTAGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAACAATAACGCCTGtt                                                                                                                                        <  2:138499/251‑1 (MQ=255)                                                                                            aaTCAGGAGGTTAAACGTAATTTTTAGTTGTGGGACTGGAGATATCGTTGATCTATCCTCTCTGGACTTACGTAATGTCGATTTAGATTATTATGATTTCACAGATAAACATATGGCTAATACTATTTTAAATCCTTTTAAATTGAATTCAGCAAA‑tttt                                                                                                                                                                                                                                  <  1:314479/160‑1 (MQ=255)                                                                                            aaTCAGGAGGTTAAACGTAATTTTTAGTTGTGGGACTGGAGATATCGTTGATCTATCCTCTCTGGACTTACGTAATGTCGATTTAGATTATTATGATTTCACAGATAAACATATGGCTAATACTATTTTAAATCCTTTTAAATTGAATTCAGCAAA‑tttt                                                                                                                                                                                                                                  >  2:314479/1‑160 (MQ=255)                                                                                                         aaCGTAATTTTTAGTTGTGGGACTGGAGATATCGTTGATCTATCCTCTCTGGACTTACGTAATGTCGATTTAGATTATTATGATTTCACAGATAAACATATGGCTAATACTATTTTAAATCCTTTTAAATTGAATTCAGCAAA‑TTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTAGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCATTAg                                                                                                                          <  1:249652/251‑1 (MQ=255)                                                                                                                 ttttAGTTGTGGGACTGGAGATATCGTTGATCTATCCTCTCTGGACTTACGTAATGTCGATTTAGATTATTATGATTTCACAGATAAACATATGGCTAATACTATTTTAAATCCTTTTAAATTGAATTCAGCAAA‑TTTTACTAAGGCCAACATGTTTCAGGTTAATTTTGTTAGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCATTAGCGATATg                                                                                                                   <  1:192953/250‑1 (MQ=255)                                                                                                                                 ggAGATATCGTTGATCTATCCTCTCTGGACTTACGTAATGTCGATTTAGATTATTATGATTTCACAGATAAACATATGGCTAATACTATTTTAAATCCTTTTAAATTGAATTCAGCAAA‑TTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTAGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCATTAGCGATATGTATTCTGAAGAAAAAAt                                                                                                  >  1:99711/1‑251 (MQ=255)                                                                                                                                          gTTGATCTATCCTCTCTGGACTTACGTAATGTCGATTTAGATTATTATGATTTCACAGATAAACATATGGCTAATACTATTTTAAATCCTTTTAAATTGAATTCAGCAAA‑TTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTAGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTa                                                                                                                                                       >  1:228422/1‑189 (MQ=255)                                                                                                                                          gTTGATCTATCCTCTCTGGACTTACGTAATGTCGATTTAGATTATTATGATTTCACAGATAAACATATGGCTAATACTATTTTAAATCCTTTTAAATTGAATTCAGCAAA‑TTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTAGTTCAACACAAAACGCCACAAGCTCCTGGGATTATTTACTa                                                                                                                                                       <  2:228422/189‑1 (MQ=255)                                                                                                                                                                                aGATTATTATGATTTCACAGATAAACATATGGCTAATACTATTTTAAATCCTTTTAAATTGAATTCAGCAAA‑TTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTAGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCATTAGCGATATGTATTCTGAAGAAAAAATCAAGTTTGTCGAAAGTTGTTTAAATGAGCCTGGAGACATTACCgaag                                                   <  1:346945/251‑1 (MQ=255)                                                                                                                                                                                    tattatGATTTCACAGATAAACATATGGCTAATACTATTTTAAATCCTTTTAAATTGAATTCAGCAAA‑TTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTAGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCATTATCGATATGTATTCTGAAGAAAAAATCAAGTTTGTCGAAAGTTGTTTAAATGAGCCTGGAGACATTACCgaagacc                                                >  1:293334/1‑248 (MQ=255)                                                                                                                                                                                                acaGATAAACATATGGCTAATACTATTTTAAATCCTTTTAAATTGAATTCAGCAAA‑TTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTAGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCATTAGCGATATGTATTCTGAAGAAAAAATCAAGTTTGTCGAAAGTTGTTTAAATGAGCCTGGAGACATTACCGAAGAACAATTAAAAATTAt                                   >  1:50572/1‑251 (MQ=255)                                                                                                                                                                                                         cATATGGCTAATACTATTTTAAATCCTTTTAAATTGAATTCAGCAAA‑TTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTAGTTCAACACAAAACGCAACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCATTAGCGATATGTATTCTGAAGAAAAAATCAAGTTTGTCGAAAGTTGTTTAAATGAGCCTGGAGACATTACCGAAGAACAATTAAAAATTATGAGATTTGc                          >  2:19674/1‑251 (MQ=255)                                                                                                                                                                                                                          tttAAATCCTTTTAAATTGAATTCAGCAAA‑TTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTAGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCATTAGCGATATGTATTCTGAAGAAAAAATCAAGTTTGTCGAAAGTTGTTTAAATGAGCCTGGAGACATTACCGAAGAACAATTAAAAATTATGAGATTTGCAATTATAAAATCTATAc         >  1:319643/1‑251 (MQ=255)                                                                                                                                                                                                                                 ccTTTTAAATTGAATTCAGCAAA‑TTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTAGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCATTAGCGATATGTATTCTGAAGAAAAAATCAAGTTTGTCGAAAGTTGTTTAAATGAGCCTGGAGACATTACCGAAGAACAATTAAAAATTATGAGATTTGCAATTATAAAATCTATACCAAGGGc  <  1:81457/251‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                          GTGAACTTACCATTAATGAACCGTTAAAAAACCAGGAAATATCCACAGGTAATATCAATATTAATTGCCCATTAAAAGATATTTACAATGAAATCAGGAGGTTAAACGTAATTTTTAGTTGTGGGACTGGAGATATCGTTGATCTATCCTCTCTGGACTTACGTAATGTCGATTTAGATTATTATGATTTCACAGATAAACATATGGCTAATACTATTTTAAATCCTTTTAAATTGAATTCAACAAA‑TTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTAGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCATTAGCGATATGTATTCTGAAGAAAAAATCAAGTTTGTCGAAAGTTGTTTAAATGAGCCTGGAGACATTACCGAAGAACAATTAAAAATTATGAGATTTGCAATTATAAAATCTATACCAAGGGC  >  NC_000913/3268728‑3269202

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