Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	2,955,260	C→T	8.9%	V761I (GTA→ATA)	ptr ←	protease III

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	2,955,260	0	C	T	8.9%	126.7 / 2.6	45	V761I (GTA→ATA)	ptr	protease III
Reads supporting (aligned to +/- strand):  ref base C (23/18);  new base T (2/2);  total (25/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.97e-01

ACGATCTGCCCCAACAGAGAGCTATAGGCTGAGCTGGTGTATTCATCGTAGCCAGTCGGTACAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCAC  >  W3110S.gb/2955193‑2955329                                                                    |                                                                      aCGATCTGCCCCAACAGAGAGCTATAGGCTGAGCTGGTGTATTCATCGTAGCCAGTCGGTACAAATACCGc                                                                    >  1:464403/1‑71 (MQ=255)  cGATCTGCCCCAACAGAGAGCTATAGGCTGAGCTGGTGTATTCATCGTAGCCAGTCGGTACAAATACCg                                                                     >  1:564916/1‑69 (MQ=255)  cGATCTGCCCCAACAGAGAGCTATAGGCTGAGCTGGTGTATTCATCGTAGCCAGTCGGTACAAATACCGCt                                                                   >  1:3819397/1‑71 (MQ=255)  cGATCTGCCCCAACAGAGAGCTATAGGCTGAGCTGGTGTATTCATCGTAGCCAGTCGGTACAAATACCGCt                                                                   >  1:1131812/1‑71 (MQ=255)  cGATCTGCCCCAACAGAGAGCTATAGGCTGAGCTGGTGTATTCATCGTAGCCAGTCGGTACAAATACCGCt                                                                   >  1:2523620/1‑71 (MQ=255)  cGATCTGCCCCAACAGAGAGCTATAGGCTGAGCTGGTGTATTCATCGTAGCCAGTCGGTACAAATACCGCt                                                                   >  1:2385157/1‑71 (MQ=255)         ccccAACAGAGAGCTATAGGCTGAGCTGGTGTATTCATCGTAGCCAGTCGGTACAAATACCGCTGc                                                                 <  1:511190/66‑1 (MQ=255)           ccAACAGAGAGCTATAGGCTGAGCTGGTGTATTCATCGTAGCCAGTCGGTACAAATACCGCTGCCAGTGCg                                                          <  1:2141905/71‑1 (MQ=255)            cAACAGAGAGCTATAGGCTGAGCTGGTGTATTCATCGTAGCCAGTCGGTACAAATACCGCTGCCAGTGCgg                                                         <  1:2047325/71‑1 (MQ=255)             aaCAGAGAGCTATAGGCTGAGCTGGTGTATTCATCGTAGCCAGTCGGTACAAATACCGCTGCCAGTGCGGa                                                        <  1:1895409/71‑1 (MQ=255)                agagagCTATAGGCTGAGCTGGTGTATTCATCGTAGCCAGTCGGTACAAATACCGCTGCCAGTGCGGAGTc                                                     <  1:235474/71‑1 (MQ=255)                 gagagCTATAGGCTGAGCTGGTGTATTCATCGTAGCCAGTCGGTACAAATACCGCTGCCAGTGCGGAGTCg                                                    >  1:3410565/1‑71 (MQ=255)                 gagagCTATAGGCTGAGCTGGTGTATTCATCGTAGCCAGTCGGTACAAATACCGCTGCCAGTGCGGAGTCg                                                    >  1:3415063/1‑71 (MQ=255)                 gagagCTATAGGCTGAGCTGGTGTATTCATCGTAGCCAGTCGGTACAAATACCGCTGCCAGTGCGGAGTCg                                                    >  1:1202827/1‑71 (MQ=255)                 gagagCTATAGGCTGAGCTGGTGTATTCATCGTAGCCAGTCGGTACAAATACCGCTGCCAGTGCGGAGTCg                                                    >  1:1838332/1‑71 (MQ=255)                   tagCTATAGGCTGAGCTGGTGTATTCATCGTAGCCAGTCGGTACAAATACCGCTGCCAGTGCGGAGTCGGt                                                  >  1:180910/2‑71 (MQ=255)                    agCTATAGGCTGAGCTGGTGTATTCATCGTAGCCAGTCGGTACAAATACCGCTGCCAGTGCGGAGTCGGTg                                                 >  1:2458069/1‑71 (MQ=255)                      cTATAGGCTGAGCTGGTGTATTCATCGTAGCCAGTCGGTACAAATACCGCTGCAAGTGCGGAGTCGGTGCt                                               <  1:252779/71‑1 (MQ=255)                          aGGCTGAGCTGGTGTATTCATCGTAGCCAGTCGGTACAAATACCGCTGCCAGTGCGGAGTCGGTGCTGtt                                            <  1:612298/70‑1 (MQ=255)                          aGGCTGAGCTGGTGTATTCATCGTAGCAAGTCGGTACAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTa                                           >  1:2829054/1‑71 (MQ=255)                           ggCTGAGCTGGTGTATTCATCGTAGCCAGTCGGTACAAATACCGCTGCCa                                                               >  1:1108462/1‑50 (MQ=255)                           ggCTGAGCTGGTGTATTCATCGTAGCCAGTCGGTACAAATACCGCTGCCAg                                                              >  1:2772441/1‑51 (MQ=255)                           ggCTGAGCTGGTGTATTCATCGTAGCCAGTCGGTACAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTAc                                          >  1:1279925/1‑71 (MQ=255)                                    ggTGTATTCATCGTAGCCAGTCGGTACAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCttttt                                 >  1:2634083/1‑71 (MQ=255)                                        tATTCATCGTAGCCAGTCGGTACAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCaaa                             <  1:3589634/71‑1 (MQ=255)                                          ttCATCGTAGCCAGTCGGTACAAATACCGCTGCCa                                                               >  1:473747/1‑35 (MQ=255)                                          ttCATCGTAGCCAGTCGGTACAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCaaa                             <  1:3496313/69‑1 (MQ=255)                                            cATCGTAGCCAGTCGGTACAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGAt                          <  1:585794/70‑1 (MQ=255)                                            cATCGTAGCCAGTCGGTACAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGAt                          <  1:729317/70‑1 (MQ=255)                                            cATCGTAGCCAGTCGGTACAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGAt                          <  1:1054435/70‑1 (MQ=255)                                            cATCGTAGCCAGTCGGTACAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGAt                          <  1:3342925/70‑1 (MQ=255)                                            cATCGTAGCCAGTCGGTACAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGAt                          <  1:2812270/70‑1 (MQ=255)                                             aTCGTAGCCAGTCGGTACAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGATg                         <  1:265535/70‑1 (MQ=255)                                             aTCGTAGCCAGTCGGTACAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGATGa                        >  1:286390/1‑71 (MQ=255)                                              tCGTAGCCAGTCGGTACAAATACCGCTGCCAGTGCGGAGTCgg                                                   <  1:1278867/43‑1 (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:3371704/1‑71 (MQ=255)                                                      aGTCGGTACAAATATCGCTGCCAGTGCGGAGTCGGTGCTGTTACCAGCTTTTTCAAAGATGACGGATTGtt               >  1:505570/1‑71 (MQ=255)                                                        tCGGTACAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGtttt             >  1:3444058/1‑71 (MQ=255)                                                              cAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATc        >  1:3149756/1‑70 (MQ=255)                                                              cAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATc        <  1:1020292/70‑1 (MQ=255)                                                              cAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCg       <  1:717820/71‑1 (MQ=255)                                                                    ccGCTGCCAGTGCGGAGTCTGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCAc  <  1:865504/70‑1 (MQ=255)                                                                    ccGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCAc  <  1:285899/70‑1 (MQ=255)                                                                    ccGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCAc  >  1:2167982/1‑70 (MQ=255)                                                                    |                                                                      ACGATCTGCCCCAACAGAGAGCTATAGGCTGAGCTGGTGTATTCATCGTAGCCAGTCGGTACAAATACCGCTGCCAGTGCGGAGTCGGTGCTGTTACCGGCTTTTTCAAAGATGACGGATTGTTTTTTATCGACCAC  >  W3110S.gb/2955193‑2955329

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: ‑