Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 297,667 C→T 15.8% R110C (CGT→TGT)  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
*minE297,6670CT15.8% 144.6 / 17.9 57R110C (CGT→TGT) mdlAfused predicted multidrug transporter subunits and ATP‑binding components of ABC superfamily
Reads supporting (aligned to +/- strand):  ref base C (19/29);  new base T (6/3);  total (25/32)
Fisher's exact test for biased strand distribution p-value = 1.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

AACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCC  >  minE/297599‑297732
                                                                    |                                                                 
aaCTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGt                                                                 <  1:560159/71‑1 (MQ=255)
 aCTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGTGTc                                                                <  1:2196591/71‑1 (MQ=255)
 aCTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTc                                                                <  1:486918/71‑1 (MQ=255)
 aCTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTc                                                                <  1:3329687/71‑1 (MQ=255)
  cTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCa                                                               <  1:1325912/71‑1 (MQ=255)
  cTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCa                                                               <  1:3862720/71‑1 (MQ=255)
  cTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCa                                                               <  1:146100/71‑1 (MQ=255)
  cTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCa                                                               <  1:2584838/71‑1 (MQ=255)
  cTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCa                                                               <  1:1957394/71‑1 (MQ=255)
     gCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCAt                                                              <  1:2527946/69‑1 (MQ=255)
      cTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTc                                                                <  1:1013974/66‑1 (MQ=255)
      cTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGc                                                           <  1:2530256/71‑1 (MQ=255)
        gTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCAc                                                         >  1:1094273/1‑71 (MQ=255)
        gTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCAc                                                         >  1:901631/1‑71 (MQ=255)
        gTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCAc                                                         >  1:2286773/1‑71 (MQ=255)
        gTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCAc                                                         >  1:193534/1‑71 (MQ=255)
         ttGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCAcc                                                        <  1:1926936/71‑1 (MQ=255)
         ttGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCAcc                                                        <  1:3070388/71‑1 (MQ=255)
         ttGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCAcc                                                        <  1:691202/71‑1 (MQ=255)
         ttGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCAcc                                                        <  1:2844772/71‑1 (MQ=255)
          tGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGTGTCATCGCACCg                                                       >  1:3845195/1‑71 (MQ=255)
          tGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGTGTCATCGCACCg                                                       >  1:3519494/1‑71 (MQ=255)
          tGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGTGTCATCGCACCg                                                       >  1:2532622/1‑71 (MQ=255)
           gAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCgg                                                      >  1:3473139/1‑71 (MQ=255)
            aaCTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGt                                                     <  1:331015/71‑1 (MQ=255)
            aaCTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGt                                                     <  1:551212/71‑1 (MQ=255)
            aaCTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGt                                                     <  1:2905320/71‑1 (MQ=255)
            aaCTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGt                                                     <  1:1516770/71‑1 (MQ=255)
             aCTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGt                                                     <  1:2777866/70‑1 (MQ=255)
              cTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGTGTCATCGCACCGGTg                                                    <  1:3062186/70‑1 (MQ=255)
                     aaGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATg                                            >  1:3695406/1‑71 (MQ=255)
                                    aGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAAt                             <  1:850028/71‑1 (MQ=255)
                                    aGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAAt                             <  1:1951483/71‑1 (MQ=255)
                                      cTAAGCCGGCAGCATCCTGAGTTTTACCTGTGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGa                           <  1:3539429/71‑1 (MQ=255)
                                          gCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAAt                             <  1:3633436/65‑1 (MQ=255)
                                              gcagcaTCCTGAGTTTTACCTGTGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGAt                    >  1:1533423/1‑70 (MQ=255)
                                              gcagcaTCCTGAGTTTTACCTGTGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATc                   >  1:1473785/1‑71 (MQ=255)
                                                   atacTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGcgtc              >  1:1677628/4‑71 (MQ=255)
                                                   aTCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCg                                       >  1:2821918/1‑46 (MQ=255)
                                                   aTCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCg                                       >  1:2213107/1‑46 (MQ=255)
                                                   aTCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCg                                       >  1:1464008/1‑46 (MQ=255)
                                                   aTCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGcgtc              >  1:2358287/1‑71 (MQ=255)
                                                   aTCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGcgtc              >  1:1842132/1‑71 (MQ=255)
                                                    tCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGAc                                 <  1:1447252/51‑1 (MQ=255)
                                                    tCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGAc                                 <  1:1193009/51‑1 (MQ=255)
                                                    tCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGAc                                 <  1:227928/51‑1 (MQ=255)
                                                      cTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATg                                            >  1:1344157/1‑38 (MQ=255)
                                                      cTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGAt                    >  1:3699166/1‑62 (MQ=255)
                                                      cTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGAt                    >  1:539088/1‑62 (MQ=255)
                                                      cTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGAt                    >  1:1201554/1‑62 (MQ=255)
                                                      cTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGAt                    >  1:670170/1‑62 (MQ=255)
                                                      cTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCgtg           >  1:700888/1‑71 (MQ=255)
                                                      cTGAGTTTTACATGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGAt                    >  1:1637804/1‑62 (MQ=255)
                                                           ttttACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGAcc                                <  1:3563669/45‑1 (MQ=255)
                                                           ttttACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGAcc                                <  1:3256329/45‑1 (MQ=255)
                                                           ttttACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGAcc                                <  1:2512436/45‑1 (MQ=255)
                                                               aCCTGTGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTgccgcc  >  1:2017491/1‑71 (MQ=255)
                                                                    |                                                                 
AACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCC  >  minE/297599‑297732

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: