Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,810,979 C→T 70.3% G589G (GGG→GGA uvrA ← ATPase and DNA damage recognition protein of nucleotide excision repair excinuclease UvrABC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,810,9790CT70.3% 50.1 / 23.8 57G589G (GGG→GGAuvrAATPase and DNA damage recognition protein of nucleotide excision repair excinuclease UvrABC
Reads supporting (aligned to +/- strand):  ref base C (5/12);  new base T (20/20);  total (25/32)
Fisher's exact test for biased strand distribution p-value = 2.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.51e-01

CTTTTTCCGGATTCGCCGGAACGCGTTTCTTCGGCACTTCAATCTTGCGTTTGCCGCTCATGTACTGCCCGGTCAACGACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTGCGACCACTTC  >  minE/2810912‑2811038
                                                                   |                                                           
cTTTTTCCGGATTCGCCGGAACGCGTTTCTTCGGCACTTCAATCTTGCGTTTGCCGCTCATGTACTGcc                                                            >  1:2505904/1‑69 (MQ=255)
cTTTTTCCGGATTCGCCGGAACGCGTTTCTTCGGCACTTCAATCTTGCGTTTGCCGCTCATGTACTGcc                                                            >  1:3294641/1‑69 (MQ=255)
cTTTTTCCGGATTCGCCGGAACGCGTTTCTTCGGCACTTCAATCTTGCGTTTGCCGCTCATGTACTGcc                                                            >  1:2936369/1‑69 (MQ=255)
cTTTTTCCGGATTCGCCGGAACGCGTTTCTTCGGCACTTCAATCTTGCGTTTGCCGCTCATGTACTGcc                                                            >  1:2015231/1‑69 (MQ=255)
cTTTTTCCGGATTCGCCGGAACGCGTTTCTTCGGCACTTCAATCTTGCGTTTGCCGCTCATGTACTGcc                                                            >  1:517587/1‑69 (MQ=255)
cTTTTTCCGGATTCGCCGGAACGCGTTTCTTCGGCACCTCAATCTTGCGTTTGCCGCTCATGTACTGcc                                                            >  1:36573/1‑69 (MQ=255)
  ttttCCGGATTCGCCGGAACGCGTTTCTTCGGCACTTCAATCTTGCGTTTGCCGCTCATGTACTGccc                                                           >  1:20837/1‑68 (MQ=255)
  ttttCCGGATTCGCCGGAACGCGTTTCTTCGGCACTTCAATCTTGCGTTTGCCGCTCATGTACTGCCCg                                                          >  1:1669718/1‑69 (MQ=255)
  ttttCCGGATTCGCCGGAACGCGTTTCTTCGGCACTTCAATCTTGCGTTTGCCGCTCATGTACTGCCCg                                                          >  1:3210032/1‑69 (MQ=255)
  ttttCCAGATTCGCCGGAACGCGTTTCTTCGGCACTTCAATCTTGCGTTTGCCGCTCATGTACTGCCCg                                                          >  1:442122/1‑69 (MQ=255)
           ttCGCCGGAACGCGTTTCTTCGGCACTTCAATCTTGCGTTTGCCGCTCATGTACTGCCCGGTCAACGAc                                                 <  1:2222149/69‑1 (MQ=255)
           ttCGCCGGAACGCGTTTCTTCGGCACTTCAATCTTGCGTTTGCCGCTCATGTACTGCCCGGTCAACGAc                                                 <  1:2289775/69‑1 (MQ=255)
             cGCCGGAACGCGTTTCTTCGGCACTTCAATCTTGCGTTTACCGCTCATGTACTGTCCGGt                                                        >  1:1813050/1‑60 (MQ=255)
             cGCCGGAACGCGTTTCTTCGGCACTTCAATCTTGCGTTTACCGCTCATGTACTGTCCGGTCAACGAc                                                 >  1:1866302/1‑67 (MQ=255)
             cGCCGGAACGCGTTTCTTCGGCACTTCAATCTTGCGTTTACCGCTCATGTACTGTCCGGTCAACGACt                                                >  1:328551/1‑68 (MQ=255)
             cGCCGGAACGCGTTTCTTCGGCACTTCAATCTTGCGTTTACCGCTCATGTACTGTCCGGTCAACGACt                                                >  1:2490292/1‑68 (MQ=255)
             cGCCGGAACGCGTTTCTTCGGCACTTCAATCTTGCGTTTACCGCTCATGTACTGTCCGGTCAACGACt                                                >  1:626741/1‑68 (MQ=255)
             cGCCGGAACGCGTTTCTTCGGCACTTCAATCTTGCGTTTACCGCTCATGTACTGTCCGGTCAACGACt                                                >  1:3211667/1‑68 (MQ=255)
             cGCCGGAACGCGTTTCTTCGGCACTTCAATCTTGCGTTTACCGCTCATGTACTGTCCGGTCAACGACt                                                >  1:2633303/1‑68 (MQ=255)
             cGCCGGAACGCGTTTCTTCGGCACTTCAATCTTGCGTTTACCGCTCATGTACTGTCCGGTCAACGACt                                                >  1:2745165/1‑68 (MQ=255)
             cGCCGGAACGCGTTTCTTCGGCACTTCAATCTTGCGTTTACCGCTCATGTACTGTCCGGTCAACGACt                                                >  1:3187332/1‑68 (MQ=255)
             cGCCGGAACGCGTTTCTTCGGCACTTCAATCTTGCGTTTACCGCTCATGTACTGTCCGGTCAACGACt                                                >  1:3014532/1‑68 (MQ=255)
             cGCCGGAACGCGTTTCTTCGGCACTTCAATCTTGCGTTTACCGCTCATGTACTGTCCGGTCAACGACt                                                >  1:3023610/1‑68 (MQ=255)
             cGCCGGAACGCGTTTCTTCGGCACTTCAATCTTGCGTTTACCGCTCATGTACTGTCCGGTCAACGACt                                                >  1:309366/1‑68 (MQ=255)
             cGCCGGAACGCGTTTCTTCGGCACTTCAATCTTGCGTTTACCGCTCATGTACTGTCCGGTCAACGACt                                                >  1:2348782/1‑68 (MQ=255)
             cGCCGGAACGCGTTTCTTCGGCACTTCAATCTTGCGTTTACCGCTCATGTACTGTCCGGTCAACGACt                                                >  1:2334716/1‑68 (MQ=255)
             cGCCGGAACGCGTTTCTTCGGCACTTCAATCTTGCGTTTACCGCTCATGTACTGTCCGGTCAACGACt                                                >  1:2153196/1‑68 (MQ=255)
             cGCCGGAACGCGTTTCTTCGGCACTTCAATCTTGCGTTTACCGCTCATGTACTGTCCGGTCAACGACt                                                >  1:1837961/1‑68 (MQ=255)
             cGCCGGAACGCGTTTCTTCGGCACTTCAATCTTGCGTTTACCGCTCATGTACTGTCCGGTCAACGACt                                                >  1:170079/1‑68 (MQ=255)
             cGCCGGAACGCGTTTCTTCGGCACTTCAATCTTGCGTTTACCGCTCATGTACTGTCCGGTCAACGACt                                                >  1:1618300/1‑68 (MQ=255)
             cGCCGGAACGCGTTTCTTCGGCACTTCAATCTTGCGTTTACCGCTCATGTACTGTCCGGTCAACGACt                                                >  1:626395/1‑68 (MQ=255)
             cGCCGGAACGCGTTTCTTCGGCACTTCAATCTTGCGTTTACCGCTCATGTACTGTCCGGTCAACGACt                                                >  1:935953/1‑68 (MQ=255)
                          ttcttcGGCACTTCAATCTTGCGTTTACCGCTCATGTACTGTCCGGTCAACGACTCCGGCAcc                                        <  1:3289077/63‑1 (MQ=255)
                          ttcttcGGCACTTCAATCTTGCGTTTACCGCTCATGTACTGTCCGGTCAACGACTCCGGCAcc                                        <  1:344988/63‑1 (MQ=255)
                          ttcttcGGCACTTCAATCTTGCGTTTACCGCTCATGTACTGTCCGGTCAACGACTCCGGCAcc                                        <  1:595102/63‑1 (MQ=255)
                          ttcttcGGCACTTCAATCTTGCGTTTACCGCTCATGTACTGTCCGGTCAACGACTCCGGCAcc                                        <  1:596259/63‑1 (MQ=255)
                          ttcttcGGCACTTCAATCTTGCGTTTACCGCTCATGTACTGTCCGGTCAACGACTCCGGCAcc                                        <  1:3158556/63‑1 (MQ=255)
                          ttcttcGGCACTTCAATCTTGCGTTTACCGCTCATGTACTGTCCGGTCAACGACTCCGGCAcc                                        <  1:639900/63‑1 (MQ=255)
                          ttcttcGGCACTTCAATCTTGCGTTTACCGCTCATGTACTGTCCGGTCAACGACTCCGGCAcc                                        <  1:2926608/63‑1 (MQ=255)
                          ttcttcGGCACTTCAATCTTGCGTTTACCGCTCATGTACTGTCCGGTCAACGACTCCGGCAcc                                        <  1:2699019/63‑1 (MQ=255)
                          ttcttcGGCACTTCAATCTTGCGTTTACCGCTCATGTACTGTCCGGTCAACGACTCCGGCAcc                                        <  1:1872706/63‑1 (MQ=255)
                          ttcttcGGCACTTCAATCTTGCGTTTACCGCTCATGTACTGTCCGGTCAACGACTCCGGCAcc                                        <  1:2044899/63‑1 (MQ=255)
                          ttcttcGGCACTTCAATCTTGCGTTTACCGCTCATGTACTGTCCGGTCAACGACTCCGGCAcc                                        <  1:3140062/63‑1 (MQ=255)
                           tcttcGGCACTTCAATCTTGCGTTTACCGCTCATGTACTGTCCGGTCAACGACTCCGGCAcc                                        <  1:1278469/62‑1 (MQ=255)
                           tcttcGGCACTTCAATCTTGCGTTTACCGCTCATGTACTGTCCGGTCAACGACTCCGGCAcc                                        <  1:900880/62‑1 (MQ=255)
                           tcttcGGCACTTCAATCTTGCGTTTACCGCTCATGTACTGTCCGGTCAACGACTCCGGCAcc                                        <  1:927466/62‑1 (MQ=255)
                                      tCAATCTTGCGTTTGCCGCTCATGTACTGCCCGGTCAACGACTCCgg                                            <  1:2086104/47‑1 (MQ=255)
                                      tCAATCTTGCGTTGGCCGCTCATGTACTGCCCGGTCAACGACTCCgg                                            <  1:3306774/47‑1 (MQ=255)
                                         aTCTTGCGTTTGCCGCTCATGTACTGCCCGGTCAACGACTCCGGCACCGCCATAATCGCTTCCAGCgg                    <  1:3354421/68‑1 (MQ=255)
                                           cTTGCGTTTGCCGCTCATGTACTGCCCGGTCAACGACTc                                               >  1:814356/1‑39 (MQ=255)
                                           cTTGCGTTTGCCGCTCATGTACTGCCCGGTCAACGACTCCGGCACCGCCATAATCGCTTCCAGCGGAc                  >  1:1936213/1‑68 (MQ=255)
                                                 tttGCCGCTCATGTACTGCCCGGTCAACGACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTGc           <  1:2576438/69‑1 (MQ=255)
                                                     ccGCTCATGTACTGCCCGGTCAACGACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTGCGAc        <  1:3327521/68‑1 (MQ=255)
                                                       gCTCATGTACTGCCCGGTCAACGACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTGCGACCAc     <  1:1092708/69‑1 (MQ=255)
                                                        cTCATGTACTGCCCGGTCAACGACTCCGGCACCGCCATAATCGCTTCCAGCgg                    <  1:2504129/53‑1 (MQ=255)
                                                         tCATGTACTGCCCGGTCAACGACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTGCGACCACtt   <  1:2318147/69‑1 (MQ=255)
                                                         tCATGTACTGCCCGGTCAACGACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTGCGACCACtt   <  1:3287930/69‑1 (MQ=255)
                                                          cATGTACTGTCCGGTCAACGACTCCGGCACCGGCATAATTGCTTCCAGTGGACCTTCCGCGACCACTTc  <  1:6413/69‑1 (MQ=255)
                                                          cATGTACTGTCCGGTCAACGACTCCGGCACCGCCATAATTGCTTCCAGTGGACCTTCCGCGACCACTTc  <  1:113150/69‑1 (MQ=255)
                                                          cATGTACTGTCCGGTCAACGACTCCGGCACCGCCATAATTGCTTCCAGTGGACCTTCCGCGACCACTTc  <  1:2548668/69‑1 (MQ=255)
                                                          cATGTACTGTCCGGTCAACGACTCCGGCACCGCCATAATTGCTTCCAGTGGACCTTCCGCGACCACTTc  <  1:2266888/69‑1 (MQ=255)
                                                          cATGTACTGTCCGGTCAACGACTCCGGCACCGCCATAATTGCTTCCAGTGGACCTTCCGCGACCACTTc  <  1:2265202/69‑1 (MQ=255)
                                                          cATGTACTGTCCGGTCAACGACTCCGGCACCGCCATAATTGCTTCCAGTGGACCTTCCGCGACCACTTc  <  1:1179268/69‑1 (MQ=255)
                                                          cATGTACTGCCCGGTCAACGACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTGCGACCACtt   <  1:283982/68‑1 (MQ=255)
                                                                   |                                                           
CTTTTTCCGGATTCGCCGGAACGCGTTTCTTCGGCACTTCAATCTTGCGTTTGCCGCTCATGTACTGCCCGGTCAACGACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTGCGACCACTTC  >  minE/2810912‑2811038

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: