Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,916,897 A→G 10.6% N293S (AAC→AGC)  glyS → glycine tRNA synthetase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,916,8970AG10.6% 119.1 / 6.9 47N293S (AAC→AGC) glySglycine tRNA synthetase, beta subunit
Reads supporting (aligned to +/- strand):  ref base A (22/20);  new base G (3/2);  total (25/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.22e-01

GGCTGAAGCGCTGGTTTACACCATGAAAGGTGACCAGAAATACTTCCCGGTGTATGCGAACGACGGCAAACTGCTGCCGAACTTTATCTTCGTTGCCAACATCGAATCGAAAGATCCGCAGCAGATTATC  >  W3110S.gb/3916838‑3916967
                                                           |                                                                      
ggCTGAAGCGCTGGTTTACACCATGAAAGGTGACCAGAAATACTTCCCGGTGTATGCGAACGACGGCAAAc                                                             <  1:1427340/71‑1 (MQ=255)
 gCTGAAGCGCTGGTTTACACCATGAAAGGTGACCAGAAATACTTCCCGGTGTATGCGAACGACGGCAAAct                                                            >  1:4589604/1‑71 (MQ=255)
        cgcTGGTTTACACCATGAAAGGTGACCAGAAATACTTCCCGGTGTATGCGAACGACGGCAAACTGCTGCCg                                                     <  1:1224836/71‑1 (MQ=255)
             gTTTACACCATGAAAGGTGACCAGAAATACTTCCCGGTGTATGCGAACGACgg                                                                  <  1:3039313/53‑1 (MQ=255)
             gTTTACACCATGAAAGGTGACCAGAAATACTTCCCGGTGTATGCGAACGACGGCa                                                                <  1:2105702/55‑1 (MQ=255)
                  cacCATGAAAGGTGACCAGAAATACTTCCCGGTGTATGCGAac                                                                       <  1:2378437/43‑1 (MQ=255)
                  cacCATGAAAGGTGACCAGAAATACTTCCCGGTGTATGCGAACGACGGCaaa                                                              >  1:4093399/1‑52 (MQ=255)
                     cATGAAAGGTGACCAGAAATACTTCCCGGTGTATGCGAGCGACGGCAAACTGCTGCCGAACTTTATCTTCg                                        >  1:355050/1‑71 (MQ=255)
                      aTGAAAGGTGACCAGAAATACTTCCCGGTGTATGCGAACGACGGCAAACTGCTGCCGAACTTTATCTTCGt                                       <  1:2334317/71‑1 (MQ=255)
                        gAAAGGTGACCAGAAATACTTCCCGGTGTATGCGAACGACGGCAAACTGCTGCCGAACt                                                 <  1:764811/59‑1 (MQ=255)
                         aaaGGTGACCAGAAATACTTCCCGGTGTATGCGAACGACGGCAAACTGCTGCCGAACTTTATCTTCGTTGc                                    >  1:3385880/1‑71 (MQ=255)
                          aaGGTGACCAGAAATACTTCCCGGTGTATGCGAACGACGGCAAACTGCTGCCGAACTTTATCTTCGTTGcc                                   <  1:894419/71‑1 (MQ=255)
                           aGGTGACCAGAAATACTTCCCGGTGTATGCGAGCGACGGCAAACTGCTGCCGaa                                                   >  1:2862486/1‑54 (MQ=255)
                            ggTGACCAGAAATACTTCCCGGTGTATGCGAACGACGGCAAACTGCTGCCGAACTTTATCTTCGTTGCCaa                                 <  1:4498732/71‑1 (MQ=255)
                              tGACCAGAAATACTTCCCGGTGTATGCGAACGACGGCAAACTGCTGCCGAACTTTAt                                             >  1:294560/1‑57 (MQ=255)
                                 ccAGAAATACTTCCCGGTGTATGCGAACGACGGCAAACTGCTGCCGAACTTTATCTTCGTTGCCAACatcg                            <  1:2783459/71‑1 (MQ=255)
                                  cAGAAATACTTCCCGGTGTATGCGAGCGACGGCAAACTGCTGCCGAACTTTATCTTCGATGCCAACatcg                            >  1:3179096/1‑70 (MQ=255)
                                   aGAAATACTTCCCGGTGTATGCGAACGACGGCAAACTGCTGCCGAACTTTATCTTCGTTGCCAACatcgaa                          <  1:886757/71‑1 (MQ=255)
                                    gAAATACTTCCCGGTGTATGCGAACGACGGCAAACTGCTGCCGAACttt                                               <  1:885515/49‑1 (MQ=255)
                                     aaaTACTTCCCGGTGTATGCGAGCGACGGCAAACTGCTGCCGAACTTTATCTTCGTTGCCAACatcgaatc                        <  1:1873795/71‑1 (MQ=255)
                                       aTACTTCCCGGTGTATGCGAACGACGGCAAACTGCTGCCGAACTTTATc                                            >  1:1881349/1‑49 (MQ=255)
                                       aTACTTCCCGGTGTATGCGAACGACGGCAAACTGCTGCCGAACTTTATCTTCGTTg                                     >  1:548243/1‑56 (MQ=255)
                                       aTACTTCCCGGTGTATGCGAACGACGGCAAACTGCTGCCGAACTTTATCTTCGTTg                                     >  1:3621000/1‑56 (MQ=255)
                                       aTACTTCCCGGTGTATGCGAACGACGGCAAACTGCTGCCGAACTTTATCTTCGTTGCCAACatcgaatcg                       >  1:3179464/1‑70 (MQ=255)
                                       aTACTTCCCGGTGTATGCGAACGACGGCAAACTGCTGCCGAACTTTATCTTCGTTGCCAACatcaaatcg                       >  1:374760/1‑70 (MQ=255)
                                         aCTTCCCGGTGTATGCGAACGACGGCAAACTGCTGCCGAACTTTAt                                             >  1:1870719/1‑46 (MQ=255)
                                         aCTTCCCGGTGTATGCGAACGACGGCAAACTGCTGCCGAACTTTATCTTCGTTGCCaa                                 >  1:1131393/1‑58 (MQ=255)
                                          cTTCCCGGTGTATGCGAACGACGGCAAACTGCTGCCGAACttt                                               >  1:485402/1‑43 (MQ=255)
                                          cTTCCCGGTGTATGCGAACGACGGCAAACTGCTGCCGAACTTTATCTTCGTTg                                     >  1:638392/1‑53 (MQ=255)
                                          cTTCCCGGTGTATGCGAACGACGGCAAACTGCTGCCGAACTTTATCTTCGTTGCCAACatc                             >  1:4598429/1‑61 (MQ=255)
                                          cTTCCCGGTGTATGCGAACGACGGCAAACTGCTGCCGAACTTTATCTTCGTTGCCAACATCGAATCGAAAg                   <  1:861804/71‑1 (MQ=255)
                                              ccGGTGTATGCGAACGACGGCAAACTGCTGCCGAACTTTATCTTCGTTGCCAACATCGAATCGAAAGATcc               <  1:1333003/71‑1 (MQ=255)
                                               cGGTGTATGCGAACGACGGCAAACTGCTGCCGAACTTTATCTTCGTTGcc                                   <  1:4154400/50‑1 (MQ=255)
                                               cGGTGTATGCGAACGACGGCAAACTGCTGCCGAACTTTATCTTCGTTGCCAACATCGAATCGAAAGATCCg              <  1:3401079/71‑1 (MQ=255)
                                                  tgtATGCGAGCGACGGCAAACTGCTGCCGAACTTTATc                                            <  1:3037097/38‑1 (MQ=255)
                                                   gtATGCGAACGACGGCAAACTGCTGCCGAACTTTATCTTCGTTCCAACATCGAATCGaaa                    >  1:4108387/1‑60 (MQ=255)
                                                    tATGCGAACGACGGCAAACTGCTGCCGAACTTTATCTTCGTTGCCAACatcgaatcga                      <  1:1255470/58‑1 (MQ=255)
                                                      tGCGAACGACGGCAAACTGCTGCCGAACTTTATCTTCGTTGCCAACatcgaa                          <  1:343488/52‑1 (MQ=255)
                                                      tGCGAACGACGGCAAACTGCTGCCGAACTTTATCTTCGTTGCCAACATCGAATCGaaa                    >  1:299513/1‑58 (MQ=255)
                                                      tGCGAACGACGGCAAACTGCTGCCGAACTTTATCTTCGTTGCCAACATCGAATCGAAAGATCCgcag           >  1:957418/1‑67 (MQ=255)
                                                      tGCGAACGACGGCAAACTGCTGCCGAACTTTATCTTCGTTGCCAACATCGAATCGAAAGATCCgcag           >  1:2670966/1‑67 (MQ=255)
                                                      tGCGAACGACGGCAAACTGCTGCCGAACTTTATCTTCGTTGCCAACATCGAATCGAAAGATCCgca            >  1:4089226/1‑66 (MQ=255)
                                                       gCGAACGACGGCAAACTGCTGCCGAACTTTATCTTCGTTGCCAACatcgaa                          >  1:2205509/1‑51 (MQ=255)
                                                        cGAACGACGGCAAACTGCTGCCGATCTTTATCTTCGTTGCCAACATCGAATCGAAAGATCCGCAGCAGAtt     <  1:2369436/71‑1 (MQ=255)
                                                        cGAACGACGGCAAACTGCTGCCGAACTTTATCTTCGTTGCCAACATCGAATCGAAAGATcc               >  1:1334039/1‑61 (MQ=255)
                                                        cGAACGACGGCAAACTGCTGCCGAACTTTATCTTCGTTGCCAACATCGAATCGAAAGATCCGCAGCAGAtt     <  1:3522389/71‑1 (MQ=255)
                                                        cGAACGACGGCAAACTGCTGCCGAACTTTATCTTCGTTGCCAACATCGAATCGAAAGATCCGCAGCAGAtt     <  1:2137570/71‑1 (MQ=255)
                                                        cGAACGACGGCAAACTGCTGCCGAACTTTATCTTCGTTGCCAACATCGAATCGAAAGATCCGCAGCAGAt      >  1:404096/1‑70 (MQ=255)
                                                          aaCGACGGCAAACTGCTGCCGAACTTTATCTTCGTTGCCAACATCGAATCGAAAGATCCGCAGCAGATTAt   <  1:853113/71‑1 (MQ=255)
                                                           acgacgGCAAACTGCTGCCGAACTTTATCTTCGTTGCCAACATCGAATCGAAAGATCCGCAGCAGATTAtc  >  1:3225766/1‑71 (MQ=255)
                                                           |                                                                      
GGCTGAAGCGCTGGTTTACACCATGAAAGGTGACCAGAAATACTTCCCGGTGTATGCGAACGACGGCAAACTGCTGCCGAACTTTATCTTCGTTGCCAACATCGAATCGAAAGATCCGCAGCAGATTATC  >  W3110S.gb/3916838‑3916967

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: