Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 707,271 C→T 14.2% 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,2710CT14.2% 92.6 / 6.2 35H253Y (CAC→TAC) glnSglutamyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base C (18/12);  new base T (2/3);  total (20/15)
Fisher's exact test for biased strand distribution p-value = 6.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

CGCTTGAGTTCCAGGACAACCGTCGTCTGTACGACTGGGTACTGGACAACATCACGATTCCTGTTCACCCGCGCCAGTATGAGTTCTCGCGCCTGAATCTGGAATACACCGTGATGTCCAAGCGTAAGTTGA  >  W3110S.gb/707206‑707337
                                                                 |                                                                  
cGCTTGAGTTCCAGGACAACCGTCGTCTGTACGACTGGGTACTGGACAACATCACGATTCCTGTTCACCcg                                                               >  1:2418529/1‑71 (MQ=255)
          ccAGGACAACCGTCGTCTGTACGACTGGGTACTGGACAACATCACGATTCCTGTTCACCCGCGCCAGTAt                                                      >  1:1428400/1‑70 (MQ=255)
           cAGGACAACCGTCGTCTGTACGACTGGGTACTGGACAACATCACGATTCCTGTTTACCCGCGCCAGTATg                                                     <  1:488424/70‑1 (MQ=255)
           cAGGACAACCGTCGTCTGTACGACTGGGTACTGGACAACATCACGATTCCTGTTCACCCGCGCCAGTATg                                                     <  1:1605744/70‑1 (MQ=255)
           cAGGACAACCGTCGTCTGTACGACTGGGTACTGGACAACATCACGATTCCTGTTCACCCGCGCCAGTATg                                                     <  1:2207233/70‑1 (MQ=255)
             ggACAACCGTCGTCTGTACGACTGGGTACTGGACAACATCACGATTCCTGTTCACCCGCGCCAGTa                                                       >  1:2197079/1‑66 (MQ=255)
                cAACGGTCGTCTGTACGACTGGGTACTGGACAACATCACGATTCCTGTTCACCCGCGCCAGTATGAGTtct                                               <  1:1757718/71‑1 (MQ=255)
                    cgtcgtCTGTACGACTGGGTACTGGACAACATCACGATTCCTGTTCACCCGCGCCAGTATGAGTTCTcgcg                                           <  1:1036217/71‑1 (MQ=255)
                    cgtcgtCTGTACGACTGGGTACTGGACAACATCACGATTCCTGTTCACCCGCGCCAGTATGAGTTCTcgcg                                           <  1:1039423/71‑1 (MQ=255)
                      tcgtcTGTACGACTGGGTACTGGACAACATCACGATTCCTGTTTACCCGCGCCAGTATGAg                                                   <  1:2460753/61‑1 (MQ=255)
                        gtcTGTACGACTGGGTACTGGACAACATCACGATTCCTGTTTACCCGCGCCAGTATGAGtt                                                 >  1:492308/1‑61 (MQ=255)
                        gtcTGTACGACTGGGTACTGGACAACATCACGATTCCTGTTCAc                                                                  >  1:257214/1‑44 (MQ=255)
                          cTGTACGACTGGGTACTGGACAACATCACGATTCCTGTTCACCCGCGCCAGTATGAGTTCTCGCGCCTGaa                                     >  1:1430545/1‑71 (MQ=255)
                           tGTACGACTGGGTACTGGACAACATCACGATTCCTGTTCACCCGCGCCAGTATGAGTTCTCGCGCCTGaa                                     >  1:2048521/1‑70 (MQ=255)
                             tACGACTGGGTACTGGACAACATCACGATTCCTGTTCACCCGCGCCAGTATGAGTTCTCGCGCCTGAATc                                   <  1:27586/70‑1 (MQ=255)
                             tACGACTGGGTACTGGACAACATCACGATTCCTGTTCACCCGCGCCAGTATGAGTTCTCGCGCCTGAATc                                   <  1:1067853/70‑1 (MQ=255)
                              aCGACTGGGTACTGGACAACATCACGATTCCTGTTCACCCGCGCCAGTATGAGTTCTCGCGCCTGAATCTg                                 >  1:1038297/1‑71 (MQ=255)
                               cGACTGGGTACTGGACAACATCACGATTCCTGTTCACCCGCGCCAGTa                                                       >  1:1287517/1‑48 (MQ=255)
                               cGACTGGGTACTGGACAACATCACGATTCCTGTTCACCCGCGCCAGTa                                                       >  1:2655486/1‑48 (MQ=255)
                                     ggTACTGGACAACATCACGATTCCTGTTCACCCGCGCCAGTATGAGTTCTCGCGCCTGAATCTGGAATaca                          <  1:2525603/71‑1 (MQ=255)
                                          tGGACAACATCACGATTCCTGTTTACCCGCGCCAGTATGAGTTCTCGCGCCTGAATCTGGaa                              <  1:1060413/62‑1 (MQ=255)
                                            gACAACATCACGATTCCTGTTCACCCGCGCCAGTATGAGTTCTCGCGCCTGAATCTGGAATACACCGTg                     <  1:966285/69‑1 (MQ=255)
                                            gACAACATCACGATTCCTGTTCACCCGCGCCAGTATGAGTTCTCGCGCCTGAATCTGGAATACACCGTGAt                   >  1:415367/1‑71 (MQ=255)
                                            gACAACATCACGATTCCTGTTCACCCGCGCCAGTATGAGTTCTCGCGCCTGAATCTGGAATACACCGTGAt                   >  1:412084/1‑71 (MQ=255)
                                            gACAACATCACGATTCCTGTTCACCCGCGCCAGTATGAGTTCTCGCGCCTGAATCTGGAATACACCGTGAt                   >  1:408953/1‑71 (MQ=255)
                                            gACAAAATCACGATTCCTGTTCACCCGCGCCAGTATGAGTTCTCGCGCCTGAATCTGGAATACACCGTGa                    >  1:2045043/1‑70 (MQ=255)
                                               aacaTCACGATTCCTGTTCACCCGCGCCAGTATGAGTTCTCGCGCCTGaa                                     >  1:712843/1‑50 (MQ=255)
                                               aacaTCACGATTCCTGTTCACCCGCGCCAGTATGAGTTCTCGCGCCTGaa                                     >  1:669844/1‑50 (MQ=255)
                                                  aTCACGATTCCTGTTTACCCGCGCCAGTATGATTTCTCGCGCCTGAATCTGGaa                              >  1:372658/1‑54 (MQ=255)
                                                  aTCACGATTCCTGTTCACCCGCGCCAGTATGAGTTCTCGCGCCTGAATCTGGAATACACCGTGATGTcc               >  1:1911837/1‑69 (MQ=255)
                                                   tCACGATTCCTGTTCACCCGCGCCAGTATGAGTTCTCGCGCCTGAATCTGGaa                              <  1:1383212/53‑1 (MQ=255)
                                                    cACGATTCCTGTTCACCCGCGCCAGTATGAGTTCTCGCGCCTGAATCTGGAATaca                          >  1:40102/1‑56 (MQ=255)
                                                     aCGATTCCTGTTCACCCGCGCCAGTATGAGTTCTCGCGCCTGAATCTGGAATACACCGTGATGTCCAAGCg          >  1:1483136/1‑71 (MQ=255)
                                                            cTGTTCACCCGCGCCAGTATGAGTTCTCGCGCCTGAATCTGGAATACACCGTGATGTCCAAGCGTAAGTTg   <  1:183839/71‑1 (MQ=255)
                                                               ttCACCCGCGCCAGTATGAGTTCTCGCGCCTGAATCTGGAATACACCGTGATGTCCAAGCGTAAGTTGa  <  1:1107121/69‑1 (MQ=255)
                                                                 cACCCGCGCCAGTATGAGTTCTCGCGCCTGAATCTGGAATACACCGTGa                    >  1:1241132/1‑49 (MQ=255)
                                                                 cACCCGCGCCAGTATGAGTTCTCGCGCCTGAATCTGGAATACACCGTGa                    >  1:1507453/1‑49 (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: