Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,605,790 T→C 18.2% D127D (GAT→GAC amiA → N‑acetylmuramoyl‑l‑alanine amidase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,605,7900TC18.2% 185.5 / 27.8 120D127D (GAT→GACamiAN‑acetylmuramoyl‑l‑alanine amidase I
Reads supporting (aligned to +/- strand):  ref base T (45/53);  new base C (15/7);  total (60/60)
Fisher's exact test for biased strand distribution p-value = 9.72e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 121 total aligned reads displayed.

TTCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCTTCGGTATTTGCCCT  >  minE/1605727‑1605858
                                                               |                                                                    
ttCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                                                                 <  1:3033629/69‑1 (MQ=255)
ttCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                                                                 <  1:3353829/69‑1 (MQ=255)
ttCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                                                                 <  1:3334399/69‑1 (MQ=255)
ttCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                                                                 <  1:1748517/69‑1 (MQ=255)
ttCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                                                                 <  1:1787308/69‑1 (MQ=255)
ttCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                                                                 >  1:3224015/1‑69 (MQ=255)
ttCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                                                                 <  1:3214388/69‑1 (MQ=255)
ttCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                                                                 <  1:3210424/69‑1 (MQ=255)
ttCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                                                                 >  1:1932207/1‑69 (MQ=255)
ttCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                                                                 <  1:958944/69‑1 (MQ=255)
ttCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                                                                 <  1:2900045/69‑1 (MQ=255)
ttCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                                                                 >  1:2023591/1‑69 (MQ=255)
ttCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                                                                 <  1:2654930/69‑1 (MQ=255)
ttCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                                                                 <  1:2083035/69‑1 (MQ=255)
ttCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                                                                 >  1:2482056/1‑69 (MQ=255)
ttCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                                                                 <  1:245217/69‑1 (MQ=255)
ttCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                                                                 <  1:2356228/69‑1 (MQ=255)
ttCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                                                                 <  1:2274470/69‑1 (MQ=255)
ttCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                                                                 <  1:1542277/69‑1 (MQ=255)
ttCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                                                                 <  1:3421363/69‑1 (MQ=255)
ttCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                                                                 <  1:476515/69‑1 (MQ=255)
ttCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                                                                 <  1:534992/69‑1 (MQ=255)
ttCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                                                                 >  1:1330926/1‑69 (MQ=255)
ttCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                                                                 <  1:568769/69‑1 (MQ=255)
ttCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                                                                 <  1:628110/69‑1 (MQ=255)
ttCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                                                                 >  1:1247453/1‑69 (MQ=255)
ttCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                                                                 >  1:1089582/1‑69 (MQ=255)
ttCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                                                                 <  1:750743/69‑1 (MQ=255)
ttCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                                                                 <  1:705450/69‑1 (MQ=255)
 tCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCATATCTGtt                                                                 <  1:2659415/68‑1 (MQ=255)
 tCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                                                                 <  1:901190/68‑1 (MQ=255)
 tCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                                                                 <  1:2754789/68‑1 (MQ=255)
             tttATCCCACTGTACGATCGCGTGGAAATCGCCCATAAACACGGCGCAGACCTGTTTATGTCAATTCac                                                    >  1:2353333/1‑68 (MQ=255)
             tttATCCCACTGTACGATCGCGTGGAAATCGCCCATAAACACGGCGCAGACCTGTTTATGTCAATTCac                                                    >  1:977726/1‑68 (MQ=255)
             tttATCCCACTGTACGATCGCGTGGAAATCGCCCATAAACACGGCGCAGACCTGTTTATGTCAATTCac                                                    >  1:701412/1‑68 (MQ=255)
             tttATCCCACTGTACGATCGCGTGGAAATCGCCCATAAACACGGCGCAGACCTGTTTATGTCAATTCac                                                    >  1:67483/1‑68 (MQ=255)
             tttATCCCACTGTACGATCGCGTGGAAATCGCCCATAAACACGGCGCAGACCTGTTTATGTCAATTCac                                                    >  1:2796265/1‑68 (MQ=255)
             tttATCCCACTGTACGATCGCGTGGAAATCGCCCATAAACACGGCGCAGACCTGTTTATGTCAATTCac                                                    >  1:498261/1‑68 (MQ=255)
             tttATCCCACTGTACGATCGCGTGGAAATCGCCCATAAACACGGCGCAGACCTGTTTATGTCAATTCac                                                    >  1:462733/1‑68 (MQ=255)
             tttATCCCACTGTACGATCGCGTGGAAATCGCCCATAAACACGGCGCAGACCTGTTTATGTCAATTCac                                                    >  1:1488004/1‑68 (MQ=255)
             tttATCCCACTGTACGATCGCGTGGAAATCGCCCATAAACACGGCGCAGACCTGTTTATGTCAATTCac                                                    >  1:1587072/1‑68 (MQ=255)
             tttATCCCACTGTACGATCGCGTGGAAATCGCCCATAAACACGGCGCAGACCTGTTTATGTCAATTCac                                                    >  1:1932714/1‑68 (MQ=255)
             tttATCCCACTGTACGATCGCGTGGAAATCGCCCATAAACACGGCGCAGACCTGTTTATGTCAATTCac                                                    >  1:1408212/1‑68 (MQ=255)
             tttATCCCACTGTACGATCGCGTGGAAATCGCCCATAAACACGGCGCAGACCTGTTTATGTCAATTCac                                                    >  1:1047646/1‑68 (MQ=255)
               taacccACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGTTTATGTCAATTCATg                                                   >  1:2555619/4‑68 (MQ=255)
               tATCCTACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGTTTATGTCaa                                                         <  1:2499350/62‑1 (MQ=255)
               tATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                                                                 >  1:1614128/1‑54 (MQ=255)
               tATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                                                                 >  1:1319736/1‑54 (MQ=255)
               tATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                                                                 >  1:1543411/1‑54 (MQ=255)
               tATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                                                                 >  1:2658287/1‑54 (MQ=255)
               tATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGt                                                                  >  1:2919258/1‑53 (MQ=255)
               tATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGTTTATGTCaa                                                         <  1:2447629/62‑1 (MQ=255)
               tATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGTTTATGTCaa                                                         <  1:1831966/62‑1 (MQ=255)
               tATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGTTTATGTCAATTCATg                                                   >  1:2743256/1‑68 (MQ=255)
               tATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGTTTATGTCAATTCATg                                                   >  1:482195/1‑68 (MQ=255)
               tATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGTTTATGTCAATTCATg                                                   >  1:2070538/1‑68 (MQ=255)
               tATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGTTTATGTCAATTCATg                                                   >  1:475762/1‑68 (MQ=255)
               tATCCCACTTTACGATCGCCTTGAAATCGCCCATAAACATGGCGCAGATc                                                                     >  1:1427936/1‑50 (MQ=255)
                 tCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGTTTATGTCAATTCATGcc                                                 <  1:1811357/68‑1 (MQ=255)
                 tCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGTTTATGTCAATTCATGcc                                                 <  1:2054120/68‑1 (MQ=255)
                 tCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGTTTATGTCAATTCATGcc                                                 <  1:453587/68‑1 (MQ=255)
                 tCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGTTTATGTCAATTCATGcc                                                 <  1:326306/68‑1 (MQ=255)
                        ttACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                                                                 <  1:927004/45‑1 (MQ=255)
                        ttACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                                                                 <  1:1210332/45‑1 (MQ=255)
                                     gAAATCGCCCATAAACATGGCGCAGATCTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGaaa                            >  1:1997699/1‑69 (MQ=255)
                                     gAAATCGCCCATAAACATGGCGCAGATCTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGaaa                            >  1:3176543/1‑69 (MQ=255)
                                     gAAATCGCCCATAAACATGGCGCAGATCTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGaaa                            >  1:318909/1‑69 (MQ=255)
                                      aaaTCGCCCATAAACACGGCGCAGACCTGTTTATGTCAATTCACGCCGACGGTTTTACCAAc                                  <  1:429010/62‑1 (MQ=255)
                                      aaaTCGCCCATAAACACGGCGCAGACCTGTTTATGTCAATTCACGCCGACGGTTTTACCAAc                                  <  1:2760816/62‑1 (MQ=255)
                                      aaaTCGCCCATAAACACGGCGCAGACCTGTTTATGTCAATTCACGCCGACGGTTTTACCAAc                                  <  1:1242439/62‑1 (MQ=255)
                                      aaaTCGCCCATAAACACGGCGCAGACCTGTTTATGTCAATTCACGCCGACGGTTTTACCAAc                                  <  1:3276798/62‑1 (MQ=255)
                                      aaaTCGCCCATAAACACGGCGCAGACCTGTTTATGTCAATTCACGCCGACGGTTTTACCAAc                                  <  1:223479/62‑1 (MQ=255)
                                                  aaCATGGCGCAGATCTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCtt               >  1:46101/1‑69 (MQ=255)
                                                  aaCATGGCGCAGATCTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCtt               >  1:1514579/1‑69 (MQ=255)
                                                  aaCATGGCGCAGATCTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCtt               >  1:364407/1‑69 (MQ=255)
                                                  aaCATGGCGCAGATCTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCtt               >  1:2283325/1‑69 (MQ=255)
                                                  aaCATGGCGCAGATCTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCtt               >  1:78794/1‑69 (MQ=255)
                                                  aaCATGGCGCAGATCTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCtt               >  1:1911441/1‑69 (MQ=255)
                                                  aaCATGGCGCAGATCTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCtt               >  1:3077431/1‑69 (MQ=255)
                                                  aaCATGGCGCAGATCTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCtt               >  1:3008078/1‑69 (MQ=255)
                                                  aaCATGGCGCAGATCTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCtt               >  1:284447/1‑69 (MQ=255)
                                                  aaCATGGCGCAGATCTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCtt               >  1:2746564/1‑69 (MQ=255)
                                                  aaCATGGCGCAGATCTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCtt               >  1:263846/1‑69 (MQ=255)
                                                  aaCATGGCGCAGATCTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCt                >  1:2244456/1‑68 (MQ=255)
                                                  aaCATGGCGCAGATCTGTTTATGTCAATTCATACCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCtt               >  1:3068604/1‑69 (MQ=255)
                                                    cATGGCGCAGATCTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCTTc              >  1:1101430/1‑68 (MQ=255)
                                                    cATGGCGCAGATCTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCTTc              >  1:2019128/1‑68 (MQ=255)
                                                    cATGGCGCAGATCTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCTTc              >  1:2637484/1‑68 (MQ=255)
                                                    cATGGCGCAGATCTGTTTATGTCAATTCATCCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCTTc              >  1:539958/1‑68 (MQ=255)
                                                      tGGCGCAGATCTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCTTCGGt           <  1:3178286/69‑1 (MQ=255)
                                                      tGGCGCAGATCTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCTTCGGt           <  1:2365550/69‑1 (MQ=255)
                                                      tGGCGCAGATCTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCTTCGGt           <  1:1198488/69‑1 (MQ=255)
                                                      tGGCGCAGATCTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCTTCGGt           <  1:2007635/69‑1 (MQ=255)
                                                           cAGATCTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCTTCGGTATTTg      >  1:128318/1‑69 (MQ=255)
                                                             gATCTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCTTCGGTATTTGcc    <  1:2933912/69‑1 (MQ=255)
                                                             gATCTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCTTCGGTATTTGcc    <  1:2968248/69‑1 (MQ=255)
                                                             gATCTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCTTCGGTATTTGcc    <  1:2252118/69‑1 (MQ=255)
                                                             gATCTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCTTCGGTATTTGcc    <  1:2314080/69‑1 (MQ=255)
                                                             gATCTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCTTCGGTATTTGcc    <  1:1037867/69‑1 (MQ=255)
                                                               tCTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCTTCGGTATTTGCCct  >  1:1040309/1‑69 (MQ=255)
                                                               |                                                                    
TTCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCTTCGGTATTTGCCCT  >  minE/1605727‑1605858

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