Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,973,916 T→C 31.1% R104R (CGT→CGC nadR → bifunctional DNA‑binding transcriptional repressor and NMN adenylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,973,9160TC31.1% 79.8 / 46.2 93R104R (CGT→CGCnadRbifunctional DNA‑binding transcriptional repressor and NMN adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base T (20/44);  new base C (14/15);  total (34/59)
Fisher's exact test for biased strand distribution p-value = 1.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

TACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCGCAGCAGCCGACCGTGCCGGATCGTCTGCGTTGGT  >  minE/2973848‑2973983
                                                                    |                                                                   
tACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGt                                                                     >  1:1792875/1‑69 (MQ=255)
tACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGt                                                                     >  1:2141262/1‑69 (MQ=255)
tACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGt                                                                     >  1:1315836/1‑69 (MQ=255)
         cAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGCGACcgcgcg                                                            <  1:3232797/69‑1 (MQ=255)
         cAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGCGACcgcgcg                                                            <  1:1975762/69‑1 (MQ=255)
         cAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGCGACcgcgcg                                                            <  1:2244697/69‑1 (MQ=255)
         cAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGCGACcgcgcg                                                            <  1:2499952/69‑1 (MQ=255)
         cAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGCGACcgcgcg                                                            <  1:2733873/69‑1 (MQ=255)
         cAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGCGACcgcgcg                                                            <  1:1602092/69‑1 (MQ=255)
         cAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGCGACcgcgcg                                                            <  1:1576596/69‑1 (MQ=255)
         cAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGCGACcgcgcg                                                            <  1:400885/69‑1 (MQ=255)
         cAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGCGACcgcgcg                                                            <  1:472235/69‑1 (MQ=255)
         cAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGCGACcgcgcg                                                            <  1:756211/69‑1 (MQ=255)
         cAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGCGACcgcgcg                                                            <  1:819974/69‑1 (MQ=255)
         cAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGCGACcgcgcg                                                            <  1:1215358/69‑1 (MQ=255)
         cAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGCGACcgcgcg                                                            <  1:89748/69‑1 (MQ=255)
         cAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGCGACcgcgcg                                                            <  1:965855/69‑1 (MQ=255)
          aGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGCGACcgcgcg                                                            <  1:1207122/68‑1 (MQ=255)
                  tGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGaa                                                   <  1:2954913/69‑1 (MQ=255)
                   gTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGCGACCGCGCgttgtttg                                                     >  1:1771240/1‑64 (MQ=255)
                   gTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGCGACCGCGCgttgtttg                                                     >  1:669204/1‑64 (MQ=255)
                   gTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGCGACCGCGCgttgtttg                                                     >  1:534143/1‑64 (MQ=255)
                   gTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGCGACCGCGCgttgtttg                                                     >  1:2526043/1‑64 (MQ=255)
                   gTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGCGACCGCGCgttgtttg                                                     >  1:3387385/1‑64 (MQ=255)
                   gTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGCGACCGCGCgttgtttg                                                     >  1:1705739/1‑64 (MQ=255)
                   gTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGCGACCGCGCgttgtttg                                                     >  1:2791842/1‑64 (MQ=255)
                   gTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGCGACCGCGCgttgtttg                                                     >  1:2177256/1‑64 (MQ=255)
                   gTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGCGACCGCGCgttgtttg                                                     >  1:2156196/1‑64 (MQ=255)
                   gTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGCGACCGCGCgttgtttg                                                     >  1:2867340/1‑64 (MQ=255)
                      tccAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGa                                                 <  1:3189779/66‑1 (MQ=255)
                      gCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGa                                                 <  1:2470990/67‑1 (MQ=255)
                      gCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGa                                                 <  1:247839/67‑1 (MQ=255)
                      gCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGa                                                 <  1:1697892/67‑1 (MQ=255)
                      gCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGa                                                 <  1:2680510/67‑1 (MQ=255)
                      gCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGa                                                 <  1:1707346/67‑1 (MQ=255)
                      gCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGa                                                 <  1:213863/67‑1 (MQ=255)
                      gCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGa                                                 <  1:2098440/67‑1 (MQ=255)
                      gCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGa                                                 <  1:2097324/67‑1 (MQ=255)
                      gCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGa                                                 <  1:3175578/67‑1 (MQ=255)
                      gCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGa                                                 <  1:325934/67‑1 (MQ=255)
                      gCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGa                                                 <  1:1458312/67‑1 (MQ=255)
                      gCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGa                                                 <  1:1278294/67‑1 (MQ=255)
                      gCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGAc                                                <  1:2595549/68‑1 (MQ=255)
                      gCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGAc                                                <  1:778125/68‑1 (MQ=255)
                       ccAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGa                                                 <  1:2580237/66‑1 (MQ=255)
                        cAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGAc                                                <  1:2372746/66‑1 (MQ=255)
                          ggTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGc                                           <  1:623469/69‑1 (MQ=255)
                             tGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCAt                                        >  1:241201/1‑69 (MQ=255)
                             tGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCAt                                        >  1:3126611/1‑69 (MQ=255)
                             tGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCAt                                        <  1:2030476/69‑1 (MQ=255)
                             tGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCAt                                        >  1:2013207/1‑69 (MQ=255)
                             tGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCAt                                        >  1:378751/1‑69 (MQ=255)
                             tGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCAt                                        <  1:2180976/69‑1 (MQ=255)
                             tGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCAt                                        >  1:924539/1‑69 (MQ=255)
                             tGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCAt                                        >  1:1463674/1‑69 (MQ=255)
                             tGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCAt                                        <  1:2199429/69‑1 (MQ=255)
                             tGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCAt                                        <  1:884341/69‑1 (MQ=255)
                             tGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCAt                                        >  1:2374246/1‑69 (MQ=255)
                              gACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATg                                       >  1:264620/1‑69 (MQ=255)
                                    cTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCgc                                   >  1:1939453/1‑67 (MQ=255)
                                    cTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCgc                                   >  1:562842/1‑67 (MQ=255)
                                    cTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCgc                                   >  1:2564059/1‑67 (MQ=255)
                                    cTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCgc                                   >  1:1435918/1‑67 (MQ=255)
                                    cTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCgc                                   >  1:960191/1‑67 (MQ=255)
                                      gCATATCATTATGGGTCTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCgcagca                               >  1:2393215/1‑69 (MQ=255)
                                         tatCATTATGGGTTTTGACTATACCCGTGACCGCGCGTTGTTCGAAGACAGTGc                                           <  1:3146761/54‑1 (MQ=255)
                                         tatCATTATGGGTTTTGACGATACGCGTGACCGCGCGTTGTTCGAAGACAGTGc                                           <  1:1207906/54‑1 (MQ=255)
                                         tatCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGc                                           <  1:2654159/54‑1 (MQ=255)
                                         tatCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGc                                           <  1:677750/54‑1 (MQ=255)
                                         tatCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGc                                           <  1:1866826/54‑1 (MQ=255)
                                         tatCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGc                                           <  1:1880901/54‑1 (MQ=255)
                                         tatCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGc                                           <  1:2217895/54‑1 (MQ=255)
                                         tatCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGc                                           <  1:562215/54‑1 (MQ=255)
                                         tatCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGc                                           <  1:2608807/54‑1 (MQ=255)
                                         tatCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGc                                           <  1:3105538/54‑1 (MQ=255)
                                         tatCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGc                                           <  1:171514/54‑1 (MQ=255)
                                         tatCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGc                                           <  1:1615319/54‑1 (MQ=255)
                                         tatCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGc                                           <  1:2813345/54‑1 (MQ=255)
                                         tatCATTATCGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGc                                           <  1:2294639/54‑1 (MQ=255)
                                            cATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCGCAGCAGCCGa                          >  1:2073900/1‑68 (MQ=255)
                                            cATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCGCAGCAGCCGAc                         >  1:1487673/1‑69 (MQ=255)
                                                            gATACCCGCGACCGCGCGTTGTTTGAAGACAGCGCCATGTCGCAGCAGCCCACGGTGCCGGATCGTCt          >  1:191601/1‑68 (MQ=255)
                                                            gATACCCGCGACCGCGCGTTGTTTGAAGACAGCGCCATGTCGCAGCAGCCCACGGTGCCGGATCGTCt          >  1:2026209/1‑68 (MQ=255)
                                                            gATACCCGCGACCGCGCGTTGTTTGAAGACAGCGCCATGTCGCAGCAGCCCACGGTGCCGGATCGTCTg         >  1:497485/1‑69 (MQ=255)
                                                            gATACCCGCGACCGCGCGTTGTTTGAAGACAGCGCCATGTCGCAGCAGCCCACGGTGCCGGATCGTCTg         >  1:2153758/1‑69 (MQ=255)
                                                                  cGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCGCAGCAGCCGACCGTGCCGGATCGTCTGCGTTg    >  1:981808/1‑68 (MQ=255)
                                                                  cGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCGCAGCAGCCGACCGTGCCGGATCGTCTGCGTTg    >  1:1014677/1‑68 (MQ=255)
                                                                  cGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCGCAGCAGCCGACCGTGCCGGATCGTCTGCGTTg    >  1:620182/1‑68 (MQ=255)
                                                                  cGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCGCAGCAGCCGACCGTGCCGGATCGTCTGCGTTg    >  1:493184/1‑68 (MQ=255)
                                                                   gTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCGCAGCAGCCGACCGTGCCGGATCGTCTGCGTTGGt  <  1:653608/69‑1 (MQ=255)
                                                                   gTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCGCAGCAGCCGACCGTGCCGGATCGTCTGCGTTGGt  <  1:475585/69‑1 (MQ=255)
                                                                   gTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCGCAGCAGCCGACCGTGCCGGATCGTCTGCGTTGGt  <  1:2919372/69‑1 (MQ=255)
                                                                   gTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCGCAGCAGCCGACCGTGCCGGATCGTCTGCGTTGGt  <  1:2772378/69‑1 (MQ=255)
                                                                   gTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCGCAGCAGCCGACCGTGCCGGATCGTCTGCGTTGGt  <  1:2519104/69‑1 (MQ=255)
                                                                   gTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCGCAGCAGCCGACCGTGCCGGATCGTCTGCGTTGGt  <  1:2143417/69‑1 (MQ=255)
                                                                   gTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCGCAGCAGCCGACCGTGCCGGATCGTCTGCGTTGGt  <  1:1859020/69‑1 (MQ=255)
                                                                    |                                                                   
TACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCGCAGCAGCCGACCGTGCCGGATCGTCTGCGTTGGT  >  minE/2973848‑2973983

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: