Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,426,403 G→A 29.7% G113R (GGG→AGG)  yfcG → predicted glutathione S‑transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,426,4030GA29.7% 51.9 / 26.3 37G113R (GGG→AGG) yfcGpredicted glutathione S‑transferase
Reads supporting (aligned to +/- strand):  ref base G (10/16);  new base A (3/8);  total (13/24)
Fisher's exact test for biased strand distribution p-value = 7.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.26e-01

CGTGAGCGCGCCGCCACATTACAGTGGTTATTCTGGCAGGTAGGCGGACTGGGGCCGATGCTTGGGCAAAATCATCATTTTAATCACGCAGCCCCCCAAACCATTCCTTACGCTATTGAACGTT  >  W3110S.gb/2426340‑2426463
                                                               |                                                            
cgTGAGCGCGCCGCCACATTACAGTGGTTATTCTGGCAGGTAGGCGGACTGGGGCCGATGCTTGGGCaaaa                                                       >  1:3838358/1‑71 (MQ=255)
      cgcgcCGCCACATTACAGTGGTTATTCTGGCAGGTAGGCGGACTGGGGCCGATGCTTGGGCAAAAtcatca                                                 <  1:3584794/71‑1 (MQ=255)
               acaTTACAGTGGTTATTCTGGCAGGTAGGCGGACTGTGGCCGATGCTTGGGCAAAATCATCATTTTAATCa                                        >  1:1422128/1‑71 (MQ=255)
                caTTACAGTGGTTATTCTGGCAGGTAGGCGGACTGGGGCCGATGCTTAGGCAAAatcat                                                   <  1:1212634/59‑1 (MQ=255)
                caTTACAGTGGTTATTCTGGCAGGTAGGCGGACTGGGGCCGATGCTTAGGCAAAatcat                                                   <  1:3678271/59‑1 (MQ=255)
                caTTACAGTGGTTATTCTGGCAGGTAGGCGGACTGGGGCCGATGCTTAGGCAAAatcat                                                   <  1:2147586/59‑1 (MQ=255)
                 aTTACAGTGGTTATTCTGGCAGGTAGGCGGACTGGGGCCGATGCTTGGGCAAAATCATCATTTTAATCACg                                      <  1:556879/71‑1 (MQ=255)
                  ttACAGTGGTTATTCTGGCAGGTAGGCGGACTGGGGCCGATGCTTGGGCAAAATCATCATTTTAATCACg                                      <  1:1662117/70‑1 (MQ=255)
                  ttACAGTGGTTATTCTGGCAGGTAGGCGGACTGGGGCCGATGCTTGGGCAAAATCATCATTTTAATCACg                                      <  1:2891338/70‑1 (MQ=255)
                  ttACAGTGGTTATTCTGGCAGGTAGGCGGACTGGGGCCGATGCTTGGGCAAAATCATCATTTTAATCACg                                      <  1:2186332/70‑1 (MQ=255)
                  ttACAGTGGTTATTCTGGCAGGTAGGCGGACTGGGGCCGATGCTTGGGCAAAATCATCATTTTAATCACGc                                     <  1:2775449/71‑1 (MQ=255)
                  ttACAGTGGTTATTCTGGCAGGTAGGCGGACTGGGGCCGATGCTTGGGCAAAATCATCATTTTAATCACGc                                     <  1:3813717/71‑1 (MQ=255)
                    aCAGTGGTTATTCTGGCAGGTAGGCGGACTGGGGCCGATGCTTGGGCAAAATCATCATTTTAATCAc                                       >  1:936419/1‑67 (MQ=255)
                    aCAGTGGTTATTCTGGCAGGTAGGCGGACTGGGGCCGATGCTTAGGCAAAATCATCATTTTAATCa                                        >  1:1552744/1‑66 (MQ=255)
                          gTTATTCTGGCAGGTAGGCGGACTGGGGCCGATGCTTGGGCAAAATCATCATTTTAATCACg                                      <  1:3469269/62‑1 (MQ=255)
                             aTTCTGGCAGGTAGGCGGACTGGGGCCGATGCTTGGGCAAAATCATCATTTTAATCACGCa                                    >  1:2100114/1‑61 (MQ=255)
                             aTTCTGGCAGGTAGGCGGACTGGGGCCGATGCTTGGGCAAAATCATCATTTTAATCACGCa                                    >  1:1761570/1‑61 (MQ=255)
                                cTGGCAGGTAGGCGGACTGGGGCCGATGCTTGGGCAAAATCATCATTTGATCACGCAGCCCCCCAAACCa                       <  1:2337333/70‑1 (MQ=255)
                                     aGGTAGGCGGACTGGGGCCGATGCTTAGGCAAAATCATCATTTTAATCACGCAGCCCCCCAAACCATTCCt                  <  1:2289056/71‑1 (MQ=255)
                                     aGGTAGGCGGACTGGGGCCGATGCTTAGGCAAAATCATCATTTTAATCACGCAGCCCCCCAAACCATTCCt                  <  1:803593/71‑1 (MQ=255)
                                     aGGTAGGCGGACTGGGGCCGATGCTTAGGCAAAATCATCATTTTAATCACGCAGCCCCCCAAACCATTCCt                  <  1:4059040/71‑1 (MQ=255)
                                     aGGTAGGCGGACTGGGGCCGATGCTTAGGCAAAATCATCATTTTAATCACGCAGCCCCCCAAACCATTCCt                  <  1:1732027/71‑1 (MQ=255)
                                       gTAGGCGGACTGGGGCCGATGCTTGGGCAAAATCATCAtttt                                             <  1:2080367/42‑1 (MQ=255)
                                       gTAGGCGGACTGGGGCCGATGCTTGGGCAAAATCATCAtttt                                             <  1:1753479/42‑1 (MQ=255)
                                       gTAGGCGGACTGGGGCCGATGCTTGGGCAAAATCATCAtttt                                             <  1:1998405/42‑1 (MQ=255)
                                        tAGGCGGACTGGGGCCGATGCTTAGGCAAAATCATCATTTTAATCACGCAGCCCCCCAAACCATTCCTTAc               >  1:2463685/1‑71 (MQ=255)
                                         aGGCGGACTGGGGCCGATGCTTGGGCAAAATCATCATTTTAATCACg                                      >  1:1185060/1‑47 (MQ=255)
                                         aGGCGGACTGGGGCCGATGCTTGGGCAAAATCATCATTTTAATCACGCAGCCCCCCAAACCATTCCTTAc               >  1:3760741/1‑70 (MQ=255)
                                            cGGACTGGGGCCGATGCTTAGGCAAAATCATCATTTTAATCACGCAGCCCCCCAAACCATTcc                   >  1:577548/1‑63 (MQ=255)
                                             ggACTGGGGCCGATGCTTGGGCAAAATCATCTTTTTAATCACGCAGCCCCCAAAACCATTCCTTACGCTAt          >  1:3228164/1‑71 (MQ=255)
                                             ggACTGGGGCCGATGCTTGGGCAAAATCATCATTTTAATCACGCAGCCCCCCAAACCATTCCTTACGCTAt          >  1:3845832/1‑71 (MQ=255)
                                             ggACTGGGGCCGATGCTTGGGCAAAATCATCATTTTAATCACGCAGCCCCCCAAACCATTCCTTACGCTAt          >  1:1386256/1‑71 (MQ=255)
                                                    cgCCGATGCTTGGGCAAAATCATCATTTTAATCACGCAGCCCCCCAAACCATTCCTTACGCTATTGAACGt   <  1:3932028/70‑1 (MQ=255)
                                                     gCCGATGCTTGGGCAAAATCATCATTTTAATCACGCAGCCCCCCAAACCATTCCTTACGCTATTGAACGtt  <  1:1751547/71‑1 (MQ=255)
                                                     gCCGATGCTTGGGCAAAATCATCATTTTAATCACGCAGCCCCCCAAACCATTCCTTACGCTATTGAACGtt  <  1:3692163/71‑1 (MQ=255)
                                                     gCCGATGCTTGGGCAAAATCATCATTTTAATCACGCAGCCCCCCAAACCATTCCTTACGCTATTGAACGtt  <  1:502717/71‑1 (MQ=255)
                                                     gCCGATGCTTAGGCAAAATCATGATTTTAATCACGCAGCCCACCAAACCATTCCTTACGCTATTGAACGtt  <  1:3972325/71‑1 (MQ=255)
                                                               |                                                            
CGTGAGCGCGCCGCCACATTACAGTGGTTATTCTGGCAGGTAGGCGGACTGGGGCCGATGCTTGGGCAAAATCATCATTTTAATCACGCAGCCCCCCAAACCATTCCTTACGCTATTGAACGTT  >  W3110S.gb/2426340‑2426463

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: