Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,369,034 A→G 67.9% R127R (CGT→CGC dusC ← tRNA‑dihydrouridine synthase C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,369,0340AG67.9% 50.3 / 27.3 59R127R (CGT→CGCdusCtRNA‑dihydrouridine synthase C
Reads supporting (aligned to +/- strand):  ref base A (9/10);  new base G (26/14);  total (35/24)
Fisher's exact test for biased strand distribution p-value = 2.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.83e-01

AGCCCAGACGCACTTTCACGCTGACGGGCAAATGCGCCGGTACAGCTTCACGCATCGCTTTTGCACCCTGGTAGATGAGTTCAGGATCTTTGAGTAACGTCGCCCCGCCG  >  minE/1368985‑1369094
                                                 |                                                            
agcccagACGCACTTTTACGCTGACAGGCAAATGCGCCGGTACAGCTTCGCGCATCGCTTTTGCACCCt                                           >  1:105164/1‑69 (MQ=255)
agcccagACGCACTTTTACGCTGACAGGCAAATGCGCCGGTACAGCTTCGCGCATCGCTTTTGCACCCt                                           >  1:1514656/1‑69 (MQ=255)
agcccagACGCACTTTTACGCTGACAGGCAAATGCGCCGGTACAGCTTCGCGCATCGCTTTTGCACCCt                                           >  1:86066/1‑69 (MQ=255)
agcccagACGCACTTTTACGCTGACAGGCAAATGCGCCGGTACAGCTTCGCGCATCGCTTTTGCACCCt                                           >  1:1289747/1‑69 (MQ=255)
agcccagACGCACTTTTACGCTGACAGGCAAATGCGCCGGTACAGCTTCGCGCATCGCTTTTGCACCCt                                           >  1:892749/1‑69 (MQ=255)
agcccagACGCACTTTTACGCTGACAGGCAAATGCGCCGGTACAGCTTCGCGCATCGCTTTTGCACCCt                                           >  1:123344/1‑69 (MQ=255)
agcccagACGCACTTTTACGCTGACAGGCAAATGCGCCGGTACAGCTTCGCGCATCGCTTTTGCACCCt                                           >  1:1130756/1‑69 (MQ=255)
    cagACGCACTTTTACGCTGACAGGCAAATGCGCCGGTACAGCTTCGCGCATCGCTTTTGCACCCTGGTa                                       <  1:876853/69‑1 (MQ=255)
    cagACGCACTTTTACGCTGACAGGCAAATGCGCCGGTACAGCTTCGCGCATCGCTTTTGCACCCTGGTa                                       <  1:2543595/69‑1 (MQ=255)
    cagACGCACTTTTACGCTGACAGGCAAATGCGCCGGTACAGCTTCGCGCATCGCTTTTGCACCCTGGTa                                       <  1:320277/69‑1 (MQ=255)
    cagACGCACTTTTACGCTGACAGGCAAATGCGCCGGTACAGCTTCGCGCATCGCTTTTGCACCCTGGTa                                       <  1:2239519/69‑1 (MQ=255)
    cagACGCACTTTTACGCTGACAGGCAAATGCGCCGGTACAGCTTCGCGCATCGCTTTTGCACCCTGGTa                                       <  1:2582926/69‑1 (MQ=255)
    cagACGCACTTTTACGCTGACAGGCAAATGCGCCGGTACAGCTTCGCGCATCGCTTTTGCACCCTGGTa                                       <  1:2584032/69‑1 (MQ=255)
    cagACGCACTTTTACGCTGACAGGCAAATGCGCCGGTACAGCTTCGCGCATCGCTTTTGCACCCTGGTa                                       <  1:3185724/69‑1 (MQ=255)
    cagACGCACTTTTACGCTGACAGGCAAATGCGCCGGTACAGCTTCGCGCATCGCTTTTGCACCCTGGTa                                       <  1:2919236/69‑1 (MQ=255)
    cagACGCACTTTTACGCTGACAGGCAAATGCGCCGGTACAGCTTCGCGCATCGCTTTTGCACCCTGGTa                                       <  1:592353/69‑1 (MQ=255)
    cagACGCACTTTTACGCTGACAGGCAAATGCGCCGGTACAGCTTCGCGCATCGCTTTTGCACCCTGGTa                                       <  1:678009/69‑1 (MQ=255)
    cagACGCACTTTTACGCTGACAGGCAAATGCGCCGGTACAGCTTCGCGCATCGCTTTTGCACCCTGGTa                                       <  1:3092281/69‑1 (MQ=255)
    cagACGCACTTTTACGCTGACAGGCAAATGCGCCGGTACAGCTTCGCGCATCGCTTTTGCACCCTGGTa                                       <  1:1283950/69‑1 (MQ=255)
    cagACGCACTTTTACGCTGACAGGCAAATGCGCCAGTACAGCTTCGCGCATCGCTTTTGCACCCTGGTa                                       <  1:2982485/69‑1 (MQ=255)
    cagACGCACTTTCACGCTGACGGGCAAATGCGCCGGTACAGCTTCACGCATCGCTTTTGCACCCTGGTa                                       <  1:257481/69‑1 (MQ=255)
     agACGCACTTTTACTCTGACAGGCAAATGCGCCGGTACAGCTTCGCGCATCGCTTTTGCACCCTGGTa                                       <  1:2792646/68‑1 (MQ=255)
     agACGCACTTTCACGCTGACGGGCAAATGCGCCGGTACAGCTTCACGCATCGCTTTTGCACCCTGGTa                                       <  1:502713/68‑1 (MQ=255)
     agACGCACTTTCACGCTGACGGGCAAATGCGCCGGTACAGCTTCACGCATCGCTTTTGCACCCTGGTa                                       <  1:477441/68‑1 (MQ=255)
                       aCGGGCAAATGCGCCGGTACAGCTTCACGCATCGCTTTTGCACCCTgg                                         <  1:298602/48‑1 (MQ=255)
                       aCGGGCAAATGCGCCGGTACAGCTTCACGCATCGCTTTTGCACCCTgg                                         <  1:3139921/48‑1 (MQ=255)
                       aCGGGCAAATGCGCCGGTACAGCTTCACGCATCGCTTTTGCACCCTgg                                         <  1:2062419/48‑1 (MQ=255)
                       aCGGGCAAATGCGCCGGTACAGCTTCACGCATCGCTTTTGCACCCTgg                                         <  1:1906090/48‑1 (MQ=255)
                       aCGGGCAAATGCGCCGGTACAGCTTCACGCATCGCTTTTGCACCCTgg                                         <  1:1493310/48‑1 (MQ=255)
                       aCGGGCAAATGCGCCGGTACAGCTTCACGCATCGCTTTTGCACCCTgg                                         <  1:1450205/48‑1 (MQ=255)
                       aCGGGCAAATGCGCCGGTACAGCTTCACGCATCGCTTTTGCACCCTgg                                         <  1:1431246/48‑1 (MQ=255)
                                         aCAGCTTCGCGCATCGCTTTTGCACCCTGGTAGATGAGTTCAGGATCTTTCAGTAACGTTGCCccacca  >  1:3271070/1‑68 (MQ=255)
                                         aCAGCTTCGCGCATCGCTTTTGCACCCTGGTAGATGAGTTCAGGATCTTTCAGTAACGTTGCCccacca  >  1:3344022/1‑68 (MQ=255)
                                         aCAGCTTCGCGCATCGCTTTTGCACCCTGGTAGATGAGTTCAGGATCTTTCAGTAACGTTGCCccacca  >  1:3200418/1‑68 (MQ=255)
                                         aCAGCTTCGCGCATCGCTTTTGCACCCTGGTAGATGAGTTCAGGATCTTTCAGTAACGTTGCCccacca  >  1:349819/1‑68 (MQ=255)
                                         aCAGCTTCGCGCATCGCTTTTGCACCCTGGTAGATGAGTTCAGGATCTTTCAGTAACGTTGCCccacca  >  1:3179600/1‑68 (MQ=255)
                                         aCAGCTTCGCGCATCGCTTTTGCACCCTGGTAGATGAGTTCAGGATCTTTCAGTAACGTTGCCccacca  >  1:480773/1‑68 (MQ=255)
                                         aCAGCTTCGCGCATCGCTTTTGCACCCTGGTAGATGAGTTCAGGATCTTTCAGTAACGTTGCCccacca  >  1:1112175/1‑68 (MQ=255)
                                         aCAGCTTCGCGCATCGCTTTTGCACCCTGGTAGATGAGTTCAGGATCTTTCAGTAACGTTGCCccacca  >  1:2932095/1‑68 (MQ=255)
                                         aCAGCTTCGCGCATCGCTTTTGCACCCTGGTAGATGAGTTCAGGATCTTTCAGTAACGTTGCCccacca  >  1:1234415/1‑68 (MQ=255)
                                         aCAGCTTCGCGCATCGCTTTTGCACCCTGGTAGATGAGTTCAGGATCTTTCAGTAACGTTGCCccacca  >  1:1458404/1‑68 (MQ=255)
                                         aCAGCTTCGCGCATCGCTTTTGCACCCTGGTAGATGAGTTCAGGATCTTTCAGTAACGTTGCCccacca  >  1:1523044/1‑68 (MQ=255)
                                         aCAGCTTCGCGCATCGCTTTTGCACCCTGGTAGATGAGTTCAGGATCTTTCAGTAACGTTGCCccacca  >  1:1659943/1‑68 (MQ=255)
                                         aCAGCTTCGCGCATCGCTTTTGCACCCTGGTAGATGAGTTCAGGATCTTTCAGTAACGTTGCCccacca  >  1:1683271/1‑68 (MQ=255)
                                         aCAGCTTCGCGCATCGCTTTTGCACCCTGGTAGATGAGTTCAGGATCTTTCAGTAACGTTGCCccacca  >  1:1722110/1‑68 (MQ=255)
                                         aCAGCTTCGCGCATCGCTTTTGCACCCTGGTAGATGAGTTCAGGATCTTTCAGTAACGTTGCCccacca  >  1:2045253/1‑68 (MQ=255)
                                         aCAGCTTCGCGCATCGCTTTTGCACCCTGGTAGATGAGTTCAGGATCTTTCAGTAACGTTGCCccacca  >  1:2236776/1‑68 (MQ=255)
                                         aCAGCTTCGCGCATCGCTTTTGCACCCTGGTAGATGAGTTCAGGATCTTTCAGTAACGTTGCCccacca  >  1:2750661/1‑68 (MQ=255)
                                         aCAGCTTCGCGCATCGCTTTTGCACCCTGGTAGATGAGTTCAGGATCTTTCAGTAACGTTGCCccacc   >  1:2345718/1‑68 (MQ=255)
                                         aCAGCTTCGCGCATCGCTTTTGCACCCTGGTAGATGAGTTCAGGATCTTTCAGTAACGTTGCCccacc   >  1:3154599/1‑68 (MQ=255)
                                         aCAGCTTCACGCATCGCTTTTGCACCCTGGTAGATGAGTTCAGGATCTTTGAGTAACGTCGCCcctccg  >  1:3270951/1‑69 (MQ=255)
                                         aCAGCTTCACGCATCGCTTTTGCACCCTGGTAGATGAGTTCAGGATCTTTGAGTAACGTCGCCccgccg  >  1:2246331/1‑69 (MQ=255)
                                         aCAGCTTCACGCATCGCTTTTGCACCCTGGTAGATGAGTTCAGGATCTTTGAGTAACGTCGCCccgccg  >  1:340112/1‑69 (MQ=255)
                                         aCAGCTTCACGCATCGCTTTTGCACCCTGGTAGATGAGTTCAGGATCTTTGAGTAACGTCGCCccgccg  >  1:2590416/1‑69 (MQ=255)
                                         aCAGCTTCACGCATCGCTTTTGCACCCTGGTAGATGAGTTCAGGATCTTTGAGTAACGTCGCCccgccg  >  1:1347403/1‑69 (MQ=255)
                                         aCAGCTTCACGCATCGCTTTTGCACCCTGGTAGATGAGTTCAGGATCTTTGAGTAACGTCGCCccgccg  >  1:880216/1‑69 (MQ=255)
                                         aCAGCTTCACGCATCGCTTTTGCACCCTGGTAGATGAGTTCAGGATCTTTGAGTAACGTCGCCccgccg  >  1:2954874/1‑69 (MQ=255)
                                         aCAGCTTCACGCATCGCTTTTGCACCCTGGTAGATGAGTTCAGGATCTTTGAGTAACGTCGCCccgcc   >  1:2820399/1‑68 (MQ=255)
                                         aCAGCTTCACGCATCGCTTTTGCACCCTGGTAGATGAGTTCAGGATCTTTGAGTAACGTCGCCccgcc   >  1:1259940/1‑68 (MQ=255)
                                                 |                                                            
AGCCCAGACGCACTTTCACGCTGACGGGCAAATGCGCCGGTACAGCTTCACGCATCGCTTTTGCACCCTGGTAGATGAGTTCAGGATCTTTGAGTAACGTCGCCCCGCCG  >  minE/1368985‑1369094

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: