Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	1,849,589	C→T	16.3%	G45D (GGT→GAT)	selD ←	selenophosphate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	1,849,589	0	C	T	16.3%	105.2 / 12.7	43	G45D (GGT→GAT)	selD	selenophosphate synthase
Reads supporting (aligned to +/- strand):  ref base C (20/16);  new base T (3/4);  total (23/20)
Fisher's exact test for biased strand distribution p-value = 6.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.98e-01

ACTGATAACGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGG  >  W3110S.gb/1849528‑1849657                                                              |                                                                     aCTGATAACGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTATCCACAAGc                                                              >  1:3469784/1‑70 (MQ=255)  cTGATAACGCTGTTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCa                                                             <  1:2611223/70‑1 (MQ=255)  cTGATAACGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCa                                                             <  1:1324259/70‑1 (MQ=255)    gATAACGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAAt                                                          >  1:3735907/1‑71 (MQ=255)    gATAACGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAAt                                                          >  1:2599952/1‑71 (MQ=255)    gATAACGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAAt                                                          >  1:1870473/1‑71 (MQ=255)    gATAACGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAAt                                                          >  1:2110388/1‑71 (MQ=255)     aTAACGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAAtt                                                         >  1:1357323/1‑71 (MQ=255)        aCGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCg                                                       >  1:3776721/1‑70 (MQ=255)        aCGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCg                                                       >  1:2245858/1‑70 (MQ=255)        aCGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCg                                                       >  1:2191422/1‑70 (MQ=255)               tGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTAc                                                                      <  1:870472/48‑1 (MQ=255)               tGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTAc                                                                      <  1:1226646/48‑1 (MQ=255)               tGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATcaaca                                               <  1:2168034/71‑1 (MQ=255)                   aTTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACt                                           >  1:3308196/1‑71 (MQ=255)                      gCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTc                                         <  1:1885327/70‑1 (MQ=255)                      gCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTc                                         <  1:1160148/70‑1 (MQ=255)                      gCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCg                                        >  1:3758998/1‑71 (MQ=255)                      gCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCg                                        >  1:2682527/1‑71 (MQ=255)                      gCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCg                                        >  1:4046393/1‑71 (MQ=255)                       cccAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCg                                        <  1:2830617/70‑1 (MQ=255)                       cccAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGc                                       <  1:2536083/71‑1 (MQ=255)                            aTCGTACACCGCCGCATCGTCGCGGGTTTCATTATCCACAAGCAAATTCGGATCAACAAACTTc                                         <  1:1304714/64‑1 (MQ=255)                               gTACACCGCCGCATCTTCGCGGTTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCAc                               <  1:3790211/71‑1 (MQ=255)                               gTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCAc                               <  1:2495817/71‑1 (MQ=255)                               gTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCAc                               <  1:3943419/71‑1 (MQ=255)                               gTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCAc                               <  1:1604533/71‑1 (MQ=255)                                  cacCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTAt                            >  1:2409242/1‑71 (MQ=255)                                      gccgcATCGTCGCGGGTTTCATTATCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAg                        <  1:2769388/71‑1 (MQ=255)                                      gccgcATCGTCGCGGGTTTCATTATCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAg                        <  1:2404175/71‑1 (MQ=255)                                      gccgcATCGTCGCGGGTTTCATTATCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAg                        <  1:332956/71‑1 (MQ=255)                                      gccgcATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAg                        <  1:2363722/71‑1 (MQ=255)                                         gcATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGAt                     >  1:1727151/1‑71 (MQ=255)                                              gtcgCGGGTTTCATTATCCACAAGCAAATTCGGATCAACAAACTTCGcc                                      >  1:3735128/1‑49 (MQ=255)                                              gtcgCGGGTTTCATTATCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGa                      >  1:3989290/1‑65 (MQ=255)                                              gtcgCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACt                                           >  1:2276514/1‑44 (MQ=255)                                              gtcgCGCGTTTCATTACCCACAAGCAAATTCGGATCAACAAACtt                                          >  1:692/1‑45 (MQ=255)                                                    ggTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAg                        <  1:3939472/57‑1 (MQ=255)                                                        tCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACttt      >  1:470848/1‑71 (MQ=255)                                                        tCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACttt      >  1:2522941/1‑71 (MQ=255)                                                          aTTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTg    <  1:3711593/71‑1 (MQ=255)                                                          aTTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTg    <  1:3136074/71‑1 (MQ=255)                                                           ttACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCAc                               <  1:3415718/43‑1 (MQ=255)                                                           ttACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTa                             >  1:1568106/1‑45 (MQ=255)                                                             aCCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTggg  >  1:528483/1‑70 (MQ=255)                                                              |                                                                     ACTGATAACGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGG  >  W3110S.gb/1849528‑1849657

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