Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,586,071 T→G 10.2% I402I (ATA→ATC mdtF ← multidrug transporter, RpoS‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,586,0710TG10.2% 297.3 / 28.3 118I402I (ATA→ATCmdtFmultidrug transporter, RpoS‑dependent
Reads supporting (aligned to +/- strand):  ref base T (42/64);  new base G (7/5);  total (49/69)
Fisher's exact test for biased strand distribution p-value = 2.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.56e-01

Only 100 of 119 total aligned reads displayed.

TAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGACAAGATCGCAAACG  >  minE/2586002‑2586141
                                                                     |                                                                      
tAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTa                                                                       <  1:290303/71‑1 (MQ=255)
tAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTa                                                                       <  1:3595756/71‑1 (MQ=255)
 aGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTAt                                                                      >  1:3489502/1‑71 (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)
      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)
           ccGCAATGACACTCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCa                                                              <  1:1871226/69‑1 (MQ=255)
           ccGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCa                                                              <  1:3080823/69‑1 (MQ=255)
            cGCAATTACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCAcc                                                            <  1:2019901/70‑1 (MQ=255)
            cGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCAcc                                                            <  1:59734/70‑1 (MQ=255)
            cGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCa                                                           <  1:3776617/71‑1 (MQ=255)
             gCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTAt                                                                      >  1:179953/1‑59 (MQ=255)
               aaTGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTa                                                                       <  1:1279990/56‑1 (MQ=255)
               aaTGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATcc                                                        >  1:253318/1‑71 (MQ=255)
               aaTGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATcc                                                        >  1:3364156/1‑71 (MQ=255)
                 tGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATc                                                                     >  1:3748009/1‑56 (MQ=255)
                 tGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATc                                                                     >  1:1467828/1‑56 (MQ=255)
                 tGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATccc                                                       >  1:416224/1‑70 (MQ=255)
                  gACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCAc                                                             <  1:786780/63‑1 (MQ=255)
                   acgcGCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGaa                                                    >  1:3790593/4‑71 (MQ=2)
                       gCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCa                                                              <  1:1200652/57‑1 (MQ=255)
                       gCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCa                                                              <  1:148099/57‑1 (MQ=255)
                       gCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCa                                                              <  1:3711481/57‑1 (MQ=255)
                       gCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACa                                                  >  1:944151/1‑69 (MQ=255)
                       gCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACa                                                  >  1:2023796/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)
                          cGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAAc                                                   <  1:413724/65‑1 (MQ=255)
                           gACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAACATCGtta                                            <  1:2084177/71‑3 (MQ=2)
                           gACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAACATCGtta                                            <  1:85629/71‑3 (MQ=2)
                           gACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGt                                              <  1:2651753/69‑1 (MQ=255)
                           gACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTc                                             <  1:2405175/70‑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:2033255/71‑1 (MQ=255)
                              gTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACg                                         >  1:1140588/1‑71 (MQ=255)
                              gTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACg                                         <  1:2210775/71‑1 (MQ=255)
                              gTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACg                                         <  1:2441673/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: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)
                                     accaccacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGAt                                  >  1:3481322/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)
                                            cGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACgtgt                                      >  1:1537788/1‑60 (MQ=255)
                                            cGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGaaa                           <  1:484202/71‑1 (MQ=255)
                                            cGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGaaa                           <  1:2915695/71‑1 (MQ=255)
                                             gATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTc                                             >  1:124831/1‑52 (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)
                                                    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:1729638/1‑70 (MQ=255)
                                                      aTCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGAccgcc                  <  1:1155688/70‑1 (MQ=255)
                                                      aTCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGAccgcc                  >  1:1610756/1‑70 (MQ=255)
                                                         caccaGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGAc               <  1:923957/70‑1 (MQ=255)
                                                         caccaGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGAc               <  1:1942710/70‑1 (MQ=255)
                                                          accaGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGAccgc                   <  1:625019/65‑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:417796/1‑55 (MQ=255)
                                                            caGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGaaa                           >  1:522119/1‑55 (MQ=255)
                                                            caGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGACAAGa           <  1:1782093/71‑1 (MQ=255)
                                                            caGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGACAAGa           <  1:2857097/71‑1 (MQ=255)
                                                            caGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGACAAGa           <  1:2619898/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)
                                                                 aCCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGACAAGATCGc       <  1:356836/70‑1 (MQ=255)
                                                                 aCCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGACAAGATCGCa      <  1:3838294/71‑1 (MQ=255)
                                                                 aCCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGACAAGATCGCa      <  1:3146710/71‑1 (MQ=255)
                                                                    cTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGa                             <  1:1406112/45‑1 (MQ=255)
                                                                     tATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGACAAGATCGCAAACg  >  1:1318356/1‑71 (MQ=255)
                                                                     |                                                                      
TAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGACAAGATCGCAAACG  >  minE/2586002‑2586141

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: