Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,620,805 A→G 31.6% intergenic (‑67/+7) zntA ← / ← yhhN zinc, cobalt and lead efflux system/conserved inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,620,8050AG31.6% 20.7 / 11.2 19intergenic (‑67/+7)zntA/yhhNzinc, cobalt and lead efflux system/conserved inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (8/5);  new base G (3/3);  total (11/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTTTTCTCATTAACCGAAGGATACACTCTGGAGTCGACTCCAGAGTCAAGTTTTATCAGAGATACAGCGAGCGGACGATCAGGAAGTGACCGGCAAAGTAGCAGG  >  minE/2620756‑2620861
                                                 |                                                        
ttttttCTCATTAACCGAAGGATACACTCTGGAGTCGACTCCAGAGTCAAGTTTTATCagag                                              <  1:57050/62‑1 (MQ=255)
   tttCTCATTAACCGAAGGATACACTCTGGAGTCGACTCCAGAGTCAGGTTTTATCagag                                              >  1:417312/1‑59 (MQ=255)
   tttCTCATTAACCGAAGGATACACTCTGGAGTCGACTCCAGAGTCAGGTTTTATCAgaga                                             >  1:526814/1‑60 (MQ=255)
   tttCTCATTAACCGAAGGATACACTCTGGAGTCGACTCCAGAGTCAGGTTTTATCAgaga                                             >  1:513614/1‑60 (MQ=255)
           tAACCGAAGGATACACTCTGGAGTCGACTCCAGAGTCAAGTTTTATCAGAGATACagcgag                                    <  1:245552/61‑1 (MQ=255)
           tAACCGAAGGATACACTCTGGAGTCGACTCCAGAGTCAAGTTTTATCAGAGATACagcgag                                    <  1:279634/61‑1 (MQ=255)
                  aGGATACACTCTGGAGTCGACTCCAGAGTCAAGTTTTATCag                                                >  1:92417/1‑42 (MQ=255)
                  aGGATACACTCTGGAGTCGACTCCAGAGTCAAGTTTTATCag                                                >  1:231789/1‑42 (MQ=255)
                  aGGATACACTCTGGAGTCGACTCCAGAGTCAAGTTTTATCag                                                >  1:271956/1‑42 (MQ=255)
                  aGGATACACTCTGGAGTCGACTCCAGAGTCAAGTTTTATCag                                                >  1:273385/1‑42 (MQ=255)
                  aGGATACACTCTGGAGTCGACTCCAGAGTCAAGTTTTATCag                                                >  1:478752/1‑42 (MQ=255)
                      tACACTCTGGAGTCGACTCCAGAGTCAAGTTTTATCAGAGATACAGCGAGCGGACg                              <  1:306360/56‑1 (MQ=255)
                           tctGGAGTCGACTCCAGAGTCAAGTTTTATCAGAGATACAGCGAGCGGACGATCAGGAAGTg                   >  1:278606/1‑62 (MQ=255)
                           tctGGAGTCGACTCCAGAGTCAAGTTTTATCAGAGATACAGCGACCGGACGATCAg                         >  1:447777/1‑56 (MQ=255)
                                   cGACTCCAGAGTCAAGTTTTATCAGAGCTACAGCGAGCGGACGATCAGGAAGTGACCGGc             >  1:256310/1‑60 (MQ=255)
                                           gagTCAAGTTTTATCAGAGATACAGCGAGCGGACg                              <  1:189641/35‑1 (MQ=255)
                                            agTCAGGTTTTATCAGAGATACAGCGAGCGGACGATCAGGAAGTGACCGGCAAAGTAGCAgg  <  1:41732/62‑1 (MQ=255)
                                            agTCAGGTTTTATCAGAGATACAGCGAGCGGACGATCAGGAAGTGACCGGCAAAGTAGCAgg  <  1:33824/62‑1 (MQ=255)
                                            agTCAGGTTTTATCAGAGATACAGCGAGCGGACGATCAGGAAGTGACCGGCAAAGTAGCAgg  <  1:200259/62‑1 (MQ=255)
                                                 |                                                        
TTTTTTCTCATTAACCGAAGGATACACTCTGGAGTCGACTCCAGAGTCAAGTTTTATCAGAGATACAGCGAGCGGACGATCAGGAAGTGACCGGCAAAGTAGCAGG  >  minE/2620756‑2620861

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: