Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 531,823 G→A 100% I247I (ATC→ATT ybhR ← predicted transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE531,8230GA79.6% 98.5 / 20.4 49I247I (ATC→ATTybhRpredicted transporter subunit
Reads supporting (aligned to +/- strand):  ref base G (0/10);  new base A (39/0);  total (39/10)
Fisher's exact test for biased strand distribution p-value = 1.22e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATACGCCCAGATACCAATCGCCAGCACAATGGTGGC  >  minE/531787‑531864
                                    |                                         
cGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATACGCCCAGATACCAATCGCCAGc             <  1:875240/67‑1 (MQ=255)
cGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATACGCCCAGATACCAATCGCCAGc             <  1:207359/67‑1 (MQ=255)
cGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATACGCCCAGATACCAATCGCCAGc             <  1:691709/67‑1 (MQ=255)
cGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATACGCCCAGATACCAATCGCCAGc             <  1:2190269/67‑1 (MQ=255)
cGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATACGCCCAGATACCAATCGCCAGc             <  1:608723/67‑1 (MQ=255)
cGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATACGCCCAGATACCAATCGCCAGc             <  1:1839356/67‑1 (MQ=255)
cGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATACGCCCAGATACCAATCGCCAGc             <  1:508378/67‑1 (MQ=255)
cGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATACGCCCAGATACCAATCGCCAGc             <  1:2247640/67‑1 (MQ=255)
cGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATACGCCCAGATACCAATCGCCAGc             <  1:1611486/67‑1 (MQ=255)
 gTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATACGCCCAGATACCAATCGCCAGc             <  1:172740/66‑1 (MQ=255)
                            gCGAAGGGAATTTGATACGCCCAGATACCAATCGCCAGCACAATGGTGGc  >  1:289221/1‑50 (MQ=255)
                            gCGAAGGGAATTTGATACGCCCAGATACCAATCGCCAGCACAATGGTGGc  >  1:2128217/1‑50 (MQ=255)
                            gCGAAGGGAATTTGATACGCCCAGATACCAATCGCCAGCACAATGGTGGc  >  1:2164518/1‑50 (MQ=255)
                            gCGAAGGGAATTTGATACGCCCAGATACCAATCGCCAGCACAATGGTGGc  >  1:2200237/1‑50 (MQ=255)
                            gCGAAGGGAATTTGATACGCCCAGATACCAATCGCCAGCACAATGGTGGc  >  1:2220985/1‑50 (MQ=255)
                            gCGAAGGGAATTTGATACGCCCAGATACCAATCGCCAGCACAATGGTGGc  >  1:2227308/1‑50 (MQ=255)
                            gCGAAGGGAATTTGATACGCCCAGATACCAATCGCCAGCACAATGGTGGc  >  1:2280098/1‑50 (MQ=255)
                            gCGAAGGGAATTTGATACGCCCAGATACCAATCGCCAGCACAATGGTGGc  >  1:287337/1‑50 (MQ=255)
                            gCGAAGGGAATTTGATACGCCCAGATACCAATCGCCAGCACAATGGTGGc  >  1:1936637/1‑50 (MQ=255)
                            gCGAAGGGAATTTGATACGCCCAGATACCAATCGCCAGCACAATGGTGGc  >  1:331561/1‑50 (MQ=255)
                            gCGAAGGGAATTTGATACGCCCAGATACCAATCGCCAGCACAATGGTGGc  >  1:489504/1‑50 (MQ=255)
                            gCGAAGGGAATTTGATACGCCCAGATACCAATCGCCAGCACAATGGTGGc  >  1:54682/1‑50 (MQ=255)
                            gCGAAGGGAATTTGATACGCCCAGATACCAATCGCCAGCACAATGGTGGc  >  1:586576/1‑50 (MQ=255)
                            gCGAAGGGAATTTGATACGCCCAGATACCAATCGCCAGCACAATGGTGGc  >  1:637032/1‑50 (MQ=255)
                            gCGAAGGGAATTTGATACGCCCAGATACCAATCGCCAGCACAATGGTGGc  >  1:691318/1‑50 (MQ=255)
                            gCGAAGGGAATTTGATACGCCCAGATACCAATCGCCAGCACAATGGTGGc  >  1:792702/1‑50 (MQ=255)
                            gCGAAGGGAATTTGATACGCCCAGATACCAATCGCCAGCACAATGGTGGc  >  1:877170/1‑50 (MQ=255)
                            gCGAAGGGAATTTGATACGCCCAGATACCAATCGCCAGCACAATGGTGGc  >  1:1480232/1‑50 (MQ=255)
                            gCGAAGGGAATTTGATACGCCCAGATACCAATCGCCAGCACAATGGTGGc  >  1:1010878/1‑50 (MQ=255)
                            gCGAAGGGAATTTGATACGCCCAGATACCAATCGCCAGCACAATGGTGGc  >  1:1023113/1‑50 (MQ=255)
                            gCGAAGGGAATTTGATACGCCCAGATACCAATCGCCAGCACAATGGTGGc  >  1:1056706/1‑50 (MQ=255)
                            gCGAAGGGAATTTGATACGCCCAGATACCAATCGCCAGCACAATGGTGGc  >  1:1091361/1‑50 (MQ=255)
                            gCGAAGGGAATTTGATACGCCCAGATACCAATCGCCAGCACAATGGTGGc  >  1:1156542/1‑50 (MQ=255)
                            gCGAAGGGAATTTGATACGCCCAGATACCAATCGCCAGCACAATGGTGGc  >  1:122105/1‑50 (MQ=255)
                            gCGAAGGGAATTTGATACGCCCAGATACCAATCGCCAGCACAATGGTGGc  >  1:1334185/1‑50 (MQ=255)
                            gCGAAGGGAATTTGATACGCCCAGATACCAATCGCCAGCACAATGGTGGc  >  1:1363828/1‑50 (MQ=255)
                            gCGAAGGGAATTTGATACGCCCAGATACCAATCGCCAGCACAATGGTGGc  >  1:1404023/1‑50 (MQ=255)
                            gCGAAGGGAATTTGATACGCCCAGATACCAATCGCCAGCACAATGGTGGc  >  1:2080651/1‑50 (MQ=255)
                            gCGAAGGGAATTTGATACGCCCAGATACCAATCGCCAGCACAATGGTGGc  >  1:1491069/1‑50 (MQ=255)
                            gCGAAGGGAATTTGATACGCCCAGATACCAATCGCCAGCACAATGGTGGc  >  1:1605761/1‑50 (MQ=255)
                            gCGAAGGGAATTTGATACGCCCAGATACCAATCGCCAGCACAATGGTGGc  >  1:1695385/1‑50 (MQ=255)
                            gCGAAGGGAATTTGATACGCCCAGATACCAATCGCCAGCACAATGGTGGc  >  1:1699799/1‑50 (MQ=255)
                            gCGAAGGGAATTTGATACGCCCAGATACCAATCGCCAGCACAATGGTGGc  >  1:1744031/1‑50 (MQ=255)
                            gCGAAGGGAATTTGATACGCCCAGATACCAATCGCCAGCACAATGGTGGc  >  1:1795233/1‑50 (MQ=255)
                            gCGAAGGGAATTTGATACGCCCAGATACCAATCGCCAGCACAATGGTGGc  >  1:183164/1‑50 (MQ=255)
                            gCGAAGGGAATTTGATACGCCCAGATACCAATCGCCAGCACAATGGTGGc  >  1:1885542/1‑50 (MQ=255)
                            gCGAAGGGAATTTGATACGCCCAGATACCAATCGCCAGCACAATGGTGGc  >  1:1904525/1‑50 (MQ=255)
                            gCGAAGGGAATTTGATACGCCCAGATACCAATCGCCAGCACAATGGTGGc  >  1:10048/1‑50 (MQ=255)
                            gCGAAGGGAATTTGATACGCCCAGATACAAATCGCCAGCACAAtggt     >  1:1450790/1‑47 (MQ=39)
                                    |                                         
CGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATACGCCCAGATACCAATCGCCAGCACAATGGTGGC  >  minE/531787‑531864

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