Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,671,075 A→G 37.5% D88D (GAT→GAC yfgL ← protein assembly complex, lipoprotein component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,671,0750AG37.5% 2.9 / 5.9 16D88D (GAT→GACyfgLprotein assembly complex, lipoprotein component
Reads supporting (aligned to +/- strand):  ref base A (5/5);  new base G (2/4);  total (7/9)
Fisher's exact test for biased strand distribution p-value = 6.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.28e-01

AATGCAGGCTCTTTAGAGAACCAGCCATCTTTCTCGGCCAGGCTGACAGACCAGATTTCTTTGCCATCATCCGCATTCAGCGCTTTTACTAAACCAGCGCGGTCCGCTGCATA  >  minE/1671010‑1671122
                                                                 |                                               
aaTGCAGGCTCTTTAGAGAACCAGCCATCTTTCTCGGCCAGGCTGACAGACCAGATTTCTTTGCcatca                                              <  1:661657/69‑1 (MQ=255)
aaTGCAGGCTCTTTAGAGAACCAGCCATCTTTCTCGGCCAGGCTGACAGACCAGATTTCTTTGCcatca                                              <  1:2776689/69‑1 (MQ=255)
aaTGCAGGCTCTTTAGAGAACCAGCCATCTTTCTCGGCCAGGCTGACAGACCAGATTTCTTTGCcatca                                              <  1:2363322/69‑1 (MQ=255)
 aTGCAGGCTCTTTAGAGAACCAGCCATCTTTCTCGGCCAGGCTGACAGACCAGATTTCTTTGCcatca                                              <  1:1526458/68‑1 (MQ=255)
      ggCTCTTTAGAGAACCAGCCATCTTTCTCGGCCAGGCTGACAGACCAGATTTCTTTGCCATCATCCGCa                                        >  1:1632287/1‑69 (MQ=255)
      ggCTCTTTAGAGAACCAGCCATCTTTCTCGGCCAGGCTGACAGACCAGATTTCTTTGCCATCATCCGCa                                        >  1:2269285/1‑69 (MQ=255)
                agaACCAGCCATCTTTCTCGGCCAGGCTGACAGACCAGATTTCTTTGCCATCATCCGCATTCAGCGCt                               <  1:662928/68‑1 (MQ=255)
                agaACCAGCCATCTTTCTCGGCCAGGCTGACAGACCAGATTTCTTTGCCATCATCCGCATTCAGCGCt                               <  1:2586411/68‑1 (MQ=255)
                 gaACCAGCCATCTTTCTCGGCCAGGCTGACAGACCAGATTTCTTTGCCATCATCCGCATTCAGCGCt                               <  1:1950265/67‑1 (MQ=255)
                             tttCTCGGCCAGGCTGACAGACCAGATTTCTTTGCCATCATCCGCATTCAGCGCTTTTACTAAACCAgc                 >  1:1146074/1‑69 (MQ=255)
                             tttCTCGGCCAGGCTGACAGACCAGATTTCTTTGCCATCATCCGCATTCAGCGCTTTTACTAAACCAgc                 >  1:2844672/1‑69 (MQ=255)
                             tttCTCGGCCAGGCTGACAGACCAGATTTCTTTGCCATCATCCGCATTCAGCGCTTTTACTAAACCAgc                 >  1:1678955/1‑69 (MQ=255)
                                          cTGACAGACCAGATTTCTTTGCCATCATCCGCATTCAGCGCTTTTACTAAACCAGCGCGGTCCGCTGCa    <  1:803473/69‑1 (MQ=255)
                                           tGACAGACCAGATTTCTTTGCCGTCATCCGCATTCAGCGCTTTTACTAAACCAGCGCGGTCCGCTg      >  1:2551009/1‑66 (MQ=255)
                                           tGACAGACCAGATTTCTTTGCCGTCATCCGCATTCAGCGCTTTTACTAAACCAGCGCGGTCCGCTg      >  1:615905/1‑66 (MQ=255)
                                            gACAGACCAGATTTCTTTGCCGTCATCCGCATTCAGCGCTTTTACTAAACCAGCGCGGTCCGCTGCATa  <  1:139852/69‑1 (MQ=255)
                                            gACAGACCAGATTTCTTTGCCGTCATCCGCATTCAGCGCTTTTACTAAACCAGCGCGGTCCGCTGCATa  <  1:3050173/69‑1 (MQ=255)
                                            gACAGACCAGATTTCTTTGCCGTCATCCGCATTCAGCGCTTTTACTAAACCAGCGCGGTCCGCTGCATa  <  1:1154665/69‑1 (MQ=255)
                                            gACAGACCAGATTTCTTTGCCGTCATCCGCATTCAGCGCTTTTACTAAACCAGCGCGGTCCGCTGCATa  <  1:1106595/69‑1 (MQ=255)
                                                       tttCTTTGACATCATCCGCATTCAGCGCTTTTACTAAACCAGCGCGGTCCGCTGCATa  <  1:2026112/58‑1 (MQ=255)
                                                                 |                                               
AATGCAGGCTCTTTAGAGAACCAGCCATCTTTCTCGGCCAGGCTGACAGACCAGATTTCTTTGCCATCATCCGCATTCAGCGCTTTTACTAAACCAGCGCGGTCCGCTGCATA  >  minE/1671010‑1671122

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: