Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,390,648 C→T 36.6% L260L (CTC→CTT nfo → endonuclease IV with intrinsic 3'‑5' exonuclease activity

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,390,6480CT36.6% 58.2 / 53.4 90L260L (CTC→CTTnfoendonuclease IV with intrinsic 3'‑5' exonuclease activity
Reads supporting (aligned to +/- strand):  ref base C (24/33);  new base T (14/19);  total (38/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.29e-01

TATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAAGAGATCGCCTGGCTGAAAGCGCAACAA  >  minE/1390582‑1390705
                                                                  |                                                         
tatCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGGATTCCGCTGATCCTCGa                                                         <  1:3178853/69‑1 (MQ=255)
 atCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTTGaa                                                        <  1:1720886/69‑1 (MQ=255)
 atCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTTGaa                                                        <  1:87606/69‑1 (MQ=255)
 atCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTTGaa                                                        <  1:2352895/69‑1 (MQ=255)
 atCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTTGaa                                                        <  1:2375357/69‑1 (MQ=255)
 atCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTTGaa                                                        <  1:2093202/69‑1 (MQ=255)
 atCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTTGaa                                                        <  1:206836/69‑1 (MQ=255)
 atCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTTGaa                                                        <  1:2518815/69‑1 (MQ=255)
 atCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTTGaa                                                        <  1:1947427/69‑1 (MQ=255)
 atCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTTGaa                                                        <  1:2667012/69‑1 (MQ=255)
 atCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTTGaa                                                        <  1:1820829/69‑1 (MQ=255)
 atCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTTGaa                                                        <  1:1775744/69‑1 (MQ=255)
 atCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTTGaa                                                        <  1:3099730/69‑1 (MQ=255)
 atCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTTGaa                                                        <  1:1403373/69‑1 (MQ=255)
 atCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTTGaa                                                        <  1:3425898/69‑1 (MQ=255)
 atCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTTGaa                                                        <  1:3177025/69‑1 (MQ=255)
 atCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGa                                                         <  1:3111002/68‑1 (MQ=255)
 atCCGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTTGaa                                                        <  1:417422/69‑1 (MQ=255)
  tCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGaa                                                        >  1:2997190/1‑68 (MQ=255)
  tCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGaa                                                        >  1:2461240/1‑68 (MQ=255)
  tCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGaa                                                        >  1:1985528/1‑68 (MQ=255)
             gCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCg                                            <  1:1956191/69‑1 (MQ=255)
             gCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCg                                            <  1:320417/69‑1 (MQ=255)
             gCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCg                                            <  1:1203382/69‑1 (MQ=255)
             gCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCg                                            <  1:2924243/69‑1 (MQ=255)
             gCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCg                                            <  1:1928741/69‑1 (MQ=255)
             gCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCg                                            <  1:1735084/69‑1 (MQ=255)
             gCGTTCAGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCg                                            <  1:1133912/69‑1 (MQ=255)
              cGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCAGATCCTCGAAACCATCAACCCg                                            >  1:2917577/1‑68 (MQ=255)
                 tCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTTGaa                                                        <  1:2048927/53‑1 (MQ=255)
                  ccGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGatat                                       >  1:2716616/1‑69 (MQ=255)
                  ccGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGatat                                       >  1:1200291/1‑69 (MQ=255)
                  ccGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGatat                                       >  1:1647533/1‑69 (MQ=255)
                    gCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATc                                      <  1:492224/68‑1 (MQ=255)
                    gCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATc                                      <  1:3129091/68‑1 (MQ=255)
                    gCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATc                                      <  1:776025/68‑1 (MQ=255)
                    gCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATc                                      <  1:777913/68‑1 (MQ=255)
                    gCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATc                                      <  1:1356256/68‑1 (MQ=255)
                    gCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATc                                      <  1:1936034/68‑1 (MQ=255)
                    gCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATc                                      <  1:1467917/68‑1 (MQ=255)
                          tCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGc                                 >  1:832425/1‑67 (MQ=255)
                            aTGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTTGaa                                                        <  1:856396/42‑1 (MQ=255)
                               ctaGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAagag                          <  1:46833/67‑1 (MQ=255)
                               cAAGACGACCGTTTCGACGGCATTCCGCTGATCCTTGaa                                                        <  1:810947/39‑1 (MQ=255)
                               cAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAagag                          <  1:2303453/69‑1 (MQ=255)
                               cAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAagag                          <  1:1031683/69‑1 (MQ=255)
                               cAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAagag                          <  1:1002243/69‑1 (MQ=255)
                               cAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAagag                          <  1:3256440/69‑1 (MQ=255)
                               cAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAagag                          <  1:985693/69‑1 (MQ=255)
                               cAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAagag                          <  1:975239/69‑1 (MQ=255)
                               cAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAagag                          <  1:2281445/69‑1 (MQ=255)
                               cAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAagag                          <  1:559389/69‑1 (MQ=255)
                               cAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAagag                          <  1:890314/69‑1 (MQ=255)
                               cAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAagag                          <  1:1895662/69‑1 (MQ=255)
                               cAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAagag                          <  1:247643/69‑1 (MQ=255)
                               cAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAagag                          <  1:2060035/69‑1 (MQ=255)
                               cAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAagag                          <  1:2149466/69‑1 (MQ=255)
                                   acgacCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCa                                                 >  1:2700974/1‑42 (MQ=255)
                                   acgacCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAAGAGATc                       >  1:1192603/1‑68 (MQ=255)
                                   acgacCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAAGAGATCg                      >  1:1622850/1‑69 (MQ=255)
                                   acgacCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAAGAGATCg                      >  1:477550/1‑69 (MQ=255)
                                   acgacCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAAGAGATCg                      >  1:1275796/1‑69 (MQ=255)
                                   acgacCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAAGAGATCg                      >  1:2126590/1‑69 (MQ=255)
                                   acgacCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAAGAGATCg                      >  1:428858/1‑69 (MQ=255)
                                   acgacCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAAGAGATCg                      >  1:1062696/1‑69 (MQ=255)
                                        cGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAagag                          <  1:422347/60‑1 (MQ=255)
                                            tCGACGGCATTCCGCTGATCCTTGAAACCATCAACCCGGatat                                       >  1:1418956/1‑43 (MQ=255)
                                            tCGACGGCATTCCGCTGATCCTTGAAACCATCAACCCGGATATCTGGGCAGAAGAGATCGCCTGGCt               >  1:678145/1‑67 (MQ=255)
                                            tCGACGGCATTCCGCTGATCCTTGAAACCATCAACCCGGATATCTGGGCAGAAGAGATCGCCTGGCt               >  1:763742/1‑67 (MQ=255)
                                            tCGACGGCATTCCGCTGATCCTTGAAACCATCAACCCGGATATCTGGGCAGAAGAGATCGCCTGGCt               >  1:971665/1‑67 (MQ=255)
                                            tCGACGGCATTCCGCTGATCCTTGAAACCATCAACCCGGATATCTGGGCAGAAGAGATCGCCTGGCt               >  1:778409/1‑67 (MQ=255)
                                            tCGACGGCATTCCGCTGATCCTTGAAACCATCAACCCGGATATCTGGGCAGAAGAGATCGCCTGGCt               >  1:88409/1‑67 (MQ=255)
                                            tCGACGGCATTCCGCTGATCCTTGAAACCATCAACCCGGATATCTGGGCAGAAGAGATCGCCTGGCt               >  1:991324/1‑67 (MQ=255)
                                            tCGACGGCATTCCGCTGATCCTTGAAACCATCAACCCGGATATCTGGGCAGAAGAGATCGCCTGGCt               >  1:126468/1‑67 (MQ=255)
                                            tCGACGGCATTCCGCTGATCCTTGAAACCATCAACCCGGATATCTGGGCAGAAGAGATCGCCTGGCt               >  1:2924747/1‑67 (MQ=255)
                                            tCGACGGCATTCCGCTGATCCTTGAAACCATCAACCCGGATATCTGGGCAGAAGAGATCGCCTGGCt               >  1:2499878/1‑67 (MQ=255)
                                            tCGACGGCATTCCGCTGATCCTTGAAACCATCAACCCGGATATCTGGGCAGAAGAGATCGCCTGGCt               >  1:2416738/1‑67 (MQ=255)
                                            tCGACGGCATTCCGCTGATCCTTGAAACCATCAACCCGGATATCTGGGCAGAAGAGATCGCCTGGCt               >  1:2240032/1‑67 (MQ=255)
                                            tCGACGGCATTCCGCTGATCCTTGAAACCATCAACCCGGATATCTGGGCAGAAGAGATCGCCTGGCt               >  1:1998527/1‑67 (MQ=255)
                                            tCGACGGCATTCCGCTGATCCTTGAAACCATCAACCCGGATATCTGGGCAGAAGAGATCACCTGGCt               >  1:3363824/1‑67 (MQ=255)
                                            tCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAAGAGATCGCCTGGCTGa             <  1:2516321/69‑1 (MQ=255)
                                                  gCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAAGAGATCGCCTGGCTGAAAgcgc       >  1:2850813/1‑69 (MQ=255)
                                                       ccGCTGATCCTCGAAACCATCAACCCGGATATCTTGGCAGAAGAGATCGCCTGGCTGAAAGCGcaacaa  >  1:1292515/1‑69 (MQ=255)
                                                       ccGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAAGAGATCGCCTGGCTGAAAGCGcaacaa  >  1:1138284/1‑69 (MQ=255)
                                                       ccGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAAGAGATCGCCTGGCTGAAAGCGcaacaa  >  1:443864/1‑69 (MQ=255)
                                                       ccGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAAGAGATCGCCTGGCTGAAAGCGcaacaa  >  1:1571974/1‑69 (MQ=255)
                                                       ccGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAAGAGATCGCCTGGCTGAAAGCGcaacaa  >  1:1807673/1‑69 (MQ=255)
                                                       ccGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAAGAGATCGCCTGGCTGAAAGCGcaacaa  >  1:3406712/1‑69 (MQ=255)
                                                       ccGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAAGAGATCGCCTGGCTGAAAGCGcaacaa  >  1:379212/1‑69 (MQ=255)
                                                            gATCCTCGAAACCATCAACCCGGATATCTGGGCAGAagag                          <  1:984449/40‑1 (MQ=255)
                                                                  |                                                         
TATCGGTCATGATGCGTTCCGCTGGATCATGCAAGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAAGAGATCGCCTGGCTGAAAGCGCAACAA  >  minE/1390582‑1390705

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: