Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	2,955,291	G→A	9.5%	A750A (GCC→GCT)	ptr ←	protease III

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	2,955,291	0	G	A	9.5%	115.5 / 2.8	42	A750A (GCC→GCT)	ptr	protease III
Reads supporting (aligned to +/- strand):  ref base G (13/25);  new base A (2/2);  total (15/27)
Fisher's exact test for biased strand distribution p-value = 6.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

GGTGTATTCATCGTAGCCAGTCGGTACAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCACTACATCTTTGTTTCGACACCACTCTGAA  >  W3110S.gb/2955228‑2955357                                                                |                                                                   ggTGTATTCATCGTAGCCAGTCGGTACAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCttttt                                                             >  1:2634083/1‑71 (MQ=255)     tATTCATCGTAGCCAGTCGGTACAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCaaa                                                         <  1:3589634/71‑1 (MQ=255)       ttCATCGTAGCCAGTCGGTACAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCaaa                                                         <  1:3496313/69‑1 (MQ=255)         cATCGTAGCCAGTCGGTACAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGAt                                                      <  1:1054435/70‑1 (MQ=255)         cATCGTAGCCAGTCGGTACAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGAt                                                      <  1:2812270/70‑1 (MQ=255)         cATCGTAGCCAGTCGGTACAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGAt                                                      <  1:729317/70‑1 (MQ=255)         cATCGTAGCCAGTCGGTACAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGAt                                                      <  1:3342925/70‑1 (MQ=255)         cATCGTAGCCAGTCGGTACAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGAt                                                      <  1:585794/70‑1 (MQ=255)          aTCGTAGCCAGTCGGTACAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGATg                                                     <  1:265535/70‑1 (MQ=255)          aTCGTAGCCAGTCGGTACAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGATGa                                                    >  1:286390/1‑71 (MQ=255)            cGTAGCCAGTCGGTACAAATATCGCTGCCAGTGCGGAGTCGGTGCTGTTACCAGCTTTTTCAAAGATg                                                     <  1:3731385/68‑1 (MQ=255)            cGTAGCCAGTCGGTACAAATATCGCTGCCAGTGCGGAGTCGGTGCTGTTACCAGCTTTTTCAAAGATg                                                     <  1:3169414/68‑1 (MQ=255)             gTAGCCAGTCGGTACAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGATGa                                                    >  1:2407925/1‑68 (MQ=255)                  cAGTCGGTACAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGt                                            >  1:1116246/1‑71 (MQ=255)                   aGTCGGTACAAATATCGCTGCCAGTGCGGAGTCGGTGCTGTTACCAGCTTTTTCAAAGATGACGGATTGtt                                           >  1:505570/1‑71 (MQ=255)                   aGTCGGTACAAATATCGCTGCCAGTGCGGAGTCGGTGCTGTTACCAGCTTTTTCAAAGATGACGGATTGtt                                           >  1:3371704/1‑71 (MQ=255)                     tCGGTACAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGtttt                                         >  1:3444058/1‑71 (MQ=255)                           cAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATc                                    <  1:1020292/70‑1 (MQ=255)                           cAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATc                                    >  1:3149756/1‑70 (MQ=255)                           cAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCg                                   <  1:717820/71‑1 (MQ=255)                                 ccGCTGCCAGTGCGGAGTCTGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCAc                              <  1:865504/70‑1 (MQ=255)                                 ccGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCAc                              >  1:2167982/1‑70 (MQ=255)                                 ccGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCAc                              <  1:285899/70‑1 (MQ=255)                                     tGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCACTa                            <  1:3199137/68‑1 (MQ=255)                                     tGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCACTa                            <  1:930608/68‑1 (MQ=255)                                     tGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCACTa                            <  1:2803858/68‑1 (MQ=255)                                       ccAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGtttttt                                       <  1:1280645/55‑1 (MQ=255)                                        cAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGAt                                                      <  1:2510762/39‑1 (MQ=255)                                        cAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCACTACATCtt                      <  1:2032910/71‑1 (MQ=255)                                              ggAGTCGGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCACTACATCTTTGTTTc                <  1:2225264/71‑1 (MQ=255)                                              ggAGTCGGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCACTACATCTTTGTTTc                <  1:1777816/71‑1 (MQ=255)                                                 gTCGGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCACt                             <  1:1929187/55‑1 (MQ=255)                                                  tCGGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGAc                                 <  1:1356620/50‑1 (MQ=255)                                                    ggTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCACTACATCTTTGTTTCGacac           >  1:1569699/1‑70 (MQ=255)                                                       gCTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCACTACATc                        >  1:287742/1‑54 (MQ=255)                                                        cTGTTACCGGCTTTTTCAAAGATGACGGATTGttttt                                        >  1:2130675/1‑37 (MQ=255)                                                          gTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATc                                    >  1:1788958/1‑39 (MQ=255)                                                           ttACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCACTACATCTTTGtt                  <  1:3727523/56‑1 (MQ=255)                                                           ttACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCACTACATCTTTGTTTCGACACCACTCTg    <  1:2589661/70‑1 (MQ=255)                                                            tACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCACTACATCTTTGTTTCGACACCACTCTGaa  >  1:2300334/1‑71 (MQ=255)                                                             aCCGGCTTTTTCAAAGATGACGGATTGTTTTTTATc                                    >  1:2804625/1‑36 (MQ=255)                                                               cGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCACTACATCTTTGTTTc                <  1:1403911/54‑1 (MQ=255)                                                                ggCTTTTTCAAAGATGACGGATTGTTTTTTATCGAc                                 <  1:3441592/36‑1 (MQ=255)                                                                |                                                                   GGTGTATTCATCGTAGCCAGTCGGTACAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCACTACATCTTTGTTTCGACACCACTCTGAA  >  W3110S.gb/2955228‑2955357

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