Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,784,666 T→C 43.9% A62A (GCT→GCC recN → recombination and repair protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,784,6660TC43.9% 4.6 / 30.0 41A62A (GCT→GCCrecNrecombination and repair protein
Reads supporting (aligned to +/- strand):  ref base T (8/15);  new base C (12/6);  total (20/21)
Fisher's exact test for biased strand distribution p-value = 6.16e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.54e-01

CTTTGTCTCGGTGGTCGCGCTGAAGCCGACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGTTTTTCTCTGAAAGATACGCCAGCGGCTCTGCGCTGGCTGGAAGAAAACCAG  >  minE/1784610‑1784732
                                                        |                                                                  
cTTTGTCTCGGTGGTCGCGCTGAAGCCGACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGcc                                                        >  1:1210108/1‑69 (MQ=255)
   tGTCTCGGTGGTCGCGCTGAAGCCGACATGGTGCGTAGCGGCGCTGCTCGCGCTGACCTGTGCGCCCGt                                                     >  1:2236047/1‑69 (MQ=255)
   tGTCTCGGTGGTCGCGCTGAAGCCGACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGt                                                     >  1:3191433/1‑69 (MQ=255)
   tGTCTCGGTGGTCGCGCTGAAGCCGACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGt                                                     >  1:3115030/1‑69 (MQ=255)
   tGTCTCGGTGGTCGCGCTGAAGCCGACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGt                                                     >  1:1204141/1‑69 (MQ=255)
   tGTCTCGGTGGTCGCGCTGAAGCCGACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGt                                                     >  1:2842824/1‑69 (MQ=255)
   tGTCTCGGTGGTCGCGCTGAAGCCGACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGt                                                     >  1:787604/1‑69 (MQ=255)
      ctcGGTGGTCGCGCTGAAGCCGACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGtttt                                                  >  1:1350142/1‑69 (MQ=255)
           tggtCGCGCTGAAGCCGACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGTTTTtct                                               <  1:700271/67‑1 (MQ=255)
           tggtCGCGCTGAAGCCGACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGTTTTtct                                               <  1:2687242/67‑1 (MQ=255)
           tggtCGCGCTGAAGCCGACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGTTTTtct                                               <  1:599181/67‑1 (MQ=255)
           tggtCGCGCTGAAGCCGACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGTTTTtct                                               <  1:1254086/67‑1 (MQ=255)
                  gcTGAAGCCGACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGTTTTTCTCTGAAAGa                                       <  1:1447034/68‑1 (MQ=255)
                  gcTGAAGCCGACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGTTTTTCTCTGAAAGa                                       <  1:3106304/68‑1 (MQ=255)
                  gcTGAAGCCGACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGTTTTTCTCTGAAAGa                                       <  1:2334057/68‑1 (MQ=255)
                  gcTGAAGCCGACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGTTTTTCTCTGAAAGa                                       <  1:1117874/68‑1 (MQ=255)
                  gcTGAAGCCGACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGTTTTTCTCTGAAAGa                                       <  1:270239/68‑1 (MQ=255)
                  gcTGAAGCAGACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGTTTTTCTCTGAAAGa                                       <  1:979544/68‑1 (MQ=255)
                      aaGCCGACATGGTGCGTACCGGCGCTGCTCGTGCCGACCTGTGCGCCCGTTTTTCTCTg                                            <  1:1930049/59‑1 (MQ=255)
                      aaGCCGACATGGTGCGTACCGGCGCTGCTCGTGCCGACCTGTGCGCCCGTTTTTCTCTg                                            <  1:2061531/59‑1 (MQ=255)
                      aaGCCGACATGGTGCGTACCGGCGCTGCTCGTGCCGACCTGTGCGCCCGTTTTTCTCTg                                            <  1:809984/59‑1 (MQ=255)
                      aaGCCGACATGGTGCGTACCGGCGCTGCTCGTGCCGACCTGTGCGCCCGTTTTTCTCTg                                            <  1:558919/59‑1 (MQ=255)
                      aaGCCGACATGGTGCGTACCGGCGCTCCTCGTGCCGACCTGTGCGCCCGTTTTTCTCTg                                            <  1:991031/59‑1 (MQ=255)
                                  tGCGTACCGGCGCTGCTCGTGCCGACCTGTGCGCCCGTTTTTCTCTg                                            <  1:534887/47‑1 (MQ=39)
                                   gCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGTTTTTCTCTGAAAGATACGCCAGCGGCTCTgc                      <  1:2357288/68‑1 (MQ=255)
                                   gCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGTTTTTCTCTGAAAGATACGCCAGCGGCTCTgc                      <  1:1190964/68‑1 (MQ=255)
                                    cGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGTTTTTCTCTGAAAGATACGCCAGCGGCTCTgc                      <  1:2282749/67‑1 (MQ=255)
                                      tACCGGCGCTGCTCGTGCCGACCTGTGCGCCCGTTTTTCTCTGAAAGATACGCCAGc                              >  1:1986173/1‑57 (MQ=255)
                                      tACCGGCGCTGCTCGTGCCGACCTGTGCGCCCGTTTTTCTCTGAAAGATACGCCAGCGGCCCTgcgtt                   >  1:2630208/1‑66 (MQ=255)
                                      tACCGGCGCTGCTCGTGCCGACCTGTGCGCCCGTTTTTCTCTGAAAGATACGCCAGCGGCCCTgcgtt                   >  1:1973951/1‑66 (MQ=255)
                                      tACCGGCGCTGCTCGTGCCGACCTGTGCGCCCGTTTTTCTCTGAAAGATACGCCAGCGGCCCTgcgtt                   >  1:1157619/1‑66 (MQ=255)
                                      tACCGGCGCTGCTCGTGCCGACCTGTGCGCCCGTTTTTCTCTGAAAGATACGCCAGCGGCCCTgcgt                    >  1:1395840/1‑66 (MQ=255)
                                      tACCGGCGCTGCTCGTGCCGACCTGTGCGCCCGTTTTTCTCTGAAAGATACGCCAGCGGCCCTGCgttg                  >  1:3425969/1‑69 (MQ=255)
                                      tACCGGCGCTGCTCGTGCCGACCTGTGCGCCCGTTTTTCTCTGAAAGATACGCCAGCGGCCCTGCgttg                  >  1:401723/1‑69 (MQ=255)
                                      tACCGGCGCTGCTCGTGCCGACCTGTGCGCCCGTTTTTCTCTGAAAGATACGCCAGCGGCCCTGCgttg                  >  1:503065/1‑69 (MQ=255)
                                      tACCGGCGCTGCTCGTGCCGACCTGTGCGCCCGTTTTTCTCTGAAAGATACGCCAGCGGCCCTGCgttg                  >  1:3034394/1‑69 (MQ=255)
                                      tACCGGCGCTGCTCGTGCCGACCTGTGCGCCCGTTTTTCTCTGAAAGATACGCCAGCGGCCCTGCgttg                  >  1:298527/1‑69 (MQ=255)
                                      tACCGGCGCTGCTCGTGCCGACCTGTGCGCCCGTTTTTCTCTGAAAGATACGCCAGCGGCCCTGCgttg                  >  1:1743723/1‑69 (MQ=255)
                                      tACCGGCGCTGCTCGTGCCGACCTGTGCGCCCGTTTTTCTCTGAAAGATACGCCAGCGGCCCTGCgttg                  >  1:1500421/1‑69 (MQ=255)
                                                      gcTGACCTGTGCGCCCGTTTTTCTCTGAAAGATACGCCAGCGGCTCTGCGCTGGCTGGAAGAAAACCAg  <  1:760517/69‑1 (MQ=255)
                                                      gcTGACCTGTGCGCCCGTTTTTCTCTGAAAGATACGCCAGCGGCTCTGCGCTGGCTGGAAGAAAACCAg  <  1:1139590/69‑1 (MQ=255)
                                                        |                                                                  
CTTTGTCTCGGTGGTCGCGCTGAAGCCGACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGTTTTTCTCTGAAAGATACGCCAGCGGCTCTGCGCTGGCTGGAAGAAAACCAG  >  minE/1784610‑1784732

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: