Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 667,181 T→C 36.5% S173S (AGT→AGC uup → fused predicted transporter subunits and ATP‑binding components of ABC superfamily

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE667,1810TC36.5% 43.8 / 37.5 63S173S (AGT→AGCuupfused predicted transporter subunits and ATP‑binding components of ABC superfamily
Reads supporting (aligned to +/- strand):  ref base T (12/28);  new base C (14/9);  total (26/37)
Fisher's exact test for biased strand distribution p-value = 3.23e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGTCGTCGCTTTCCGGCGGCTGGTTGCGTAAAGCGGCATTAGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAACCACCTGGATATTGAAACCATCGACT  >  minE/667122‑667242
                                                           |                                                             
cTGTCGTCGCTTTCCGGCGGCTGGTTGCGTAAAGCGGCATTAGGACGCGCGCTGGTGAGTAATCcgcgc                                                      >  1:2596288/1‑69 (MQ=255)
cTGTCGTCGCTTTCCGGCGGCTGGTTGAGTAAATCGGCATTAGGACGCGCGCTGGTGAGTAATCcgcgc                                                      >  1:775733/1‑69 (MQ=255)
                             tAAAGCGGCATTAGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAacca                         <  1:754027/69‑1 (MQ=255)
                             tAAAGCGGCATTAGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAacca                         <  1:2383042/69‑1 (MQ=255)
                             tAAAGCGGCATTAGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAacca                         <  1:2140134/69‑1 (MQ=255)
                             tAAAGCGGCATTAGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAacca                         <  1:213087/69‑1 (MQ=255)
                             tAAAGCGGCATTAGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAacca                         <  1:93299/69‑1 (MQ=255)
                             tAAAGCGGCATTAGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAacca                         <  1:148241/69‑1 (MQ=255)
                                   ggCATTAGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAACCACCTGGa                   <  1:1612659/69‑1 (MQ=255)
                                   ggCATTAGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAACCACCTGGa                   <  1:2278230/69‑1 (MQ=255)
                                   ggCATTAGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAACCACCTGGa                   <  1:539736/69‑1 (MQ=255)
                                   ggCATTAGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAACCACCTGGa                   <  1:2218883/69‑1 (MQ=255)
                                   ggCATTAGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAACCACCTGGa                   <  1:432259/69‑1 (MQ=255)
                                   ggCATTAGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAACCACCTGGa                   <  1:2085175/69‑1 (MQ=255)
                                   ggCATTAGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAACCACCTGGa                   <  1:2082279/69‑1 (MQ=255)
                                   ggCATTAGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAACCACCTGGa                   <  1:2078241/69‑1 (MQ=255)
                                   ggCATTAGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAACCACCTGGa                   <  1:2910455/69‑1 (MQ=255)
                                   ggCATTAGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAACCACCTGGa                   <  1:1454586/69‑1 (MQ=255)
                                   ggCATTAGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAACCACCTGGa                   <  1:1901771/69‑1 (MQ=255)
                                   ggCATTAGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAACCACCTGGa                   <  1:3066584/69‑1 (MQ=255)
                                   ggCATTAGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAACCACCTGGa                   <  1:1587939/69‑1 (MQ=255)
                                         aGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAACCACCTGGATATTGa             >  1:278094/1‑69 (MQ=255)
                                         aGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAACCACCTGGATATTGa             >  1:479633/1‑69 (MQ=255)
                                         aGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAACCACCTGGATATTGa             >  1:2727131/1‑69 (MQ=255)
                                         aGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAACCACCTGGATATTGa             >  1:3176883/1‑69 (MQ=255)
                                         aGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAACCACCTGGATATTGa             >  1:2468840/1‑69 (MQ=255)
                                         aGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAACCACCTGGATATTGa             >  1:2422555/1‑69 (MQ=255)
                                         aGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAACCACCTGGATATTGa             >  1:1337583/1‑69 (MQ=255)
                                         aGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAACCACCTGGATATTGa             >  1:2177903/1‑69 (MQ=255)
                                         aGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAACCACCTGGATATTGa             >  1:2055440/1‑69 (MQ=255)
                                             cgcgcgCTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAACCACCTGGATATTGAAAcc         >  1:3154142/1‑69 (MQ=255)
                                                gcgcTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAACCACCTGGATATTGAAACCATc      <  1:961651/69‑1 (MQ=255)
                                                gcgcTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAACCACCTGGATATTGAAACCATc      <  1:700485/69‑1 (MQ=255)
                                                gcgcTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAACCACCTGGATATTGAAACCATc      <  1:696481/69‑1 (MQ=255)
                                                gcgcTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAACCACCTGGATATTGAAACCATc      <  1:1853036/69‑1 (MQ=255)
                                                gcgcTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAACCACCTGGATATTGAAACCATc      <  1:298459/69‑1 (MQ=255)
                                                gcgcTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAACCACCTGGATATTGAAACCATc      <  1:1987580/69‑1 (MQ=255)
                                                gcgcTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAACCACCTGGATATTGAAACCATc      <  1:2728768/69‑1 (MQ=255)
                                                gcgcTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAACCACCTGGATATTGAAACCATc      <  1:2499306/69‑1 (MQ=255)
                                                  gcTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAACCACCTGGATATTGAAACCATCGa    <  1:1053652/69‑1 (MQ=255)
                                                  gcTGGTGAGCAATCCGCGCGTATTGCTGCTCGACGAACCGACTAACCACCTGGATATTGAAACCATCGa    >  1:420757/1‑69 (MQ=18)
                                                  gcTGGTGAGCAATCCGCGCGTATTGCTGCTCGACGAACCGACTAACCACCTGGATATTGAAACCATCGa    >  1:2399708/1‑69 (MQ=18)
                                                  gcTGGTGAGCAATCCGCGCGTATTGCTGCTCGACGAACCGACTAACCACCTGGATATTGAAACCATCGa    >  1:636333/1‑69 (MQ=18)
                                                  gcTGGTGAGCAATCCGCGCGTATTGCTGCTCGACGAACCGACTAACCACCTGGATATTGAAACCATCGa    >  1:563901/1‑69 (MQ=18)
                                                  gcTGGTGAGCAATCCGCGCGTATTGCTGCTCGACGAACCGACTAACCACCTGGATATTGAAACCATCGa    >  1:556943/1‑69 (MQ=18)
                                                  gcTGGTGAGCAATCCGCGCGTATTGCTGCTCGACGAACCGACTAACCACCTGGATATTGAAACCATCGa    >  1:10772/1‑69 (MQ=18)
                                                  gcTGGTGAGCAATCCGCGCGTATTGCTGCTCGACGAACCGACTAACCACCTGGATATTGAAACCATCGa    >  1:1207253/1‑69 (MQ=18)
                                                  gcTGGTGAGCAATCCGCGCGTATTGCTGCTCGACGAACCGACTAACCACCTGGATATTGAAACCATCGa    >  1:2220447/1‑69 (MQ=18)
                                                  gcTGGTGAGCAATCCGCGCGTATTGCTGCTCGACGAACCGACTAACCACCTGGATATTGAAACCATCGa    >  1:3380215/1‑69 (MQ=18)
                                                  gcTGGTGAGCAATCCGCGCGTATTGCTGCTCGACGAACCGACTAACCACCTGGATATTGAAACCATCGa    >  1:3232348/1‑69 (MQ=18)
                                                  gcTGGTGAGCAATCCGCGCGTATTGCTGCTCGACGAACCGACTAACCACCTGGATATTGAAACCATCGa    >  1:1697239/1‑69 (MQ=18)
                                                  gcTGGTGAGCAATCCGCGCGTATTGCTGCTCGACGAACCGACTAACCACCTGGATATTGAAACCATCGa    >  1:1814099/1‑69 (MQ=18)
                                                  gcTGGTGAGCAATCCGCGCGTATTGCTGCTCGACGAACCGACTAACCACCTGGATATTGAAACCATCGa    >  1:2659817/1‑69 (MQ=18)
                                                  gcTGGTGAGCAATCCGCGCGTATTGCTGCTCGACGAACCGACTAACCACCTGGATATTGAAACCATCGa    >  1:2935919/1‑69 (MQ=18)
                                                    tGGTGAGCAATCCGCGCGTATTGCTGCTCGACGAACCGACTAACCACCTGGATATTGAAACCATCGAct  <  1:923040/69‑1 (MQ=18)
                                                    tGGTGAGCAATCCGCGCGTATTGCTGCTCGACGAACCGACTAACCACCTGGATATTGAAACCATCGAct  <  1:825662/69‑1 (MQ=18)
                                                    tGGTGAGCAATCCGCGCGTATTGCTGCTCGACGAACCGACTAACCACCTGGATATTGAAACCATCGAct  <  1:2655021/69‑1 (MQ=18)
                                                    tGGTGAGCAATCCGCGCGTATTGCTGCTCGACGAACCGACTAACCACCTGGATATTGAAACCATCGAct  <  1:734046/69‑1 (MQ=18)
                                                    tGGTGAGCAATCCGCGCGTATTGCTGCTCGACGAACCGACTAACCACCTGGATATTGAAACCATCGAct  <  1:3170519/69‑1 (MQ=18)
                                                    tGGTGAGCAATCCGCGCGTATTGCTGCTCGACGAACCGACTAACCACCTGGATATTGAAACCATCGAct  <  1:3018352/69‑1 (MQ=18)
                                                    tGGTGAGCAATCCGCGCGTATTGCTGCTCGACGAACCGACTAACCACCTGGATATTGAAACCATCGAct  <  1:1912247/69‑1 (MQ=18)
                                                    tGGTGAGCAATCCGCGCGTATTGCTGCTCGACGAACCGACTAACCACCTGGATATTGAAACCATCGAct  <  1:1065511/69‑1 (MQ=18)
                                                       tGAGCAATCCGCGCGTTTTGCTGCTCGACGAACCGACTAACCACCTGGATATTGAAACCATCGAct  <  1:2833382/66‑1 (MQ=18)
                                                           |                                                             
CTGTCGTCGCTTTCCGGCGGCTGGTTGCGTAAAGCGGCATTAGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAACCACCTGGATATTGAAACCATCGACT  >  minE/667122‑667242

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: