Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 25,217 G→A 52.3% E83K (GAA→AAA)  rihC → ribonucleoside hydrolase 3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE25,2170GA52.3% ‑1.2 / 56.3 65E83K (GAA→AAA) rihCribonucleoside hydrolase 3
Reads supporting (aligned to +/- strand):  ref base G (9/22);  new base A (17/17);  total (26/39)
Fisher's exact test for biased strand distribution p-value = 1.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.14e-01

CCAAGGGGCCGCTGTGCCACTGGTACGCGCACCGCGTGATGCGGCATCTGTGCACGGCGAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCGCTCGGGATACCGGCGT  >  minE/25159‑25280
                                                          |                                                               
ccACGGGGCCGCTGTGCCACTGGTACGCGCACCGCGTGATGCGGCATCTGTGCACGGCGAATCGGGAAt                                                       <  1:3187982/69‑1 (MQ=255)
ccAAGGGGCCGCTGTGCCACTGGTACGCGCACCGCGTGATGCGGCATCTGTGCACGGCGAATCGGGAAt                                                       <  1:1005002/69‑1 (MQ=255)
ccAAGGGGCCGCTGTGCCACTGGTACGCGCACCGCGTGATGCGGCATCTGTGCACGGCGAATCGGGAAt                                                       <  1:2595168/69‑1 (MQ=255)
ccAAGGGGCCGCTGTGCCACTGGTACGCGCACCGCGTGATGCGGCATCTGTGCACGGCGAATCGGGAAt                                                       <  1:1698799/69‑1 (MQ=255)
      ggCCGCTGTGCCACTGGTACGCGCACCGCGTGATGCGGCATCTGTGCACGGCGAATCGGGAAtggctg                                                  >  1:862950/1‑68 (MQ=255)
      ggCCGCTGTGCCACTGGTACGCGCACCGCGTGATGCGGCATCTGTGCACGGCGAATCGGGAAtggctg                                                  >  1:3231207/1‑68 (MQ=255)
                    tGGTACGCGCACCGCGTGATGCGGCATCTGTGCACGGCAAATCGGGAATGGCTGGCTACGACTttgtt                                    <  1:180640/68‑1 (MQ=255)
                    tGGTACGCGCACCGCGTGATGCGGCATCTGTGCACGGCAAATCGGGAATGGCTGGCTACGACTttgtt                                    <  1:1046992/68‑1 (MQ=255)
                    tGGTACGCGCACCGCGTGATGCGGCATCTGTGCACGGCAAATCGGGAATGGCTGGCTACGACTttgtt                                    <  1:2309003/68‑1 (MQ=255)
                    tGGTACGCGCACCGCGTGATGCGGCATCTGTGCACGGCAAATCGGGAATGGCTGGCTACGACTttgtt                                    <  1:2661859/68‑1 (MQ=255)
                    tGGTACGCGCACCGCGTGATGCGGCATCTGTGCACGGCAAATCGGGAATGGCTGGCTACGACTttgtt                                    <  1:2751259/68‑1 (MQ=255)
                    tGGTACGCGCACCGCGTGATGCGGCATCTGTGCACGGCAAATCGGGAATGGCTGGCTACGACTttgtt                                    <  1:2803350/68‑1 (MQ=255)
                    tGGTACGCGCACCGCGTGATGCGGCATCTGTGCACGGCAAATCGGGAATGGCTGGCTACGACTttgtt                                    <  1:2966723/68‑1 (MQ=255)
                    tGGTACGCGCACCGCGTGATGCGGCATCTGTGCACGGCAAATCGGGAATGGCTGGCTACGACTttgtt                                    <  1:2977338/68‑1 (MQ=255)
                    tGGTACGCGCACCGCGTGATGCGGCATCTGTGCACGGCAAATCGGGAATGGCTGGCTACGACTttgtt                                    <  1:3019928/68‑1 (MQ=255)
                    tGGTACGCGCACCGCGTGATGCGGCATCTGTGCACGGCAAATCGGGAATGGCTGGCTACGACTttgtt                                    <  1:627330/68‑1 (MQ=255)
                    tGGTACGCGCACCGCGTGATGCGGCATCTGTGCACGGCAAATCGGGAATGGCTGGCTACGACTttgtt                                    <  1:761745/68‑1 (MQ=255)
                    tGGTACGCGCACCGCGTGATGCGGCATCTGTGCACGGCAAATCGGGAATGGCTGGCTACGACTttgtt                                    <  1:1787128/68‑1 (MQ=255)
                    tGGTACGCGCACCGCGTGATGCGGCATCTGTGCACGGCAAATCGGGAATGGCTGGCTACGACTttgtt                                    <  1:1538995/68‑1 (MQ=255)
                    tGGTACGCGCACCGCGTGATGCGGCATCTGTGCACGGCAAATCGGGAATGGCTGGCTACGACTttgtt                                    <  1:1640758/68‑1 (MQ=255)
                    tGGTACGCGCACCGCGTGATGCGGCATCTGTGCACGGCAAATCGGGAATGGCTGGCTACGACTttctt                                    <  1:1462106/68‑1 (MQ=255)
                    tGGTACGCGCACCGCGTGATGCGGCATCTGTGCACGGCAAATCAGGAATGGCTGGCTACGACTttgtt                                    <  1:3051650/68‑1 (MQ=255)
                     ggTACGCGCACCGCGTGATGCGGCATCTGTGCACGGCGAATCGGGAATGGCTGGCTACGACTttgttg                                   >  1:3385366/1‑68 (MQ=255)
                     ggTACGCGCACCGCGTGATGCGGCATCTGTGCACGGCGAATCGGGAATGGCTGGCTACGACTttgttg                                   >  1:1508272/1‑68 (MQ=255)
                                    tGATGCGGCATCTGTGCACGGCAAATCGGGAATGGCTGGCTACGACTttgtt                                    <  1:2451733/52‑1 (MQ=255)
                                       tGCGGCATCTGTGCACGGCAAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCGc                 >  1:3324432/1‑68 (MQ=255)
                                       tGCGGCATCTGTGCACGGCAAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCGc                 >  1:3262355/1‑68 (MQ=255)
                                       tGCGGCATCTGTGCACGGCAAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCGc                 >  1:623781/1‑68 (MQ=255)
                                       tGCGGCATCTGTGCACGGCAAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCGc                 >  1:3104016/1‑68 (MQ=255)
                                       tGCGGCATCTGTGCACGGCAAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCGc                 >  1:697446/1‑68 (MQ=255)
                                       tGCGGCATCTGTGCACGGCAAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCGc                 >  1:702972/1‑68 (MQ=255)
                                       tGCGGCATCTGTGCACGGCAAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCGc                 >  1:2865523/1‑68 (MQ=255)
                                       tGCGGCATCTGTGCACGGCAAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCGc                 >  1:2829172/1‑68 (MQ=255)
                                       tGCGGCATCTGTGCACGGCAAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCGc                 >  1:1194181/1‑68 (MQ=255)
                                       tGCGGCATCTGTGCACGGCAAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCGc                 >  1:2642596/1‑68 (MQ=255)
                                       tGCGGCATCTGTGCACGGCAAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCGc                 >  1:2052222/1‑68 (MQ=255)
                                       tGCGGCATCTGTGCACGGCAAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCGc                 >  1:1803872/1‑68 (MQ=255)
                                       tGCGGCATCTGTGCACGGCAAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCGc                 >  1:2418380/1‑68 (MQ=255)
                                       tGCGGCATCTGTGCACGGCAAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCGc                 >  1:2066310/1‑68 (MQ=255)
                                       tGCGGCATCTGTGCACGGCAAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCGc                 >  1:2347838/1‑68 (MQ=255)
                                       tGCGGCATCTGTGCACGGCAAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCGc                 >  1:225721/1‑68 (MQ=255)
                                       tGCGGCATCTGTGCACGGCAAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCGc                 >  1:268086/1‑68 (MQ=255)
                                           gCATCTGTGCACGGCGAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCGCTCgg             <  1:1304067/68‑1 (MQ=255)
                                              tCTGTGCACGGCGAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCGc                 >  1:898333/1‑61 (MQ=255)
                                              tCTGTGCACGGCGAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCGCTCGGGa           >  1:2276542/1‑67 (MQ=255)
                                              tCTGTGCACGGCGAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCGCTCGGGa           >  1:2960549/1‑67 (MQ=255)
                                                tgtgCACGGCGAATCGGGAATGGCTGGCTACGACTttg                                      <  1:1270908/38‑1 (MQ=255)
                                                tgtgCACGGCGAATCGGGAATGGCTGGCTACGACTttg                                      <  1:66541/38‑1 (MQ=255)
                                                tgtgCACGGCGAATCGGGAATGGCTGGCTACGACTttg                                      <  1:2442067/38‑1 (MQ=255)
                                                tgtgCACGGCGAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCGCTCGGGATAcc       >  1:1766931/1‑69 (MQ=255)
                                                tgtgCACGGCGAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCGCTCGGGATAc        >  1:2344162/1‑68 (MQ=255)
                                                 gtgCACGGCGAATCGGGAATGGCTGGCTACGACTttg                                      <  1:3403412/37‑1 (MQ=255)
                                                      cGGCGAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAACGCCGCTCGGGATACCGGCGt  <  1:1900193/68‑1 (MQ=255)
                                                      cGGCGAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCGCTCGGGATACCGGCGt  <  1:3296369/68‑1 (MQ=255)
                                                      cGGCGAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCGCTCGGGATACCGGCGt  <  1:3211712/68‑1 (MQ=255)
                                                      cGGCGAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCGCTCGGGATACCGGCGt  <  1:3202590/68‑1 (MQ=255)
                                                      cGGCGAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCGCTCGGGATACCGGCGt  <  1:3167589/68‑1 (MQ=255)
                                                      cGGCGAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCGCTCGGGATACCGGCGt  <  1:2977980/68‑1 (MQ=255)
                                                      cGGCGAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCGCTCGGGATACCGGCGt  <  1:1523469/68‑1 (MQ=255)
                                                      cGGCGAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCGCTCGGGATACCGGCGt  <  1:694321/68‑1 (MQ=255)
                                                      cGGCGAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCGCTCGGGATACCGGCGt  <  1:1534256/68‑1 (MQ=255)
                                                      cGGCGAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCGCTCGGGATACCGGCGt  <  1:288640/68‑1 (MQ=255)
                                                      cGGCGAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCGCTCGGGATACCGGCGt  <  1:1205260/68‑1 (MQ=255)
                                                      cGGCGAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCGCTCGGGATACCGGCGt  <  1:1719236/68‑1 (MQ=255)
                                                      cGGCGAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCGCGCGGGATACCGGCGt  <  1:2289979/68‑1 (MQ=255)
                                                          |                                                               
CCAAGGGGCCGCTGTGCCACTGGTACGCGCACCGCGTGATGCGGCATCTGTGCACGGCGAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCGCTCGGGATACCGGCGT  >  minE/25159‑25280

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