Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 410,370 A→G 60.9% G662G (GGT→GGC leuS ← leucyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE410,3700AG60.9% 21.1 / 27.4 46G662G (GGT→GGCleuSleucyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base A (10/8);  new base G (20/8);  total (30/16)
Fisher's exact test for biased strand distribution p-value = 3.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.37e-01

GTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGCCCCTTCCACACCGGATTCCTGCCATTCGAGAGTCATATCAGCCGGAGAAGCAAACATCATAAACA  >  minE/410314‑410425
                                                        |                                                       
gtgtgCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGCCCCTTCCACACCGGATTCCTGc                                             <  1:322762/69‑1 (MQ=255)
   tgCTCGTAAACCAGTTTCCAGACACGTTTCAGGAGGCGGTTAGCCCCTTCCACACCGGATTCCTGCCAt                                          >  1:2182725/1‑69 (MQ=255)
   tgCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGCCCCTTCCACACCGGATTCCTGCCAt                                          >  1:1386713/1‑69 (MQ=255)
   tgCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGCCCCTTCCACACCGGATTCCTGCCAt                                          >  1:2942226/1‑69 (MQ=255)
   tgCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGCCCCTTCCACACCGGATTCCTGCCAt                                          >  1:909814/1‑69 (MQ=255)
                   tCCAGACACGTTTCAGGAAGCGGTTAGCCCCTTCCACACCGGAt                                                   >  1:3251271/1‑44 (MQ=255)
                   tCCAGACACGTTTCAGGAAGCGGTTAGCCCCTTCCACACCGGAt                                                   >  1:805770/1‑44 (MQ=255)
                   tCCAGACACGTTTCAGGAAGCGGTTAGCCCCTTCCACACCGGAt                                                   <  1:800046/44‑1 (MQ=255)
                   tCCAGACACGTTTCAGGAAGCGGTTAGCCCCTTCCACACCGGAt                                                   <  1:1787496/44‑1 (MQ=255)
                   tCCAGACACGTTTCAGGAAGCGGTTAGCCCCTTCCACACCGGAt                                                   <  1:1882417/44‑1 (MQ=255)
                   tCCAGACACGTTTCAGGAAGCGGTTAGCCCCTTCCACACCGGAt                                                   >  1:500270/1‑44 (MQ=255)
                       gACACGTTTCAGGAAGCGGTTAGCCCCTTCCACGCCGGATTCCTGCCATTCGAGAGTCATATCAGCCgg                      <  1:2273883/69‑1 (MQ=255)
                       gACACGTTTCAGGAAGCGGTTAGCCCCTTCCACGCCGGATTCCTGCCATTCGAGAGTCATATCAGCCgg                      <  1:1001407/69‑1 (MQ=255)
                       gACACGTTTCAGGAAGCGGTTAGCCCCTTCCACGCCGGATTCCTGCCATTCGAGAGTCATATCAGCCgg                      <  1:630420/69‑1 (MQ=255)
                       gACACGTTTCAGGAAGCGGTTAGCCCCTTCCACGCCGGATTCCTGCCATTCGAGAGTCATATCAGCCgg                      <  1:776119/69‑1 (MQ=255)
                       gACACGTTTCAGGAAGCGGTTAGCCCCTTCCACGCCGGATTCCTGCCATTCGAGAGTCATATCAGCCgg                      <  1:2367564/69‑1 (MQ=255)
                       gACACGTTTCAGGAAGCGGTTAGCCCCTTCCACGCCGGATTCCTGCCATTCGAGAGTCATATCAGCCgg                      <  1:1068114/69‑1 (MQ=255)
                       gACACGTTTCAGGAAGCGGTTAGCCCCTTCCACGCCGGATTCCTGCCATTCGAGAGTCATATCAGCCgg                      <  1:1886955/69‑1 (MQ=255)
                       gACACGTTTCAGGAAGCGGTTAGCCCCTTCCACGCCGGATTCCTGCCATTCGAGAGTCATATCAGCCgg                      <  1:1281444/69‑1 (MQ=255)
                              ttCAGGAAGCGGTTAGCCCCTTCCACGCCGGATTCCTGCCATTCGAGAGTCATATCAGCCGGAGAAGc                >  1:2814417/1‑68 (MQ=255)
                              ttCAGGAAGCGGTTAGCCCCTTCCACGCCGGATTCCTGCCATTCGAGAGTCATATCAGCCGGAGAAGc                >  1:989278/1‑68 (MQ=255)
                              ttCAGGAAGCGGTTAGCCCCTTCCACGCCGGATTCCTGCCATTCGAGAGTCATATCAGCCGGAGAAGCa               >  1:1035246/1‑69 (MQ=255)
                              ttCAGGAAGCGGTTAGCCCCTTCCACGCCGGATTCCTGCCATTCGAGAGTCATATCAGCCGGAGAAGCa               >  1:605115/1‑69 (MQ=255)
                              ttCAGGAAGCGGTTAGCCCCTTCCACGCCGGATTCCTGCCATTCGAGAGTCATATCAGCCGGAGAAGCa               >  1:499764/1‑69 (MQ=255)
                              ttCAGGAAGCGGTTAGCCCCTTCCACGCCGGATTCCTGCCATTCGAGAGTCATATCAGCCGGAGAAGCa               >  1:949381/1‑69 (MQ=255)
                              ttCAGGAAGCGGTTAGCCCCTTCCACGCCGGATTCCTGCCATTCGAGAGTCATATCAGCCGGAGAAGCa               >  1:3381766/1‑69 (MQ=255)
                              ttCAGGAAGCGGTTAGCCCCTTCCACGCCGGATTCCTGCCATTCGAGAGTCATATCAGCCGGAGAAGCa               >  1:332800/1‑69 (MQ=255)
                              ttCAGGAAGCGGTTAGCCCCTTCCACGCCGGATTCCTGCCATTCGAGAGTCATATCAGCCGGAGAAGCa               >  1:2776140/1‑69 (MQ=255)
                              ttCAGGAAGCGGTTAGCCCCTTCCACGCCGGATTCCTGCCATTCGAGAGTCATATCAGCCGGAGAAGCa               >  1:1273758/1‑69 (MQ=255)
                              ttCAGGAAGCGGTTAGCCCCTTCCACGCCGGATTCCTGCCATTCGAGAGTCATATCAGCCGGAGAAGCa               >  1:1401039/1‑69 (MQ=255)
                              ttCAGGAAGCGGTTAGCCCCTTCCACGCCGGATTCCTGCCATTCGAGAGTCATATCAGCCGGAGAAGCa               >  1:1633610/1‑69 (MQ=255)
                              ttCAGGAAGCGGTTAGCCCCTTCCACGCCGGATTCCTGCCATTCGAGAGTCATATCAGCCGGAGAAGCa               >  1:1734053/1‑69 (MQ=255)
                              ttCAGGAAGCGGTTAGCCCCTTCCACGCCGGATTCCTGCCATTCGAGAGTCATATCAGCCGGAGAAGCa               >  1:247641/1‑69 (MQ=255)
                              ttCAGGAAGCGGTTAGCCCCTTCCACGCCGGATTCCTGCCATTCGAGAGTCATATCAGCCGGAGAAGCa               >  1:2696075/1‑69 (MQ=255)
                              ttCAGGAAGCGGTTAGCCCCTTCCACGCCGGATTCCTGCCATTCGAGAGTCATATCAGCCGGAGAAGCa               >  1:3160128/1‑69 (MQ=255)
                              ttCAGGAAGCGGTTAGCCCCTTCCACGCCGGATTCCTGCCATTCGAGAGTCATATCAGCCGGAGAAGCa               >  1:279146/1‑69 (MQ=255)
                              ttCAGGAAGCGGTTAGCCCCTTCCACGCCGGATTCCTGCCATTCGAGAGTCATATCAGCCGGAGAAGCa               >  1:3065416/1‑69 (MQ=255)
                              ttCAGGAAGCGGTTAGCCCCTTCCACGCCGGATTCCTGCCATTCGAGAGTCATATAAGCCGGAGAAGCa               >  1:485440/1‑69 (MQ=255)
                              atCAGGAAGCGGTTAGCCCCTTCAACGCCGGATTCCTGCCATTAGAGAGTCATATCAGCCGgaga                   >  1:1393398/2‑65 (MQ=255)
                                     aGCGGTTAGCCCCTTCCACACCGGATTCCTGCCATTCGAGAGTCATATCAGCCGGAGAAGCAAAcatc         >  1:3051133/1‑68 (MQ=255)
                                           tAGCCCCTTCCACACCGGATTCCTGCCATTCGAGAGTCATATCAGCCGGAGAAGCAAACATCATAAACa  <  1:2814055/69‑1 (MQ=255)
                                           tAGCCCCTTCCACACCGGATTCCTGCCATTCGAGAGTCATATCAGCCGGAGAAGCAAACATCATAAACa  <  1:3049389/69‑1 (MQ=255)
                                           tAGCCCCTTCCACACCGGATTCCTGCCATTCGAGAGTCATATCAGCCGGAGAAGCAAACATCATAAACa  <  1:809209/69‑1 (MQ=255)
                                            aGCCCCTTCCACACCGGATTCCTGCCATTCGAGAGTCATATCAGCCGGAGAAGCAAACATCATAAACa  <  1:3137525/68‑1 (MQ=255)
                                                ccTTCCACACCGGATTCCTGCCATTCGAGAGTCATATCAGCCGGAGAAGCAAACATCATa      >  1:716907/1‑60 (MQ=255)
                                                ccTTCCACACCGGATTCCTGCCATTCGAGAGTCATATCAGCCGGAGAAGCAAACATCATa      >  1:3353661/1‑60 (MQ=255)
                                                        |                                                       
GTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGCCCCTTCCACACCGGATTCCTGCCATTCGAGAGTCATATCAGCCGGAGAAGCAAACATCATAAACA  >  minE/410314‑410425

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: