Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,408,722 G→A 37.1% A331T (GCG→ACG)  yejF → fused predicted oligopeptide 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
*minE1,408,7220GA37.1% 30.3 / 33.9 35A331T (GCG→ACG) yejFfused predicted oligopeptide transporter subunits and ATP‑binding components of ABC superfamily
Reads supporting (aligned to +/- strand):  ref base G (6/16);  new base A (3/10);  total (9/26)
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.67e-01

CTACGAGCGGGTGAAACACTGGGTTTAGTGGGCGAGTCCGGTTCCGGGAAAAGTACGACGGGACTGGCGCTGCTGCGACTGATTAATTCTCAGGGCAGCATCATCTTTGACGGTCAGCCACTGCAAAA  >  minE/1408656‑1408783
                                                                  |                                                             
cTACGAGCGGGTGAAACACTGGGTTTAGTGGGCGAGTCCGGTTCCGGGAAAAGTACGACGGGACTGGCgct                                                           <  1:3661995/71‑1 (MQ=255)
cTACGAGCGGGTGAAACACTGGGTTTAGTGGGCGAGTCCGGTTCCGGGAAAAGTACGACGGGACTGACgct                                                           <  1:2217939/71‑1 (MQ=255)
     aGCGGGTGAAACACTGGGTTTAGTGGGCGAGTCCGGTTCCGGGAAAAGTACGACGGGACTGGCGCtgctgc                                                      <  1:2820745/71‑1 (MQ=255)
     aGCGGGTGAAACACTGGGTTTAGTGGGCGAGTCCGGTTCCGGGAAAAGTACGACGGGACTGGCGCtgctgc                                                      <  1:2534652/71‑1 (MQ=255)
     aGCGGGTGAAACACTGGGTTTAGTGGGCGAGTCCGGTTCCGGGAAAAGTACGACGGGACTGGCGCtgctgc                                                      <  1:1892244/71‑1 (MQ=255)
                            tGGGCGAGTCCGGTTCCGGGAAAAGTACGACGGGACTGGCGCtgctgc                                                      <  1:2550218/48‑1 (MQ=255)
                            tGGGCGAGTCCGGTTCCGGGAAAAGTACGACGGGACTGGCGCtgctgc                                                      <  1:94885/48‑1 (MQ=255)
                              ggCGAGTCCGGTTCCGGGAAAAGTACGACGGGACTGGCGCTGCTGCGACTGATTAATTCTCAGGGCAGcat                             >  1:3772550/1‑71 (MQ=255)
                               gCGAGTCCGGTTCCGGGAAAAGTACGACGGGACTGGCgct                                                           <  1:2086039/40‑1 (MQ=255)
                               gCGAGTCCGGTTCCGGGAAAAGTACGACGGGACTGGCgct                                                           <  1:618386/40‑1 (MQ=255)
                               gCGAGTCCGGTTCCGGGAAAAGTACGACGGGACTGGCgct                                                           <  1:3349550/40‑1 (MQ=255)
                               gCGAGTCCGGTTCCGGGAAAAGTACGACGGGACTGGCgct                                                           <  1:3084938/40‑1 (MQ=255)
                               gCGAGTCCGGTTCCGGGAAAAGTACGACGGGACTGGCGCtgctgc                                                      <  1:1620873/45‑1 (MQ=255)
                               gCGAGTCCGGTTCCGGGAAAAGTACGACGGGACTGGCGCTGCTGCGACTGATTAATTCTCAGGGCAGcatc                            <  1:1578167/71‑1 (MQ=255)
                               gCGAGTCCGGTTCCGGGAAAAGTACGACGGGACTGGAGCTGCTGCGACTGATTAATTCTCAGGGCAGcatc                            <  1:2506846/71‑1 (MQ=255)
                                   gTCCGGTTCCGGGAAAAGTACGACGGGACTGACGCTGCTGCGACTGATTAATtc                                         <  1:2506459/54‑1 (MQ=255)
                                   gTCCGGTTCCGGGAAAAGTACGACGGGACTGACGCTGCTGCGACTGATTAATtc                                         <  1:3841280/54‑1 (MQ=255)
                                   gTCCGGTTCCGGGAAAAGTACGACGGGACTGACGCTGCTGCGACTGATTAATtc                                         <  1:2973838/54‑1 (MQ=255)
                                   gTCCGGTTCCGGGAAAAGTACGACGGGACTGACGCTGCTGCGACTGATTAATtc                                         <  1:3546493/54‑1 (MQ=255)
                                    tCCGGTTCCGGGAAAAGTACGACGGGACTGGCGCTGCTGCGACTGATTAATTCTCAGGGCAGCATCATCtt                       >  1:1212143/1‑71 (MQ=255)
                                    tCCGGTTCCGGGAAAAGTACGACGGGACTGGCGCTGCTGCGACTGATTAATTCTCAGGGCAGCATCATCtt                       >  1:3445884/1‑71 (MQ=255)
                                      cGGTTCCGGGAAAAGTACGACGGGACTGACGCTGCTGCGACTGATTAATTCTCAGGGCAGCATCATCttt                      <  1:139627/70‑1 (MQ=255)
                                          tCCGGGAAAAGTACGACGGGACTGGCGCTGCTGCGACTGATTAATTCTCAGGGCAGCATCATCTTTGACgg                 >  1:2230134/1‑71 (MQ=255)
                                          tCCGGGAAAAGTACGACGGGACTGGCGCTGCTGCGACTGATTAATTCTCAGGGCAGCATCATCTTTGACg                  >  1:4030512/1‑70 (MQ=255)
                                              ggAAAAGTACGACGGGACTGACGCTGCTGCGACTGATTAATTCTCAGGGCAGCATCATCTTTGACGGTCa              >  1:2086416/1‑70 (MQ=255)
                                              ggAAAAGTACGACGGGACTGACGCTGCTGCGACTGATTAATTCTCAGGGCAGCATCATCTTTGACGGTCAg             >  1:107640/1‑71 (MQ=255)
                                               gAAAAGTACGACGGGACTGGCGCTGCTGCGACTGATTAATTCTCAGGGCAGCATCATCTTTGACGGTCAGc            <  1:238480/71‑1 (MQ=255)
                                                       cgacgGGACTGACGCTGCTGCGACTGATTAATTCTCAGGGCAGCATCATCTTTGACGGTCAGCCACTGCa     >  1:3158738/1‑70 (MQ=255)
                                                          cgGGACTGGCGCTGCTGCGACTGATTAATTGTCAGGGCAGCATCATCTTTGACGGTCAGCCACTGCaaaa  <  1:2278862/70‑1 (MQ=255)
                                                          cgGGACTGGCGCTGCTGCGACTGATTAATTCTCAGGGCAGCATCATCTTTGACGGTCAGCCACTGCaaaa  <  1:2640787/70‑1 (MQ=255)
                                                              aCTGGCGCTGCTGCGACTGATTAATTCTCAGGgcag                                >  1:1448079/1‑36 (MQ=255)
                                                               cTGACGCTGCTGCGACTGATTAATTCTCAGGGCAGCATCATCTTTGa                    <  1:333818/47‑1 (MQ=255)
                                                               cTGACGCTGCTGCGACTGATTAATTCTCAGGGCAGCATCATCTTTGa                    <  1:4087993/47‑1 (MQ=255)
                                                               cTGACGCTGCTGCGACTGATTAATTCTCAGGGCAGCATCATCTTTGa                    <  1:1835315/47‑1 (MQ=255)
                                                               cTGACGCTGCTGCGACTGATTAATTCTCAGGGCAGCATCATCTTTGa                    <  1:63777/47‑1 (MQ=255)
                                                                  |                                                             
CTACGAGCGGGTGAAACACTGGGTTTAGTGGGCGAGTCCGGTTCCGGGAAAAGTACGACGGGACTGGCGCTGCTGCGACTGATTAATTCTCAGGGCAGCATCATCTTTGACGGTCAGCCACTGCAAAA  >  minE/1408656‑1408783

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: