Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,816,284 T→C 10.0% K217K (AAA→AAG 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,2840TC10.0% 175.8 / 9.5 60K217K (AAA→AAGrecADNA strand exchange and recombination protein with protease and nuclease activity
Reads supporting (aligned to +/- strand):  ref base T (26/28);  new base C (3/3);  total (29/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.34e-01

ACCCACCACGTTTTCGCCCTCTTTCACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATA  >  minE/1816215‑1816344
                                                                     |                                                            
aCCCACCACGTTTTCGCCCTCTTTCACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAAtt                                                             >  1:3090448/1‑71 (MQ=255)
 cccACCACGTTTTCGCCCTCTTTCACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAAttt                                                            <  1:2371690/71‑1 (MQ=255)
   caccacGTTTTCGCCCTCTTTCACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAAttt                                                            >  1:2306463/1‑69 (MQ=255)
   caccacGTTTTCGCCCTCTTTCACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAAttt                                                            >  1:397485/1‑69 (MQ=255)
   caccacGTTTTCACCCTCTTTCACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAAttt                                                            >  1:1908758/1‑69 (MQ=255)
   caccacGTTTTCACCCTCTTTCACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAAtt                                                             >  1:1912736/1‑68 (MQ=255)
    accacGTTTTCGCCCTCTTTCACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAg                                                         >  1:1110647/1‑71 (MQ=255)
     ccacGTTTTCGCCCTCTTTCACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTc                                                           <  1:282270/68‑1 (MQ=255)
      cacGTTTTCGCCCTCTTTCACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAgcg                                                       <  1:3342749/71‑1 (MQ=255)
      cacGTTTTCGCCCTCTTTCACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAgcg                                                       <  1:2982892/71‑1 (MQ=255)
       acGTTTTCGCCCTCTTTCACCGCGCGGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAgcgc                                                      >  1:1667655/1‑71 (MQ=255)
        cGTTTTCGCCCTCTTTCACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGcgcg                                                     <  1:1747289/71‑1 (MQ=255)
        cGTTTTCGCCCTCTTTCACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGcgcg                                                     <  1:1942852/71‑1 (MQ=255)
        cGTTTTCGCCCTCTTTCACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGcgcg                                                     <  1:2930665/71‑1 (MQ=255)
        cGTTTTCGCCCTCTTTCACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGcgcg                                                     <  1:1524326/71‑1 (MQ=255)
         gTTTTCGCCCTCTTTCACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTc                                                           >  1:1281405/1‑64 (MQ=255)
         gTTTTCGCCCTCTTTCACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTc                                                           >  1:1511858/1‑64 (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:1859678/1‑71 (MQ=255)
                cccTCTTTCACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTaccacc                                             >  1:571632/1‑71 (MQ=255)
                      ttCACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGt                                       >  1:1852267/1‑71 (MQ=255)
                        cACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCa                                                          >  1:2861118/1‑50 (MQ=255)
                        cACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTgg                                     >  1:2309427/1‑71 (MQ=255)
                         aCCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGt                                    <  1:2613020/71‑1 (MQ=255)
                           cgcgCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTac                                                 >  1:598213/1‑56 (MQ=255)
                            gcgcCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAAttt                                                            >  1:2835945/1‑44 (MQ=255)
                            gcgcCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAAttt                                                            >  1:801481/1‑44 (MQ=255)
                            gcgcCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAAttt                                                            >  1:944781/1‑44 (MQ=255)
                               ccGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCgg                              <  1:3023153/71‑1 (MQ=255)
                               ccGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCgg                              <  1:1646926/71‑1 (MQ=255)
                                      gacgGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCg                       <  1:3101420/71‑1 (MQ=255)
                                      gacgGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCg                       <  1:749540/71‑1 (MQ=255)
                                      gacgGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCg                       <  1:695987/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:116627/50‑1 (MQ=255)
                                       acgGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCgg                                           <  1:692674/50‑1 (MQ=255)
                                       acgGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCgg                                           <  1:3015873/50‑1 (MQ=255)
                                       acgGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCgg                                           <  1:198053/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:1114898/1‑71 (MQ=255)
                                             gTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATca                >  1:2874121/1‑71 (MQ=255)
                                                  gaCGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACAcc            >  1:2948023/1‑70 (MQ=255)
                                                  gaCGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACAcc            >  1:1993333/1‑70 (MQ=255)
                                                      aaCAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAAttt       <  1:149718/71‑1 (MQ=255)
                                                      aaCAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAAttt       <  1:71285/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:1430072/1‑70 (MQ=255)
                                                          gagGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCa    >  1:1184975/1‑70 (MQ=255)
                                                          gagGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCa    >  1:2545682/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:2716624/1‑70 (MQ=255)
                                                           agGCGTAGAACTTCAGTGCGTTACCGCCGGTGGTGGTTTCCGGGTTGCCGAACATGACGCCGATCTTCAtg  >  1:1114137/1‑70 (MQ=255)
                                                                   aaTTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTAcc                        >  1:2488896/1‑41 (MQ=255)
                                                                   aaTTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACAt                  <  1:1817197/47‑1 (MQ=255)
                                                                   aaTTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACAt                  <  1:640349/47‑1 (MQ=255)
                                                                   aaTTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACAt                  <  1:2765711/47‑1 (MQ=255)
                                                                     |                                                            
ACCCACCACGTTTTCGCCCTCTTTCACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATA  >  minE/1816215‑1816344

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: