Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,605,763 T→G 34.7% V118V (GTT→GTG 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,7630TG34.7% 109.2 / 121.6 147V118V (GTT→GTGamiAN‑acetylmuramoyl‑l‑alanine amidase I
Reads supporting (aligned to +/- strand):  ref base T (31/65);  new base G (15/36);  total (46/101)
Fisher's exact test for biased strand distribution p-value = 8.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 147 total aligned reads displayed.

CATGGGATTGATGCGCGTTTAACGCGTTCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGTTTATGTCAATTCATGCCGATGGCTTTA  >  minE/1605701‑1605821
                                                              |                                                          
cATGGCATTGATGCGCGTTTAACGCGTTCTGGCGATACGTTTATCCCACTGTACGATCGCGTGGAAATc                                                      <  1:1296455/69‑1 (MQ=255)
cATGGCATTGATGCGCGTTTAACGCGTTCTGGCGATACGTTTATCCCACTGTACGATCGCGTGGAAATc                                                      <  1:1577198/69‑1 (MQ=255)
cATGGCATTGATGCGCGTTTAACGCGTTCTGGCGATACGTTTATCCCACTGTACGATCGCGTGGAAATc                                                      <  1:2353143/69‑1 (MQ=255)
cATGGCATTGATGCGCGTTTAACGCGTTCTGGCGATACGTTTATCCCACTGTACGATCGCGTGGAAATc                                                      <  1:254728/69‑1 (MQ=255)
cATGGCATTGATGCGCGTTTAACGCGTTCTGGCGATACGTTTATCCCACTGTACGATCGCGTGGAAATc                                                      <  1:2556141/69‑1 (MQ=255)
cATGGCATTGATGCGCGTTTAACGCGTTCTGGCGATACGTTTATCCCACTGTACGATCGCGTGGAAATc                                                      <  1:2338581/69‑1 (MQ=255)
cATGGCATTGATGCGCGTTTAACGCGTTCTGGCGATACGTTTATCCCACTGTACGATCGCGTGGAAATc                                                      <  1:742394/69‑1 (MQ=255)
cATGGCATTGATGCGCGTTTAACGCGTTCTGGCGATACGTTTATCCCACTGTACGATCGCGTGGAAATc                                                      <  1:1407930/69‑1 (MQ=255)
cATGGCATTGATGCGCGTTTAACGCGTTCTGGCGATACGTTTATCCCACTGTACGATCGCGTGGAAATc                                                      <  1:3425678/69‑1 (MQ=255)
cATGGCATTGATGCGCGTTTAACGCGTTCTGGCGATACGTTTATCCCACTGTACGATCGCGTGGAAATc                                                      <  1:2742745/69‑1 (MQ=255)
cATGGCATTGATGCGCGTTTAACGCGTTCTGGCGATACGTTTATCCCACTGTACGATCGCGTGGAAATc                                                      <  1:130293/69‑1 (MQ=255)
cATGGCATTGATGCGCGTTTAACGCGTTCTGGCGATACGTTTATCCCACTGTACGATCGCGTGGAAATc                                                      <  1:286459/69‑1 (MQ=255)
cATGGCATTGATGCGCGTTTAACGCGTTCTGGCGATACGTTTATCCCACTGTACGATCGCGTGGAAATc                                                      <  1:1290072/69‑1 (MQ=255)
cATGGCATTGATGCGCGTTTAACGCGTTCTGGCGATACGTTTATCCCACTGTACGATCGCGTGGAAATc                                                      <  1:3265350/69‑1 (MQ=255)
cATGGCATTGATGCGCGTTTAACGCGTTCTGGCGATACGTTTATCCCACTGTACGATCGCGTGGAAATc                                                      <  1:2951836/69‑1 (MQ=255)
cATGGCATTGATGCGCGTTTAACGCGTTCTGGCGATACGTTTATCCCACTGTACGATCGCGTGGAAATc                                                      <  1:2292436/69‑1 (MQ=255)
cATGGCATTGATGCGCGTTTAACGCGTTCTGGCGATACGTTTATCCCACTGTACGATCGCGTGGAAATc                                                      <  1:1158317/69‑1 (MQ=255)
cATGGCATTGATGCGCGTTTAACGCGTTCTGGCGATACGTTTATCCCACTGTACGATCGCGTGGAAATc                                                      <  1:3068621/69‑1 (MQ=255)
cATGGCATTGATGCGCGTTTAACGCGTTCTGGCGATACGTTTATCCCACTGTACGATCGCGTGGAAATc                                                      <  1:1796158/69‑1 (MQ=255)
cATGGCATTGATGCGCGTTTAACGCGTTCTGGCGATACGTTTATCCCACTGTACGATCGCGTGGAAATc                                                      <  1:505425/69‑1 (MQ=255)
      aTTGATGCGCGTTTAACGCGTTCTGGCGATACGTTTATCCCACTGTACGATCGCGTGGAAATCGCCCa                                                 <  1:512797/68‑1 (MQ=255)
      aTTGATGCGCGTTTAACGCGTTCTGGCGATACGTTTATCCCACTGTACGATCGCGTGGAAATCGCCCa                                                 <  1:2191593/68‑1 (MQ=255)
      aTTGATGCGCGTTTAACGCGTTCTGGCGATACGTTTATCCCACTGTACGATCGCGTGGAAATCGCCCa                                                 <  1:1167114/68‑1 (MQ=255)
      aTTGATGCGCGTTTAACGCGTTCTGGCGATACGTTTATCCCACTGTACGATCGCGTGGAAATCGCCCa                                                 <  1:3123911/68‑1 (MQ=255)
      aTTGATGCGCGTTTAACGCGTTCTGGCGATACGTTTATCCCACTGTACGATCGCGTGGAAATCGCCCa                                                 <  1:1010726/68‑1 (MQ=255)
        tGATGCGCGTTTAACGCGTTCTGGCGATACGTTTATCCCACTGTACGATCGCGTGGAAATCGCCCa                                                 <  1:3185810/66‑1 (MQ=255)
         gATGCGCGTTTAACGCGTTCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATaaa                                             <  1:2526348/69‑1 (MQ=255)
         gATGCGCGTTTAACGCGTTCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATaaa                                             <  1:2629389/69‑1 (MQ=255)
         gATGCGCGTTTAACGCGTTCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATaaa                                             <  1:2268792/69‑1 (MQ=255)
         gATGCGCGTTTAACGCGTTCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATaaa                                             <  1:2139864/69‑1 (MQ=255)
         gATGCGCGTTTAACGCGTTCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATaaa                                             <  1:902707/69‑1 (MQ=255)
         gATGCGCGTTTAACGCGTTCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATaaa                                             <  1:614278/69‑1 (MQ=255)
         gATGCGCGTTTAACGCGTTCTGGCGATACGTTTATCCCACTGTACGATCGCGTGGAAATCGCCCa                                                 <  1:2644099/65‑1 (MQ=255)
           tGCGCGTTTAACGCGTTCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATaaa                                             <  1:3218251/67‑1 (MQ=255)
           tGCGCGTTTAACGCGTTCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATaaa                                             <  1:870916/67‑1 (MQ=255)
           tGCGCGTTTAACGCGTTCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATaaa                                             <  1:1273074/67‑1 (MQ=255)
           tGCGCGTTTAACGCGTTCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATaaa                                             <  1:1769632/67‑1 (MQ=255)
           tGCGCGTTTAACGCGTTCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATaaa                                             <  1:3073556/67‑1 (MQ=255)
           tGCGCGTTTAACGCGTTCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATaaa                                             <  1:1087127/67‑1 (MQ=255)
           tGCGCGTTTAACGCGTTCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATaaa                                             <  1:2134177/67‑1 (MQ=255)
                gTTTAACGCGTTCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGgcg                                      <  1:1597699/69‑1 (MQ=255)
                   tAACGCGTTCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATgg                                        >  1:1090221/1‑64 (MQ=255)
                   tAACGCGTTCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATgg                                        >  1:451891/1‑64 (MQ=255)
                   tAACGCGTTCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATgg                                        >  1:2765635/1‑64 (MQ=255)
                   tAACGCGTTCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATgg                                        >  1:2436349/1‑64 (MQ=255)
                   tAACGCGTTCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATgg                                        >  1:2351415/1‑64 (MQ=255)
                   tAACGCGTTCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATgg                                        >  1:2303292/1‑64 (MQ=255)
                   tAACGCGTTCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAg                                   >  1:3150670/1‑69 (MQ=255)
                          ttCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                            <  1:3334399/69‑1 (MQ=255)
                          ttCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                            <  1:3353829/69‑1 (MQ=255)
                          ttCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                            >  1:1247453/1‑69 (MQ=255)
                          ttCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                            <  1:3421363/69‑1 (MQ=255)
                          ttCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                            <  1:3214388/69‑1 (MQ=255)
                          ttCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                            <  1:958944/69‑1 (MQ=255)
                          ttCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                            >  1:1330926/1‑69 (MQ=255)
                          ttCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                            <  1:1496762/69‑1 (MQ=255)
                          ttCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                            <  1:2356228/69‑1 (MQ=255)
                          ttCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                            <  1:568769/69‑1 (MQ=255)
                          ttCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                            >  1:2023591/1‑69 (MQ=255)
                          ttCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                            <  1:48867/69‑1 (MQ=255)
                          ttCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                            <  1:1787308/69‑1 (MQ=255)
                          ttCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                            <  1:2274470/69‑1 (MQ=255)
                          ttCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                            <  1:1748517/69‑1 (MQ=255)
                          ttCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                            <  1:705450/69‑1 (MQ=255)
                          ttCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                            <  1:245217/69‑1 (MQ=255)
                          ttCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                            >  1:2482056/1‑69 (MQ=255)
                          ttCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                            <  1:65162/69‑1 (MQ=255)
                          ttCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                            >  1:1932207/1‑69 (MQ=255)
                           tCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCATATCTGtt                            <  1:2659415/68‑1 (MQ=255)
                                       tttATCCCACTGTACGATCGCGTGGAAATCGCCCATAAACACGGCGCAGACCTGTTTATGTCAATTCac               >  1:462733/1‑68 (MQ=255)
                                       tttATCCCACTGTACGATCGCGTGGAAATCGCCCATAAACACGGCGCAGACCTGTTTATGTCAATTCac               >  1:977726/1‑68 (MQ=255)
                                       tttATCCCACTGTACGATCGCGTGGAAATCGCCCATAAACACGGCGCAGACCTGTTTATGTCAATTCac               >  1:498261/1‑68 (MQ=255)
                                       tttATCCCACTGTACGATCGCGTGGAAATCGCCCATAAACACGGCGCAGACCTGTTTATGTCAATTCac               >  1:67483/1‑68 (MQ=255)
                                       tttATCCCACTGTACGATCGCGTGGAAATCGCCCATAAACACGGCGCAGACCTGTTTATGTCAATTCac               >  1:1932714/1‑68 (MQ=255)
                                       tttATCCCACTGTACGATCGCGTGGAAATCGCCCATAAACACGGCGCAGACCTGTTTATGTCAATTCac               >  1:1047646/1‑68 (MQ=255)
                                       tttATCCCACTGTACGATCGCGTGGAAATCGCCCATAAACACGGCGCAGACCTGTTTATGTCAATTCac               >  1:1408212/1‑68 (MQ=255)
                                       tttATCCCACTGTACGATCGCGTGGAAATCGCCCATAAACACGGCGCAGACCTGTTTATGTCAATTCac               >  1:1488004/1‑68 (MQ=255)
                                       tttATCCCACTGTACGATCGCGTGGAAATCGCCCATAAACACGGCGCAGACCTGTTTATGTCAATTCac               >  1:1587072/1‑68 (MQ=255)
                                       tttATCCCACTGTACGATCGCGTGGAAATCGCCCATAAACACGGCGCAGACCTGTTTATGTCAATTCac               >  1:2566684/1‑68 (MQ=255)
                                       tttATCCCACTGTACGATCGCGTGGAAATCGCCCATAAACACGGCGCAGACCTGTTTATGTCAATTCac               >  1:2841310/1‑68 (MQ=255)
                                       tttATCCCACTGTACGATCGCGTGGAAATCGCCCATAAACACGGCGCAGACCTGTTTATGTCAATTCac               >  1:2796265/1‑68 (MQ=255)
                                         tATCCTACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGTTTATGTCaa                    <  1:2499350/62‑1 (MQ=255)
                                         tATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                            >  1:1111890/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:1831966/62‑1 (MQ=255)
                                         tATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGTTTATGTCAATTCATg              >  1:2743256/1‑68 (MQ=255)
                                         tATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGTTTATGTCAATTCATg              >  1:1290256/1‑68 (MQ=255)
                                         tATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGTTTATGTCAATTCATg              >  1:2070538/1‑68 (MQ=255)
                                         tATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGTTTATGTCAATTCATg              >  1:482195/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:326306/68‑1 (MQ=255)
                                           tCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGTTTATGTCAATTCATGcc            <  1:2054120/68‑1 (MQ=255)
                                           tCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGTTTATGTCAATTCATGcc            <  1:453587/68‑1 (MQ=255)
                                                  ttACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGtt                            <  1:927004/45‑1 (MQ=255)
                                                     cGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGTTTATGTCAATTCATGCCGATGGCTTTa  <  1:861832/68‑1 (MQ=255)
                                                              |                                                          
CATGGGATTGATGCGCGTTTAACGCGTTCTGGCGATACGTTTATCCCACTTTACGATCGCGTTGAAATCGCCCATAAACATGGCGCAGATCTGTTTATGTCAATTCATGCCGATGGCTTTA  >  minE/1605701‑1605821

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: