Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,487,895 T→C 21.4% K209E (AAA→GAA)  valS ← valyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,487,8950TC21.4% 63.7 / 14.5 34K209E (AAA→GAA) valSvalyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base T (14/12);  new base C (2/5);  total (17/17)
Fisher's exact test for biased strand distribution p-value = 3.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCGCCCAGCAGGGTTTCTGGACGGGTAGTCGCGACCACCAGATAATCTTTACCGTCTGCGGTTTTCGCACCGTCAGCCAGCGGATAGCGGATGTGCCACATCGAACCTTTCGATTCGCGGTTTTCGACTTCCAG  >  W3110S.gb/4487830‑4487964
                                                                 |                                                                     
aTCGCCCAGCAGGGTTTCTTTACGGGTAGTCGCGACCACCAGATAATCTTTACCGTCTGCGGTTTTCg                                                                     <  1:3995806/68‑1 (MQ=255)
 tCGCCCAGCAGGGTTTCTGGACGGGTAGTCGCGACCACCAGATAATCTTTACCGTCTGCGGTTTTCGCAcc                                                                 <  1:2768884/71‑1 (MQ=255)
 tCGCCCAGCAGGGTTTCTGGACGGGTAGTCGCGACCACCAGATAATCTTTACCGTCTGCGGTTTTCGCAcc                                                                 >  1:3233939/1‑71 (MQ=255)
 tCGCCCAGCAGGGTTTCTGGACGGGTAGTCGCGACCACCAGATAATCTTTACCGTCTGCGGTTTTCGCAcc                                                                 >  1:3798846/1‑71 (MQ=255)
 tCGCCCAGCAGGGTTTCTGGACGGGTAGTCGCGACCACCAGATAATCTTTACCGTCTGCGGTTTTCGCAcc                                                                 >  1:3805335/1‑71 (MQ=255)
 tCGCCCAGCAGGGTTTCTGGACGGGTAGTCGCGACCACAAGATAATCTTTACCGTCTGCGGTTTACGCa                                                                   >  1:2575791/1‑69 (MQ=255)
   gCCCAGCAGGGTTTCTGGACGGGTAGTCGCGACCACCAGATAATCTTTACCGTCTGCGGTTTTCGCACCGt                                                               >  1:297980/1‑71 (MQ=255)
          agGGTTTCTGGACGGGTAGTCGCGACCACCAGATAATCTTTACCGTCTGCGGTTTTCGCACCGTcagccag                                                        <  1:547382/71‑1 (MQ=255)
          agGGTTTCTGGACGGGTAGTCGCGACCACCAGATAATCTTTACCGTCTGCGGTTTTCGCACCGTcagccag                                                        <  1:2647450/71‑1 (MQ=255)
            ggTTTCTGGACGGGTAGTCGCGACCACCAGATAATCTTTACCGTCTGCGGTTTCCGCACCGTCAGCCagcg                                                      <  1:2273470/71‑1 (MQ=255)
                     aCGGGTAGTCGCGACCACCAGATAATCTTTACCGTCTGCGGTTTTCGCACCGTCAGCCagcggatagcgga                                             <  1:2393032/71‑1 (MQ=255)
                     aCGGGTAGTCGCGACCACCAGATAATCTTTACCGTCTGCGGTTTCCGCACCGTCAGCCagcggatagcgga                                             <  1:670717/71‑1 (MQ=255)
                               gcgACCACCAGATAATCTTTACCGTCTGCGGTTTTCGCACCGTCAGCCAGCGGATAGCGGATGTGCCACAt                                   >  1:3412710/1‑71 (MQ=255)
                               gcgACCACCAGATAATCTTTACCGTCTGCGGTTTTCGCACCGTCAGCCAGCGGATAGCGGATGTGCCACAt                                   >  1:3892977/1‑71 (MQ=255)
                               gcgACCACCAGATAATCTTTACCGTCTGCGGTTTTCGCACCGTCAGCCAGCGGATAGCGGATGTGCCACAt                                   >  1:1837181/1‑71 (MQ=255)
                                  accaccAGATAATCTTTACCGTCTGCGGTTTTCGCACCGTCAGCCAGCGGATAGCGGATGTGCCACATCGa                                <  1:2554648/71‑1 (MQ=255)
                                  accaccAGATAATCTTTACCGTCTGCGGTTTTCGCACCGTCAGCCAGCGGATAGCGGATGTGCCACATCGa                                <  1:2096142/71‑1 (MQ=255)
                                  accaccAGATAATCTTTACCGTCTGCGGTTTTCGCACCGTCAGCCAGCGGATAGCGGATGTGCCACATCGa                                <  1:1321940/71‑1 (MQ=255)
                                        aGATAATCTTTACCGTCTGCGGTTTTCGCACCGTCAGCCagcg                                                      <  1:335979/43‑1 (MQ=255)
                                         gATAATCTTTACCGTCTGCGGTTTTCGCACCGTCAGCCAGCGGATAGCGGATGTGCCACATCGAACCTTTc                         >  1:3856069/1‑71 (MQ=255)
                                           tAATCTTTACCGTCTGCGGTTTTCGCACCGTCAGCCAGCGGATAGCGGATGTGCCACATCGAACCTTTCGa                       >  1:1439521/1‑71 (MQ=255)
                                            aaTCTTTACCGTCTGCGGTTTTCGCACCGTCAGCCAGCGGATAGCGGATGTGCCACATCGAACCTTTCGAt                      >  1:2958/1‑71 (MQ=255)
                                            aaTCTTTACCGTCTGCGGTTTTCGCACCGTCAGCCAGCGGATAGCGGATGTGCCACATCGAACCTTTCGAt                      >  1:1638602/1‑71 (MQ=255)
                                            aaTCTTTACCGTCTGCGGTTTTCGCACCGTCAGCCAGCGGATAGCGGATGTGCCACATCGAACCTTTCGAt                      >  1:3144657/1‑71 (MQ=255)
                                             aTCTTTACCGTCTGCGGTTTCCGCACCGTCAGCCagcg                                                      <  1:2309330/38‑1 (MQ=255)
                                             aTCTTTACCGTCTGCGGTTTCCGCACCGTCAGCCagcg                                                      <  1:219171/38‑1 (MQ=255)
                                             aTCTTTACCGTCTGCGGTTTCCGCACCGTCAGCCagcg                                                      <  1:1801744/38‑1 (MQ=255)
                                                  tACCGTCTGCGGTTTCCGCACCGTCAGCCAGCGGATAGCGGATGTGCCACATCGaa                               >  1:3730792/1‑56 (MQ=255)
                                                  tACCGTCTGCGGTTTCCGCACCGTCAGCCAGCGGATAGCGGATGTGCCACATCGAACCTTTCGATTCGCgg                >  1:1979343/1‑71 (MQ=255)
                                                       tCTGCGGTTTTCGCACCGTCAGCCAGCGGATAGCGGATGTGCCACATCGAACCTTTCGATTCGCgg                <  1:3063094/66‑1 (MQ=255)
                                                       tCTGCGGTTTTCGCACCGTCAGCCAGCGGATAGCGGATGTGCCACATCGAACCTTTCGATTCGCGGTTTTc           <  1:2626641/71‑1 (MQ=255)
                                                       tCTGCGGTTTTCGCACCGTCAGCCAGCGGATAGCGGATGTGCCACATCGAACCTTTCGATTCGCGGTTTTc           <  1:2505510/71‑1 (MQ=255)
                                                         tGCGGTTTTCGCACCGTCAGCCAGCGGATAGCGGATGTGCCACATCGaa                               >  1:2561531/1‑49 (MQ=255)
                                                         tGCGGTTTTCGCACCGTCAGCCAGCGGATAGCGGATGTGCCACATCGaa                               >  1:2812598/1‑49 (MQ=255)
                                                                ttCGCACCGTCAGCCAGCGGATAGCGGATGTGCCACATCGAACCTTTCGATTCGCGGTTTTCGACTTCCAg  >  1:2693874/1‑71 (MQ=255)
                                                                ttCGCACCGTCAGCCAGCGGATAGCGGATGTGCCACATCGAACCTTTCGATTCGCGGTTTTCGACTTCCAg  >  1:2688952/1‑71 (MQ=255)
                                                                ttCGCACCGTCAGCCAGCGGATAGCGGATGTGCCACATCGAACCTTTCGATTCGCGGTTTTCGACTTCCAg  >  1:3649357/1‑71 (MQ=255)
                                                                ttCGCACCGTCAGCCAGCGGATAGCGGATGTGCCACATCGAACCTTTCGATTCGCGGTTTTCGACTTCCAg  >  1:2536670/1‑71 (MQ=255)
                                                                ttCGCACCGTCAGCCAGCGGATAGCGGATGTGCCACATCGAACCTTTCGATTCGCGGTTTTCGACTTCCAg  >  1:561183/1‑71 (MQ=255)
                                                                 |                                                                     
ATCGCCCAGCAGGGTTTCTGGACGGGTAGTCGCGACCACCAGATAATCTTTACCGTCTGCGGTTTTCGCACCGTCAGCCAGCGGATAGCGGATGTGCCACATCGAACCTTTCGATTCGCGGTTTTCGACTTCCAG  >  W3110S.gb/4487830‑4487964

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: