Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,816,290 C→T 10.0% A215A (GCG→GCA recA ← DNA strand exchange and recombination protein with protease and nuclease activity

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,816,2900CT10.0% 176.1 / 9.5 60A215A (GCG→GCArecADNA strand exchange and recombination protein with protease and nuclease activity
Reads supporting (aligned to +/- strand):  ref base C (27/27);  new base T (3/3);  total (30/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.12e-01

CACGTTTTCGCCCTCTTTCACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGAT  >  minE/1816221‑1816358
                                                                     |                                                                    
cacGTTTTCGCCCTCTTTCACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAgcg                                                                     <  1:2982892/71‑1 (MQ=255)
cacGTTTTCGCCCTCTTTCACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAgcg                                                                     <  1:3342749/71‑1 (MQ=255)
 acGTTTTCGCCCTCTTTCACCGCGCGGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAgcgc                                                                    >  1:1667655/1‑71 (MQ=255)
  cGTTTTCGCCCTCTTTCACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGcgcg                                                                   <  1:2930665/71‑1 (MQ=255)
  cGTTTTCGCCCTCTTTCACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGcgcg                                                                   <  1:1942852/71‑1 (MQ=255)
  cGTTTTCGCCCTCTTTCACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGcgcg                                                                   <  1:1747289/71‑1 (MQ=255)
  cGTTTTCGCCCTCTTTCACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGcgcg                                                                   <  1:1524326/71‑1 (MQ=255)
   gTTTTCGCCCTCTTTCACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGcgcg                                                                   >  1:2463947/1‑70 (MQ=255)
     tttCGCCCTCTTTCACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTa                                                                >  1:1539961/1‑71 (MQ=255)
      ttCGCCCTCTTTCACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTac                                                               >  1:1655465/1‑71 (MQ=255)
         gCCCTCTTTCACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTaccac                                                            <  1:647448/71‑1 (MQ=255)
          cccTCTTTCACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTaccacc                                                           >  1:1768035/1‑71 (MQ=255)
          cccTCTTTCACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTaccacc                                                           >  1:571632/1‑71 (MQ=255)
          cccTCTTTCACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTaccacc                                                           >  1:1859678/1‑71 (MQ=255)
                ttCACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGt                                                     >  1:1852267/1‑71 (MQ=255)
                  cACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTgg                                                   >  1:2309427/1‑71 (MQ=255)
                   aCCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGt                                                  <  1:2613020/71‑1 (MQ=255)
                     cgcgCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTac                                                               >  1:598213/1‑56 (MQ=255)
                         ccGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCgg                                            <  1:3023153/71‑1 (MQ=255)
                         ccGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCgg                                            <  1:1646926/71‑1 (MQ=255)
                                gacgGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCg                                     <  1:695987/71‑1 (MQ=255)
                                gacgGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCg                                     <  1:3101420/71‑1 (MQ=255)
                                gacgGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCg                                     <  1:749540/71‑1 (MQ=255)
                                gacgGATGTCCATGCGCACGGAGGCGTAGAACTTCAGTGCGtt                                                                 <  1:485267/43‑1 (MQ=255)
                                gacgGATGTCCATGCGCACGGAGGCGTAGAACTTCAGTGCGtt                                                                 <  1:1210086/43‑1 (MQ=255)
                                gacgGATGTCCATGCGCACGGAGGCGTAGAACTTCAGTGCGtt                                                                 <  1:1557634/43‑1 (MQ=255)
                                 acgGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCgg                                                         <  1:3015873/50‑1 (MQ=255)
                                 acgGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCgg                                                         <  1:198053/50‑1 (MQ=255)
                                 acgGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCgg                                                         <  1:692674/50‑1 (MQ=255)
                                 acgGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCgg                                                         <  1:116627/50‑1 (MQ=255)
                                  cgGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGaa                                   >  1:151706/1‑71 (MQ=255)
                                  cgGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGaa                                   >  1:2463824/1‑71 (MQ=255)
                                   ggATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCgg                                                         <  1:407249/48‑1 (MQ=255)
                                   ggATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACAACCGGTa                                                       >  1:2031346/1‑50 (MQ=255)
                                       gTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATca                              >  1:2874121/1‑71 (MQ=255)
                                       gTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATca                              >  1:1114898/1‑71 (MQ=255)
                                            gaCGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACAcc                          >  1:1993333/1‑70 (MQ=255)
                                            gaCGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACAcc                          >  1:2948023/1‑70 (MQ=255)
                                                aaCAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAAttt                     <  1:71285/71‑1 (MQ=255)
                                                aaCAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAAttt                     <  1:149718/71‑1 (MQ=255)
                                                 aCAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAAtttt                    <  1:369028/71‑1 (MQ=255)
                                                 aCAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAAtttt                    <  1:3294925/71‑1 (MQ=255)
                                                   agagGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCa                  <  1:1879185/71‑1 (MQ=255)
                                                    gagGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCa                  >  1:2545682/1‑70 (MQ=255)
                                                    gagGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCa                  >  1:1430072/1‑70 (MQ=255)
                                                    gagGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCa                  >  1:1184975/1‑70 (MQ=255)
                                                    gagGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCa                  >  1:104092/1‑70 (MQ=255)
                                                     agGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATa                <  1:1644776/71‑1 (MQ=255)
                                                     agGCGTAGAACTTCAGTGCGTTACCGCCGGTGGTGGTTTCCGGGTTGCCGAACATGACGCCGATCTTCAtg                >  1:510541/1‑70 (MQ=255)
                                                     agGCGTAGAACTTCAGTGCGTTACCGCCGGTGGTGGTTTCCGGGTTGCCGAACATGACGCCGATCTTCAtg                >  1:1114137/1‑70 (MQ=255)
                                                     agGCGTAGAACTTCAGTGCGTTACCGCCGGTGGTGGTTTCCGGGTTGCCGAACATGACGCCGATCTTCAtg                >  1:2716624/1‑70 (MQ=255)
                                                             aaTTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTAcc                                      >  1:2488896/1‑41 (MQ=255)
                                                             aaTTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACAt                                <  1:2765711/47‑1 (MQ=255)
                                                             aaTTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACAt                                <  1:1817197/47‑1 (MQ=255)
                                                             aaTTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACAt                                <  1:640349/47‑1 (MQ=255)
                                                                 tCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGtt     <  1:3396454/70‑1 (MQ=255)
                                                                 tCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGtt     >  1:879352/1‑70 (MQ=255)
                                                                 tCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGtt     <  1:1354823/70‑1 (MQ=255)
                                                                 tCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGtt     <  1:1516998/70‑1 (MQ=255)
                                                                 tCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGtt     >  1:2300310/1‑70 (MQ=255)
                                                                 tCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGtt     <  1:398839/70‑1 (MQ=255)
                                                                 tCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACAAATTTTCATACGGAt           >  1:1756692/1‑64 (MQ=255)
                                                                  cAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAAc                                  <  1:612668/40‑1 (MQ=255)
                                                                  cAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTtg    <  1:3179918/70‑1 (MQ=255)
                                                                   aGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTtgat  >  1:2016725/1‑71 (MQ=255)
                                                                   aGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTtgat  >  1:2082225/1‑71 (MQ=255)
                                                                   aGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTtgat  >  1:1463126/1‑71 (MQ=255)
                                                                   aGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTtgat  >  1:846545/1‑71 (MQ=255)
                                                                     |                                                                    
CACGTTTTCGCCCTCTTTCACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGAT  >  minE/1816221‑1816358

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: