Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 707,271 C→T 16.1% H253Y (CAC→TAC)  glnS → glutamyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb707,2710CT16.1% 75.6 / 7.5 31H253Y (CAC→TAC) glnSglutamyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base C (15/11);  new base T (2/3);  total (17/14)
Fisher's exact test for biased strand distribution p-value = 6.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.20e-01

CGCTTGAGTTCCAGGACAACCGTCGTCTGTACGACTGGGTACTGGACAACATCACGATTCCTGTTCACCCGCGCCAGTATGAGTTCTCGCGCCTGAATCTGGAATACACCGTGATGTCCAAGCGTAAGTTGA  >  W3110S.gb/707206‑707337
                                                                 |                                                                  
cGCTTGAGTTCCAGGACAACCGTCGTCTGTACGACTGGGTACTGGACAACATCACGATTCCTGTTTACCcg                                                               >  1:2582194/1‑71 (MQ=255)
    tGAGTTCCAGGACAACCGTCGTCTGTACGACTGGGTACTGGACAACATCACGATTCCTGTTCACCCgcgc                                                            >  1:2460343/1‑70 (MQ=255)
     gAGTTCCAGGACAACCGTCGTCTGTACGACTGGGTACTGGACAACATCACGATTCCTGTTCACCCgcgc                                                            >  1:2015777/1‑69 (MQ=255)
          ccAGGACAACCGTCGTCTGTACGACTGGGTACTGGACAACATCACGATTCCTGTTCACCCGCGCCAGTAt                                                      <  1:2247576/70‑1 (MQ=255)
                cAACCGTCGTCTGTACGACTGGGTACTGGACAACATCACGATTCCTGTTCACCCGCGCCAGTATGAGTtct                                               <  1:3446587/71‑1 (MQ=255)
                         tcTGTACGACTGGGTACTGGACAACATCACGATTCCTGTTCACCCGCGCCAg                                                         <  1:1744271/52‑1 (MQ=255)
                           tGTACGACTGGGTACTGGACAACATCACGATTCCTGTTTACCCGCGCCAGTATGAGTTCTCGCGCCTGAAt                                    >  1:2745979/1‑71 (MQ=255)
                             tACGACTGGGTACTGGACAACATCACGATTCCTGTTCACCCCGCCAGTATGAGTTCTCGCGCCTGAATCTg                                 <  1:1568536/71‑1 (MQ=255)
                               cGACTGGGTACTGGACAACATCACGATTCCTGTTTACCCGCGCCAGTATGAGTTCTCGCGCCTGAATCTgg                                <  1:2795622/71‑1 (MQ=255)
                               cGACTGGGTACTGGACAACATCACGATTCCTGTTCACCCGCGCCAGTCTGAGTTCTCGCGCCTGAATCTgg                                <  1:1776912/71‑1 (MQ=255)
                               cGACTGGGTACTGGACAACATCACGATTCCTGTTCACCCGCGCCAGTATGAGTTCTCGCGCCTGAATCTg                                 >  1:1501932/1‑70 (MQ=255)
                               cGACTGGGTACTGGACAACATCACGATTCCTGTTCACCCGCGCCAGTATGAGTTCTCGCGCCTGAATCTg                                 <  1:3961870/70‑1 (MQ=255)
                                         cTGGACAACATCACGATTCCTGTTCACCCGCGCCAGTATGAGTTCTCGCGCCTGAATCt                                  <  1:926077/59‑1 (MQ=255)
                                            gACAACATCACGATTCCTGTTTACCCGCGCCAGTATGAGTTCTCGCGCCTGAATCTGGAATACACCGTg                     <  1:2195495/69‑1 (MQ=255)
                                               aacaTCACGATTCCTGTTCACCCGCGCCAGTATGAGTTCTCGCGCCTGAATCTGGAATACACCGTGATGTc                >  1:257712/1‑71 (MQ=255)
                                                   tCACGATTCCTGTTCACCCGCGCCAGTATGAGTtct                                               <  1:982708/36‑1 (MQ=255)
                                                    cACGATTCCTGTTCACCCGCGCCAGTATGAGTTCTCGCGCCTGAATCTGGaa                              >  1:1525600/1‑52 (MQ=255)
                                                    cACGATTCCTGTTCACCCGCGCCAGTATGAGTTCTCGCGCCTGAATCTGGaa                              >  1:2138853/1‑52 (MQ=255)
                                                    cACGATTCCTGTTCACCCGCGCCAGTATGAGTTCTCGCGCCTGAATCTGGaa                              >  1:1807399/1‑52 (MQ=255)
                                                     aCGATTCCTGTTCACCCGCGCCAGTATGAGTTCTcg                                             >  1:2475251/1‑36 (MQ=255)
                                                     aCGATTCCTGTTCACCCGCGCCAGTATGAGTTCTcg                                             >  1:2783769/1‑36 (MQ=255)
                                                     aCGATTCCTGTTCACCCGCGCCAGTATGAGTTCTcg                                             >  1:3953978/1‑36 (MQ=255)
                                                     aCGATTCCTGTTCACCCGCGCCAGTATGAGTTCTCGCGCCTGAATCTg                                 >  1:357599/1‑48 (MQ=255)
                                                     aCGATTCCTGTTCACCCGCGCCAGTATGAGTTCTCGCGCCTGAATCTGGAATACACCGTGATGTCCAAGCg          >  1:2768422/1‑71 (MQ=255)
                                                     aCGATTCCTGTTCACCCGCGCCAGTATGAGTTCTCGCGCCTGAATCTGGAATACACCGTGATGTCCAAGCg          >  1:2174777/1‑71 (MQ=255)
                                                     aCGATTCCTGTTCACCCGCGCCAGTATGAGTTCTCGCGCCTGAATCTGGAATACACCGTGATGTCCAAGCg          >  1:1895360/1‑71 (MQ=255)
                                                     aCGATTCCTGTTCACCCGCGCCAGTATGAGTTCTCGCGCCTGAATCTGGAATACACCGTGATGTCCAAGCg          >  1:928023/1‑71 (MQ=255)
                                                            cTGTTTACCCGCGCCAGTATGAGTTCTCGCGCCTGAATCTGGAATACACCGTGATGTCCAAGCGTAAGTTg   <  1:2441523/71‑1 (MQ=255)
                                                            cTGTTCACCCGCGCCAGTATGAGTTCTCGCGCCTGAATCTGGAATACACCGTGATGTCCAAGCGTAAGTTg   <  1:3802789/71‑1 (MQ=255)
                                                               ttCACCCGCGCCAGTATGAGTTCTCGCGCCTGAATCTGGAATACACCGTGATGTCCAAGCGTAAGTTGa  <  1:2314447/69‑1 (MQ=255)
                                                               ttCACCCGCGCCAGTATGAGTTCTCGCGCCTGAATCTGGAATACACCGTGATGTCCAAGCGTAAGTTGa  <  1:2810635/69‑1 (MQ=255)
                                                                 |                                                                  
CGCTTGAGTTCCAGGACAACCGTCGTCTGTACGACTGGGTACTGGACAACATCACGATTCCTGTTCACCCGCGCCAGTATGAGTTCTCGCGCCTGAATCTGGAATACACCGTGATGTCCAAGCGTAAGTTGA  >  W3110S.gb/707206‑707337

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