Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,268,835	G→T	V141L (GTA→TTA)	yhaC →	uncharacterized protein YhaC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,268,835	0	G	T	95.2%	64.5 / NA	21	V141L (GTA→TTA)	yhaC	uncharacterized protein YhaC
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (8/12);  total (8/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 = 1.00e+00

TGATTTTCAAAAAGAATTAGCGTACTTACAAAATGTAATTCTTTATAATGTCTCGTCTCTTCGTCTGGATTTTTTAGGTTATAACGCCCAAATTATTCAACGATCGGACAATACTTGTGAACTTACCATTAATGAACCGTTAAAAAACCAGGAAATATCCACAGGTAATATCAATATTAATTGCCCATTAAAAGATATTTACAATGAAATCAGGAGGTTAAACGTAATTTTTAGTTGTGGGACTGGAGATATCGTTGATCTATCCTCTCTGGACTTACGTAATGTCGATTTAGATTATTATGATTTCACAGATAAACATATGGCTAATACTATTTTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTAGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCATTAGCG  >  NC_000913/3268612‑3269084                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                          tGATTTTCAAAAAGAATTAGCGTACTTACAAAATGTAATTCTTTATAATGTCTCGTCTCTTCGTCTGGATTTTTTAGGTTATAACGCCCAAATTATTCAACGATCGGACAATACTTGTGAACTTACCATTAATGAACCGTTAAAAAACCAGGAAATATCCACAGGTAATATCAATATTAATTGCCCATTAAAAGATATTTACAATGAAATCAGGAGGTTAAACTTAATTTTTAGTTGTGGGACTGGAGat                                                                                                                                                                                                                                 >  2:573949/1‑250 (MQ=255)               aaaaaGCGTACTTACAAAATGTAATTCTTTATAATGTCTCGTCTCTTCGTCTGGATTTTTTAGGTTATAACGCCCAAATTATTCAACGATCGGACAATACTTGTGAACTTACCATTAATGAACCGTTAAAAAACCAGGAAATATCCACAGGTAATATCAATATTAATTGCCCATTAAAAGATATTTACAATGAAATCAGGAGGTTAAACTTAATTTTTAGTTGTGGGATTGGAGATATCGTTGATCTATcc                                                                                                                                                                                                                  >  2:123990/5‑251 (MQ=255)                  tAGCGTACTTACAAAATGTAATTCTTTATAATGTCTCGTCTCTTCGTCTGGATTTTTTAGGTTATAACGCCCAAATTATTCAACGATCGGACAATACTTGTGAACTTACCATTAATGAACCGTTAAAAAACCAGGAAATATCCACAGGTAATATCAATATTAATTGCCCATTAAAAGATATTTACAATGAAATCAGGAGGTTAAACTTAATTTTTAGTTGTGGGACTGGAGATATCGTTGATCTATCctc                                                                                                                                                                                                                <  1:8892/250‑1 (MQ=255)                    gCGTACTTACAAAATGTAATTCTTTATAATGTCTCGTCTCTTCGTCTGGATTTTTTAGGTTATAACGCCCAAATTATTCAACGATCGGACAATACTTGTGAACTTACCATTAATGAACCGTTAAAAAACCAGGAAATATCCACAGGTAATATCAATATTAATTGCCCATTAAAAGATATTTACAATGAAATCAGGAGGTTAAACTTAATTTTTAGTTGTGGGACTGGAGATATCGTTGATCTATCctctct                                                                                                                                                                                                             <  2:225645/251‑1 (MQ=255)                                   gtaaatccttATAATGGCTCGTCTCTTCGTCTGGATTGGTTGGGTTATAACGCCCAAATTATTCAACGATCGGACAATACTTGTGAACGTACCATTAATGAACCGGTAAAAAACCAGGAAAGATCCACAGGGAATATCAATATTAATTGCCAAGGAAAAGATATTTACAATGAAATCAGGTGGTTAAACTTAATTATTATTTGTGGGACAGGAGATATCGATGATCTAGCCTCTCTGGACTTACGTAATGt                                                                                                                                                                                              <  2:121673/243‑1 (MQ=255)                                   gTAATTCTTTATAATGTCTCGTCTCTTCGTCTGGATTTTTTAGGTTATAACGCCCAAATTATTCAACGATCGGACAATACTTGTGAACTTACCATTAATGAACCGTTAAAAAACCAGGAAATATCCACAGGTAATATCAATAGTAATTGCCCATTAAAAGATATTTATAATGAAGTCAGGAGTTTATACTTAATTTTTAGTTGTGGGACTGGATATATCGTTGATCTATCCTCTCTGGACTTACGTAATGt                                                                                                                                                                                              <  1:345845/251‑1 (MQ=255)                                     aaTTCTTTATAATGTCTCGTCTCTTCGTCTGGATTTTTTAGGTTATAACGCCCAAATTATTCAACGATCGGACAATACTTGTGAACTTACCATTAATGAACCGTTAAAAAACCAGGAAATATCCACAGGTAATATCAATATTAATTGCCCATTAAAAGATATTTACAATGAAATCAGGAGGTTAAACTTAATTTTTAGTTGTGGGACTGGAGATATCGTTGATCTATCCTCTCTGGACTTACGTAATGTCg                                                                                                                                                                                            <  1:591494/251‑1 (MQ=255)                                                                   ggATTTTTTAGGTTATAACGCCCAAATTATTCAACGATCGGCCAATACTTGTGAACTTACCATTAATGAACCGTTAAAAAACCAGGACATATCCACAGGTAATATCAATATTAATTGCCCATTAAAAGATATTTACAATGAAATCAGGAGGTTACACTTAATTTTTAGTTGTGGGACTGGAGATATCGTTGATCTATCCTCTCTGGACTTACGTAATGTCGATTTAGATTATTATGATTTCACAGATaaa                                                                                                                                                               >  2:269597/1‑250 (MQ=255)                                                                              gTTATAACGCCCAAATTATTCAACGATCGGACAATACTTGTGAACTTACCATTAATGAACCGTTAAAAAACCAGGAAATAGCCACAGGTAATAACAATATTAATTGCCCATTAAAAGATATTTACAATGAAATCAGGAGGTTAAACTTAATTTTTAGTTGTGGTACTGGAGATATCGTTGATCTATCCTCTCTGGACTTGCGTTATGTCGATGTAGATTATTATGATTTCACAGATAAACATATGGCTAAt                                                                                                                                                   <  1:618495/251‑1 (MQ=255)                                                                                               attCAACGATCGGACAATACTTGTGAACTTACCATTAATGAACCGTTAAAAAAACAGGAAATATCAACAGGTAATCTCAATATTAATTGCCCATTAAAAGATATTTACAATGAAATCAGGAGGTTAAACTTAATTTTTAGTTGTGGGACTGGAGATATCGTTGATCTATCCTCTCTGGACTTACGTAATGTCGATTTAGATTATTATGATTTCACAGATAAACATATGGCTAATACTATTTTAAATCCttt                                                                                                                                  >  1:281125/1‑251 (MQ=255)                                                                                                                     gtgAACTTACCATTAATGAACCGTTAAAAAACCAGGAAATATCCACAGGTAATAGCAATATTAATTGCCCATTAAAAGATATTTACAATGAAATCAGGAGGTTAAACTTAAGTTTTAGTTGTGGGACTTGAGATATCGTTGATCTATCCTCTCTGGACGTACGTAATTTCGATTTAGATTATTATGATTTCACAGATAATCATATGGCTAATACTATTTTAAAGCCTTTTAAATTGAATTCAACAAAttt                                                                                                             <  2:123841/250‑1 (MQ=255)                                                                                                                                    aTGAACCGTTAAAAAACCAGGAACTATCCACAGGTAATATCAATATTAATTGCCCATTAAAAGATATTTACAATGAAATCAGGAGGTTAAACTTAATTTTTAGTTGTGGGACTGGAGATATCGTTGATCTATCCTCTCTGGACTTACGTAATGTCGATTTAGATTATTATGATTTCACAGATAAACATATGGCTAATACTATTTTAACTCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATg                                                                                             >  2:333557/1‑251 (MQ=255)                                                                                                                                              aaaaaaCCAGGAAATATCCACAGGTAATATCAATATTAATTGCCCATTAAAAGATATTTACAATGAAATCAGGAGGTTAAACTTAATTTTTAGTTGTGGGACTGGAGATATCGTTGATCTATCCTCTCTGGACTTACGTAATGTCGATTTAGATTATTATGATTTCACAGATAAACATATGGCTAATACTATTTTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGtt                                                                                    >  2:176443/1‑250 (MQ=255)                                                                                                                                                             atCCACAGGTAATATCAATATTAATTGCCCATTAAAAGATATTTACAATGAAATCAGGAGGTTAAACTTAATTTTTAGTTGTGGGACTGGAGATATCGTTGATCTATCCTCTCTGGACTTACGTAATGTCGATTTAGATTATTATGATTTCACAGATAAACATATGGCTAATACTATTTTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTAGTTCa                                                                     >  2:208890/1‑250 (MQ=255)                                                                                                                                                                     gTAATATCAATATTAATTGCCCATTAAAAGATATTTACAATTAAATCAGGAGGTTAAACTTAATTTTTAGTTGTGAGACTGGAGATATCGTTGATCTATCCTCTCTTGACTTTCGTAATTGCGATTTAGATTATTATGATTTCACAGATAAACATATGGCTAATACTATGTTAAATCCTTTTAAAGTGAATTCAACAAATTTTACTAATGCCAACATGTGTCAGGTTCATTTTGTTAGTTCATCACAAAAc                                                            <  2:579656/251‑1 (MQ=255)                                                                                                                                                                                       gCCCATTAAAAGATATTTACAATGAAATCAGGAGGTTAAACTTAATTTTTAGTTGTGGGACTGGAGATATCGTTGATCTATCCTCTCTGGACTTACGTAATGTCGATTTAGATTATGATGATTTCACAGATAAACATATGGCTAATACTATTTTAAATCCGTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCATGTTAATTGTGTTAGTTCAACACAAAACGCCACAATCTCCTGGGat                                          <  1:307890/251‑1 (MQ=255)                                                                                                                                                                                                 aGATATTTACAATGAAATCAGGAGGTTAAACTTAATTTTTAGTTGTGGGACTGGAGATATCGTTGATCTATCCTCTCTGGACTTACGTAATGTCGATTTAGATTATTATGATTTCACAGATAAACATATGGCTAATACTATTTTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTAGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTa                                 >  1:507772/1‑250 (MQ=255)                                                                                                                                                                                                     atTTACAATGAAATCAGGAGGTTAAACTTAATTTTTAGTTGTGGGACTGGAGATATCGTTGATCTATCCTCTCTGGACTTACGTAATGTCGATTTAGATTATTATGATTTCACAGATAAACATATGGCTAATACTATTTTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTAGTGCAACACAAAACGCCACAATCTCCTGGGATTATTTACTaaaaa                             <  2:552258/250‑1 (MQ=255)                                                                                                                                                                                                         aCAATGAAATCAGGAGGTTAAACTTAATTTTTAGTTGTGGGACTGGAGATATCGTTGATCTATCCTCTCTGGACTTACGTAATGTCGATTTAGATTATTATGATTTCACAGATAAACATATGGCTAATACTATTTTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTGAATTTTGTTAGTTCAACACAAAGCGCCACAATCTCCTGGGATTATTTACTAAAAATAAc                         <  2:593269/250‑1 (MQ=255)                                                                                                                                                                                                                     ggaggTTAAACTTAATTTTTAGTTGTGGGACTGGAGATATCGTTGATCTATCCTCTCTGGACTTACGTAATGTCGATTTAGATTATTATGATTTCACAGATAAACATATGGCTAATACTATTTTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTAGTTCAACACAAAACGCCGCAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTtaata            <  1:76196/251‑1 (MQ=255)                                                                                                                                                                                                                               cGTAATTTTTAGTTGTGGGACTGGAGATATCGTTGATCTATCCTCTCTGGACTTACGTAATGTCGATTTAGATTATTATGATTTCACAGATAAACATATGGCTAATACTATTTTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTAGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCATTAGCg  <  1:405401/251‑1 (MQ=255)                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                          TGATTTTCAAAAAGAATTAGCGTACTTACAAAATGTAATTCTTTATAATGTCTCGTCTCTTCGTCTGGATTTTTTAGGTTATAACGCCCAAATTATTCAACGATCGGACAATACTTGTGAACTTACCATTAATGAACCGTTAAAAAACCAGGAAATATCCACAGGTAATATCAATATTAATTGCCCATTAAAAGATATTTACAATGAAATCAGGAGGTTAAACGTAATTTTTAGTTGTGGGACTGGAGATATCGTTGATCTATCCTCTCTGGACTTACGTAATGTCGATTTAGATTATTATGATTTCACAGATAAACATATGGCTAATACTATTTTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTAGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCATTAGCG  >  NC_000913/3268612‑3269084

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