Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,586,059 C→G 7.6% V406V (GTG→GTC mdtF ← multidrug transporter, RpoS‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,586,0590CG 7.6% 327.3 / 11.6 105V406V (GTG→GTCmdtFmultidrug transporter, RpoS‑dependent
Reads supporting (aligned to +/- strand):  ref base C (48/49);  new base G (4/4);  total (52/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-01

Only 100 of 121 total aligned reads displayed.

TTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGA  >  minE/2585990‑2586127
                                                                     |                                                                    
ttccttcGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATcca                                                                     >  1:2393417/1‑71 (MQ=255)
  ccttcGGCGGTATCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATccac                                                                    <  1:442400/70‑1 (MQ=255)
  ccttcGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATccac                                                                    >  1:3671983/1‑70 (MQ=255)
     ttggcggGAGCTTATCCTCCATCATCACGCGCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAg                                                                 >  1:3930902/3‑70 (MQ=12)
       ggcggGAGCTTATCCTCCATCATCACGCGCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAgcag                                                              >  1:825352/1‑68 (MQ=9)
       ggcggGAGCTTATCCTCCATCATCACGCGCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAgcag                                                              >  1:1710291/1‑68 (MQ=9)
       ggcggGAGCTTATCCTCCATCATCACGCGCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAgcag                                                              >  1:1569165/1‑68 (MQ=9)
        gcggTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAAc                                                             <  1:3972463/71‑1 (MQ=255)
        gcggTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAAc                                                             <  1:3136815/71‑1 (MQ=255)
        gcggTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAAc                                                             <  1:3005287/71‑1 (MQ=255)
        gcggGAGCTTATCCTCCATCATCACGCGCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAgcag                                                              <  1:3445634/70‑4 (MQ=9)
          ggTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAAccc                                                           >  1:933768/1‑71 (MQ=255)
           gTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATccac                                                                    <  1:444999/61‑1 (MQ=255)
           gTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATccac                                                                    <  1:1864993/61‑1 (MQ=255)
           gTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATccac                                                                    <  1:3232997/61‑1 (MQ=255)
             aGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATccacca                                                                  >  1:4102139/1‑61 (MQ=255)
             aGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTAt                                                        >  1:5560486/1‑71 (MQ=255)
              gCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGAAACCCTATc                                                       >  1:2644426/1‑71 (MQ=255)
               cTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCg                                                      <  1:1158534/71‑1 (MQ=255)
                ttATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATcca                                                                     >  1:5611922/1‑55 (MQ=255)
                 tATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATccacca                                                                  >  1:612211/1‑57 (MQ=255)
                  aTCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGcc                                                    <  1:1981326/70‑1 (MQ=255)
                  aTCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGcc                                                    <  1:5123431/70‑1 (MQ=255)
                  aTCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGcc                                                    <  1:3455893/70‑1 (MQ=255)
                     ttCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCa                                                >  1:2154697/1‑71 (MQ=255)
                       ccGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAg                                                  <  1:2701471/67‑1 (MQ=255)
                       ccGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAg                                                  >  1:5306042/1‑67 (MQ=255)
                       ccGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCAc                                               >  1:36240/1‑70 (MQ=255)
                       ccGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCAc                                               <  1:3462431/70‑1 (MQ=255)
                       ccGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCAc                                               >  1:727284/1‑70 (MQ=255)
                       ccGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCAc                                               >  1:4482433/1‑70 (MQ=255)
                       ccGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCAc                                               >  1:1496372/1‑70 (MQ=255)
                       ccGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCAc                                               >  1:1000977/1‑70 (MQ=255)
                        cGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATccacca                                                                  >  1:5595408/1‑50 (MQ=255)
                        cGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATccacca                                                                  >  1:13833/1‑50 (MQ=255)
                           aaTGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAAcc                                                            >  1:487722/1‑53 (MQ=255)
                           aaTGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCCTc                                           >  1:2445586/1‑70 (MQ=255)
                             tGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCg                                        <  1:239552/71‑1 (MQ=255)
                              gACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCg                                        >  1:3710229/1‑70 (MQ=255)
                              gACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCg                                        >  1:513850/1‑70 (MQ=255)
                              gACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCg                                        >  1:4479163/1‑70 (MQ=255)
                               acacGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTa                                                         <  1:3211752/52‑1 (MQ=255)
                               acacGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCa                                             >  1:4387819/1‑64 (MQ=255)
                               acacGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCa                                             >  1:2238589/1‑64 (MQ=255)
                               acacGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGaa                                      <  1:2038152/71‑1 (MQ=255)
                               acacGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGaa                                      <  1:902688/71‑1 (MQ=255)
                               acacGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGaa                                      >  1:1026468/1‑71 (MQ=255)
                               acacGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGa                                       <  1:4655363/70‑1 (MQ=255)
                               acacGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGa                                       <  1:1874588/70‑1 (MQ=255)
                                 gcGCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAACa                                    <  1:159455/70‑1 (MQ=2)
                                 acGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATc                                                       >  1:243287/1‑52 (MQ=255)
                                 acGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACa                                    <  1:3496881/71‑1 (MQ=255)
                                 acGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACa                                    <  1:3254561/71‑1 (MQ=255)
                                 acGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACa                                    <  1:5507608/71‑1 (MQ=255)
                                 acGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACa                                    >  1:610043/1‑71 (MQ=255)
                                   gCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACAt                                   <  1:3537085/70‑1 (MQ=255)
                                    cTCGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAACAtcg                                 <  1:602302/71‑3 (MQ=2)
                                    cTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAg                                 >  1:3221352/1‑71 (MQ=255)
                                      cGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGt                                >  1:5515794/1‑70 (MQ=255)
                                       gACGTTCTCCACCACCACGATGGCGTCATCCACCAGTATCCCTATCGCCAGCACCTTCCCg                                        <  1:3733546/61‑1 (MQ=255)
                                       gACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTc                               <  1:5561870/70‑1 (MQ=255)
                                       gACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTc                               >  1:1944879/1‑70 (MQ=255)
                                       gACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTc                               <  1:849761/70‑1 (MQ=255)
                                       gACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTc                               >  1:5289962/1‑70 (MQ=255)
                                       gACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTc                               >  1:5310107/1‑70 (MQ=255)
                                       gACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCa                              <  1:5507873/71‑1 (MQ=255)
                                       gACGTTCTCCACCACCACGATGGCATCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCa                              <  1:3569061/71‑1 (MQ=255)
                                        aCGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTc                               <  1:146853/69‑1 (MQ=255)
                                            tctcCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACgt                          <  1:2751877/70‑1 (MQ=255)
                                             ctcCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCaa                             >  1:4981563/1‑66 (MQ=255)
                                               ccaccaccaCGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGtt                       >  1:1307566/1‑70 (MQ=255)
                                               ccaccaccaCGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGtt                       <  1:4429117/70‑1 (MQ=255)
                                               ccaccaccaCGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGtt                       <  1:11034/70‑1 (MQ=255)
                                                 accaccacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGAt                    >  1:5508989/1‑71 (MQ=255)
                                                  ccaccacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCAcc                                              >  1:2825773/1‑44 (MQ=255)
                                                  ccaccacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCAcc                                              >  1:3903108/1‑44 (MQ=255)
                                                   caccacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATcc                                          >  1:2570334/1‑47 (MQ=255)
                                                   caccacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATcc                                          >  1:1357341/1‑47 (MQ=255)
                                                   caccacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAAc                                     <  1:2948496/52‑1 (MQ=255)
                                                    accacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGt                 <  1:1407205/71‑1 (MQ=255)
                                                    accacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGt                 <  1:1467012/71‑1 (MQ=255)
                                                    accacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGt                 <  1:914571/71‑1 (MQ=255)
                                                     ccacGATGGCGTCATCCACCAGTAACCCTATCGcc                                                    <  1:5219577/35‑1 (MQ=255)
                                                     ccacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCAt                                            >  1:5176160/1‑43 (MQ=255)
                                                     ccacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGa                                       <  1:2856327/48‑1 (MQ=255)
                                                     ccacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACa                                    >  1:3947346/1‑51 (MQ=255)
                                                     ccacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACa                                    >  1:1259390/1‑51 (MQ=255)
                                                     ccacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACAt                                   >  1:4855156/1‑52 (MQ=255)
                                                       acGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGa               >  1:5439424/1‑70 (MQ=255)
                                                          aTGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAacc           <  1:5213378/71‑1 (MQ=255)
                                                           tGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGa               <  1:4520277/66‑1 (MQ=255)
                                                           tGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAacc           >  1:2316960/1‑70 (MQ=255)
                                                             gCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAaccga         <  1:266197/70‑1 (MQ=255)
                                                                 cATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAaccga         >  1:1072437/1‑66 (MQ=255)
                                                                  aTCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGt                                <  1:4471326/42‑1 (MQ=255)
                                                                  aTCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACgtgt                        >  1:4056023/1‑50 (MQ=255)
                                                                  aTCCACCAGTAACCCTATCGCCAGCACCATCACGAACATAGTCAACgtgt                        >  1:5002358/1‑50 (MQ=255)
                                                                   tCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGa  <  1:1654899/71‑1 (MQ=255)
                                                                    ccaccaGTAACCCTATCGCCAGCACCATCCCGAACATAGTCa                              <  1:5265999/42‑1 (MQ=255)
                                                                     caccaGTAACCCTATCGCCAGCACCATCCCGAACATa                                  >  1:5619273/1‑37 (MQ=255)
                                                                     |                                                                    
TTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCGACCGCCGA  >  minE/2585990‑2586127

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: