Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,586,077 C→A 9.7% L400L (CTG→CTT mdtF ← multidrug transporter, RpoS‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,586,0770CA 9.7% 359.5 / 13.9 113L400L (CTG→CTTmdtFmultidrug transporter, RpoS‑dependent
Reads supporting (aligned to +/- strand):  ref base C (41/61);  new base A (6/5);  total (47/66)
Fisher's exact test for biased strand distribution p-value = 5.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 115 total aligned reads displayed.

ATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGACAAGATCGCAAACGTCCCG  >  minE/2586008‑2586146
                                                                     |                                                                     
aTCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGcc                                                                       <  1:3554745/70‑1 (MQ=255)
 tCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCa                                                                      <  1:2884450/70‑1 (MQ=255)
   ttCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATACACCAGTAACCCTATCGCCAGCa                                                                   >  1:3488866/1‑71 (MQ=255)
    tCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCAc                                                                  >  1:521560/1‑71 (MQ=255)
     ccGCAATGACACTCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCa                                                                   <  1:1871226/69‑1 (MQ=255)
     ccGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCa                                                                   <  1:3080823/69‑1 (MQ=255)
      cGCAATTACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCAcc                                                                 <  1:2019901/70‑1 (MQ=255)
      cGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCa                                                                <  1:3776617/71‑1 (MQ=255)
         aaTGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATcc                                                             >  1:253318/1‑71 (MQ=255)
           tGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATccc                                                            >  1:416224/1‑70 (MQ=255)
            gACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCAc                                                                  <  1:786780/63‑1 (MQ=255)
             acgcGCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGaa                                                         >  1:1284789/4‑71 (MQ=2)
             acgcGCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGaa                                                         >  1:3790593/4‑71 (MQ=2)
             acgcGCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGaa                                                         >  1:1833409/4‑71 (MQ=2)
                cGCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAACa                                                       >  1:1706856/1‑70 (MQ=2)
                cGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACAt                                                      <  1:3795494/71‑1 (MQ=255)
                 gCTCGACGTTCTCCACCACCGCGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCg                                                           >  1:3353477/1‑65 (MQ=255)
                 gCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCa                                                                   <  1:148099/57‑1 (MQ=255)
                 gCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCa                                                                   <  1:3711481/57‑1 (MQ=255)
                 gCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCa                                                                   <  1:1200652/57‑1 (MQ=255)
                 gCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACa                                                       >  1:2023796/1‑69 (MQ=255)
                 gCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACa                                                       >  1:944151/1‑69 (MQ=255)
                 gCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACa                                                       >  1:1872219/1‑69 (MQ=255)
                   tCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGt                                                   >  1:187113/1‑71 (MQ=255)
                    cGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAACATCGt                                                   >  1:3416428/1‑70 (MQ=2)
                    cGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAAc                                                        <  1:928251/65‑1 (MQ=255)
                     gACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAACATCGtta                                                 <  1:85629/71‑3 (MQ=2)
                     gACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAACATCGtta                                                 <  1:2084177/71‑3 (MQ=2)
                     gACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGt                                                   <  1:2651753/69‑1 (MQ=255)
                     gACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCa                                                 <  1:2225889/71‑1 (MQ=255)
                      aCGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAg                                                    >  1:2130414/1‑67 (MQ=255)
                       cGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAAc                                                        >  1:2417964/1‑62 (MQ=255)
                       cGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAAc                                               >  1:417058/1‑71 (MQ=255)
                        gTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACg                                              >  1:1140588/1‑71 (MQ=255)
                        gTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACg                                              <  1:2441673/71‑1 (MQ=255)
                        gTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACg                                              <  1:2033255/71‑1 (MQ=255)
                        gTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACg                                              <  1:2210775/71‑1 (MQ=255)
                         ttCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAAc                                               <  1:1396246/69‑1 (MQ=255)
                         ttCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACgt                                             >  1:1738287/1‑71 (MQ=255)
                          tctcCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAg                                                                     <  1:3877786/46‑1 (MQ=255)
                          tctcCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACgtg                                            >  1:142341/1‑71 (MQ=255)
                          tctcCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACgt                                             <  1:526951/70‑1 (MQ=255)
                               accaccacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACAt                                                      <  1:1956328/56‑1 (MQ=255)
                               accaccacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGAt                                       >  1:3481322/1‑71 (MQ=255)
                               accaccacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGAt                                       <  1:916688/71‑1 (MQ=255)
                               accaccacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGAt                                       >  1:3757272/1‑71 (MQ=255)
                               accaccacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGAt                                       >  1:2685255/1‑71 (MQ=255)
                                 caccacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACgt                                             <  1:1421996/63‑1 (MQ=255)
                                    cacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGa                                  <  1:2792971/71‑1 (MQ=255)
                                     acGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGaa                                 <  1:2680586/71‑1 (MQ=255)
                                     acGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGaa                                 <  1:1450063/71‑1 (MQ=255)
                                     acGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGaa                                 <  1:40117/71‑1 (MQ=255)
                                     acGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGaa                                 <  1:3900925/71‑1 (MQ=255)
                                      cGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGaaa                                <  1:2915695/71‑1 (MQ=255)
                                      cGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGaaa                                <  1:484202/71‑1 (MQ=255)
                                       gATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTc                                                  >  1:124831/1‑52 (MQ=255)
                                        aTGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAacc                              <  1:3099023/71‑1 (MQ=255)
                                        aTGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAacc                              <  1:3423260/71‑1 (MQ=255)
                                         tGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAaccg                             <  1:3864657/71‑1 (MQ=255)
                                           gCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAaccga                            <  1:1753339/70‑1 (MQ=255)
                                             gTCATCCACCAGTAACCCTATCGCCAGCACCATCTCGAACATAGTCAAc                                               <  1:3740599/49‑1 (MQ=255)
                                              tCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATg                                      <  1:2577487/57‑1 (MQ=255)
                                              tCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGAccgc                        <  1:3925589/71‑1 (MQ=255)
                                               cATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACgt                                             >  1:3763548/1‑49 (MQ=255)
                                               cATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGAccgcc                       >  1:2434099/1‑71 (MQ=255)
                                                aTCCACCAGTAACCCTATCGCCAGCACCATCCCGAACAt                                                      >  1:178616/1‑39 (MQ=255)
                                                aTCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAac                               >  1:212543/1‑62 (MQ=255)
                                                aTCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGAccgcc                       <  1:1155688/70‑1 (MQ=255)
                                                aTCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGAccgcc                       >  1:1610756/1‑70 (MQ=255)
                                                aTCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGAccgcc                       >  1:1729638/1‑70 (MQ=255)
                                                 tCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGa                     >  1:3519300/1‑71 (MQ=255)
                                                   caccaGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGAc                    <  1:923957/70‑1 (MQ=255)
                                                    accaGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGAccgc                        <  1:625019/65‑1 (MQ=255)
                                                    accaGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGACaa                  <  1:673907/71‑1 (MQ=255)
                                                     ccaGTAACCCTATCGCCAGCACCATCCCGAACATAGTc                                                  <  1:2271318/38‑1 (MQ=255)
                                                     ccaGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGAc                    >  1:2024760/1‑68 (MQ=255)
                                                     ccaGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGACAAg                 <  1:694486/71‑1 (MQ=255)
                                                      caGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGa                                        <  1:738470/47‑1 (MQ=255)
                                                      caGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGaaa                                >  1:522119/1‑55 (MQ=255)
                                                      caGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGaaa                                >  1:417796/1‑55 (MQ=255)
                                                      caGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGACAAGa                <  1:2857097/71‑1 (MQ=255)
                                                      caGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGACAAGa                <  1:3929530/71‑1 (MQ=255)
                                                      caGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGACAAGa                <  1:852538/71‑1 (MQ=255)
                                                       aGTAACCCTATCGCCAGCACCATCCCGAACATAGTCaa                                                >  1:2509047/1‑38 (MQ=255)
                                                       aGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTg                                         >  1:146125/1‑45 (MQ=255)
                                                       aGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTg                                         >  1:664354/1‑45 (MQ=255)
                                                        gTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGACAAGAt               <  1:690295/70‑1 (MQ=255)
                                                           aCCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGACAAGATCGc            <  1:356836/70‑1 (MQ=255)
                                                           aCCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGACAAGATCGCa           <  1:3146710/71‑1 (MQ=255)
                                                             ccTATCGCCTGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGACAAGATCGCaa          <  1:2559852/70‑1 (MQ=255)
                                                               tATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGACAAGATCGCAAACg       >  1:1318356/1‑71 (MQ=255)
                                                                aTCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGACAAGATCGCAAACGt      >  1:1824727/1‑71 (MQ=255)
                                                                 tCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATTGTGAAACCGACCGCCGACAAGATCGCAAACGt      <  1:1570086/70‑1 (MQ=255)
                                                                 tCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGaa                                 >  1:3305288/1‑43 (MQ=255)
                                                                 tCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGaa                                 >  1:1935721/1‑43 (MQ=255)
                                                                 tCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGaa                                 >  1:2267649/1‑43 (MQ=255)
                                                                 tCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGa                                  <  1:1491114/42‑1 (MQ=255)
                                                                 tCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGACAAGATCGCAAACGt      <  1:2014738/70‑1 (MQ=255)
                                                                 tCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGACAAGATCGCAAACGTc     >  1:1436613/1‑71 (MQ=255)
                                                                    ccAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGACAAGATCGCAAACGTCCCg  <  1:3149276/71‑1 (MQ=255)
                                                                     |                                                                     
ATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGACAAGATCGCAAACGTCCCG  >  minE/2586008‑2586146

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: