Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,237,256 C→T 48.1% T267T (ACC→ACT purD → phosphoribosylglycinamide synthetase phosphoribosylamine‑glycine ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,237,2560CT48.1% 4.7 / 68.1 79T267T (ACC→ACTpurDphosphoribosylglycinamide synthetase phosphoribosylamine‑glycine ligase
Reads supporting (aligned to +/- strand):  ref base C (26/15);  new base T (9/29);  total (35/44)
Fisher's exact test for biased strand distribution p-value = 5.97e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.17e-01

CCATGGAACGTATTATCTGGCCAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACCGGTTTTCTCTACGCGGGCCTGATGATCGACAAACAGGGCAATCCGAAGGTTATTGAATTTA  >  minE/2237195‑2237317
                                                             |                                                             
cgaTGGAACTTATCATCTGGCCAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACTGGTTTtc                                                        <  1:3413730/67‑1 (MQ=255)
cgaTGGAACGTATCATCTTGCCAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACTGGTTTtc                                                        <  1:1767025/67‑1 (MQ=255)
cgaTGGAACGTATCATCTGGCCAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACTGGTTTtc                                                        <  1:316764/67‑1 (MQ=255)
cgaTGGAACGTATCATCTGGCCAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACTGGTTTtc                                                        <  1:2296628/67‑1 (MQ=255)
cgaTGGAACGTATCATCTGGCCAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACTGGTTTtc                                                        <  1:625602/67‑1 (MQ=255)
cgaTGGAACGTATCATCTGGCCAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACTGGTTTtc                                                        <  1:3009562/67‑1 (MQ=255)
cgaTGGAACGTATCATCTGGCCAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACTGGTTTtc                                                        <  1:1543870/67‑1 (MQ=255)
cgaTGGAACGTATCATCTGGCCAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACTGGTTTtc                                                        <  1:1539981/67‑1 (MQ=255)
cgaTGGAACGTATCATCTGGCCAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACTGGTTTtc                                                        <  1:1466786/67‑1 (MQ=255)
cgaTGGAACGTATCATCTGGCCAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACTGGTTTtc                                                        <  1:2381697/67‑1 (MQ=255)
cgaTGGAACGTATCATCTGGCCAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACTGGTTTtc                                                        <  1:2747145/67‑1 (MQ=255)
cgaTGGAACGTATCATCTGGCCAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACTGGTTTtc                                                        <  1:1392844/67‑1 (MQ=255)
cgaTGGAACGTATCATCTGGCCAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACTGGTTTtc                                                        <  1:3169456/67‑1 (MQ=255)
cgaTGGAACGTATCATCTGGCCAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACTGGTTTtc                                                        <  1:1300613/67‑1 (MQ=255)
cgaTGGAACGTATCATCTGGCCAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACTGGTTTtc                                                        <  1:2724498/67‑1 (MQ=255)
ccATGGAACGTATTATCTGGCCAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACCGGTTTtc                                                        <  1:2328584/69‑1 (MQ=255)
ccATGGAACGTATTATCTGGCCAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACCGGTTTtc                                                        <  1:1742077/69‑1 (MQ=255)
 gaTGGAACGTATCATCTGGCCAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACTGGTTTtc                                                        <  1:2374016/67‑1 (MQ=255)
 gaTGGAACGTATCATCTGGCCAACCGGGAAAGGCATGGCGGCGGAAGGCAACACCTACACTGGTTTtc                                                        <  1:2691266/67‑1 (MQ=255)
  aTGGAACGTATTATCTGGCCAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACCGGTTTtctc                                                      <  1:2224677/69‑1 (MQ=255)
  aTGGAACGTATTATCTGGCCAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACCGGTTTtctc                                                      <  1:2516787/69‑1 (MQ=255)
   tGGAACGTATTATCTGGCCAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACCGGTTTtctc                                                      <  1:803330/68‑1 (MQ=255)
    ggAACGTATTATCTTGCCAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACCGGTTTTCTCTa                                                    <  1:1573690/69‑1 (MQ=255)
    ggAACGTATTATCTGGCCAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACCGGTTTTCTCTa                                                    <  1:3174764/69‑1 (MQ=255)
              atCTGGCCAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACCGGTTTTCTCTACGCGGGCCtg                                          >  1:177560/1‑69 (MQ=255)
              atCTGGCCAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACCGGTTTTCTCTACGCGGGCCtg                                          >  1:3138806/1‑69 (MQ=255)
              atCTGGCCAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACCGGTTTTCTCTACGCGGGCCtg                                          >  1:1821432/1‑69 (MQ=255)
              atCTGGCCAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACCGGTTTTCTCTACGCGGGCCtg                                          >  1:402404/1‑69 (MQ=255)
              atCTGGCCAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACCGGTTTTCTCTACGCGGGCCtg                                          >  1:697330/1‑69 (MQ=255)
              atCTGGCCAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACCGGTTTTCTCTACGCGGGCCtg                                          >  1:1365191/1‑69 (MQ=255)
              atCTGGCCAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACCGGTTTTCTCTACGCGGGCCtg                                          >  1:922085/1‑69 (MQ=255)
              atCTGGCCAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACCGGTTTTCTCTACGCGGGCCt                                           <  1:1195699/68‑1 (MQ=255)
              atCTGGCCAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACAACGGTTTTCTCTACGCGGGCCtg                                          >  1:612252/1‑69 (MQ=255)
                     cAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACTGGTTTTCTCTACGCTGGCCTGATGATc                                    <  1:727317/68‑1 (MQ=255)
                     cAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACTGGTTTTCTCTACGCGGGCCTGATGATc                                    <  1:3161769/68‑1 (MQ=255)
                     cAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACTGGTTTTCTCTACGCGGGCCTGATGATc                                    <  1:3116437/68‑1 (MQ=255)
                     cAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACTGGTTTTCTCTACGCGGGCCTGATGATc                                    <  1:97953/68‑1 (MQ=255)
                     cAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACTGGTTTTCTCTACGCGGGCCTGATGATc                                    <  1:3767/68‑1 (MQ=255)
                     cAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACTGGTTTTCTCTACGCGGGCCTGATGATc                                    <  1:2989194/68‑1 (MQ=255)
                     cAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACTGGTTTTCTCTACGCGGGCCTGATGATc                                    <  1:72294/68‑1 (MQ=255)
                     cAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACTGGTTTTCTCTACGCGGGCCTGATGATc                                    <  1:1737031/68‑1 (MQ=255)
                     cAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACTGGTTTTCTCTACGCGGGCCTGATGATc                                    <  1:232747/68‑1 (MQ=255)
                     cAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACTGGTTTTCTCTACGCGGGCCTGATGATc                                    <  1:2087831/68‑1 (MQ=255)
                     cAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACTGGTTTTCTCTACGCGGGCCTGATGATc                                    <  1:2029630/68‑1 (MQ=255)
                     cAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACTGGTTTTCTCTACGCGGGCCTGATGATc                                    <  1:2012588/68‑1 (MQ=255)
                      aaCCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACCGGTTTTCTCTACGCGGGCCTGATGATCg                                   >  1:1453894/1‑68 (MQ=255)
                      aaCCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACCGGTTTTCTCTACGCGGGCCTGATGATCg                                   >  1:3331799/1‑68 (MQ=255)
                          gTGAAAGGCATGGCGGCGGAAGGCAACACCTACACCGGTTTTCTCTACGCGGGCCTGATGATCGACaaa                              <  1:920067/69‑1 (MQ=255)
                          gTGAAAGGCATGGCGGCGGAAGGCAACACCTACACCGGTTTTCTCTACGCGGGCCTGATGATCGACaaa                              <  1:2686332/69‑1 (MQ=255)
                          gTGAAAGGCATGGCGGCGGAAGGCAACACCTACACCGGTTTTCTCTACGCGGGCCTGATGATCGACaaa                              <  1:3278146/69‑1 (MQ=255)
                          gTGAAAGGCATGGCGGCGGAAGGCAACACCTACACCGGTTTTCTCTACGCGGGCCTGATGATCGACaaa                              <  1:585971/69‑1 (MQ=255)
                                ggCATGGCGGCGGAAGGCAACACCTACACCGGTTTTCTCTACGCGGGCCtgatg                                       >  1:945704/1‑54 (MQ=255)
                                ggCATGGCGGCGGAAGGCAACACCTACACCGGTTTTCTCTACGCGGGCCtgatg                                       >  1:1938698/1‑54 (MQ=255)
                                     ggcggcGGAAGGCAACACCTACACTGGTTTTCTCTACGCGGGCCTGATGATCGACAAACAGGGCAAt                     >  1:3433055/1‑67 (MQ=255)
                                     ggcggcGGAAGGCAACACCTACACTGGTTTTCTCTACGCGGGCCTGATGATCGACAAACAGGGCAATc                    >  1:385004/1‑68 (MQ=255)
                                     ggcggcGGAAGGCAACACCTACACTGGTTTTCTCTACGCGGGCCTGATGATCGACAAACAGGGCAATc                    >  1:1609250/1‑68 (MQ=255)
                                     ggcggcGGAAGGCAACACCTACACTGGTTTTCTCTACGCGGGCCTGATGATCGACAAACAGGGCAATc                    >  1:1607183/1‑68 (MQ=255)
                                     ggcggcGGAAGGCAACACCTACACTGGTTTTCTCTACGCGGGCCTGATGATCGACAAACAGGGCAATc                    >  1:1602593/1‑68 (MQ=255)
                                     ggcggcGGAAGGCAACACCTACACTGGTTTTCTCTACGCGGGCCTGATGATCGACAAACAGGGCAATc                    >  1:2948436/1‑68 (MQ=255)
                                     ggcggcGGAAGGCAACACCTACACTGGTTTTCTCTACGCGGGCCTGATGATCGACAAACAGGGCAATc                    >  1:2845646/1‑68 (MQ=255)
                                     ggcggcGGAAGGCAACACCTACACTGGTTTTCTCTACGCGGGCCTGATGATCGACAAACAGGGCAATc                    >  1:1238327/1‑68 (MQ=255)
                                     ggcggcGGAAGGCAACACCTACACTGGTTTTCTCTACGCGGGCCTGATGATCGACAAACAGGGCAATc                    >  1:410654/1‑68 (MQ=255)
                                           ggAAGGCAACACCTACACCGGTTTTCTCTACGCGGGCCTGATGATCGACAAACAGGGCAATCCGAAGGt             <  1:2083197/69‑1 (MQ=255)
                                              aGGCAACACCTACACCGTTTTTCTCTACGCGGGCCTGATTATCGACaaa                              <  1:109093/49‑1 (MQ=255)
                                               ggCAACACCTACACCGGTTTTCTCTACGCGGGCCTGATGATCGACAAACAGGGCAATCCGAAGGTTAt          >  1:457574/1‑68 (MQ=255)
                                                 cAACACCTACACCGGTTTTCTCTACGCGGGCCTGATGATCGACa                                >  1:629919/1‑44 (MQ=255)
                                                 cAACACCTACACCGGTTTTCTCTACGCGGGCCTGATGATCGACa                                >  1:1780586/1‑44 (MQ=255)
                                                 cAACACCTACACCGGTTTTCTCTACGCGGGCCTGATGATCGACa                                >  1:1621190/1‑44 (MQ=255)
                                                 cAACACCTACACCGGTTTTCTCTACGCGGGCCTGATGATCGACa                                >  1:362252/1‑44 (MQ=255)
                                                 cAACACCTACACCGGTTTTCTCTACGCGGGCCTGATGATCGACa                                >  1:2528163/1‑44 (MQ=255)
                                                 cAACACCTACACCGGTTTTCTCTACGCGGGCCTGATGATCGACa                                >  1:2704454/1‑44 (MQ=255)
                                                     acCTACACCGGTTTTCTCTACGCGGGCCTGATGATCGACAAACAGGGCAATCCGAAGGTTATTGAAttt   >  1:2074734/1‑69 (MQ=255)
                                                     acCTACACCGGTTTTCTCTACGCGGGCCTGATGATCGACAAACAGGGCAATCCGAAGGTTATTGAAttt   >  1:3076071/1‑69 (MQ=255)
                                                     acCTACACCGGTTTTCTCTACGCGGGCCTGATGATCGACAAACAGGGCAATCCGAAGGTTATTGAAttt   >  1:1305898/1‑69 (MQ=255)
                                                     acCTACACCGGTTTTCTCTACGCGGGCCTGATGATCGACAAACAGGGCAATCCGAAGGTTATTGAAttt   <  1:2766999/69‑1 (MQ=255)
                                                     acCTACACCGGTTTTCTCTACGCGGGCCTGATGATCGACAAACAGGGCAATCCGAAGGTTATTGAAttt   >  1:2458426/1‑69 (MQ=255)
                                                     acCTACACCGGTTTTCTCTACGCGGGCCTGATGATCGACAAACAGGGCAATCCGAAGGTTATTGAAtt    >  1:203328/1‑68 (MQ=255)
                                                      ccTACACCGGTTTTCTCTACGCGGGCCTGATGATCGACAAACAGGGCAATCCGAAGGTTATTGAATTTa  >  1:3011119/1‑69 (MQ=255)
                                                      ccTACACCGGTTTTCTCTACGCGGGCCTGATGATCGACAAACAGGGCAATCCGAAGGTTATTGAATTTa  >  1:1407553/1‑69 (MQ=255)
                                                             |                                                             
CCATGGAACGTATTATCTGGCCAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACCGGTTTTCTCTACGCGGGCCTGATGATCGACAAACAGGGCAATCCGAAGGTTATTGAATTTA  >  minE/2237195‑2237317

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: