Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 667,205 T→C 31.9% D181D (GAT→GAC 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,2050TC31.9% 65.7 / 35.6 72D181D (GAT→GACuupfused predicted transporter subunits and ATP‑binding components of ABC superfamily
Reads supporting (aligned to +/- strand):  ref base T (13/36);  new base C (14/9);  total (27/45)
Fisher's exact test for biased strand distribution p-value = 8.45e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.34e-01

TAAAGCGGCATTAGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAACCACCTGGATATTGAAACCATCGACTGGCTGGAAGGGTTTTTGAA  >  minE/667151‑667261
                                                      |                                                        
tAAAGCGGCATTAGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAacca                                            <  1:2140134/69‑1 (MQ=255)
tAAAGCGGCATTAGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAacca                                            <  1:213087/69‑1 (MQ=255)
tAAAGCGGCATTAGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAacca                                            <  1:2383042/69‑1 (MQ=255)
tAAAGCGGCATTAGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAacca                                            <  1:148241/69‑1 (MQ=255)
tAAAGCGGCATTAGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAacca                                            <  1:754027/69‑1 (MQ=255)
tAAAGCGGCATTAGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAacca                                            <  1:93299/69‑1 (MQ=255)
      ggCATTAGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAACCACCTGGa                                      <  1:2278230/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:2218883/69‑1 (MQ=255)
      ggCATTAGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAACCACCTGGa                                      <  1:2910455/69‑1 (MQ=255)
      ggCATTAGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAACCACCTGGa                                      <  1:1901771/69‑1 (MQ=255)
      ggCATTAGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAACCACCTGGa                                      <  1:3066584/69‑1 (MQ=255)
      ggCATTAGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAACCACCTGGa                                      <  1:539736/69‑1 (MQ=255)
      ggCATTAGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAACCACCTGGa                                      <  1:1612659/69‑1 (MQ=255)
      ggCATTAGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAACCACCTGGa                                      <  1:1587939/69‑1 (MQ=255)
      ggCATTAGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAACCACCTGGa                                      <  1:1454586/69‑1 (MQ=255)
      ggCATTAGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAACCACCTGGa                                      <  1:432259/69‑1 (MQ=255)
            aGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAACCACCTGGATATTGa                                >  1:278094/1‑69 (MQ=255)
            aGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAACCACCTGGATATTGa                                >  1:2055440/1‑69 (MQ=255)
            aGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAACCACCTGGATATTGa                                >  1:2177903/1‑69 (MQ=255)
            aGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAACCACCTGGATATTGa                                >  1:2422555/1‑69 (MQ=255)
            aGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAACCACCTGGATATTGa                                >  1:2468840/1‑69 (MQ=255)
            aGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAACCACCTGGATATTGa                                >  1:2727131/1‑69 (MQ=255)
            aGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAACCACCTGGATATTGa                                >  1:1337583/1‑69 (MQ=255)
            aGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAACCACCTGGATATTGa                                >  1:3176883/1‑69 (MQ=255)
            aGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAACCACCTGGATATTGa                                >  1:479633/1‑69 (MQ=255)
                cgcgcgCTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAACCACCTGGATATTGAAAcc                            >  1:3154142/1‑69 (MQ=255)
                   gcgcTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAACCACCTGGATATTGAAACCATc                         <  1:1987580/69‑1 (MQ=255)
                   gcgcTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAACCACCTGGATATTGAAACCATc                         <  1:298459/69‑1 (MQ=255)
                   gcgcTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAACCACCTGGATATTGAAACCATc                         <  1:1853036/69‑1 (MQ=255)
                   gcgcTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAACCACCTGGATATTGAAACCATc                         <  1:700485/69‑1 (MQ=255)
                   gcgcTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAACCACCTGGATATTGAAACCATc                         <  1:961651/69‑1 (MQ=255)
                   gcgcTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAACCACCTGGATATTGAAACCATc                         <  1:696481/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:636333/1‑69 (MQ=18)
                     gcTGGTGAGCAATCCGCGCGTATTGCTGCTCGACGAACCGACTAACCACCTGGATATTGAAACCATCGa                       >  1:2399708/1‑69 (MQ=18)
                     gcTGGTGAGCAATCCGCGCGTATTGCTGCTCGACGAACCGACTAACCACCTGGATATTGAAACCATCGa                       >  1:2659817/1‑69 (MQ=18)
                     gcTGGTGAGCAATCCGCGCGTATTGCTGCTCGACGAACCGACTAACCACCTGGATATTGAAACCATCGa                       >  1:563901/1‑69 (MQ=18)
                     gcTGGTGAGCAATCCGCGCGTATTGCTGCTCGACGAACCGACTAACCACCTGGATATTGAAACCATCGa                       >  1:2220447/1‑69 (MQ=18)
                     gcTGGTGAGCAATCCGCGCGTATTGCTGCTCGACGAACCGACTAACCACCTGGATATTGAAACCATCGa                       >  1:2935919/1‑69 (MQ=18)
                     gcTGGTGAGCAATCCGCGCGTATTGCTGCTCGACGAACCGACTAACCACCTGGATATTGAAACCATCGa                       >  1:10772/1‑69 (MQ=18)
                     gcTGGTGAGCAATCCGCGCGTATTGCTGCTCGACGAACCGACTAACCACCTGGATATTGAAACCATCGa                       >  1:556943/1‑69 (MQ=18)
                     gcTGGTGAGCAATCCGCGCGTATTGCTGCTCGACGAACCGACTAACCACCTGGATATTGAAACCATCGa                       >  1:420757/1‑69 (MQ=18)
                     gcTGGTGAGCAATCCGCGCGTATTGCTGCTCGACGAACCGACTAACCACCTGGATATTGAAACCATCGa                       >  1:1814099/1‑69 (MQ=18)
                     gcTGGTGAGCAATCCGCGCGTATTGCTGCTCGACGAACCGACTAACCACCTGGATATTGAAACCATCGa                       >  1:1697239/1‑69 (MQ=18)
                     gcTGGTGAGCAATCCGCGCGTATTGCTGCTCGACGAACCGACTAACCACCTGGATATTGAAACCATCGa                       >  1:1207253/1‑69 (MQ=18)
                     gcTGGTGAGCAATCCGCGCGTATTGCTGCTCGACGAACCGACTAACCACCTGGATATTGAAACCATCGa                       >  1:3232348/1‑69 (MQ=18)
                     gcTGGTGAGCAATCCGCGCGTATTGCTGCTCGACGAACCGACTAACCACCTGGATATTGAAACCATCGa                       >  1:3380215/1‑69 (MQ=18)
                       tGGTGAGCAATCCGCGCGTATTGCTGCTCGACGAACCGACTAACCACCTGGATATTGAAACCATCGAct                     <  1:1065511/69‑1 (MQ=18)
                       tGGTGAGCAATCCGCGCGTATTGCTGCTCGACGAACCGACTAACCACCTGGATATTGAAACCATCGAct                     <  1:923040/69‑1 (MQ=18)
                       tGGTGAGCAATCCGCGCGTATTGCTGCTCGACGAACCGACTAACCACCTGGATATTGAAACCATCGAct                     <  1:825662/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:2655021/69‑1 (MQ=18)
                       tGGTGAGCAATCCGCGCGTATTGCTGCTCGACGAACCGACTAACCACCTGGATATTGAAACCATCGAct                     <  1:1912247/69‑1 (MQ=18)
                          tGAGCAATCCGCGCGTTTTGCTGCTCGACGAACCGACTAACCACCTGGATATTGAAACCATCGAct                     <  1:2833382/66‑1 (MQ=18)
                                   cgcgcgTGCTGTTGCTTGATGAACCGACAAACCACCTGGATATTGAAACCATCGACTGGCTGGAAGGGt         >  1:2012608/1‑69 (MQ=255)
                                   cgcgcgTGCTGTTGCTTGATGAACCGACAAACCACCTGGATATTGAAACCATCGACTGGCTGGAAGGGt         >  1:1469273/1‑69 (MQ=255)
                                   cgcgcgTGCTGTTGCTTGATGAACCGACAAACCACCTGGATATTGAAACCATCGACTGGCTGGAAGGGt         >  1:965127/1‑69 (MQ=255)
                                          gCTGTTGCTTGATGAACCGACAAACCACCTTGATATTGAAACCATCGACTGGCTGGAAGGGTTTTTGaa  <  1:477533/69‑1 (MQ=255)
                                          gCTGTTGCTTGATGAACCGACAAACCACCTGGATATTGAAACCATCGACTGGCTGGAAGGGTTTTTGaa  <  1:2397067/69‑1 (MQ=255)
                                          gCTGTTGCTTGATGAACCGACAAACCACCTGGATATTGAAACCATCGACTGGCTGGAAGGGTTTTTGaa  <  1:2258060/69‑1 (MQ=255)
                                          gCTGTTGCTTGATGAACCGACAAACCACCTGGATATTGAAACCATCGACTGGCTGGAAGGGTTTTTGaa  <  1:1858905/69‑1 (MQ=255)
                                          gCTGTTGCTTGATGAACCGACAAACCACCTGGATATTGAAACCATCGACTGGCTGGAAGGGTTTTTGaa  <  1:1609516/69‑1 (MQ=255)
                                          gCTGTTGCTTGATGAACCGACAAACCACCTGGATATTGAAACCAACGACTGGCTGGAAGGGTTTTTGaa  <  1:2330886/69‑1 (MQ=255)
                                          gCTGTCGCTTGATGAACCGACAAACCACCTGGATATTGAAACCATCGACTGGCTGGAAGGGTTTTTGaa  <  1:872517/69‑1 (MQ=255)
                                           cTGTTGCTTGATGAACCGACAAACCACCTGGATATTGAAACCATCGACTGGCTGGAAGGGTTTTTGaa  <  1:97191/68‑1 (MQ=255)
                                                      |                                                        
TAAAGCGGCATTAGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAACCACCTGGATATTGAAACCATCGACTGGCTGGAAGGGTTTTTGAA  >  minE/667151‑667261

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