Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,936,148 A→G 26.6% R409R (CGT→CGC recJ ← ssDNA exonuclease, 5' ‑‑> 3'‑specific

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,936,1480AG26.6% 79.9 / 28.9 75R409R (CGT→CGCrecJssDNA exonuclease, 5' ‑‑> 3'‑specific
Reads supporting (aligned to +/- strand):  ref base A (26/29);  new base G (7/13);  total (33/42)
Fisher's exact test for biased strand distribution p-value = 4.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.56e-01

TCGCATGACCGCCAAACTTCAGCATCATGCCAGGGTAGAGTGTGTCTAATCGCTCCAGCGCATCACGCATATGCAGCCCCTGAATGGAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGC  >  minE/1936084‑1936216
                                                                |                                                                    
tCGCATGGCCGCCAAACTTGAGCATCATGCCGGGGTAAAGCGTGTCTAGCCGCTCCAGCGCATCGCGCa                                                                  <  1:869978/69‑1 (MQ=255)
tCGCATGGCCGCCAAACTTGAGCATCATGCCGGGGTAAAGCGTGTCTAGCCGCTCCAGCGCATCGCGCa                                                                  <  1:1807689/69‑1 (MQ=255)
tCGCATGGCCGCCAAACTTGAGCATCATGCCGGGGTAAAGCGTGTCTAGCCGCTCCAGCGCATCGCGCa                                                                  <  1:387738/69‑1 (MQ=255)
tCGCATGGCCGCCAAACTTGAGCATCATGCCGGGGTAAAGCGTGTCTAGCCGCTCCAGCGCATCGCGCa                                                                  <  1:1326998/69‑1 (MQ=255)
tCGCATGGCCGCCAAACTTGAGCATCATGCCGGGGTAAAGCGTGTCTAGCCGCTCCAGCGCATCGCGCa                                                                  <  1:1953554/69‑1 (MQ=255)
tCGCATGGCCGCCAAACTTGAGCATCATGCCGGGGTAAAGCGTGTCTAGCCGCTCCAGCGCATCGCGCa                                                                  <  1:3104832/69‑1 (MQ=255)
tCGCATGGCCGCCAAACTTGAGCATCATGCCGGGGTAAAGCGTGTCTAGCCGCTCCAGCGCATCGCGCa                                                                  <  1:310013/69‑1 (MQ=255)
tCGCATGGCCGCCAAACTTGAGCATCATGCCGGGGTAAAGCGTGTCTAGCCGCTCCAGCGCATCGCGCa                                                                  <  1:148750/69‑1 (MQ=255)
tCGCATGGCCGCCAAACTTGAGCATCATGCCGGGGTAAAGCGTGTCTAGCCGCTCCAGCGCATCGCGCa                                                                  <  1:1957506/69‑1 (MQ=255)
tCGCATGGCCGCCAAACTTGAGCATCATGCCGGGGTAAAGCGTGTCTAGCCGCTCCAGCGCATCGCGCa                                                                  <  1:2532447/69‑1 (MQ=255)
tCGCATGGCCGCCAAACTTGAGCATCATGCCGGGGTAAAGCGTGTCTAGCCGCTCCAGCGCATCGCGCa                                                                  <  1:2662865/69‑1 (MQ=255)
tCGCACGGCCGCCAAACTTGAGCATCATGCCGGGGTAAAGCGTGTCTAGCCGCTCCAGCGCATCGCGCa                                                                  <  1:3348752/69‑1 (MQ=255)
  gCATGACCGCCAAACTTCAGCATCATGCCAGGGTAGAGTGTGTCTAATCGCTCCAGCGCATCACGCata                                                                >  1:180666/1‑69 (MQ=255)
  gCATGACCGCCAAACTTCAGCATCATGCCAGGGTAGAGTGTGTCTAATCGCTCCAGCGCATCACGCata                                                                >  1:1131505/1‑69 (MQ=255)
  gCATGACCGCCAAACTTCAGCATCATGCCAGGGTAGAGTGTGTCTAATCGCTCCAGCGCATCACGCata                                                                >  1:2894227/1‑69 (MQ=255)
  gCATGACCGCCAAACTTCAGCATCATGCCAGGGTAGAGTGTGTCTAATCGCTCCAGCGCATCACGCata                                                                >  1:2926894/1‑69 (MQ=255)
  gCATGACCGCCAAACTTCAGCATCATGCCAGGGTAGAGTGTGTCTAATCGCTCCAGCGCATCACGCata                                                                >  1:387279/1‑69 (MQ=255)
  gCATGACCGCCAAACTTCAGCATCATGCCAGGGTAGAGTGTGTCTAATCGCTCCAGCGCATCACGCata                                                                >  1:800464/1‑69 (MQ=255)
  gCATGACCGCCAAACTTCAGCATCATGCCAGGGTAGAGTGTGTCTAATCGCTCCAGCGCATCACGCata                                                                >  1:948308/1‑69 (MQ=255)
  gCATGACCGCCAAACTTCAGCATCATGCCAGGGTAGAGTGTGTCTAATCGCTCCAGCGCATCACGCata                                                                >  1:1674839/1‑69 (MQ=255)
  gCATGACCGCCAAACTTCAGCATCATGCCAGGGTAGAGTGTGTCTAATCGCTCCAGCGCATCACGCa                                                                  >  1:2564959/1‑67 (MQ=255)
  gCATGACCGCCAAACTTCAGCATCATGCCAGGGTAGAGTGTGTCTAATCGCTCCAGCGCATCACGCa                                                                  >  1:2298498/1‑67 (MQ=255)
  gCATGACCGCCAAACTTCAGAATCATGCCAGGGTAGAGTGTGTCTAATAGCTCCAGCGCATCACGCata                                                                >  1:1262290/1‑69 (MQ=255)
               aCTTCAGCATCATGCCAGGGTAGAGTGTGTCTAATCGCTCCAGCGCATCACGCATATGCAGCCCCTGaa                                                   <  1:1767031/69‑1 (MQ=255)
               aCTTCAGCATCATGCCAGGGTAGAGTGTGTCTAATCGCTCCAGCGCATCACGCATATGCAGCCCCTGaa                                                   <  1:3027917/69‑1 (MQ=255)
               aCTTCAGCATCATGCCAGGGTAGAGTGTGTCTAATCGCTCCAGCGCATCACGCATATGCAGCCCCTGaa                                                   <  1:3204769/69‑1 (MQ=255)
               aCTTCAGCATCATGCCAGGGTAGAGTGTGTCTAATCGCTCCAGCGCATCACGCATATGCAGCCCCTGaa                                                   <  1:2744237/69‑1 (MQ=255)
               aCTTCAGCATCATGCCAGGGTAGAGTGTGTCTAATCGCTCCAGCGCATCACGCATATGCAGCCCCTGaa                                                   <  1:3419415/69‑1 (MQ=255)
               aCTTCAGCATCATGCCAGGGTAGAGTGTGTCTAATCGCTCCAGCGCATCACGCATATGCAGCCCCTGaa                                                   <  1:2206679/69‑1 (MQ=255)
               aCTTCAGCATCATGCCAGGGTAGAGTGTGTCTAATCGCTCCAGCGCATCACGCATATGCAGCCCCTGaa                                                   <  1:591937/69‑1 (MQ=255)
               aCTTCAGCATCATGCCAGGGTAGAGTGTGTCTAATCGCTCCAGCGCATCACGCATATGCAGCCCCTGaa                                                   <  1:864946/69‑1 (MQ=255)
               aCTTCAGCATCATGCCAGGGTAGAGTGTGTCTAATCGCTCCAGCGCATCACGCATATGCAGCCACTGaa                                                   <  1:2436144/69‑1 (MQ=255)
                cTTCAGCATCATGCCAGGGTAGAGTGTGTCTAATCGCTCCAGCGCATCACGCATATGCAGCCCCTGaa                                                   <  1:3176982/68‑1 (MQ=255)
                        tcatGCCAGGGTAGAGTGTGTCTAATCGCTCCAGCGCATCACGCa                                                                  >  1:1992312/1‑45 (MQ=255)
                        tcatGCCAGGGTAGAGTGTGTCTAATCGCTCCAGCGCATCACGCa                                                                  >  1:529705/1‑45 (MQ=255)
                        tcatGCCAGGGTAGAGTGTGTCTAATCGCTCCAGCGCATCACGCa                                                                  >  1:576026/1‑45 (MQ=255)
                        tcatGCCAGGGTAGAGTGTGTCTAATCGCTCCAGCGCATCACGCa                                                                  >  1:1908372/1‑45 (MQ=255)
                          atGCCGGGGTAAAGCGTGTCTAGCCGCTCCAGCGCATCGCGCa                                                                  <  1:136600/43‑1 (MQ=16)
                             ccAGGGTAGAGTGTGTCTAATCGCTCCAGCGCATCACGCATATGCAGCCCCTGAATGGAGCGACCGGaa                                     <  1:2926042/69‑1 (MQ=255)
                             ccAGGGTAGAGTGTGTCTAATCGCTCCAGCGCATCACGCATATGCAGCCCCTGAATGGAGCGACCGGaa                                     <  1:3176392/69‑1 (MQ=255)
                                   tAGAGTGTGTCTAATCGCTCCAGCGCATCACGCATATGCAGCCCCTGAATGGAGCGACCGGAACCTTTc                               >  1:1677255/1‑69 (MQ=255)
                                   tAGAGTGTGTCTAATCGCTCCAGCGCATCACGCATATGAAGCCCCTGAATGGAGCGACCGGAACCTTTc                               >  1:806464/1‑69 (MQ=255)
                                     gagTGTGTCTAATCGCTCCAGCGCATCACGCATATGCAGCCCCTg                                                     <  1:62570/45‑1 (MQ=255)
                                     gagTGTGTCTAATCGCTCCAGCGCATCACGCATATGCAGCCCCTg                                                     <  1:3315400/45‑1 (MQ=255)
                                     gagTGTGTCTAATCGCTCCAGCGCATCACGCATATGCAGCCCCTg                                                     <  1:3134717/45‑1 (MQ=255)
                                     gagTGTGTCTAATCGCTCCAGCGCATCACGCATATGCAGCCCCTg                                                     <  1:1458387/45‑1 (MQ=255)
                                     gagTGTGTCTAATCGCTCCAGCGCATCACGCATATGCAGCCCCTg                                                     <  1:2962859/45‑1 (MQ=255)
                                     gagTGTGTCTAATCGCTCCAGCGCATCACGCATATGCAGCCCCTg                                                     <  1:1614537/45‑1 (MQ=255)
                                     gagTGTGTCTAATCGCTCCAGCGCATCACGCATATGCAGCCCCTg                                                     <  1:1591676/45‑1 (MQ=255)
                                     gagTGTGTCTAATCGCTCCAGCGCATCACGCATATGCAGCCCCTg                                                     <  1:1595568/45‑1 (MQ=255)
                                     gagTGTGTCTAATCGCTCCAGCGCATCACGCATATGCAGCCCCTg                                                     <  1:2734312/45‑1 (MQ=255)
                                     gagTGTGTCTAATCGCTCCAGCGCATCACGCATATGCAGCCCCTGAATGGAGCGACCGGAACCTTTCAg                             <  1:3399272/69‑1 (MQ=255)
                                      agTGTGTCTAATCTCTCCAGCGCATCACGCATATGCAGCCCCTg                                                     <  1:179436/44‑1 (MQ=255)
                                      agTGTGTCTAATCGCTCCAGCGCATCACGCATATGCAGCCCCTg                                                     <  1:2983289/44‑1 (MQ=255)
                                        tgtgtCTAATCGCTCCAGCGCATCACGCATATGCAGCCCCTGAATGGAGCGACCGGAACCTTTCAGCGt                          >  1:1049379/1‑69 (MQ=255)
                                        tgtgtCTAATCGCTCCAGCGCATCACGCATATGCAGCCCCTGAATGGAGCGACCGGAACCTTTCAGCGt                          >  1:666164/1‑69 (MQ=255)
                                        tgtgtCTAATCGCTCCAGCGCATCACGCATATGCAGCCCCTGAATGGAGCGACCGGAACCTTTCAGCGt                          >  1:2055010/1‑69 (MQ=255)
                                        tgtgtCTAATCGCTCCAGCGCATCACGCATATGCAGCCCCTGAATGGAGCGACCGGAACCTTTCAGCGt                          >  1:1594388/1‑69 (MQ=255)
                                        tgtgtCTAATCGCTCCAGCGCATCACGCATATGCAGCCCCTGAATGGAGCGACCGGAACCTTTCAGCGt                          >  1:2673955/1‑69 (MQ=255)
                                        tgtgtCTAATCGCTCCAGCGCATCACGCATATGCAGCCCCTGAATGGAGCGACCGGAACCTTTCAGCGt                          >  1:1417411/1‑69 (MQ=255)
                                        tgtgtCTAATCGCTCCAGCGCATCACGCATATGCAGCCCCTGAATGGAGCGACCGGAACCTTTCAGCGt                          >  1:1555005/1‑69 (MQ=255)
                                        tgcgtCTAATCGCTCCAGCGCATCACGCATATGCAGCCCCTGAATGGAGCGACCGGAACCTTTCAGCGt                          >  1:72134/4‑69 (MQ=255)
                                              tAATCGCTCCAGCGCATCACGCATATGCAGCCCCTGAATGGAGCGACCGGAACCTTTCAGCGTACCGTc                    <  1:2132962/69‑1 (MQ=255)
                                              tAATCGCTCCAGCGCATCACGCATATGCAGCCCCTGAATGGAGCGACCGGAACCTTTCAGCGTACCGTc                    <  1:2819861/69‑1 (MQ=255)
                                              tAATCGCTCCAGCGCATCACGCATATGCAGCCCCTGAATGGAGCGACCGGAACCTTTCAGCGTACCGTc                    <  1:1911885/69‑1 (MQ=255)
                                               aaTCGCTCCAGCGCATCACGCATATGCAGCCCCTGAATGGAGCGACCGGAACCTTTCAGCGTACCGTc                    <  1:3225365/68‑1 (MQ=255)
                                                  cGCTCCAGCGCATCACGCATATGCAGCCCCTGAATGGAGCGACCGGAACCTTTCAGCGTACCGTCAcc                 <  1:1986391/68‑1 (MQ=255)
                                                   gCTCCAGCGCATCACGCATATGCAGCCCCTGAATGGAGCGACCGGAACCTTTCAGCGTACCGTCACctg               >  1:1162025/1‑69 (MQ=255)
                                                       cAGCGCATCGCGAATATGCAGCCCCTTAATGGAGCGACCGGAACCtt                                 >  1:484937/1‑47 (MQ=255)
                                                       cAGCGCATCGCGAATATGCAGCCCCTGAATGGAGCGACCGGAACCtt                                 >  1:669047/1‑47 (MQ=39)
                                                       cAGCGCATCGCGAATATGCAGCCCCTGAATGGAGCGACCGGAACCtt                                 >  1:912837/1‑47 (MQ=39)
                                                       cAGCGCATCGCGAATATGCAGCCCCTGAATGGAGCGACCGGAACCtt                                 >  1:1216753/1‑47 (MQ=39)
                                                       cAGCGCATCGCGAATATGCAGCCCCTGAATGGAGCGACCGGAACCtt                                 >  1:732176/1‑47 (MQ=39)
                                                       cAGCGCATCGCGAATATGCAGCCCCTGAATGGAGCGACCGGAACCtt                                 >  1:3386667/1‑47 (MQ=39)
                                                       cAGCGCATCGCGAATATGCAGCCCCTGAATGGAGCGACCGGAACCtt                                 >  1:3190995/1‑47 (MQ=39)
                                                                aCGCATATGCAGCCCCTGAATGGAGCGACCGGAACCTTTCAGCGTACCGTCa                   >  1:2412480/1‑52 (MQ=255)
                                                                aCGCATATGCAGCCCCTGAATGGAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAAcgc  >  1:1478072/1‑69 (MQ=255)
                                                                |                                                                    
TCGCATGACCGCCAAACTTCAGCATCATGCCAGGGTAGAGTGTGTCTAATCGCTCCAGCGCATCACGCATATGCAGCCCCTGAATGGAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGC  >  minE/1936084‑1936216

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: