Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,274,249 A→G 23.5% T64A (ACT→GCT)  rtn → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,274,2490AG23.5% 59.3 / 17.7 34T64A (ACT→GCT) rtnconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (11/15);  new base G (5/3);  total (16/18)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

GTCAAATACGACACACTGATTACCGACGTACAAAAGTATCTCGATACCTATTTTGCCGACCTGAAATCCACTACTGACCGGCTCCAGCCGCTGACCTTAGATACCTGCCAGCAAGCTAACCCCGAACTGACCGC  >  W3110S.gb/2274180‑2274313
                                                                     |                                                                
gTCAAATACGACACACTGATTACCGACGTACAAAAGTATCTCGATACCTATTTTGCCGACCTGAAATCCac                                                                 <  1:2739517/71‑1 (MQ=255)
 tCAAATACGACACACTGATTACAGACGTACAAAAGTATCTCGATACCTATTTTGCCGACCTGAAATCCac                                                                 >  1:3041412/1‑70 (MQ=255)
     aTACGACACACTGATTACCGACGTACAAAAGTATCTCGATACCTATTTTGCCGACCTGAAATCCACTACTg                                                            <  1:2789545/71‑1 (MQ=255)
              acTGATTACCGACGTACAAAAGTATCTCGATACCTATTTTGCCGACCTGAAATCCACTACTGACCGGCTc                                                    >  1:3917373/1‑70 (MQ=255)
                 gATTACCGACGTACAAAAGTATCTCGATACCTATTTTGCCGACCTGAAATCCACTACTGACCGGCTCCAg                                                 >  1:2046083/1‑70 (MQ=255)
                  aTTACCGACGTACAAAAGTATCTCGATACCTATTTTGCCGACCTGAAATCCGCTACTGACCGGCTCCAGc                                                <  1:1566261/70‑1 (MQ=255)
                  aTTACCGACGTACAAAAGTATCTCGATACCTATTTTGCCGACCTGAAATCCGCTACTGACCGGCTCCAGc                                                <  1:1610517/70‑1 (MQ=255)
                          cGTACAAAAGTATCTCGATACCTATTTTGCCGACCTGAAATCCACTACTGACCGGCTCCAGCCGCTGACCt                                       <  1:1062662/71‑1 (MQ=255)
                            tACAAAAGTATCTCGATACCTATTTTGCCGACCTGAAATCCACTACTGACCGGCTCCAGCCGCTGACCTTa                                     <  1:3820087/71‑1 (MQ=255)
                            tACAAAAGTATCTCGATACCTATTTTGCCGACCTGAAATCCACTACTGACCGGCTCCAGCCGCTGACCTTa                                     <  1:2870504/71‑1 (MQ=255)
                            tACAAAAGTATCTCGATACCTATTTTGCCGACCTGAAATCCACTACTGACCGGCTCCAGCCGCTGACCTTa                                     <  1:267159/71‑1 (MQ=255)
                              cAAAAGTATCTCGATACCTATTTTGCCGACCTGAAATCCACTACTGACCGGCTCCAGCCGCTGACCt                                       >  1:1362823/1‑67 (MQ=255)
                              cAAAAGTATCTCGATACCTATTTTGCCGACCTGAAATCCACTACTGACCGGCTCCAGCCGCTGACCt                                       >  1:2736059/1‑67 (MQ=255)
                                      tctcGATACCTATTTTGCCGACCTGAAATCCACTACTGACCGGCTCCAGCCGCTGACCTTa                                     <  1:406589/61‑1 (MQ=255)
                                              ccTATTTTGCCGACCTGAAATCCGCTACTGACCGGCTCCAGCCGCt                                            >  1:1685106/1‑46 (MQ=255)
                                              ccTATTTTGCCGACCTGAAATCCGCTACTGACCGGCTCCAGCCGCTGACCTTAGATACCTGCCAGCAAGc                    >  1:2609098/1‑70 (MQ=255)
                                              ccTATTTTGCCGACCTGAAATCCGCTACTGACCGGCTCCAGCCGCTGACCTTAGATACCTGCCAGCAAGc                    >  1:974548/1‑70 (MQ=255)
                                              ccTATTTTGCCGACCTGAAATCCGCTACTGACCGGCTCCAGCCGCTGACCTTAGATACCTGCCAGCAAGCt                   <  1:2565105/71‑1 (MQ=255)
                                                  ttttGCCGACCTGAAATCCACTACTGACCGGCTCCAGCCGCTGACCTTa                                     <  1:3088128/49‑1 (MQ=255)
                                                  ttttGCCGACCTGAAATCCACTACTGACCGGCTCCAGCCGCTGACCTTa                                     <  1:3097785/49‑1 (MQ=255)
                                                  ttttGCCGACCTGAAATCCACTACTGACCGGCTCCAGCCGCTGACCTTa                                     <  1:59909/49‑1 (MQ=255)
                                                  ttttGCCGACCTGAAATCCACTACTGACCGGCTCCAGCCGCTGACCTTa                                     <  1:2684046/49‑1 (MQ=255)
                                                  ttttGCCGACCTGAAATCCACTACTGACCGGCTCCAGCCGCTGACCTTa                                     <  1:2202900/49‑1 (MQ=255)
                                                  ttttGCCGACCTGAAATCCACTACTGACCGGCTCCAGCCGCTGACCTTa                                     <  1:1639212/49‑1 (MQ=255)
                                                  ttttGCCGACCTGAAATCCACTACTGACCGGCTCCAGCCGCTGACCTTAGa                                   >  1:3911696/1‑51 (MQ=255)
                                                  ttttGCCGACCTGAAATCCACTACTGACCGGCTCCAGCCGCTGACCTTAGa                                   >  1:3916295/1‑51 (MQ=255)
                                                  ttttGCCGACCTGAAATCCACTACTGACCGGCTCCAGCCGCTGACCTTAGATACCTGCCAGCAAGCTAAcc               >  1:735475/1‑71 (MQ=255)
                                                    ttGCCGACCTGAAATCCGCTACTGACCGGCTCCAGCCGCTGACCTTAGATACCTGCCAGCAAGCTAAcccc             >  1:3553450/1‑71 (MQ=255)
                                                    ttGCCGACCTGAAATCCGCTACTGACCGGCTCCAGCCGCTGACCTTAGATACCTGCCAGCAAGCTAAcccc             >  1:515417/1‑71 (MQ=255)
                                                    ttGCCGACCTGAAATCCACTACTGACCGGCTCCAGCCGCTGACCTTAGATACCTGCCAGCAAGCTAAcccc             >  1:2196991/1‑71 (MQ=255)
                                                    ttGCCGACCTGAAATCCACTACTGACCGGCTCCAGCCGCTGACCTTAGATACCTGCCAGCAAGCTAAccc              >  1:672348/1‑70 (MQ=255)
                                                               aaaTCCACTACTGACCGGCTCCAGCCGCTGACCTTAGATACCTGCCAGCAAGCTAACCCCGAACTGAccgc  >  1:3826733/1‑71 (MQ=255)
                                                               aaaTCCACTACTGACCGGCTCCAGCCGCTGACCTTAGATACCTGCCAGCAAGCTAACCCCGAACTGAccgc  >  1:3263463/1‑71 (MQ=255)
                                                                  tCCACTACTGACCGGCTCCAGCCGCTGACCTTAGATACCTGCCAGCAAGCt                   <  1:977187/51‑1 (MQ=255)
                                                                  tCCACTACTGACCGGCTCCAGCCGCTGACCTTAGATACCTGCCAGCAAGCt                   <  1:3639733/51‑1 (MQ=255)
                                                                  tCCACTACTGACCGGCTCCAGCCGCTGACCTTAGATACCTGCCAGCAAGCt                   <  1:2093553/51‑1 (MQ=255)
                                                                     |                                                                
GTCAAATACGACACACTGATTACCGACGTACAAAAGTATCTCGATACCTATTTTGCCGACCTGAAATCCACTACTGACCGGCTCCAGCCGCTGACCTTAGATACCTGCCAGCAAGCTAACCCCGAACTGACCGC  >  W3110S.gb/2274180‑2274313

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: