Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,605,781 T→C 18.3% H124H (CAT→CAC 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,7810TC18.3% 182.7 / 28.0 119H124H (CAT→CACamiAN‑acetylmuramoyl‑l‑alanine amidase I
Reads supporting (aligned to +/- strand):  ref base T (54/43);  new base C (15/7);  total (69/50)
Fisher's exact test for biased strand distribution p-value = 3.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 124 total aligned reads displayed.

GTTTAACGCGTTCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCTTCGGT  >  minE/1605717‑1605849
                                                                |                                                                    
gTTTAACGCGTTCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGgcg                                                                  <  1:1597699/69‑1 (MQ=255)
   tAACGCGTTCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATgg                                                                    >  1:366049/1‑64 (MQ=255)
   tAACGCGTTCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATgg                                                                    >  1:1090221/1‑64 (MQ=255)
   tAACGCGTTCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATgg                                                                    >  1:2765635/1‑64 (MQ=255)
   tAACGCGTTCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATgg                                                                    >  1:2436349/1‑64 (MQ=255)
   tAACGCGTTCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATgg                                                                    >  1:2184880/1‑64 (MQ=255)
   tAACGCGTTCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATgg                                                                    >  1:2351415/1‑64 (MQ=255)
   tAACGCGTTCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATgg                                                                    >  1:1278411/1‑64 (MQ=255)
   tAACGCGTTCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAg                                                               >  1:3150670/1‑69 (MQ=255)
          ttCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                                                        <  1:2294409/69‑1 (MQ=255)
          ttCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                                                        <  1:2900045/69‑1 (MQ=255)
          ttCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                                                        <  1:2274470/69‑1 (MQ=255)
          ttCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                                                        >  1:1089582/1‑69 (MQ=255)
          ttCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                                                        <  1:2356228/69‑1 (MQ=255)
          ttCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                                                        <  1:2083035/69‑1 (MQ=255)
          ttCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                                                        <  1:3210424/69‑1 (MQ=255)
          ttCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                                                        >  1:2023591/1‑69 (MQ=255)
          ttCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                                                        <  1:245217/69‑1 (MQ=255)
          ttCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                                                        <  1:3214388/69‑1 (MQ=255)
          ttCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                                                        >  1:3224015/1‑69 (MQ=255)
          ttCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                                                        <  1:3334399/69‑1 (MQ=255)
          ttCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                                                        <  1:3421363/69‑1 (MQ=255)
          ttCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                                                        >  1:1932207/1‑69 (MQ=255)
          ttCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                                                        >  1:2482056/1‑69 (MQ=255)
          ttCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                                                        <  1:476515/69‑1 (MQ=255)
          ttCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                                                        <  1:1748517/69‑1 (MQ=255)
          ttCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                                                        <  1:48867/69‑1 (MQ=255)
          ttCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                                                        <  1:958944/69‑1 (MQ=255)
          ttCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                                                        <  1:568769/69‑1 (MQ=255)
          ttCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                                                        <  1:1542277/69‑1 (MQ=255)
          ttCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                                                        <  1:628110/69‑1 (MQ=255)
          ttCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                                                        <  1:1496762/69‑1 (MQ=255)
          ttCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                                                        <  1:65162/69‑1 (MQ=255)
          ttCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                                                        <  1:750743/69‑1 (MQ=255)
          ttCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                                                        >  1:1247453/1‑69 (MQ=255)
          ttCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                                                        >  1:1330926/1‑69 (MQ=255)
          ttCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                                                        <  1:705450/69‑1 (MQ=255)
           tCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCATATCTGtt                                                        <  1:2659415/68‑1 (MQ=255)
           tCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                                                        <  1:901190/68‑1 (MQ=255)
                       tttATCCCACTGTACGATCGCGTGGAAATCGCCCATAAACACGGCGCAGACCTGTTTATGTCAATTCac                                           >  1:701412/1‑68 (MQ=255)
                       tttATCCCACTGTACGATCGCGTGGAAATCGCCCATAAACACGGCGCAGACCTGTTTATGTCAATTCac                                           >  1:2841310/1‑68 (MQ=255)
                       tttATCCCACTGTACGATCGCGTGGAAATCGCCCATAAACACGGCGCAGACCTGTTTATGTCAATTCac                                           >  1:3138047/1‑68 (MQ=255)
                       tttATCCCACTGTACGATCGCGTGGAAATCGCCCATAAACACGGCGCAGACCTGTTTATGTCAATTCac                                           >  1:2796265/1‑68 (MQ=255)
                       tttATCCCACTGTACGATCGCGTGGAAATCGCCCATAAACACGGCGCAGACCTGTTTATGTCAATTCac                                           >  1:498261/1‑68 (MQ=255)
                       tttATCCCACTGTACGATCGCGTGGAAATCGCCCATAAACACGGCGCAGACCTGTTTATGTCAATTCac                                           >  1:1047646/1‑68 (MQ=255)
                       tttATCCCACTGTACGATCGCGTGGAAATCGCCCATAAACACGGCGCAGACCTGTTTATGTCAATTCac                                           >  1:977726/1‑68 (MQ=255)
                       tttATCCCACTGTACGATCGCGTGGAAATCGCCCATAAACACGGCGCAGACCTGTTTATGTCAATTCac                                           >  1:462733/1‑68 (MQ=255)
                       tttATCCCACTGTACGATCGCGTGGAAATCGCCCATAAACACGGCGCAGACCTGTTTATGTCAATTCac                                           >  1:2566684/1‑68 (MQ=255)
                       tttATCCCACTGTACGATCGCGTGGAAATCGCCCATAAACACGGCGCAGACCTGTTTATGTCAATTCac                                           >  1:1408212/1‑68 (MQ=255)
                       tttATCCCACTGTACGATCGCGTGGAAATCGCCCATAAACACGGCGCAGACCTGTTTATGTCAATTCac                                           >  1:2353333/1‑68 (MQ=255)
                         taacccACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGTTTATGTCAATTCATg                                          >  1:2555619/4‑68 (MQ=255)
                         tATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                                                        >  1:1319736/1‑54 (MQ=255)
                         tATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                                                        >  1:1543411/1‑54 (MQ=255)
                         tATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                                                        >  1:1614128/1‑54 (MQ=255)
                         tATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                                                        >  1:460806/1‑54 (MQ=255)
                         tATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                                                        >  1:1111890/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:2070538/1‑68 (MQ=255)
                         tATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGTTTATGTCAATTCATg                                          >  1:482195/1‑68 (MQ=255)
                         tATCCCACTTTACGATCGCCTTGAAATCGCCCATAAACATGGCGCAGATc                                                            >  1:1427936/1‑50 (MQ=255)
                           tCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGTTTATGTCAATTCATGcc                                        <  1:1811357/68‑1 (MQ=255)
                           tCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGTTTATGTCAATTCATGcc                                        <  1:326306/68‑1 (MQ=255)
                           tCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGTTTATGTCAATTCATGcc                                        <  1:453587/68‑1 (MQ=255)
                           tCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGTTTATGTCAATTCATGcc                                        <  1:338006/68‑1 (MQ=255)
                                  ttACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                                                        <  1:927004/45‑1 (MQ=255)
                                  ttACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                                                        <  1:1210332/45‑1 (MQ=255)
                                       aTCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGTTTATGTCAATTCATGCCGATGGCTTTa                              <  1:200862/66‑1 (MQ=255)
                                               gAAATCGCCCATAAACATGGCGCAGATCTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGaaa                   >  1:318909/1‑69 (MQ=255)
                                               gAAATCGCCCATAAACATGGCGCAGATCTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGaaa                   >  1:1997699/1‑69 (MQ=255)
                                                aaaTCGCCCATAAACACGGCGCAGACCTGTTTATGTCAATTCACGCCGACGGTTTTACCAAc                         <  1:2681066/62‑1 (MQ=255)
                                                aaaTCGCCCATAAACACGGCGCAGACCTGTTTATGTCAATTCACGCCGACGGTTTTACCAAc                         <  1:1242439/62‑1 (MQ=255)
                                                aaaTCGCCCATAAACACGGCGCAGACCTGTTTATGTCAATTCACGCCGACGGTTTTACCAAc                         <  1:2760816/62‑1 (MQ=255)
                                                aaaTCGCCCATAAACACGGCGCAGACCTGTTTATGTCAATTCACGCCGACGGTTTTACCAAc                         <  1:429010/62‑1 (MQ=255)
                                                aaaTCGCCCATAAACACGGCGCAGACCTGTTTATGTCAATTCACGCCGACGGTTTTACCAAc                         <  1:3276798/62‑1 (MQ=255)
                                                 aaTCGCCCATAAACACGGCGCAGACCTGTTTATGTCAATTCACGCCGACGGTTTTACCAAc                         <  1:1737809/61‑1 (MQ=255)
                                                            aaCATGGCGCAGATCTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCtt      >  1:1911441/1‑69 (MQ=255)
                                                            aaCATGGCGCAGATCTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCtt      >  1:78794/1‑69 (MQ=255)
                                                            aaCATGGCGCAGATCTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCtt      >  1:263846/1‑69 (MQ=255)
                                                            aaCATGGCGCAGATCTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCtt      >  1:46101/1‑69 (MQ=255)
                                                            aaCATGGCGCAGATCTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCtt      >  1:3077431/1‑69 (MQ=255)
                                                            aaCATGGCGCAGATCTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCtt      >  1:364407/1‑69 (MQ=255)
                                                            aaCATGGCGCAGATCTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCtt      >  1:3008078/1‑69 (MQ=255)
                                                            aaCATGGCGCAGATCTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCtt      >  1:2964791/1‑69 (MQ=255)
                                                            aaCATGGCGCAGATCTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCtt      >  1:2283325/1‑69 (MQ=255)
                                                            aaCATGGCGCAGATCTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCtt      >  1:1514579/1‑69 (MQ=255)
                                                            aaCATGGCGCAGATCTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCtt      >  1:2746564/1‑69 (MQ=255)
                                                            aaCATGGCGCAGATCTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCt       >  1:2244456/1‑68 (MQ=255)
                                                            aaCATGGCGCAGATCTGTTTATGTCAATTCATACCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCtt      >  1:3068604/1‑69 (MQ=255)
                                                              cATGGCGCAGATCTGTTTATGTCAATTCATGCCGATGGCTTTTCCAACCCGAAAGCTGCCGGTGCTTc     >  1:2227352/1‑68 (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:2007635/69‑1 (MQ=255)
                                                                tGGCGCAGATCTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCTTCGGt  <  1:3178286/69‑1 (MQ=255)
                                                                tGGCGCAGATCTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCTTCGGt  <  1:2365550/69‑1 (MQ=255)
                                                                |                                                                    
GTTTAACGCGTTCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGTTTATGTCAATTCATGCCGATGGCTTTACCAACCCGAAAGCTGCCGGTGCTTCGGT  >  minE/1605717‑1605849

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: