Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,784,663 C→T 46.1% R61R (CGC→CGT recN → recombination and repair protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,784,6630CT46.1% 2.2 / 29.5 39R61R (CGC→CGTrecNrecombination and repair protein
Reads supporting (aligned to +/- strand):  ref base C (8/13);  new base T (12/6);  total (20/19)
Fisher's exact test for biased strand distribution p-value = 1.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATAGATGCCCTCGGTCTTTGTCTCGGTGGTCGCGCTGAAGCCGACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGTTTTTCTCTGAAAGATACGCCAGCGGCTCTGCGCTG  >  minE/1784595‑1784716
                                                                    |                                                     
aTAGATGCCCTCGGTCTTTGTCTCGGTGGTCGCGCTGAAGCCGACATGGTGCGTACCGGCGCTGCTcgc                                                       >  1:3408462/1‑69 (MQ=255)
 tAGATGCCCTCGGTCTTTGTCTCGGTGGTCGCGCTGAAGCCGACATGGTGCGTACCGGCGCTGCTcgcg                                                      >  1:3189758/1‑69 (MQ=255)
               cTTTGTCTCGGTGGTCGCGCTGAAGCCGACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGcc                                        >  1:1210108/1‑69 (MQ=255)
                  tGTCTCGGTGGTCGCGCTGAAGCCGACATGGTGCGTAGCGGCGCTGCTCGCGCTGACCTGTGCGCCCGt                                     >  1:2236047/1‑69 (MQ=255)
                  tGTCTCGGTGGTCGCGCTGAAGCCGACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGt                                     >  1:3115030/1‑69 (MQ=255)
                  tGTCTCGGTGGTCGCGCTGAAGCCGACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGt                                     >  1:3191433/1‑69 (MQ=255)
                  tGTCTCGGTGGTCGCGCTGAAGCCGACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGt                                     >  1:787604/1‑69 (MQ=255)
                  tGTCTCGGTGGTCGCGCTGAAGCCGACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGt                                     >  1:1204141/1‑69 (MQ=255)
                  tGTCTCGGTGGTCGCGCTGAAGCCGACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGt                                     >  1:2842824/1‑69 (MQ=255)
                     ctcGGTGGTCGCGCTGAAGCCGACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGtttt                                  >  1:1350142/1‑69 (MQ=255)
                          tggtCGCGCTGAAGCCGACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGTTTTtct                               <  1:2687242/67‑1 (MQ=255)
                          tggtCGCGCTGAAGCCGACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGTTTTtct                               <  1:700271/67‑1 (MQ=255)
                          tggtCGCGCTGAAGCCGACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGTTTTtct                               <  1:1254086/67‑1 (MQ=255)
                          tggtCGCGCTGAAGCCGACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGTTTTtct                               <  1:599181/67‑1 (MQ=255)
                                 gcTGAAGCCGACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGTTTTTCTCTGAAAGa                       <  1:270239/68‑1 (MQ=255)
                                 gcTGAAGCCGACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGTTTTTCTCTGAAAGa                       <  1:2334057/68‑1 (MQ=255)
                                 gcTGAAGCCGACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGTTTTTCTCTGAAAGa                       <  1:3106304/68‑1 (MQ=255)
                                 gcTGAAGCCGACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGTTTTTCTCTGAAAGa                       <  1:1447034/68‑1 (MQ=255)
                                 gcTGAAGCCGACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGTTTTTCTCTGAAAGa                       <  1:1117874/68‑1 (MQ=255)
                                 gcTGAAGCAGACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGTTTTTCTCTGAAAGa                       <  1:979544/68‑1 (MQ=255)
                                     aaGCCGACATGGTGCGTACCGGCGCTGCTCGTGCCGACCTGTGCGCCCGTTTTTCTCTg                            <  1:2061531/59‑1 (MQ=255)
                                     aaGCCGACATGGTGCGTACCGGCGCTGCTCGTGCCGACCTGTGCGCCCGTTTTTCTCTg                            <  1:558919/59‑1 (MQ=255)
                                     aaGCCGACATGGTGCGTACCGGCGCTGCTCGTGCCGACCTGTGCGCCCGTTTTTCTCTg                            <  1:1930049/59‑1 (MQ=255)
                                     aaGCCGACATGGTGCGTACCGGCGCTGCTCGTGCCGACCTGTGCGCCCGTTTTTCTCTg                            <  1:809984/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:3034394/1‑69 (MQ=255)
                                                     tACCGGCGCTGCTCGTGCCGACCTGTGCGCCCGTTTTTCTCTGAAAGATACGCCAGCGGCCCTGCgttg  >  1:298527/1‑69 (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:1743723/1‑69 (MQ=255)
                                                     tACCGGCGCTGCTCGTGCCGACCTGTGCGCCCGTTTTTCTCTGAAAGATACGCCAGCGGCCCTGCgttg  >  1:1500421/1‑69 (MQ=255)
                                                                    |                                                     
ATAGATGCCCTCGGTCTTTGTCTCGGTGGTCGCGCTGAAGCCGACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGTTTTTCTCTGAAAGATACGCCAGCGGCTCTGCGCTG  >  minE/1784595‑1784716

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: