Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,915,241 A→G 57.9% E48G (GAG→GGG)  glyQ → glycine tRNA synthetase, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,915,2410AG57.9% 6.3 / 19.3 19E48G (GAG→GGG) glyQglycine tRNA synthetase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (3/5);  new base G (8/3);  total (11/8)
Fisher's exact test for biased strand distribution p-value = 1.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-01

TTGTTCAACCATTGGACATGGAAGTCGGCGCGGGAACCTCTCACCCAATGACCTGTCTGCGCGAGCTGGGGCCAGAACCGATGGCGGCTGCTTATGTTCAGCCTTCTCGTCGCCCGACCGATGGTCGCTACG  >  W3110S.gb/3915178‑3915309
                                                               |                                                                    
ttGTTCAACCATTGGACATGGAAGTCGGCGCGGGAACCTCTCACCCAATGACCTGTCTGCGCGGGCTgggg                                                               >  1:1629382/1‑71 (MQ=255)
ttGTTCAACCATAGGACATGGAAGTCGGCGCGGGAACCTCTCACCCAATGACCTGTCTGCGCGAGCTg                                                                  <  1:1566620/68‑1 (MQ=255)
     cAACCATTGGACATGGAAGTCGGCGCGGGAACCTCTCACCCAATGACCTGTCTGCGCGGGCTGGGGCCAg                                                           >  1:540437/1‑70 (MQ=255)
     cAACCATTGGACATGGAAGTCGGCGCGGGAACCTCTCACCCAATGACCTGTCTGCGCGAGCTg                                                                  <  1:599050/63‑1 (MQ=255)
                cATGGAAGTCGGCGCGGGAACCTCTCACCCAATGACCTGTCTGCGCGGGCTGGGGCCAGAACCGATggcg                                                >  1:1475591/1‑70 (MQ=255)
                          ggCGCGGGAACCTCTCACCCAATGACCTGTCTGCGCGGGCTGGGGCCAGAACCGATGGCGGCTGCTTa                                        >  1:462976/1‑68 (MQ=255)
                             gcgGGAACCTCTCACCCAATGACCTGTCTGCGCGGGCTGGGGCCAGAACCGATGGCGGCTGCTTATGTTCa                                  <  1:723107/71‑1 (MQ=255)
                                ggAACCTCTCACCCAATGACCTGTCTGCGCGGGCTGGGGCCAGAACCGATGGCGGCTGCTTATGTTCAGcc                               <  1:323246/71‑1 (MQ=255)
                                     ctctCACCCAATGACCTGTCTGCGCGGGCTGGGGCCAGAAcc                                                       <  1:1134483/42‑1 (MQ=255)
                                      tctcACCCAATGACCTGTCTGCGCGAGCTGGGGCCAg                                                           >  1:1166689/1‑37 (MQ=255)
                                      tctcACCCAATGACCTGTCTGCGCGAGCTGGGGCCAGAACCGATGGCGGCTGCTTATGTTCAGCCt                              <  1:692478/66‑1 (MQ=255)
                                        tcACCCAATGACCTGTCTGCGCGAGCTGGGGCCAGAACCGATGGCGGCTGCTTATGTTCAGCCTTCtcg                         <  1:1385039/69‑1 (MQ=255)
                                           cccAAAGACATGTCTGCGCGGGCTGGGGCCAGAACCGATGGCGGCTGCTTATGTTCAGCCTTCTCGTCGcc                    >  1:1387539/1‑71 (MQ=255)
                                            ccAATGACCTGTCTGCGCGGGCTGGGGCCAGAACCGATGGCGGCTGCTTATGTTCAGCCTTCTCGTCGccc                   >  1:428524/1‑71 (MQ=255)
                                               aTGACCTGTCTGCGCGGGCTGGGGCCAGAACCGATgg                                                  >  1:1371584/1‑37 (MQ=255)
                                               aTGACCTGTCTGCGCGGGCTGGGGCCAGAACCGATgg                                                  >  1:994972/1‑37 (MQ=255)
                                                       tCTGCGCGATCTGGGGCCAGAACCGATGGCGGCTGCTTATGTTCAGCCTTCTCGTCGCCCGACCGATgg          <  1:907659/69‑1 (MQ=255)
                                                             cgAGCTGGGGCCAGAACCGATGGCGGCTGCTTATGTTCAGCCTTCTCGTCGCCCGACCGATGGTCGCTACg  >  1:143952/1‑71 (MQ=255)
                                                             cgAGCTGGGGCCAGAACCGATGGCGGCTGCTTATGTTCAGCCTTCTCGTCGCCAGACCGATGGTCGCTa    >  1:1196313/1‑69 (MQ=255)
                                                               |                                                                    
TTGTTCAACCATTGGACATGGAAGTCGGCGCGGGAACCTCTCACCCAATGACCTGTCTGCGCGAGCTGGGGCCAGAACCGATGGCGGCTGCTTATGTTCAGCCTTCTCGTCGCCCGACCGATGGTCGCTACG  >  W3110S.gb/3915178‑3915309

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