Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 298,537 C→T 45.5% L400L (CTG→TTG)  mdlA → fused predicted multidrug 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
*minE298,5370CT45.5% 22.7 / 105.9 77L400L (CTG→TTG) mdlAfused predicted multidrug transporter subunits and ATP‑binding components of ABC superfamily
Reads supporting (aligned to +/- strand):  ref base C (21/21);  new base T (15/20);  total (36/41)
Fisher's exact test for biased strand distribution p-value = 6.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.23e-01

TACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCTTTTT  >  minE/298467‑298605
                                                                      |                                                                    
tACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCctc                                                                      >  1:2324663/1‑71 (MQ=255)
tACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCctc                                                                      >  1:2868097/1‑71 (MQ=255)
 aCCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCctc                                                                      >  1:2108543/1‑70 (MQ=255)
         gtCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGt                                                             >  1:3624360/1‑71 (MQ=255)
         gtCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGt                                                             >  1:1517550/1‑71 (MQ=255)
         gtCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGt                                                             >  1:636618/1‑71 (MQ=255)
          tCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTTTGACGAAGtt                                                            <  1:3623949/71‑1 (MQ=255)
          tCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTTTGACGAAGtt                                                            <  1:3408050/71‑1 (MQ=255)
            gCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTAc                                                          <  1:1176810/71‑1 (MQ=255)
            gCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTAc                                                          <  1:2717650/71‑1 (MQ=255)
            gCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTAc                                                          <  1:80751/71‑1 (MQ=255)
             cTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTTTGACGAAGTTAc                                                          >  1:819551/1‑70 (MQ=255)
               cATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTTTGACGAAGTTACAAc                                                       >  1:1944369/1‑71 (MQ=255)
               cATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTTTGACGAAGTTACAAc                                                       >  1:3452158/1‑71 (MQ=255)
                 ttCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTTTGACGAAGTTACAACTc                                                     <  1:3249827/71‑1 (MQ=255)
                 ttCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTTTGACGAAGTTACAACTc                                                     <  1:695947/71‑1 (MQ=255)
                 ttCAGCGTCATTTCGACGTCAGCGAGGGGCATATTCGCTTTCATGATATTCCTTTGACGAAGTTACAACTc                                                     <  1:3741520/71‑1 (MQ=255)
                  tCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTTTGACGAAGTTACAACTc                                                     >  1:1933893/1‑70 (MQ=255)
                        tCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTc                                                     <  1:857225/64‑1 (MQ=255)
                           tttCGACGTCAGCTAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTgg                                            >  1:478479/1‑70 (MQ=255)
                           tttCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGAAAGCTGGc                                           >  1:4045314/1‑71 (MQ=255)
                             tCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTTTGACGAAGTTACAACTCGATAGCTGGCg                                          <  1:2424737/70‑1 (MQ=255)
                             tCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACaa                                                        >  1:3610518/1‑56 (MQ=255)
                             tCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACaa                                                        >  1:2659218/1‑56 (MQ=255)
                                  gTCAGCGAGGGGGATATTCGCTTTCATGATATTCctct                                                                     <  1:887381/38‑1 (MQ=255)
                                  gTCAGCGAGGGGGATATTCGCTTTCATGATATTCctct                                                                     <  1:859594/38‑1 (MQ=255)
                                     aGCGAGGGGGATATTCGCTTTCATGATATTCCTTTGACGAAGTTACaa                                                        >  1:1818167/1‑48 (MQ=255)
                                     aGCGAGGGGGATATTCGCTTTCATGATATTCCTTTGACGAAGTTACaa                                                        >  1:2491219/1‑48 (MQ=255)
                                     aGCGAGGGGGATATTCGCTTTCATGATATTCCTTTGACGAAGTTACaa                                                        >  1:3687044/1‑48 (MQ=255)
                                     aGCGAGGGGGATATTCGCTTTCATGATATTCCTTTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCt                                 >  1:1338595/1‑71 (MQ=255)
                                     aGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCt                                 >  1:4226938/1‑71 (MQ=255)
                                      gCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGc                                               <  1:2502673/56‑1 (MQ=255)
                                      gCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGc                                               <  1:2148825/56‑1 (MQ=255)
                                      gCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGc                                               <  1:1523333/56‑1 (MQ=255)
                                       cGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGc                                               <  1:624574/55‑1 (MQ=255)
                                          gggggATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAAc                                                       <  1:2215726/44‑1 (MQ=255)
                                           ggggATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGt                           <  1:3829800/71‑1 (MQ=255)
                                            gggATATTCGCTTTCATGATATTCCTCTGACGAAGTTACa                                                         <  1:359254/40‑1 (MQ=255)
                                            gggATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTc                                                     >  1:2901733/1‑44 (MQ=255)
                                            gggATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTc                                                     >  1:554952/1‑44 (MQ=255)
                                             ggATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAg                         <  1:148108/71‑1 (MQ=255)
                                             ggATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAg                         <  1:1299568/71‑1 (MQ=255)
                                              gATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAg                         <  1:695033/70‑1 (MQ=255)
                                               atatTCGCTTTCATGATATTCCTTTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGtt                       <  1:1811960/71‑1 (MQ=255)
                                                tatTCGCTTTCATGATATTCCTTTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTa                      <  1:1309917/71‑1 (MQ=255)
                                                tatTCGCTTTCATGATATTCCTTTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTa                      <  1:2443287/71‑1 (MQ=255)
                                                tatTCGCTTTCATGATATTCCTCTGACGAAGTTACAAc                                                       <  1:1831392/38‑1 (MQ=255)
                                                tatTCGCTTTCATGATATTCCTCTGACGAAGTTACAAc                                                       <  1:666894/38‑1 (MQ=255)
                                                tatTCGCTTTCATGATATTCCTCTGACGAAGTTACAAc                                                       <  1:2690467/38‑1 (MQ=255)
                                                tatTCGCTTTCATGATATTCCTCTGACGAAGTTACAAc                                                       <  1:397957/38‑1 (MQ=255)
                                                tatTCGCTTTCATGATATTCCTCTGACGAAGTTACAAc                                                       <  1:3940520/38‑1 (MQ=255)
                                                  ttCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAgc                                      >  1:141609/1‑53 (MQ=255)
                                                  ttCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAgc                                      >  1:182529/1‑53 (MQ=255)
                                                    cGCTTTCATGATATTCCTTTGACGAAGTTACAACTCGa                                                   >  1:2533914/1‑38 (MQ=255)
                                                    cGCTTTCATGATATTCCTTTGACGAAGTTACAACTCGa                                                   >  1:376693/1‑38 (MQ=255)
                                                    cGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAg                         >  1:1132067/1‑64 (MQ=255)
                                                     gCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCa                  <  1:3462589/70‑1 (MQ=255)
                                                       tttCATGATATTCCTCTGATGAAGTTACAACTCGAt                                                  >  1:4224403/1‑36 (MQ=255)
                                                       tttCATGATATTCCTCTGACGAAGTTACAACTCGAt                                                  >  1:1755693/1‑36 (MQ=255)
                                                       tttCATGATATTCCTCTGACGAAGTTACAACTCGAt                                                  >  1:3177592/1‑36 (MQ=255)
                                                       tttCATGATATTCCTCTGACGAAGTTACAACTCGAt                                                  >  1:810007/1‑36 (MQ=255)
                                                       tttCATGATATTCCTCTGACGAAGTTACAACTCGAt                                                  >  1:3547859/1‑36 (MQ=255)
                                                          cATGATATTCCTTTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGcc                   >  1:3381213/1‑64 (MQ=255)
                                                          cATGATATTCCTTTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGcc                   >  1:3376395/1‑64 (MQ=255)
                                                          cATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGc             <  1:3040721/70‑1 (MQ=255)
                                                           aTGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCa           >  1:704990/1‑71 (MQ=255)
                                                           aTGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCa           >  1:2896178/1‑71 (MQ=255)
                                                                 ttCCTTTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCt      <  1:1375103/70‑1 (MQ=255)
                                                                 ttCCTTTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCt      <  1:778507/70‑1 (MQ=255)
                                                                 ttCCTTTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCt      <  1:2055432/70‑1 (MQ=255)
                                                                 ttCCTTTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCt      <  1:3065374/70‑1 (MQ=255)
                                                                 ttCCTTTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCt      <  1:4192618/70‑1 (MQ=255)
                                                                 ttCCTTTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCt      <  1:1754896/70‑1 (MQ=255)
                                                                 ttCCTTTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCt      <  1:3607692/70‑1 (MQ=255)
                                                                 ttCCTTTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCt      <  1:2590816/70‑1 (MQ=255)
                                                                  tCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCAt          >  1:2676071/1‑65 (MQ=255)
                                                                   ccTTTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCtttt   >  1:2546655/1‑71 (MQ=255)
                                                                   ccTTTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCtttt   >  1:1644854/1‑71 (MQ=255)
                                                                   ccTTTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCtttt   >  1:3966739/1‑71 (MQ=255)
                                                                   ccTTTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCttt    <  1:4031277/70‑1 (MQ=255)
                                                                   ccTTTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCttt    <  1:2141396/70‑1 (MQ=255)
                                                                   ccTTTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCttt    <  1:2152734/70‑1 (MQ=255)
                                                                    ctctGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCttttt  <  1:3092244/71‑1 (MQ=255)
                                                                    ctctGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCttttt  <  1:2209090/71‑1 (MQ=255)
                                                                    ctctGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCtttt   <  1:1668010/70‑1 (MQ=255)
                                                                    ctctGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCtttt   <  1:2779718/70‑1 (MQ=255)
                                                                    ctctGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCtttt   <  1:217341/70‑1 (MQ=255)
                                                                     tctGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCtttt   <  1:828361/69‑1 (MQ=255)
                                                                      |                                                                    
TACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCTTTTT  >  minE/298467‑298605

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: