Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,915,241 A→G 30.2% 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,2410AG30.2% 53.6 / 33.3 43E48G (GAG→GGG) glyQglycine tRNA synthetase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (18/12);  new base G (5/8);  total (23/20)
Fisher's exact test for biased strand distribution p-value = 3.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

TTGTTCAACCATTGGACATGGAAGTCGGCGCGGGAACCTCTCACCCAATGACCTGTCTGCGCGAGCTGGGGCCAGAACCGATGGCGGCTGCTTATGTTCAGCCTTCTCGTCGCCCGACCGATGGTCGCTACG  >  W3110S.gb/3915178‑3915309
                                                               |                                                                    
ttGTTCAACCATTGGACATGGTAGTCGGCGCGGGAACCTCTCACCCAATGACCTGTCTGCGCGAGCTgggg                                                               >  1:3201253/1‑71 (MQ=255)
ttGTTCAACCATTGGACATGGAAGTCGGCGCGGGAACCTCTCACCCAATGACCTGTCTGCGCGGGCTgg                                                                 <  1:73552/69‑1 (MQ=255)
  gTTCAACCATTGGACATGGAAGTCGGCGCGGGAACCTCTCACCCAATGACCTGTCTGCGCGAGCTGGGGc                                                              >  1:3786094/1‑70 (MQ=255)
        ccATTGGACATGGAAGTCGGCGCGGGAACCTCTCACCCAATGACCTGTCTGCGCGGGCTGGGGCCAGAAcc                                                       <  1:3790545/71‑1 (MQ=255)
           ttGGACATGGAAGTCGGCGCGGGAACCTCTCACCCAATGACCTGTCTGCGCGGGCTGGGGCCAGAACCGAt                                                    >  1:4523004/1‑71 (MQ=255)
           ttGGACATGGAAGTCGGCGCGGGAACCTCTCACCCAATGACCTGTCTGCGCGAGCTGGGGCCAGAACCGAt                                                    <  1:3808148/71‑1 (MQ=255)
           ttGGACATGGAAGTCGGCGCGGGAACCTCTCACCCAATGACCTGTCTGCGCGAGCTGGGGCCAGAACCGAt                                                    <  1:2913504/71‑1 (MQ=255)
           ttGGACACGGAAGTCGGCGCGGGAACCTCTCACCCAATGACCTGTCTGCGCGAGCTGGGGCCAGAACCGAt                                                    <  1:5014180/71‑1 (MQ=255)
               aCATGGAAGTCGGCGCGGGAACCTCTCACCCAATGACCTGTCTGCGCGAGCTGGGGCCAGAACCGa                                                     <  1:4756912/66‑1 (MQ=255)
                   ggAAGTCGGCGCGGGAACCTCTCACCCAATGACCTGTCTGCGCGGGCTGGGGCCAGAACCGATGGCGgctg                                            >  1:544761/1‑71 (MQ=255)
                   ggAAGTCGGCGCGGGAACCTCTCACCCAATGACCTGTCTGCGCGGGCTGGGGCCAGAACCGATGGCGgctg                                            >  1:4171129/1‑71 (MQ=255)
                     aaGTCGGCGCGGGAACCTCTCACCCAATGACCTGTCTGCGCGAGCTGGGGCCAGAACCGATGGCGgctgct                                          <  1:2332488/71‑1 (MQ=255)
                       gTCGGCGCGGGAACCTCTCACCCAATGACCTGTCTGCGCGAGc                                                                    <  1:1670094/43‑1 (MQ=255)
                        tCGGCGCGGGAACCTCTCACCCAATGACCTGTCTGCGCGGGCTGGGGCCAGAAcc                                                       <  1:2777322/55‑1 (MQ=255)
                        tCGGCGCGGGAACCTCTCACCCAATGACCTGTCTGCGCGGGCTGGGGCCAGAAcc                                                       <  1:2436518/55‑1 (MQ=255)
                          ggCGCGGGAACCTCTCACCCAATGACCTGTCTGCGGGGGCTGGGGCCAGAACCGATGGCGGCTGCTTATg                                      <  1:2722932/70‑1 (MQ=255)
                          ggCGCGGGAACCTCTCACCCAATGACCTGTCTGCGCGGGCTGGGGCCAGAACCGATGGCGGCTGCTTATg                                      <  1:636914/70‑1 (MQ=255)
                          ggCGCGGGAACCTCTCACCCAATGACCTGTCTGCGCGAGCTGGGGCCAGAACCGATGGCGGCTGCTTATg                                      <  1:2252079/70‑1 (MQ=255)
                             gcgGGAACCTCTCACCCAATGACCTGTCTGCGCGAGCTGGGGCCAGAACCGATGGCGGCTGCTTATGTTCa                                  <  1:4455536/71‑1 (MQ=255)
                             gcgGGAACCTCTCACCCAATGACCTGTCTGCGCGAGCTGGGGCCAGAACCGATGGCGGCTGCTTATGTTCa                                  <  1:2942134/71‑1 (MQ=255)
                                 gAACCTCTCACCCAATGACCTGTCTGCGCGGGCTGGGGCCAGAACCGATGGCGGCTGCTTATGTTCAGCCt                              >  1:2809654/1‑71 (MQ=255)
                                     ctctCACCCAATGACCTGTCTTCGCGGGCTGGGGCCAGAACCGATGGCGGCTGCTTATGTTCAGCCTtctc                          <  1:1346710/71‑1 (MQ=255)
                                     ctctCACCCAATGACCTGTCTGCGCGAGCTGGGGCCAGAACCGATGGCGGCTGCTTATg                                      <  1:4388981/59‑1 (MQ=255)
                                       ctcACCCAATGACCTGTCTGCGCGAGCTGGGGCCAGAACCGATGGCGgctg                                            >  1:2927088/1‑51 (MQ=255)
                                       ctcACCCAATGACCTGTCTGCGCGAGCTGGGGCCAGAACCGATGGCGgctg                                            >  1:1103514/1‑51 (MQ=255)
                                       ctcACCCAATGACCTGTCTGCGCGAGCTGGGGCCAGAACCGATGGCGgctg                                            >  1:226846/1‑51 (MQ=255)
                                       ctcACCCAATGACCTGTCTGCGCGAGCTGGGGCCAGAACCGATGGCGgctg                                            >  1:4940333/1‑51 (MQ=255)
                                       ctcACCCAATGACCTGTCTGCGCGAGCTGGGGCCAGAACCGATGGCGgctg                                            >  1:3090733/1‑51 (MQ=255)
                                       ctcACCCAATGACCTGTCTGCGCGAGCTGGGGCCAGAACCGATGGCGGCTGCTTATGTTCAGCCTTCtcgt                        >  1:3284374/1‑71 (MQ=255)
                                            ccAATGACCTGTCTGCGCGGGCTGGGGCCAGAACCGATGGCGGCTGCTTATGTTCAGCCTTCTCGTCGccc                   >  1:4725686/1‑71 (MQ=255)
                                            ccAATGACCTGTCTGCGCGAGCTGGGGCCAGAACCGATGGCGGCTGCTTATGTTCAGCCTTCTCGTCGccc                   >  1:4119435/1‑71 (MQ=255)
                                            ccAATGACCTGTCTGCGCGAGCTGGGGCCAGAACCGATGGCGGCTGCTTATGTTCAGCCTTCTCGTCGccc                   >  1:2892561/1‑71 (MQ=255)
                                               aTGACCTGTCTGCGCGATCTGGGGCCAGAACCGATGGCGGCTGCTTATGTTCAGc                                >  1:2979030/1‑55 (MQ=255)
                                               aTGACCTGTCTGCGCGAGCTGGGGCCAGAACCGATGGCGGCTGCTTATGTTCAGCCTtctc                          <  1:3566218/61‑1 (MQ=255)
                                               aTGACCTGTCTGCGCGAGCTGGGGCCAGAACCGATGGCGGCTGCTTATGTTCAGCCTTCTCGTCGCccgac                >  1:3618838/1‑71 (MQ=255)
                                                 gACCTGTCTGCGCGAGCTGGGGCCAGAACCGATGGCGGCTGCTTATGTTCAGCCTTCtcgt                        >  1:2982550/1‑61 (MQ=255)
                                                       tCTGCGCGGGCTGGGGCCAGAACCGATGGCGGCTGCTTAt                                       <  1:2847011/40‑1 (MQ=255)
                                                       tCTGCGCGAGCTGGGGCCAGAACCGATGGCGGCTGCTTATGTTCAGCCTTCTCGTCGCCCGACCGATGGTc        >  1:2417448/1‑71 (MQ=255)
                                                         tGCTCGATCTGGCGCCAGAACCGATGGCGGCTGCTTATGTTCAGcc                               >  1:2347602/1‑46 (MQ=38)
                                                         tGCGCGAGCTGGGGCCAGAACCGATGGCGGCTGCTTATGTTCAGCCTTCTCGTCGCCCGACCGATGGTCGc      >  1:3046975/1‑71 (MQ=255)
                                                           cgcgAGCTGGGGCCAGAACCGATGGCGGCTGCTTATGTTCAGCCTTCTCGTCGCCCGACCGATGGTCGc      <  1:2718205/69‑1 (MQ=255)
                                                             cgAGCTGGGGCCAGAACCGATGGCGGCTGCTTATGTTCAGCCTTCTCGTCGCCCGACCGATGGTCGCTACg  >  1:3672868/1‑71 (MQ=255)
                                                             cgAGCTGGGGCCAGAACCGATGGCGGCTGCTTATGTTCAGCCTTCTCGTCGCCCGACCGATGGTCGCTACg  >  1:3166786/1‑71 (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: