Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,358,227	G→A	R321C (CGT→TGT)	puuP ←	putrescine:H(+) symporter PuuP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,358,227	0	G	A	92.3%	25.1 / ‑6.0	13	R321C (CGT→TGT)	puuP	putrescine:H(+) symporter PuuP
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (4/8);  total (4/9)
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

TACGCCGCCAGAAATGATTAAACACCGACAGATTCACGAAGGTAAACGCCACCAGAGCACCGAAGTTAATTAACGCTGTGGCGGTAACTAAATCGAAGAACAGCGCCGACAACGCGACAATCCCGACCATAATGACGTTGAGTGCCGGAGTCCGCCATTTTGGGTGCACATAGCCAAACACGCGCTCCGGAAACACATTGTCACGCCCCATCACATACAGCAGACGCGACACGCTGGCATGTGAAGCCAGGCCAGAGGCTAACGTGTTCACAAACGTGGTGCAGAGGAAAATCGACTGGAAGAGCTTGCCGCCGACGTACAGCGCAATTTCAGGCAGTGCGGCATCCGGATCTTTAAAGCGGCTGATATCGGGAAAGAACAGCTGCATAAAGAACGACGCCGCGATAAAGATAACGC  >  NC_000913/1358023‑1358439                                                                                                                                                                                                             |                                                                                                                                                                                                                     tACGCCGCCAGAAATGATTAAACACCGACAGATTCACGAAGGTAAACGCCACCAGAGCACCGAATTTAATTAACCCTGTGGCGGTAACGAAATCAAAGAACAGCGCCGACAACGCGACAATCCCGACCATAATGACGTTGATGGCCGGAGGCCGCCATTTTGGGGGCACATTGCCAAACACGCCCTCCGGTAACACATTAGCACACCCCATCACAAACAGCAGACGCGACACGCTGGCATGTGAAgccag                                                                                                                                                                         <  1:293627/250‑1 (MQ=255)                 aTTAAACACCGACAGATTCACGAATGTAAACGCCACCATAGCACCGAAGTTAATTAACGCTGTGGCGGTAACTAAATCGAAGAACATCGCCGACAACGCGACAATCCCGACCATAATTACGTTGAGTGCCGGAGTCCGCCATTTTGGGTGCACATAGCCATACACGCGCTCCTGAATTTCATTGTTACACCCCATCACATACAGCAGACGCGACACGCTGGCATGTGAAGCCAGGCCAGAGGCTAACgtgt                                                                                                                                                        <  1:304564/251‑1 (MQ=255)                    aaaCACCGACAGATGCACGAAGGTAAACGCCACCAGAGCACCGAAGTTAATGAACGCTGTGGCGGTAACTAAATCGAAGAACAGCGCCGACAACGCGACAATCCCGACCATAATGACGTTGAGTGCCGGAGTCCGCCATTTTGGGTTCACATAGCCAAACACGCGCTCCGGAAACACATTGTCACACCCCATCACATACATCAGACGCAACGCGCTGGCATGTGAAGCCAGGCCAGAGGCTAACGTGttn                                                                                                                                                      <  1:62227/250‑2 (MQ=255)                          cGACAGATTCACGAAGGTAAACGCCACCAGAGCACCGAAGTTAATTAACGCTGTGGCGGTAACTAAATCGAAGAACAGCGCCGACAACGCGACAATCCCGACCATAATGACGTTGAGTGCCGGAGTCCGCCATTTTGGGTGCACATAGCCAAACACGCGCTCCGGAAACACATTGTCACACCCCATCCCATACAGCAGAc                                                                                                                                                                                                  >  1:88806/1‑200 (MQ=255)                                    aCGAAGGTAAACGCCACCAGAGCACCGAAGTTAATTAACGCTGTGGCGGTAACTAAATCGAAGAACAGCGCCGACAACGCGACAATCCCGACCATAATGACGTTGAGTGCCGGAGTCCGCCATTTTGGGTGCACATAGCCAAACACGCGCTCCGGAAACACATTGTCACACCCCATCACATACAGCAGACGCGACACGCTGGCATGTGAAGCCAGGCCAGAGGCTAACGTGTTCACAAACGTGGTGCagan                                                                                                                                     <  1:204142/251‑2 (MQ=255)                                      naaGGTAAACGCCACCAGAGCACCGAAGTTAATTAACGCTGTGGCGGTAACTAAATCGAAGAACAGCGCCGACAACGCGACAATCCCGACCATAATGACGTTGAGTGCCGGAGTCCGCCATTTTGGGTGCACATAGCCAAACACGCGCTCCGGAAACACATTGTCACACCCCATCACATACCGCAGACGCGACACGCTGGCATGTGAAGCCAGGCCAGAGGCTAACGGGTTCACaaa                                                                                                                                                 >  1:85438/2‑237 (MQ=255)                                            aaaCGCCACCAGAGCACCGAAGTTAATTAACGCTGAGGCGGAAACTAAATAAAAGAACAGCGCCGCCAACGCGACAAACCCCACCCTAATGACGATGAATGCCGGAGTCCGCCATTTTCGGTGTACCTAGCCACACAAGAGCGCCGGAAACACATTGTCACACACCATCACATACAGCAGAAGCGAACCGCTGGCCTGTGACGACAGGCCAGAGGATAACGTGTTCACAAACGAGGTGCAGAGGaaactc                                                                                                                              >  1:295564/1‑247 (MQ=255)                                                           aCCGAAGTTAATTAACGCTGTGGCGGTAACTAAATCGAAGAACAGCGCCGACAACGCGACAATCCCGACCATAATGACGTTGAGTGCCGGAGTCCGCCATTTTGGGTGCACATAGCCAAACACGCGCTCCGGAAACACATTGTCACACCCCATCGCATACAGCAGACGCGACACGCTGGCATGTGAAGCCAGGCCATAGGCTAACGTGTTCACAAACGTGGTGCAGAGGAAAATCGACTGGAAGAGCTTgc                                                                                                              <  1:147994/251‑1 (MQ=255)                                                                     attaaCGCTGTGGCGGTAACTAAATCGAAGAACAGCGCCGACAACGCGACAATCCCGACCATAATGACGTTGAGTGCCGGAGTCCGCCATTTTGGGTGCACATAGCCAAACACGCGCTCCGGAAACACATTGTCACACCCCATCACATACAGCAGACGCGACACGCTGGCATGTGAAGCCAGGCCAGAggn                                                                                                                                                                <  1:95180/191‑2 (MQ=255)                                                                                                                                             nntGNNNNNGTCNCCCATTTTGGGTGCACATAGCCAAACACGCGCTCCGGAAACACATTGTCACGCCCCATCACATACGGCAGGCACGACACGCCGGCATGTGAAGCCAGGCCAGAGGCTAACGTGTTCACAAACGGGGTGCagag                                                                                                                                     <  2:194738/144‑1 (MQ=255)                                                                                                                                                          gCCATTTTGGGTGCACATAGCCAAACACGCGCTCCGGAAACACATTGTCACACCCCATCACATACAGCAGACGCGACACGCTGGCATGTGAAGCCAGGCCAGAGGCTAACGTGTTCACAAACGTGGTGCAGAGGAAAATCGACTGGAAGAGCTTGCCGCCGACGTACAGCGCAATTTCAGGCAGTGCGGCATCCGGATCTTTAAAGc                                                           <  1:73485/207‑1 (MQ=255)                                                                                                                                                          gCCATTTTGGGTGCACATAGCCAAAANCGCGCTNNGGNAACACATTGTCACACCCCANCANNNANNNCNGACGCNACNNGCTGGCATGNGAAGCCAGGCCAGAGGCTAACGTGTTCACAAACGTGGTGCAGAGGAAAATCGANTGGAAGAGCNNGCCGCCGACGTACAGCGCAATNTCANNCNNTNCGGNAtn                                                                         >  2:73485/1‑192 (MQ=255)                                                                                                                                                                       cacaTAGCCAAACACGCGCTCCGGAAACACATTGTCACACCCCATCACATACAGCAGACGCGACACGCTGGCAGGTGAAGCCAGGCCAGAGGCTAACGTGTTCACAAACGTGGTGCAGAGGAAAATCGACGGGAAGAGCTTGCCGCCGACGTACAGCGCAATTTCAGGCATTGCGGCGTCCGGATCTTTAAAGCGGCGTATATCGGGAAGGAACAGCTGCATAAAGAACGACGCCGCGATAAAGATAACgn  <  1:36513/251‑2 (MQ=255)                                                                                                                                                                                                             |                                                                                                                                                                                                                     TACGCCGCCAGAAATGATTAAACACCGACAGATTCACGAAGGTAAACGCCACCAGAGCACCGAAGTTAATTAACGCTGTGGCGGTAACTAAATCGAAGAACAGCGCCGACAACGCGACAATCCCGACCATAATGACGTTGAGTGCCGGAGTCCGCCATTTTGGGTGCACATAGCCAAACACGCGCTCCGGAAACACATTGTCACGCCCCATCACATACAGCAGACGCGACACGCTGGCATGTGAAGCCAGGCCAGAGGCTAACGTGTTCACAAACGTGGTGCAGAGGAAAATCGACTGGAAGAGCTTGCCGCCGACGTACAGCGCAATTTCAGGCAGTGCGGCATCCGGATCTTTAAAGCGGCTGATATCGGGAAAGAACAGCTGCATAAAGAACGACGCCGCGATAAAGATAACGC  >  NC_000913/1358023‑1358439

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