Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 531,429 A→G 24.1% V115V (GTA→GTG ybhQ → predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE531,4290AG24.1% 36.7 / 7.6 33V115V (GTA→GTGybhQpredicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (14/11);  new base G (5/3);  total (19/14)
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

GAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACC  >  minE/531373‑531497
                                                        |                                                                    
gAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTGGTGTCTTTGCTg                                                          >  1:1028949/1‑69 (MQ=255)
gAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTGGTGTCTTTGCTg                                                          >  1:974048/1‑69 (MQ=255)
gAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTGGTGTCTTTGCTg                                                          >  1:1581355/1‑69 (MQ=255)
gAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTGGTGTCTTTGCTg                                                          >  1:2454876/1‑69 (MQ=255)
gAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTGGTGTCTTTGCTg                                                          >  1:2264992/1‑69 (MQ=255)
gAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTg                                                          >  1:2630725/1‑69 (MQ=255)
gAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTg                                                          >  1:2626601/1‑69 (MQ=255)
gAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTg                                                          >  1:3392786/1‑69 (MQ=255)
gAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTg                                                          >  1:2319979/1‑69 (MQ=255)
                           gggCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCt                                <  1:2744773/68‑1 (MQ=255)
                           gggCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCt                                <  1:865814/68‑1 (MQ=255)
                           gggCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCt                                <  1:599348/68‑1 (MQ=255)
                           gggCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCt                                <  1:2114870/68‑1 (MQ=255)
                              cTGGCAATCCTTGCCGGCCCGGCGGTGGTGTCTTTGCTGGCGa                                                      <  1:2922943/43‑1 (MQ=255)
                              cTGGCAATCCTTGCCGGCCCGGCGGTGGTGTCTTTGCTGGCGa                                                      <  1:2122655/43‑1 (MQ=255)
                              cTGGCAATCCTTGCCGGCCCGGCGGTGGTGTCTTTGCTGGCGa                                                      <  1:1759278/43‑1 (MQ=255)
                                   aaTCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTa                        <  1:2100855/68‑1 (MQ=255)
                                   aaTCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTa                        <  1:1833923/68‑1 (MQ=255)
                                    aTCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTaaa                      <  1:3231876/69‑1 (MQ=255)
                                    aTCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTaaa                      <  1:2647801/69‑1 (MQ=255)
                                    aTCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTaaa                      <  1:1546683/69‑1 (MQ=255)
                                      ccTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTaaa                      >  1:2087253/1‑67 (MQ=255)
                                      ccTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACa                    >  1:541957/1‑69 (MQ=255)
                                      ccTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACa                    >  1:3408632/1‑69 (MQ=255)
                                      ccTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACa                    >  1:955852/1‑69 (MQ=255)
                                      ccTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACa                    >  1:309304/1‑69 (MQ=255)
                                      ccTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACa                    >  1:1226619/1‑69 (MQ=255)
                                      ccTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACa                    >  1:2447066/1‑69 (MQ=255)
                                      ccTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACa                    >  1:1692917/1‑69 (MQ=255)
                                      ccTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACa                    >  1:1544478/1‑69 (MQ=255)
                                      ccTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCAAAACa                    >  1:2921575/1‑69 (MQ=255)
                                                 cggcggTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTAcgcc         <  1:2636438/69‑1 (MQ=255)
                                                       tAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTaaa                      <  1:592821/50‑1 (MQ=255)
                                                        aGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGAcc  >  1:2590629/1‑69 (MQ=255)
                                                        aGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGAcc  >  1:473159/1‑69 (MQ=255)
                                                        aGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGAcc  >  1:2339506/1‑69 (MQ=255)
                                                        aGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGAcc  >  1:1739442/1‑69 (MQ=255)
                                                        aGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGAc   >  1:8621/1‑68 (MQ=255)
                                                        aGTGTCTTTGCAGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGTGTACGCCGCTGAcc  >  1:1642361/1‑69 (MQ=255)
                                                        |                                                                    
GAAAACAACTTTTTGCTGGCAATTGCAGGGCTGGCAATCCTTGCCGGCCCGGCGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACC  >  minE/531373‑531497

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