Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,162,897 G→A 24.2% V133V (GTG→GTA ydjZ → conserved inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,162,8970GA24.2% 44.1 / 8.9 37V133V (GTG→GTAydjZconserved inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (11/17);  new base A (3/6);  total (14/23)
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

CTTTTTTATCGCCAGAGTGATGGGCCGCGAAGTGGTGGAAAAATTAACCGGCAAAACCGTGCTTGACAGTATGGACGGCTTTTTCACTCGCTACGGCAAACACACCATTCTGGTCTGTCGGT  >  minE/1162837‑1162958
                                                            |                                                             
cTTTTTTATCGCCAGAGTGATGGGCCGCGAAGTGGTGGAAAAATTAACCGGCAAAACCGTGCTTGACAg                                                       >  1:1475160/1‑69 (MQ=255)
     ttATCGCCAGAGTGATGGGCCGCGAAGTGGTGGAAAAATTAACCGGCAAAACCGTACTTGACAGTAtt                                                   <  1:3295404/68‑2 (MQ=255)
     ttATCGCCAGAGTGATGGGCCGCGAAGTGGTGGAAAAATTAACCGGCAAAACCGTACTTGACAGTATg                                                   <  1:1218658/68‑1 (MQ=255)
     ttATCGCCAGAGTGATGGGCCGCGAAGTGGTGGAAAAATTAACCGGCAAAACCGTACTTGACAGTATg                                                   <  1:2662806/68‑1 (MQ=255)
     ttATCGCCAGAGTGATGGGCCGCGAAGTGGTGGAAAAATTAACCGGCAAAACCGTACTTGACAGTATg                                                   <  1:169847/68‑1 (MQ=255)
     ttATCGCCAGAGTGATGGGCCGCGAAGTGGTGGAAAAATTAACCGGCAAAACCGTACTTGACAGTATg                                                   <  1:1509387/68‑1 (MQ=255)
     ttATCGCCAGAGTGATGGGCCGCGAAGTGGTGGAAAAATTAACCGGCAAAACCGTACTTGACAGTATg                                                   <  1:371712/68‑1 (MQ=255)
                gTGATGTGCCGCGAAGTGGTGGAAAAATTAACCGCAAAACCGTGCTTGACAGTATGGACGGCttttt                                        <  1:1652832/67‑1 (MQ=255)
                gTGATGGGCCGCGAAGTGGTGGAAAAATTAACCGGCAAAACCGTGCTTGACAGTa                                                     >  1:201751/1‑55 (MQ=255)
                gTGATGGGCCGCGAAGTGGTGGAAAAATTAACCGGCAAAACCGTGCTTGACAGTATg                                                   >  1:1492527/1‑57 (MQ=255)
                gTGATGGGCCGCGAAGTGGTGGAAAAATTAACCGGCAAAACCGTGCTTGACAGTATGGACGGCttttt                                        <  1:3336720/68‑1 (MQ=255)
                gTGATGGGCCGCGAAGTGGTGGAAAAATTAACCGGCAAAACCGTGCTTGACAGTATGGACGGCttttt                                        <  1:560335/68‑1 (MQ=255)
                gTGATGGGCCGCGAAGTGGTGGAAAAATTAACCGGCAAAACCGTGCTTGACAGTATGGACGGCttttt                                        <  1:755047/68‑1 (MQ=255)
                gTGATGGGCCGCGAAGTGGTGGAAAAATTAACCGGCAAAACCGTGCTTGACAGTATGGACGGCttttt                                        <  1:931122/68‑1 (MQ=255)
                gTGATGGGCCGCGAAGTGGTGGAAAAATTAACCGGCAAAACCGTGCTTGACAGTATGGACGGCttttt                                        <  1:2834130/68‑1 (MQ=255)
                gTGATGGGCCGCGAAGTGGTGGAAAAATTAACCGGCAAAACCGTGCTTGACAGTATGGACGGCttttt                                        <  1:281987/68‑1 (MQ=255)
                gTGATGGGCCGCGAAGTGGTGGAAAAATTAACCGGCAAAACCGTGCTTGACAGTATGGACGGCttttt                                        <  1:2763205/68‑1 (MQ=255)
                gTGATGGGCCGCGAAGTGGTGGAAAAATTAACCGGCAAAACCGTGCTTGACAGTATGGACGGCttttt                                        <  1:12740/68‑1 (MQ=255)
                gTGATGGGCCGCGAAGTGGTGGAAAAATTAACCGGCAAAACCGTGCTTGACAGTATGGACGGCttttt                                        >  1:2592435/1‑68 (MQ=255)
                gTGATGGGCCGCGAAGTGGTGGAAAAATTAACCGGCAAAACCGTGCTTGACAGTATGGACGGCttttt                                        >  1:1440000/1‑68 (MQ=255)
                gTGATGGGCCGCGAAGTGGTGGAAAAATTAACCGGCAAAACCGTGCTTGACAGTATGGACGGCttttt                                        <  1:1922304/68‑1 (MQ=255)
                gTGATGGGCCGCGAAGTGGTGGAAAAATTAACCGGCAAAACCGTGCTTGACAGTATGGACGGCttttt                                        <  1:1938157/68‑1 (MQ=255)
                gTGATGGGCCGCGAAGTGGTGGAAAAATTAACCGGCAAAACCGTGCTTGACAGTATGGACGGCttttt                                        <  1:2181323/68‑1 (MQ=255)
                gTGATAGGCCGCGAAGTGGTGGAAAAATTAACCGGCAAAACCGTGCTTGACAGTATGGACGGCttttt                                        <  1:2479553/68‑1 (MQ=255)
                       gCCGCGAAGTGGTGGAAAAATTAACCGGCAAAACCGTGCTTGACAGTATGGACGGCTTTTTCACTCGCt                                >  1:1293461/1‑69 (MQ=255)
                       gCCGCGAAGTGGTGGAAAAATTAACCGGCAAAACCGTGCTTGACAGTATGGACGGCTTTTTCACTCGCt                                >  1:3279035/1‑69 (MQ=255)
                          gcgAAGTGGTGGAAAAATTAACCGGCAAAACCGTGCTTGACAGTATGGACGGCTTTTTCACTCGCTACg                             <  1:2592392/69‑1 (MQ=255)
                                  gtggAAAAATTAACCGGCAAAACCGTGCTTGACAGTATGGACGGCTTTTTCACTCGCTACGGCAAacac                     >  1:1362576/1‑69 (MQ=255)
                                  gtggAAAAATTAACCGGCAAAACCGTACTTGACAGTATGGACGGCTTTTTCACTCGTTACGGCAAacac                     >  1:429153/1‑69 (MQ=255)
                                  gtggAAAAATTAACCGGCAAAACCGTACTTGACAGTATGGACGGCTTTTTCACTCGTTACGGCAAacac                     >  1:1331707/1‑69 (MQ=255)
                                  gtggAAAAATTAACCGGCAAAACCGTACTTGACAGTATGGACGGCTTTTTCACTCGTTACGGCAAacac                     >  1:830142/1‑69 (MQ=255)
                                                    aaaaCCGTGCTTGACAGTATGGACGGCTTTTTCACTCGCTACGGCAAACACAcc                  >  1:2314380/1‑54 (MQ=255)
                                                    aaaaCCGTGCTTGACAGTATGGACGGCTTTTTCACTCGCTACGGCAAACACACCATTCTGGTCTGTCg    >  1:2299559/1‑68 (MQ=255)
                                                    aaaaCCGTGCTTGACAGTATGGACGGCTTTTTCACTCGCTACGACAAACACACCATTCTGGTCTGTCg    >  1:2741583/1‑68 (MQ=255)
                                                     aaaCCGTGCTTGACAGTATGGACGGCTTTTTCACTCGCTACGGCAAACACACCATTCTGGTCTGTCGGt  <  1:1849478/69‑1 (MQ=255)
                                                     aaaCCGTGCTTGACAGTATGGACGGCTTTTTCACTCGCTACGGCAAACACACCATTCTGGTCTGTCGGt  <  1:1328898/69‑1 (MQ=255)
                                                      aaCCGTGCTTGACAGTATGGACGGCTTTTTCACTCGCTACGGCAAACACACCATTCTGGTCTGTCGGt  <  1:2364919/68‑1 (MQ=255)
                                                            |                                                             
CTTTTTTATCGCCAGAGTGATGGGCCGCGAAGTGGTGGAAAAATTAACCGGCAAAACCGTGCTTGACAGTATGGACGGCTTTTTCACTCGCTACGGCAAACACACCATTCTGGTCTGTCGGT  >  minE/1162837‑1162958

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: