Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,584,820 G→A 8.7% G819G (GGC→GGT mdtF ← multidrug transporter, RpoS‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,584,8200GA 8.7% 213.2 / 8.2 69G819G (GGC→GGTmdtFmultidrug transporter, RpoS‑dependent
Reads supporting (aligned to +/- strand):  ref base G (34/29);  new base A (2/4);  total (36/33)
Fisher's exact test for biased strand distribution p-value = 4.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.38e-01

AAATTTCATGGCGTCACCGGTACTTTTCCCGGCCGCCGCTTCACCTAAAATCTCCATTGACGGGATGCCGTTGTAGCGTTCCAGTCGCGGTGAACCATAGGTCCATTCGGTAGACGAGTAGGCAGAAAGCGGTGC  >  minE/2584754‑2584888
                                                                  |                                                                    
aaaTTTCATGGCGTCACCGGTACTTTTCCCGGCCGCCGCTTCACCTAAAATCTCCATTGACGGGATGCCGt                                                                  >  1:871210/1‑71 (MQ=255)
aaaTTTCATGGCGTCACCGGTACTTTTCCCGGCCGCCGCTTCACCTAAAATCTCCATTGACGGGATGCCGt                                                                  >  1:1591199/1‑71 (MQ=255)
 aaTTTCATGGCGTCACCGGTACTTTTCCCGGCCGCCGCTTCACCTAAAATCTCCATTGACGGGATGCCGtt                                                                 <  1:291906/71‑1 (MQ=255)
 aaTTTCATGGCGTCACCGGTACTTTTCCCGGCCGCCGCTTCACCTAAAATCTCCATTGACGGGATGCCGtt                                                                 >  1:46264/1‑71 (MQ=255)
      cATGGCGTCACCGGTACTTTTCCCGGCCGCCGCTTCACCTAAAATCTCCATTGACGGGATGCCGTTGTAGc                                                            >  1:3851903/1‑71 (MQ=255)
      cATGGCGTCACCGGTACTTTTCCCGGCCGCCGCTTCACCTAAAATCTCCATTGACGGGATGCCGTTGTAGc                                                            >  1:3856223/1‑71 (MQ=255)
          gCGTCACCGGTACTTTTCCCGGCCGCCGCTTCACCTAAAATCTCCATTGACGGGATGCCGTTGTAGCGTTc                                                        <  1:476602/71‑1 (MQ=255)
             tCACCGGTACTTTTCCCGGCCGCCGCTTCACCTAAAATCTCCATTGACGGGATGCCg                                                                   <  1:3469748/57‑1 (MQ=255)
               aCCGGTACTTTTCCCGGCCGCCGCTTCACCTAAAATCTCCATTGACGGGATGCCGTTg                                                                >  1:2064563/1‑58 (MQ=255)
               aCCGGTACTTTTCCCGGCCGCCGCTTCACCTAAAATCTCCATTGACGGGATGCCGTTGTa                                                              >  1:3543763/1‑60 (MQ=255)
                 cGGTACTTTTCCCGGCCGCCGCTTCACCTAAAATCTCCATTGACGGGATGCCGTTGTAGCGTTCCAGTcgc                                                 <  1:1214200/71‑1 (MQ=255)
                  ggTACTTTTCCCGGCCGCCGCTTCACCTAAAATCTCCATTGACGGGATGCCGTTGTAGCGTTCCAGTcgcg                                                <  1:887341/71‑1 (MQ=255)
                       ttttCCCGGCCGCCGCTTCACCTAAAATCTCCATTGACGGGATGCCGTTGt                                                               <  1:1448764/51‑1 (MQ=255)
                       ttttCCCGGCCGCCGCTTCACCTAAAATCTCCATTGACGGGATGCCGTTGTa                                                              <  1:322310/52‑1 (MQ=255)
                        tttCCCGGCCGCCGCTTCACCTAAAATCTCCATTGACGGGATGCCGTTGTAGCGTTCCAGTCGCGGTGAAc                                          >  1:2404640/1‑71 (MQ=255)
                        tttCCCGGCCGCCGCTTCACCTAAAATCTCCATTGACGGGATGCCGTTGTAGCGTTCCAGTCGCGGTGAAc                                          >  1:566124/1‑71 (MQ=255)
                         ttCCCGGCCGCCGCTTCACCTAAAATCTCCATTGATGGGATGCCGTTGTAGCGTTCCAGTCGCGGTGAAc                                          >  1:1982648/1‑70 (MQ=255)
                         ttCCCGGCCGCCGCTTCACCTAAAATCTCCATTGACGGGATGCCGTTGTAGCGTTCCAg                                                     >  1:3266592/1‑59 (MQ=255)
                         ttCCCGGCCGCCGCTTCACCTAAAATCTCCATTGACGGGATGCCGTTGTAGCGTTCCAGTCGCGGTGaa                                           <  1:2093379/69‑1 (MQ=255)
                                  gccgcTTCACCTAAAATCTCCATTGACGGGATGCCGTTGTAGCGTTCCAGTCGCGGTGAACCATAGGTCCa                                <  1:3766289/71‑1 (MQ=255)
                                      cTTCACCTAAAATCTCCATTGACGGGATGCCGTTg                                                                <  1:3826346/35‑1 (MQ=255)
                                        tCACCTAAAATCTCCATTGACGGGATGCCGTTGTAGCGTTCCAGTCGCGGt                                              >  1:1408989/1‑51 (MQ=255)
                                        tCACCTAAAATCTCCATTGACGGGATGCCGTTGTAGCGTTCCAGTCGCGGTGAAc                                          <  1:1954103/55‑1 (MQ=255)
                                          aCCTAAAATCTCCATTGACGGGATGCCGTTGTAGCGTTCCAGTCGCGGTGAACCATAGGTCCATTCGGTa                         <  1:438634/70‑1 (MQ=255)
                                            cTAAAATCTCCATTGACGGGATGCCGTTGTAGCGTTCCAGTCGCGGTGAACCATa                                      >  1:3111644/1‑55 (MQ=255)
                                             tAAAATCTCCATTGACGGGATGCCGTTGTAGCGTTCCAGTCGCGGt                                              >  1:2109649/1‑46 (MQ=255)
                                             tAAAATCTCCATTGACGGGATGCCGTTGTAGCGTTCCAGTCGCGGTGAACCa                                        <  1:2430727/52‑1 (MQ=255)
                                             tAAAATCTCCATTGACGGGATGCCGTTGTAGCGTTCCAGTCGCGGTGAACCATAGGTCCa                                <  1:1076991/60‑1 (MQ=255)
                                             tAAAATCTCCATTGACGGGATGCCGTTGTAGCGTTCCAGTCGCGGTGAACCATAGGTCCATTCGGt                          >  1:2330216/1‑66 (MQ=255)
                                              aaaaTCTCCATTGACGGGATGCCGTTGTAGCGTTCCAGTCGCGGTGAACCATAGGTCCATTCGGTAGAc                      >  1:3488731/1‑69 (MQ=255)
                                               aaaTCTCCATTGACGGGATGCCGTTGTAGCGTTCCAGTCGCGGTGAAc                                          <  1:1051946/48‑1 (MQ=255)
                                               aaaTCTCCATTGACGGGATGCCGTTGTAGCGTTCCAGTCGCGGTGAACCATAg                                     >  1:1176556/1‑53 (MQ=255)
                                               aaaTCTCCATTGACGGGATGCCGTTGTAGCGTTCCAGTCGCGGTGAACCATAGGTCCATTCGGTAGACGa                    <  1:932127/70‑1 (MQ=255)
                                               aaaTCTCCATTGACGGGATGCCGTTGTAGCGTTCCAGTCGCGGTGAACCATAGGTCCATTCGGTAGACGAg                   >  1:1196257/1‑71 (MQ=255)
                                               aaaTCTCCATTGACGGGATGCCGTTGTAGCGTTCCAGTCGCGGTGAACCATAGGTCCATTCGGTAGACGAg                   >  1:328266/1‑71 (MQ=255)
                                                aaTCTCCATTGACGGCAGACCGTTGTAGCGTTCCAGTcgcg                                                <  1:1413110/41‑1 (MQ=37)
                                                 aTCTCCATTGACGGGATGCCGTTGTAGCGTTCCAGTCGCGGTGAAc                                          >  1:316198/1‑46 (MQ=255)
                                                  tctcCATTGACGGGATGCCGTTGTAGCGTTCCAGTCGCGGTGAACCATAg                                     >  1:1243760/1‑50 (MQ=255)
                                                   ctcCATTGACGGGATGCCGTTGTAGCGTTCCAGTCGCGGTGAACCATAGGTCCATTCGGTAGACg                     <  1:2841894/65‑1 (MQ=255)
                                                    tcCATTGACGGGATGCCGTTGTAGGGTTCCAGTCGCGGTGAACCATAGGTCCATTCGGTAGACGAGTAgg               >  1:626195/1‑70 (MQ=255)
                                                    tcCATTGACGGGATGCCGTTGTAGCGTTCCAGTCGCGGTGAACCAt                                       <  1:3789759/46‑1 (MQ=255)
                                                      cATTGACGGGATGCCGTTGTAGCGTTCCAGTCGCGGTGAACCATAGGTCCATTCGGTAGAc                      <  1:3869215/61‑1 (MQ=255)
                                                      cATTGACGGCAGACCGTTGTAGCGTTCCAGTCGCGGAGAGCCATac                                     >  1:855210/1‑45 (MQ=255)
                                                       aatGACGGGATGCCGTTGAAGCGTTCCAGTCGCGGTGAACCATACGTCCATTCGGTAGACg                     >  1:908317/3‑61 (MQ=255)
                                                        ttGACGGTAGACCGTTGTAGCGTTCCAGTCGCGGAGAGCCa                                        <  1:493807/41‑1 (MQ=255)
                                                        ttGACGGGATGCCGTTGTAGCGTTCCAGTCGCGGTGAACCAt                                       <  1:3410269/42‑1 (MQ=255)
                                                        ttGACGGGATGCCGTTGTAGCGTTCCAGTCGCGGTGAACCATAGGTCCATTCGGTAGACGAGTAGGCAg            <  1:3431801/69‑1 (MQ=255)
                                                        ttGACGGGATGCCGTTGTAGCGTTCCAGTCGCGGTGAACCATAGGTCCATTCGGTAGACGAGTAGGCAg            <  1:3384800/69‑1 (MQ=255)
                                                        ttGACGGGATGCCGTTGTAGCGTTCCAGTCGCGGTGAACCATAGGTCCATTCGGTAGACGAGTAGGCAGaa          <  1:1562975/71‑1 (MQ=255)
                                                        ttGACGGCAGACCGTTGTAGCGTTCCAGTCGCGGAGAGCCa                                        <  1:2176676/41‑1 (MQ=255)
                                                        ttGACGGCAGACCGTTGTAGCGTTCCAGTCGCGGAGAGCCa                                        <  1:3305448/41‑1 (MQ=255)
                                                         tGACGGGATGCCGTTGTAGCGTTCCAGTCGCGGTGAACCATAGGTCCATTCGGTAGACGAGTAGGCAGa           >  1:13702/1‑69 (MQ=255)
                                                          gACGGGATGCCGTTGTAGCGTTCCAGTCGCGGTGAACCATa                                      >  1:1388946/1‑41 (MQ=255)
                                                          gACGGGATGCCGTTGTAGCGTTCCAGTCGCGGTGAACCATAGGTCCATTCGGTAGACGAGTAGGCAGaaa         >  1:379974/1‑70 (MQ=255)
                                                          gACGGGATACCGTTGTAGCGTTCCAGTCGCGGTGAACCATAGGTCCATTCGGTAGACGa                    >  1:3422552/1‑59 (MQ=255)
                                                           aCGGGATGCCGTTGTAGCGTTCCAGTCGCGGTGAACCATAGGTCCATTCGGTAGACGAGTAGGCAGAAAg        <  1:2937752/70‑1 (MQ=255)
                                                           aCGGGATGCCGTTGTAGCGTTCCAGTCGCGGTGAACCATAGGTCCATTCGGTAGACGAGTAGGCAGAAAg        <  1:3825307/70‑1 (MQ=255)
                                                            cGGGATGCCGTTGTAGCGTTCCAGTCGCGGTGAACCATAGGTCCATTCGGTAGACGAGTAGGCAGAAAGCg      >  1:1764679/1‑71 (MQ=255)
                                                            cGGGATGCCGTTGTAGCGTTCCAGTCGCGGTGAACCATAGGTCCATTCGGTAGACGAGTAGGCAGAAAGCg      >  1:402288/1‑71 (MQ=255)
                                                            cGGGATGCCGTTGTAGCGTTCCAGTCGCGGTGAACCATAGGTCCATTCGGTAGACGAGTAGGCAGAAAGCg      >  1:2158479/1‑71 (MQ=255)
                                                            cGGGATGCCGTTGTAGCGTTCCAGTCGCGGTGAACCATAGGTCCATTCGGTAGACGAGTAGGCAGAAAGCg      >  1:2572199/1‑71 (MQ=255)
                                                             gggATGCCGTTGTAGCGTTCCAGTCGCGGTGAACCATAGGTCCATTCGGTa                         <  1:3721744/51‑1 (MQ=255)
                                                               gATGCCGTTGTAGCGTTCCAGTCGCGGTGAACCATAGGTCCATTCGGTAGACGAGTAGGCAGAAAGCGGTg   >  1:3675228/1‑71 (MQ=255)
                                                                aTGCCGTTGTAGCGTTCCAGTCGCGGTGAACCATAGGTCCATTCGGTAGACGAGTAGGCAGAAAGCg      <  1:2459784/67‑1 (MQ=255)
                                                                aTGCCGTTGTAGCGTTCCAGTCGCGGTGAACCATAGGTCCATTCGGTAGACGAGTAGGCAGAAAGCg      <  1:2960790/67‑1 (MQ=255)
                                                                aTGCCGTTGTAGCGTTCCAGTCGCGGTGAACCATAGGTCCATTCGGTAGACGAGTAGGCAGAAAGCGGTGc  <  1:1714575/71‑1 (MQ=255)
                                                                aTGCCGTTGTAGCGTTCCAGTCGCGGTGAACCATAGGTCCATTCGGTAGACGAGTAGGCAGAAAGCGGTGc  >  1:1827325/1‑71 (MQ=255)
                                                                 tGCCGTTGTAGCGTTCCAGTCGCGGTGAACCATAGGTCCATTCGGTAGACGAGTAGGc              >  1:3750150/1‑58 (MQ=255)
                                                                 tGCCGTTGTAGCGTTCCAGTCGCGGTGAACCATAGGTCCATTCGGTAGACGAGTAGGCAGAAAGCGGTGc  >  1:2922522/1‑70 (MQ=255)
                                                                  gCCGTTGTAGCGTTCCAGTCGCGGTGAACCATAGGTCCATTCGGTAGACGAGTAGGCAGa           <  1:334843/60‑1 (MQ=255)
                                                                  gCCGTTGTAGCGTTCCAGTCGCGGTGAACCATAGGTCCATTCGGTAGACGAGTAGGCAGa           <  1:10468/60‑1 (MQ=255)
                                                                  |                                                                    
AAATTTCATGGCGTCACCGGTACTTTTCCCGGCCGCCGCTTCACCTAAAATCTCCATTGACGGGATGCCGTTGTAGCGTTCCAGTCGCGGTGAACCATAGGTCCATTCGGTAGACGAGTAGGCAGAAAGCGGTGC  >  minE/2584754‑2584888

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: