Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	120,663	C→T	63.8%	C181Y (TGC→TAC)	speE ←	spermidine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	120,663	0	C	T	63.8%	40.7 / 46.5	47	C181Y (TGC→TAC)	speE	spermidine synthase
Reads supporting (aligned to +/- strand):  ref base C (7/10);  new base T (12/18);  total (19/28)
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.76e-01

ATGGCTTCTTCCTGCTGTAAAAAGCAGACGCCGTTTTGTGCGACGAAGATACCGCCAGGATTCAGGCAACGTTTGCAGCCTTCATAAAATGCCGAAGTGAAAAGGCTTTCGCCGGGACCGATAGGATCGGT  >  Exported/120597‑120727                                                                   |                                                                 aTGGCTTCTTCCTGCTGTAAAAAGCAGACGCCGTTTTGTGCGACGAAGATACCGCCAGGATTCAGGTAACg                                                              <  1:2561606/71‑1 (MQ=255) aTGGCTTCTTCCTGCTGTAAAAAGCAGACGCCGTTTTGTGCGACGAAGATACCGCCAGGATTCAGGTAACg                                                              >  1:1785477/1‑71 (MQ=255)    gCTTCTTCCTGCTGTAAAAAGCAGACGCCGTTTTGTGCGACGAAGATACCGCCAGGATTCAGGTAACGttt                                                           >  1:2838109/1‑71 (MQ=255)    gCTTCTTCCTGCTGTAAAAAGCAGACGCCGTTTTGTGCGACGAAGATACCGCCAGGATTCAGGTAACGtt                                                            <  1:2644060/70‑1 (MQ=255)    gCTTCTTCCTGCTGTAAAAAGCAGACGCCGTTTTGTGCGACGAAGATACCGCCAGGATTCAGGCAACGttt                                                           >  1:605565/1‑71 (MQ=255)    gCTTCTTCCTGCTGTAAAAAGCAGACGCCGTTTTGTGCGACGAAGATACCGCCAGGATTCAGGCAACGtt                                                            <  1:2052688/70‑1 (MQ=255)    gCTTCTTCCTGCTGTAAAAAGCAGACGCCGTTTTGTGCGACGAAGATACCGCCAGGATTCAGGCAACGtt                                                            <  1:1296774/70‑1 (MQ=255)           ccTGCTGTAAAAAGCAGACGCCGTTTTGTGCGACGAAGATACCGCCAGGATTCAGGTAACGTTTGCAGcc                                                     <  1:1135895/70‑1 (MQ=255)              gctgTAAAAAGCAGACGCCGTTTTGTGCGACGAAGATACCGCCAGGATTCAGGTAACGTTTGCAGCCTTc                                                  >  1:3289983/1‑70 (MQ=255)                  tAAAAAGCAGACGCCGTTTTGTGCGACGAAGATACCGCCAGGATTCAGGCAACGTTTGCAGCCTTCATaa                                              >  1:546930/1‑70 (MQ=255)                   aaaaaGCAGACGCCGTTTTGTGCGACGAAGATACCGCCAGGATTCAGGTAACGTTTGCAGCCTTCATaaa                                             >  1:1955669/1‑70 (MQ=255)                   aaaaaGCAGACGCCGTTTTGTGCGACGAAGATACCG‑CAGGATTCAGGTAACGTTTGCAGCCTTCATAAAAt                                           >  1:201481/1‑71 (MQ=255)                    aaaaGCAGACGCCGTTTTGTGCGACGAAGATACCGCCAGGATTCAGGTAACGTTTGCAGCCTTCATaaaa                                            <  1:1991168/70‑1 (MQ=255)                    aaaaGCAGACGCCGTTTTGTGCGACGAAGATACCGCCAGGATTCAGGCAACGTTTGCAGc                                                      >  1:429238/1‑60 (MQ=255)                       aGCAGACGCCGTTTTGTGCGACGAAGATACCGCCAGGATTCAGGTAACGTTTGCAGCCTTCATAAAATGcc                                        >  1:1690835/1‑71 (MQ=255)                       aGCAGACGCCGTTTTGTGCGACGAAGATACCGCCAGGATTCAGGCAACGTTTGCAGCCTTCATAAAATGcc                                        <  1:2990865/71‑1 (MQ=255)                                cGTTTTGTGCGACGAAGATACCGCCAGGATTCAGGTAACGTTTGCAGCCTTCATAAAATGCCGAAGTGaaa                               <  1:2781220/71‑1 (MQ=255)                                cGTTTTGTGCGACGAAGATACCGCCAGGATTCAGGTAACGTTTGCAGCCTTCATAAAATGCCGAAGTGaaa                               <  1:1135709/71‑1 (MQ=255)                                cGTTTTGTGCGACGAAGATACCGCCAGGATTCAGGTAACGTTTGCAGCCTTCATAAAATGCCGAAGTGaaa                               <  1:2372098/71‑1 (MQ=255)                                cGTTTTGTGCGACGAAGATACCGCCAGGATTCAGGCAACGTTTGCAGCCTTCa                                                 <  1:507915/53‑1 (MQ=255)                                  ttttGTGCGACGAAGATACCGCCAGGATTCAGGCAACGTTTGCAGCCTTCATAAAATGCCGAAGTGAAAAg                             <  1:1154960/71‑1 (MQ=255)                                    ttGTGCGACGAAGATACCGCCAGGATTCAGGCAACGTTTGCAGCCTTCATAAAATGCCGa                                      <  1:1203568/60‑1 (MQ=255)                                    ttGTGCCACGAAGATACCGCCAGGATTCAGGTAACGTTTGCAGCCtt                                                   >  1:1315297/1‑47 (MQ=255)                                    ttGTGCCACGAAGATACCGCCAGGATTCAGGTAACGTTCGCAGCCtt                                                   >  1:3420379/1‑47 (MQ=255)                                       tgCGACGAAGATACCGCCAGGATTCAGGTATCGtt                                                            <  1:2037589/35‑1 (MQ=38)                                         cgacgaAGATACCGCCAGGATTCAGGTAACGTTTGCAGCCTTCATAAAATGCCGAAGTGAAAAGGCTTTc                       <  1:977004/70‑1 (MQ=255)                                         cgacgaAGATACCGCCAGGATTCAGGTAACGTTTGCAGCCTTCATAAAATGCCGAAGTGAAAAGGCTTTc                       <  1:564089/70‑1 (MQ=255)                                          gacgaAGATACCGCCAGGATTCAGGTAACGTTTGCAGCCTTCATAAAATGCCGAAGTGAAAAGGCTTTCGc                     >  1:3651902/1‑71 (MQ=255)                                           acgaAGATACCGCCAGGATTCAGGTAACGTTTGCAGCCTTCATAAAATGCCGaa                                     <  1:555297/54‑1 (MQ=255)                                           acgaAGATACCGCCAGGATTCAGGTAACGTTTGCAGCCTTCATAAAATGCCGaa                                     <  1:2330267/54‑1 (MQ=255)                                            cgaAGATACCGCCAGGATTCAGGCAACGTTTGCAGCCTTCATAAAATGCCGAAGTGAAAAGGCTTTCGCCg                   >  1:59149/1‑71 (MQ=255)                                             gaAGATACCGCCAGGATTCAGGTAACGTTTGCAGCCTTCATAAAATGCCGAAGt                                   >  1:1716872/1‑54 (MQ=255)                                              aaGATACCGCCAGGATTCAGGCAACGTTTGCAGCCTTCATAAAATGCCGAAGTGAAAAGGCTTTCGCCggg                 <  1:2032362/71‑1 (MQ=255)                                                 aTACCGCCAGGATTCAGGTAACGTTTGCAGCCTTCATAAAATGCCGAAGTGAAAAGGCTTTCGCCGGGAcc              <  1:1454399/71‑1 (MQ=255)                                                 aTACCGCCAGGATTCAGGTAACGTTTGCAGCCTTCATAAAATGCCGAAGTGAAAAGGCTTTCGCCGGGAcc              >  1:1357415/1‑71 (MQ=255)                                                 aTACCGCCAGGATTCAGGTAACGTTTGCAGCCTTCATAAAATGCCGAAGTGAAAAGGCTTTCGCCGGGAcc              <  1:1323189/71‑1 (MQ=255)                                                 aTACCGCCAGGATTCAGGCAACGTTTGCAGCCTTCATAAAATGCCGAAGTGAAAAGGCTTTCGCCGGGAcc              <  1:3444552/71‑1 (MQ=255)                                                  tACCGCCAGGATTCAGGTAACGTTTGCAGCCTTCATAAAATGCCGAAGTGAAAAGGCTTTCGCCGGGACCg             <  1:3568590/71‑1 (MQ=255)                                                   aCCGCCAGGATTCAGGTAACGTTTGCAGCCTTCATAAAATGCCg                                       <  1:638743/44‑1 (MQ=255)                                                     cGCCAGGATTCAGGCAACGTTTGCAGCCTTCATAAAATGCCGAAGTg                                  <  1:81256/47‑1 (MQ=255)                                                       ccAGGATTCAGGCAACGTTTGCAGCCTTCATAAAATGCCGAAGTGAAAAGGCTTTCGCCGGGACCGATAg         >  1:2604384/1‑70 (MQ=255)                                                          ggATTCAGGTAACGTTTGCAGCCTTCATAAAATGCCGAAGTGAAAAGGCt                          <  1:22118/50‑1 (MQ=255)                                                          ggATTCAGGCAACGTTTGCAGCCTTCATAAAATGCCGAAGTGAAAAGGCTTTCGCCGGGACCGATAGGATc     >  1:1708975/1‑71 (MQ=255)                                                           gATTCAGGCAACGTTTGCAGCCTTCATAAAATGCCGAAGTg                                  <  1:2236110/41‑1 (MQ=255)                                                            aTTCAGGTAACGTTTGCAGCCTTCATAAAATGCCGAAGTGAAAAGGCTTTCGCCGGGACCGATAGGATCgg   <  1:3373779/71‑1 (MQ=255)                                                             ttCAGGTAACGTTTGCAGCCTTCATAAAATGCCGAAGTGAAAAGGCTTTCGCCGGGACCGATAGGATCGGt  >  1:1836954/1‑71 (MQ=255)                                                               cAGGCAACGTTTGCAGCCTTCATAAAATGCCGAAGTGAAAAGGCTTTCGCCg                   >  1:2619897/1‑52 (MQ=255)                                                                   |                                                                 ATGGCTTCTTCCTGCTGTAAAAAGCAGACGCCGTTTTGTGCGACGAAGATACCGCCAGGATTCAGGCAACGTTTGCAGCCTTCATAAAATGCCGAAGTGAAAAGGCTTTCGCCGGGACCGATAGGATCGGT  >  Exported/120597‑120727

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