Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 534,192 A→G 51.3% G415G (GGT→GGC ybhF ← fused predicted transporter subunits and ATP‑binding components of ABC superfamily

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE534,1920AG51.3% 0.3 / 72.3 80G415G (GGT→GGCybhFfused predicted transporter subunits and ATP‑binding components of ABC superfamily
Reads supporting (aligned to +/- strand):  ref base A (16/23);  new base G (29/12);  total (45/35)
Fisher's exact test for biased strand distribution p-value = 1.27e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGCCCGACCGCGTAAGCCATACACACCAGAGAAAAAGCGTAAATTCTGTTCGACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAACTCTCTTTCAG  >  minE/534128‑534260
                                                                |                                                                    
gcgcCCGACCGCGTAAGCCATACACACCAGAGAAAAAGCGTAAATTCTGTTCGACCGTCAGGTTACCg                                                                   <  1:2571085/68‑1 (MQ=255)
 cgcCCGACCGCGTAAGCCATACACACCAGAGAAAAAGCGTAAATTCTGTTCGACCGTCAGGTTACCg                                                                   >  1:1032716/1‑67 (MQ=255)
   cccGACCGCGTAAGCCATACACACCAGAGAAAAAGTGTAAATTCTGTTCGACCGTCAGGTTGCCGTa                                                                 <  1:1868100/67‑1 (MQ=255)
   cccGACCGCGTAAGCCATACACACCAGAGAAAAAGCGTAAATTCTGTTCGACCGTCAGGTTGCCGTa                                                                 <  1:1073692/67‑1 (MQ=255)
   cccGACCGCGTAAGCCATACACACCAGAGAAAAAGCGTAAATTCTGTTCGACCGTCAGGTTGCCGTa                                                                 <  1:911235/67‑1 (MQ=255)
   cccGACCGCGTAAGCCATACACACCAGAGAAAAAGCGTAAATTCTGTTCGACCGTCAGGTTGCCGTa                                                                 <  1:723565/67‑1 (MQ=255)
   cccGACCGCGTAAGCCATACACACCAGAGAAAAAGCGTAAATTCTGTTCGACCGTCAGGTTGCCGTa                                                                 <  1:588648/67‑1 (MQ=255)
   cccGACCGCGTAAGCCATACACACCAGAGAAAAAGCGTAAATTCTGTTCGACCGTCAGGTTGCCGTa                                                                 <  1:3192614/67‑1 (MQ=255)
   cccGACCGCGTAAGCCATACACACCAGAGAAAAAGCGTAAATTCTGTTCGACCGTCAGGTTGCCGTa                                                                 <  1:3025563/67‑1 (MQ=255)
   cccGACCGCGTAAGCCATACACACCAGAGAAAAAGCGTAAATTCTGTTCGACCGTCAGGTTGCCGTa                                                                 <  1:2981765/67‑1 (MQ=255)
   cccGACCGCGTAAGCCATACACACCAGAGAAAAAGCGTAAATTCTGTTCGACCGTCAGGTTGCCGTa                                                                 <  1:2929101/67‑1 (MQ=255)
   cccGACCGCGTAAGCCATACACACCAGAGAAAAAGCGTAAATTCTGTTCGACCGTCAGGTTGCCGTa                                                                 <  1:2018360/67‑1 (MQ=255)
   cccGACCGCGTAAGCCATACACACCAGAGAAAAAGCGTAAATTCTGTTCGACCGTCAGGTTGCCGTa                                                                 <  1:1406937/67‑1 (MQ=255)
   cccGACCGCGTAAGCCATACACACCAGAGAAAAAGCGTAAATTCTGTTCGACCGTCAGGTTGCCGTa                                                                 <  1:1316844/67‑1 (MQ=255)
                 ccATACACACCAGAGAAAAAGCGTAAATTCTGTTCGACCGTCAGGTTACCGTAGAGCGAAAATTTTTgc                                                 >  1:2677336/1‑69 (MQ=255)
                 ccATACACACCAGAGAAAAAGCGTAAATTCTGTTCGACCGTCAGGTTACCGTAGAGCGAAAATTTTTgc                                                 >  1:2498919/1‑69 (MQ=255)
                 ccATACACACCAGAGAAAAAGCGTAAATTCTGTTCGACCGTCAGGTTACCGTAGAGCGAAAATTTTTgc                                                 >  1:2143245/1‑69 (MQ=255)
                 ccATACACACCAGAGAAAAAGCGTAAATTCTGTTCGACCGTCAGGTTACCGTAGAGCGAAAATTTTTgc                                                 >  1:3090988/1‑69 (MQ=255)
                 ccATACACACCAGAGAAAAAGCGTAAATTCTGTTCGACCGTCAGGTTACCGTAGAGCGAAAATTTTTgc                                                 >  1:309795/1‑69 (MQ=255)
                 ccATACACACCAGAGAAAAAGCGTAAATTCTGTTCGACCGTCAGGTTACCGTAGAGCGAAAATTTTTgc                                                 >  1:1922549/1‑69 (MQ=255)
                 ccATACACACCAGAGAAAAAGCGTAAATTCTGTTCGACCGTCAGGTTACCGTAGAGCGAAAATTTTTgc                                                 >  1:1958226/1‑69 (MQ=255)
                      cacacCAGAGAAAAAGCGTAAATTCTGTTCGACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCa                                             >  1:2432697/1‑68 (MQ=255)
                      cacacCAGAGAAAAAGCGTAAATTCTGTTCGACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCa                                             >  1:1184381/1‑68 (MQ=255)
                            agagAAAAAGCGTAAATTCTGTTCGACCGTCCGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGcc                                      <  1:876962/69‑1 (MQ=255)
                            agagAAAAAGCGTAAATTCTGTTCGACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGcc                                      <  1:748063/69‑1 (MQ=255)
                             gagaAAAAGCGTAAATTCTGTTCGACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGcc                                      <  1:555376/68‑1 (MQ=255)
                             gagaAAAAGCGTAAATTCTGTTCGACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGcc                                      <  1:3283215/68‑1 (MQ=255)
                              agaAAAAGCGTAAATTCTGTTCTACCGTCAGGTTGCCGTAGAGCGAAAATTTTTgcg                                                >  1:2130328/1‑57 (MQ=255)
                              agaAAAAGCGTAAATTCTGTTCGACCGTCAGGTTGCCGTAGAGCGAAAATTTTTgcg                                                >  1:2339638/1‑57 (MQ=255)
                              agaAAAAGCGTAAATTCTGTTCGACCGTCAGGTTGCCGTAGAGCGAAAATTTTTgcg                                                >  1:1143920/1‑57 (MQ=255)
                              agaAAAAGCGTAAATTCTGTTCGACCGTCAGGTTGCCGTAGAGCGAAAATTTTTgcg                                                >  1:2715885/1‑57 (MQ=255)
                              agaAAAAGCGTAAATTCTGTTCGACCGTCAGGTTGCCGTAGAGCGAAAATTTTTgcg                                                >  1:1423905/1‑57 (MQ=255)
                              agaAAAAGCGTAAATTCTGTTCGACCGTCAGGTTGCCGTAGAGCGAAAATTTTTgcg                                                >  1:2632049/1‑57 (MQ=255)
                              agaAAAAGCGTAAATTCTGTTCGACCGTCAGGTTGCCGTAGAGCGAAAATTTTTgcg                                                >  1:257398/1‑57 (MQ=255)
                              agaAAAAGCGTAAATTCTGTTCGACCGTCAGGTTGCCGTAGAGCGAAAATTTTTgcg                                                >  1:1236834/1‑57 (MQ=255)
                              agaAAAAGCGTAAATTCTGTTCGACCGTCAGGTTGCCGTAGAGCGAAAATTTTTgcg                                                >  1:692357/1‑57 (MQ=255)
                              agaAAAAGCGTAAATTCTGTTCGACCGTCAGGTTGCCGTAGAGCGAAAATTTTTgcg                                                >  1:759308/1‑57 (MQ=255)
                              agaAAAAGCGTAAATTCTGTTCGACCGTCAGGTTGCCGTAGAGCGAAAATTTTTgcg                                                >  1:1463877/1‑57 (MQ=255)
                              agaAAAAGCGTAAATTCTGTTCGACCGTCAGGTTGCCGTAGAGCGAAAATTTTTgcg                                                >  1:246635/1‑57 (MQ=255)
                              agaAAAAGCGTAAATTCTGTTCGACCGTCAGGTTGCCGTAGAGCGAAAATTTTTgcg                                                >  1:3361580/1‑57 (MQ=255)
                              agaAAAAGCGTAAATTCTGTTCGACCGTCAGGTTGCCGTAGAGCGAAAATTTTTgcg                                                >  1:1221088/1‑57 (MQ=255)
                              agaAAAAGCGTAAATTCTGTTCGACCGTCAGGTTGCCGTAGAGCGAAAATTTTTgcg                                                >  1:1934747/1‑57 (MQ=255)
                              agaAAAAGCGTAAATTCTGTTCGACCGTCAGGTTGCCGTAGAGCGAAAATTTTTgcg                                                >  1:2934359/1‑57 (MQ=255)
                              agaAAAAGCGTAAATTCTGTTCGACCGTCAGGTTGCCGTAGAGCGAAAATTTTTgcg                                                >  1:2969123/1‑57 (MQ=255)
                              agaAAAAGCGTAAATTCTGTTCGACCGTCAGGTTGCCGTAGAGCGAAAATTTTTgcg                                                >  1:3340730/1‑57 (MQ=255)
                              agaAAAAGCGTAAATTCTGTTCGACCGTCAGGTTGCCGTAGAGCGAAAATTTTTgcg                                                >  1:3035041/1‑57 (MQ=255)
                              agaAAAAGCGTAAATTCTGTTCGACCGTCAGGTTGCCGTAGAGCGAAAATTTTTgcg                                                >  1:1710011/1‑57 (MQ=255)
                              agaAAAAGCGTAAATTCTGTTCGACCGTCAGGTTGCCGTAGAGCGAAAATTTTTgcg                                                >  1:2146451/1‑57 (MQ=255)
                              agaAAAAGCGTAAATTCTGTTCGACCGTCAGGTTGCCGTAGAGCGAAAATTTTTgcg                                                >  1:1737741/1‑57 (MQ=255)
                              agaAAAAGCGTAAATTCTGTTCGACCGTCAGGTTGCCGTAGAGCGAAAATTTTTgcg                                                >  1:3337634/1‑57 (MQ=255)
                              agaAAAAGCGTAAATTCTGTTCGACCGTCAGGTTGCCGTAGAGCGAAAATTTTTgcg                                                >  1:2054596/1‑57 (MQ=255)
                              agaAAAAGCGTAAATTCTGTTCGACCGTCAGGTTGCCGTAGAGCGAAAATTTTTgcg                                                >  1:1892002/1‑57 (MQ=255)
                              agaAAAAGCGTAAATTCTGTTCGACCGTCAGGTTGCCGTAGAGCGAAAATTTTTgcg                                                >  1:1905322/1‑57 (MQ=255)
                              agaAAAAGCGTAAATTCTGTTCGACCGTCAGGTTGCCGTAGAGCGAAAATTTTTgc                                                 >  1:1146711/1‑56 (MQ=255)
                              agaAAAAGCGTAAATTATGTTCGACCGTCAGGTTGCCGTAGAGCGAAAATTTTTgcg                                                >  1:2201899/1‑57 (MQ=255)
                              agaAAAAGAGTAAATTCTGTTCGACCGTCAGGTTGCCGTAGAGCGAAAATTTTTgcg                                                >  1:199671/1‑57 (MQ=255)
                                                   cGACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCgg               <  1:2793860/69‑1 (MQ=255)
                                                   cGACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCgg               <  1:1646456/69‑1 (MQ=255)
                                                   cGACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCgg               <  1:1808596/69‑1 (MQ=255)
                                                   cGACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCgg               <  1:2958087/69‑1 (MQ=255)
                                                   cGACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCgg               <  1:878618/69‑1 (MQ=255)
                                                        gTCAGGTTACCGTAGAGCGAAAATTTTTGCGGCATATAGCCGAGATGCTGGCGCGCTTTACCGGAActc          <  1:2820723/69‑1 (MQ=255)
                                                        gTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGATATGCTGGCGCGCTTTACCGGAActc          <  1:118675/69‑1 (MQ=255)
                                                        gTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAActc          <  1:2217191/69‑1 (MQ=255)
                                                        gTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAActc          <  1:224138/69‑1 (MQ=255)
                                                        gTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAActc          <  1:1301606/69‑1 (MQ=255)
                                                        gTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAActc          <  1:1147185/69‑1 (MQ=255)
                                                        gTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAActc          <  1:2259108/69‑1 (MQ=255)
                                                        gTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAActc          <  1:3334642/69‑1 (MQ=255)
                                                        gTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAActc          <  1:2344216/69‑1 (MQ=255)
                                                        gTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAActc          <  1:2521825/69‑1 (MQ=255)
                                                        gTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAActc          <  1:1632847/69‑1 (MQ=255)
                                                        gTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAActc          <  1:518376/69‑1 (MQ=255)
                                                        gTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAActc          <  1:133726/69‑1 (MQ=255)
                                                              ttACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAActctct       >  1:1321178/1‑66 (MQ=255)
                                                              ttACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAActctct       >  1:2779417/1‑66 (MQ=255)
                                                              ttACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAActctct       >  1:2520428/1‑66 (MQ=255)
                                                              ttACCGTAGAGCGAAAATTTTT‑‑GCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAActctct       >  1:2850313/1‑64 (MQ=255)
                                                               tACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGAt                                 >  1:2185845/1‑39 (MQ=255)
                                                               tACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGAt                                 >  1:908630/1‑39 (MQ=255)
                                                                aCCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAACTCTCTTTCAg  <  1:2776532/69‑1 (MQ=255)
                                                                aCCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAACTCTCTTTCAg  <  1:2089268/69‑1 (MQ=255)
                                                                aCCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAACTCTCTTTCAg  <  1:111641/69‑1 (MQ=255)
                                                                |                                                                    
GCGCCCGACCGCGTAAGCCATACACACCAGAGAAAAAGCGTAAATTCTGTTCGACCGTCAGGTTACCGTAGAGCGAAAATTTTTGCGCCATATAGCCGAGATGCTGGCGCGCTTTACCGGAACTCTCTTTCAG  >  minE/534128‑534260

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: