Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,586,065 T→C 10.6% L404L (TTA→TTG) ‡ mdtF ← multidrug transporter, RpoS‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,586,0650TC10.6% 332.1 / 24.4 122L404L (TTA→TTG) ‡mdtFmultidrug transporter, RpoS‑dependent
Reads supporting (aligned to +/- strand):  ref base T (45/64);  new base C (8/5);  total (53/69)
Fisher's exact test for biased strand distribution p-value = 2.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 124 total aligned reads displayed.

CGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGACAAGAT  >  minE/2585996‑2586133
                                                                     |                                                                    
 ggcggTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTa                                                                     <  1:1024032/70‑1 (MQ=255)
 ggcggTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTa                                                                     <  1:1783259/70‑1 (MQ=255)
  gcggTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAAc                                                                   >  1:3736467/1‑71 (MQ=255)
     ggaGCTTATCCTCCATCATCACGCGCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGccc                                                                 >  1:1364286/3‑70 (MQ=255)
      tAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTa                                                               <  1:290303/71‑1 (MQ=255)
       aGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTAt                                                              >  1:3489502/1‑71 (MQ=255)
        gCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATc                                                             <  1:1975963/71‑1 (MQ=255)
        gCTTATCCTCCATCATCACGCGCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGCCCGAtg                                                             >  1:3295347/1‑70 (MQ=255)
          ttATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGc                                                           >  1:2878324/1‑71 (MQ=255)
           tATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGc                                                           >  1:1274685/1‑70 (MQ=255)
           tATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGc                                                           >  1:3234193/1‑70 (MQ=255)
            aTCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAAccc                                                                 <  1:1540280/63‑1 (MQ=255)
            aTCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCg                                                            <  1:2119080/68‑1 (MQ=255)
            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)
                  cGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCa                                                   <  1:3776617/71‑1 (MQ=255)
                   gCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTAt                                                              >  1:179953/1‑59 (MQ=255)
                     aaTGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTaa                                                                    >  1:2727636/1‑51 (MQ=255)
                     aaTGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTa                                                               <  1:1279990/56‑1 (MQ=255)
                     aaTGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATcc                                                >  1:3364156/1‑71 (MQ=255)
                     aaTGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATcc                                                >  1:253318/1‑71 (MQ=255)
                       tGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATc                                                             >  1:3748009/1‑56 (MQ=255)
                       tGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATccc                                               >  1:416224/1‑70 (MQ=255)
                        gACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCAc                                                     <  1:786780/63‑1 (MQ=255)
                         acgcGCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGaa                                            >  1:1833409/4‑71 (MQ=2)
                         acgcGCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGaa                                            >  1:1284789/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:3711481/57‑1 (MQ=255)
                             gCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCa                                                      <  1:148099/57‑1 (MQ=255)
                             gCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCa                                                      <  1:1200652/57‑1 (MQ=255)
                             gCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGt                                                                      <  1:1123023/41‑1 (MQ=255)
                             gCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACa                                          >  1:944151/1‑69 (MQ=255)
                             gCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACa                                          >  1:1872219/1‑69 (MQ=255)
                             gCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACa                                          >  1:2023796/1‑69 (MQ=255)
                              cTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAAccc                                                                 <  1:2353192/45‑1 (MQ=255)
                                cGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAACATCGt                                      >  1:3416428/1‑70 (MQ=2)
                                cGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAAc                                           <  1:928251/65‑1 (MQ=255)
                                cGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAAc                                           <  1:413724/65‑1 (MQ=255)
                                 gACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAACATCGtta                                    <  1:2084177/71‑3 (MQ=2)
                                 gACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGt                                      <  1:2651753/69‑1 (MQ=255)
                                 gACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTc                                     <  1:2405175/70‑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:2033255/71‑1 (MQ=255)
                                    gTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACg                                 <  1:2210775/71‑1 (MQ=255)
                                    gTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACg                                 <  1:2441673/71‑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:836523/56‑1 (MQ=255)
                                           accaccacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACAt                                         <  1:1956328/56‑1 (MQ=255)
                                           accaccacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGAt                          >  1:3757272/1‑71 (MQ=255)
                                           accaccacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGAt                          <  1:916688/71‑1 (MQ=255)
                                           accaccacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGAt                          >  1:2685255/1‑71 (MQ=255)
                                             caccacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACgt                                <  1:1421996/63‑1 (MQ=255)
                                                cacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGa                     <  1:3618446/71‑1 (MQ=255)
                                                 acGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGaa                    <  1:3900925/71‑1 (MQ=255)
                                                 acGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGaa                    <  1:40117/71‑1 (MQ=255)
                                                 acGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGaa                    <  1:2680586/71‑1 (MQ=255)
                                                 acGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGaa                    <  1:1450063/71‑1 (MQ=255)
                                                  cGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACgtgt                              >  1:1537788/1‑60 (MQ=255)
                                                  cGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGaaa                   <  1:484202/71‑1 (MQ=255)
                                                  cGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGaaa                   <  1:2915695/71‑1 (MQ=255)
                                                    aTGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAacc                 <  1:3423260/71‑1 (MQ=255)
                                                    aTGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAacc                 <  1:3099023/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)
                                                            aTCCACCAGTAACCCTATCGCCAGCACCATCCCGAACAt                                         >  1:178616/1‑39 (MQ=255)
                                                            aTCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAac                  >  1:212543/1‑62 (MQ=255)
                                                            aTCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGAccgcc          >  1:1610756/1‑70 (MQ=255)
                                                            aTCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGAccgcc          >  1:1729638/1‑70 (MQ=255)
                                                            aTCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGAccgcc          <  1:1155688/70‑1 (MQ=255)
                                                             tCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGa        >  1:3519300/1‑71 (MQ=255)
                                                               caccaGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGAc       <  1:923957/70‑1 (MQ=255)
                                                               caccaGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGAc       <  1:1942710/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)
                                                                  caGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGa                           <  1:738470/47‑1 (MQ=255)
                                                                  caGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGaaa                   >  1:522119/1‑55 (MQ=255)
                                                                  caGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGACAAGa   <  1:852538/71‑1 (MQ=255)
                                                                  caGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGACAAGa   <  1:2857097/71‑1 (MQ=255)
                                                                   aGTAACCCTATCGCCAGCACCATCCCGAACATAGTCaa                                   >  1:2509047/1‑38 (MQ=255)
                                                                   aGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTg                            >  1:664354/1‑45 (MQ=255)
                                                                   aGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTg                            >  1:146125/1‑45 (MQ=255)
                                                                    gTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGACAAGAt  <  1:690295/70‑1 (MQ=255)
                                                                     |                                                                    
CGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGACAAGAT  >  minE/2585996‑2586133

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