Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,438,289	G→A	E195K (GAA→AAA)	ytfI →	protein YtfI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,438,289	0	G	A	94.4%	57.3 / NA	18	E195K (GAA→AAA)	ytfI	protein YtfI
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (7/10);  minor base C (1/0);  total (8/10)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGCGATTGGTGACGTGTCTTTAGGTTATGCCTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGGAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAAATAATTGATGATAATTTTGCCATAGTCATAGAATCTTTGGGCTTTGCGCTTAAAATCTATGGAGGGGATGAACGTTTTCGTAATGGTTCCTCTGTTGTTTTGGAGGATGAAGATTACTCTCTTATTTATAATTTTCTGGTTAATGCAGGTTGCCAACAAGTCGAATTAGCCGAACAAGTGGATGCGATAGTAAGTGCGAATCTGGCCGCAGAT  >  NC_000913/4438042‑4438531                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   cATAACAATAAATGAATATAACTCATTCATGGATTTCAAAACCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGCGATTGGTGACGTGTCTTTAGGTTATGCCTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGGAGGTGGAATTAAAAAAt                                                                                                                                                                                                                                                 >  1:151857/1‑251 (MQ=255)           aaTGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGCGATTGGTGACGTGTCTTTAGGTTATGCCTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGGAGGTGGAATTAAAAAATGCgatggatg                                                                                                                                                                                                                                       <  2:280360/251‑1 (MQ=255)           aaTGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGCGATTGGTGACGTGTCTTTAGGTTATGCCTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGGAGGGGGAATTAAAAAATGCgatggatg                                                                                                                                                                                                                                       >  2:242536/1‑251 (MQ=255)                                gATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGCGATTGGTGACGTGTCTTTAGGTTATGCCTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGGAGGTGGAATTAAAAAATGCGATGGATGGTATATTTACGAAAAaataat                                                                                                                                                                                                                  <  1:63708/251‑1 (MQ=255)                                                 aTTTCCGTTAGAACAAAACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCCGCAGGCGATTGGTGACGTGTCTTTAGGTTATGCCTTATTGGAGGATATAGTAAGCGAATTAAAACGGGGTAAATAGGAACTCCAAATACCTGGAGGTGGAATTAAAAAATGCGATGGATGGTATATTTACGAAAAAATAATTGATGATAATTGTGCCa                                                                                                                                                                                                 <  2:39755/251‑1 (MQ=255)                                                            aacaacACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGCGATTGGTGACGTGTCTTTAGGTTATGCCTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGGAGGTGGAATTAAAAAATGCGATGGATGGTATATTTACGAAAAAATAATTGATGATAATTTTGCCATAGTCATAGa                                                                                                                                                                                       >  1:230393/1‑250 (MQ=255)                                                                                    ggaatagtCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGCGATTGGTGACGTGTTTTTAGGTTATGCCTTATTGGAGGATATAGTAAGCAAATTAAAACGAGGTAAATTTGAACTCCAAATACCAGGAGATTGAATTAAAAAAAGCGATGGATGGTATATTTACGAAAAAATAATTGATGATAATTTTGCCATAGTCATAGAATCTTTGGGCTTTGCGCTTAAAAt                                                                                                                                                               <  2:15817/243‑1 (MQ=255)                                                                                                                     cGAAAAGGATTAAAGGTTATATTCACCAGGCGATTGGTGACGTGTCTTTAGGTTATGCCTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGGAGGTGGAATTAAAAAATGCGATGGATGGTATATTTACGAAAAAATAATTGATGATAATTTTGCCATAGTCATAGAATCTTTGGGCTTTGCGCTTAAAATCTATGGAGGGGATGAACGTTTTCGTAATGGTTcc                                                                                                                             <  1:133270/251‑1 (MQ=255)                                                                                                                       aaaaTGATTAAAGGTTATATTCACCATGCAATTGGTGACGTGTCTTTAGGTTATGCCTTATTGGAAGATATAGTAAGCAAATTAAAACGGTGTAAATTTGAACTCCAAATACCAGGAATTGGAATTAAAAAATGCGATGGATGGTATATTTGCGATAAAATAATTGATGATAATTTTGCCAGAGTCATAGGATCTTTGGGCTTTGCGCTTCAAATCTATGGAGGGGAGGAACGTTGTCGTAATGGGTCctc                                                                                                                           <  2:133520/251‑1 (MQ=255)                                                                                                                                                        ggTGACGTGTCTTTAGGTTATGCCTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGGAGGTGGAATTAAAAAATGCGATGGATGGTATATTTACGAAAAAATAATTGATGATAATTTTGCCATAGTCATAGAATCTTTGGGCTTTGCGCTTAAAATCTATGGAGGGGATGAACGTTTTCGTAATGGTTCCTCTGTTGTTTTGGAGGATGAAGATTACTCTCttat                                                                                          >  2:108909/1‑251 (MQ=255)                                                                                                                                                                gtCTTTAGGTCATGACTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGGAGGTGGAATTAAAAAATGCGATGGATGGAATATTTACGAAAAAATAATTGATGATAATTTTGCCATAGTCATAGAATCTTTGGGCTTTGCGCTTAAAATCTATGGAGGGGATGAACGTTTTCGTAATGGTTCCTCTGTTGTTTTGGAGGATGAAGATTACTCTCTTATTTATNAtt                                                                                  >  2:259256/1‑251 (MQ=255)                                                                                                                                                                             gCCTTATTGGATGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGGAGGTGGAATTAAAAAATGCGATGGATGGTATATTTACGAAAAAATAATTGATGATAATTTTGCCATAGTCATAGAATCTTTGGGCTTTGCGCTTAAAATCTATGGAGGGGATGAACGTTTTCGTAATGGTTCCTCTGTTGTTTTGGAGTATGAAGATTACTCTCTTATTTATAATTTTCTGGTTAATGc                                                                     <  2:155741/251‑1 (MQ=255)                                                                                                                                                                              ccTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGGAGGTGGAATTAAAAAATGCGATGGATGGTATATTTACGAAAAAATAATTGATGATAATTTTGCCATAGTCATAGAATCTTTGGGCTTTGCGCTTAAAATCTATGGAGGGGATGAACGTTTTCGTAATGGTTCCTCTGTTGTTTTGGAGGATGAAGATTACTCTCTTATTTATAATTTTCTGGTTAATGCa                                                                    <  1:102503/251‑1 (MQ=255)                                                                                                                                                                                    tGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACAAGGAGGTGGAATTAAAAAATGCGATGGATGGTATATTTACGAAAAAATGATTGATGATAATTTTGCCATAGTCATAGAATCTTAGGGCTTTGCGCTTAAAATCTATGGAGGGGATGACCGTTTTCGTAATGGTTCCTCTGTTGTTTTGGAGGATGAAGCTTACTCTCTTATTTATAATTTTCTGGTTACTCCAGGTTGc                                                              >  2:147747/1‑251 (MQ=255)                                                                                                                                                                                     ggaggaTATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGGAGGTGGAATTAAAAAATTCGATGGATGGTATATTTACGAAAAAATAATTGATGATAATTTTGCCATAGTCATAGAATCTTTGGGCTTTGCGCTTAAAATCTATGGATGGGATGAACGTTTTCGTAATTGTTCCTCTGTTGTTTTGGAGGATGAAGATTCCTCTCTTATTTATAATTTTCCGGTTAATGCAGGTTGc                                                              <  2:66677/250‑1 (MQ=255)                                                                                                                                                                                                                                         ggaggTGGAATTAAAAAATGCGATGGATGGTATATTTACGAAAAAATAATTGATGATAATTTTGCCATAGTCATAGAATCTTTGGGCTTTGCGCTTAAAATCTATGGAGGGGATGAACGTTTTCGTAATGGTTCCTCTGTTGTTTTGGAGGATGAAGATTACTCTCTTATTTATAATTTTCTGGTTAATGCAGGTTGCCAACAAGTCGAATTAGCCGAACAAGTGGATGCGATAGt                        <  2:138753/236‑1 (MQ=255)                                                                                                                                                                                                                                         ggaggTGGAATTAAAAAATGCGATGGATGGTATATTTACGAAAAAATAATTGATGATAATTTTGCCATAGTCATAGAATCTTTGGGCTTTGCGCTTAAAATCTATGGAGGGGATGAACGTTTTCGTAATGGTTCCTCTGTTGTTTTGGAGGATGAAGATTACTCTCTTATTTATAATTTTCTGGTTAATGCAGGTTGCCAACAAGTCGAATTAGCCGAACAAGTGGATGCGATAGt                        >  1:138753/1‑236 (MQ=255)                                                                                                                                                                                                                                                 aaTTAAACAATGCGATGGATGGTATATTTACGAAAAAATAATTGATGATAATTTTGCCATAGTCATAGAATCTTTGGGCTTTGCGCTTAAAATCTATGGAGGGGATGAACGTTTTCGTAATGGTTCCTCTGTTGTTTTGGAGGATGAAGATTACTCTCTTATTTATAATTTTCTGGTTAATGCAGGTTGCCAACAAGTCGAATTAGCCGAACAAGTGGATGCGATAGTAAGTGCGAATCTGGCCGCAGAt  >  1:200443/1‑250 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   CATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGCGATTGGTGACGTGTCTTTAGGTTATGCCTTATTGGAGGATATAGTAAGCAAATTAAAACGGGGTAAATTTGAACTCCAAATACCAGGAGGTGGAATTAAAGAATGCGATGGATGGTATATTTACGAAAAAATAATTGATGATAATTTTGCCATAGTCATAGAATCTTTGGGCTTTGCGCTTAAAATCTATGGAGGGGATGAACGTTTTCGTAATGGTTCCTCTGTTGTTTTGGAGGATGAAGATTACTCTCTTATTTATAATTTTCTGGTTAATGCAGGTTGCCAACAAGTCGAATTAGCCGAACAAGTGGATGCGATAGTAAGTGCGAATCTGGCCGCAGAT  >  NC_000913/4438042‑4438531

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