Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,629,654 G→A 37.8% D31D (GAC→GAT yffH ← predicted NUDIX hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,629,6540GA37.8% 34.6 / 41.9 66D31D (GAC→GATyffHpredicted NUDIX hydrolase
Reads supporting (aligned to +/- strand):  ref base G (30/11);  new base A (8/17);  total (38/28)
Fisher's exact test for biased strand distribution p-value = 1.82e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

AGGAGGATCGTCGCACCATTGCCGCGATCGTACACTTCACGTTTATGGCGGATAACTTCGCCGTCTTTGCGCGTGAGATCGTAAGTAATGTTGTGCAGGGTGAAATAGTTATCGGAGAGAATTTTGTCTTT  >  minE/1629592‑1629722
                                                              |                                                                    
aggaggATCGTCGCACCATTACCGCGATCGTAAACTTCACGTTTATGGCGGATAACTTCGCCATCttt                                                                 >  1:838985/1‑68 (MQ=255)
aggaggATCGTCGCACCATTACCGCGATCGTAAACTTCACGTTTATGGCGGATAACTTCGCCATCTTTg                                                                >  1:1055048/1‑69 (MQ=255)
aggaggATCGTCGCACCATTACCGCGATCGTAAACTTCACGTTTATGGCGGATAACTTCGCCATCTTTg                                                                >  1:1211259/1‑69 (MQ=255)
aggaggATCGTCGCACCATTACCGCGATCGTAAACTTCACGTTTATGGCGGATAACTTCGCCATCTTTg                                                                >  1:595694/1‑69 (MQ=255)
aggaggATCGTCGCACCATTACCGCGATCGTAAACTTCACGTTTATGGCGGATAACTTCGCCATCTTTg                                                                >  1:306405/1‑69 (MQ=255)
aggaggATCGTCGCACCATTACCGCGATCGTAAACTTCACGTTTATGGCGGATAACTTCGCCATCTTTg                                                                >  1:2861500/1‑69 (MQ=255)
aggaggATCGTCGCACCATTACCGCGATCGTAAACTTCACGTTTATGGCGGATAACTTCGCCATCTTTg                                                                >  1:237746/1‑69 (MQ=255)
aggaggATCGTCGCACCATTACCGCGATCGTAAACTTCACGTTTATGGCGGATAACTTCGCCATCTTTg                                                                >  1:2023582/1‑69 (MQ=255)
               ccATTGCCGCGATCGTACACTTCACGTTTATGGCGGATAACTTCGCCGTCTTTGCGCGTGAGATCgta                                                  >  1:2722622/1‑68 (MQ=255)
                   tGCCGCGATCGTACACTTGACGTTTATGGCGGATAACTTCGCCGTCTTTGCGCGt                                                           <  1:38613/55‑1 (MQ=255)
                   tGCCGCGATCGTACACTTCACGTTTATGGCGGATAACTTCGCCGTCTTTGCGCGt                                                           <  1:1186388/55‑1 (MQ=255)
                   tGCCGCGATCGTACACTTCACGTTTATGGCGGATAACTTCGCCGTCTTTGCGCGt                                                           <  1:764681/55‑1 (MQ=255)
                   tGCCGCGATCGTACACTTCACGTTTATGGCGGATAACTTCGCCGTCTTTGCGCGt                                                           <  1:3107985/55‑1 (MQ=255)
                   tGCCGCGATCGTACACTTCACGTTTATGGCGGATAACTTCGCCGTCTTTGCGCGt                                                           <  1:2576370/55‑1 (MQ=255)
                   tGCCGCGATCGTACACTTCACGTTTATGGCGGATAACTTCGCCGTCTTTGCGCGt                                                           <  1:195090/55‑1 (MQ=255)
                   tGCCGCGATCGTACACTTCACGTTTATGGCGGATAACTTCGCCGTCTTTGCGCGt                                                           <  1:1456822/55‑1 (MQ=255)
                    gCCGCGATCGTACACTTCACGTTTATGGCGGATAACTTCGCCGTCTTTGCGCGTGAGATCGTAAGTAAt                                            >  1:533264/1‑69 (MQ=255)
                    gCCGCGATCGTACACTTCACGTTTATGGCGGATAACTTCGCCGTCTTTGCGCGTGAGATCGTAAGTAAt                                            >  1:3095670/1‑69 (MQ=255)
                    gCCGCGATCGTACACTTCACGTTTATGGCGGATAACTTCGCCGTCTTTGCGCGTGAGATCGTAAGTAAt                                            >  1:3124948/1‑69 (MQ=255)
                    gCCGCGATCGTACACTTCACGTTTATGGCGGATAACTTCGCCGTCTTTGCGCGTGAGATCGTAAGTAAt                                            >  1:2434388/1‑69 (MQ=255)
                    gCCGCGATCGTACACTTCACGTTTATGGCGGATAACTTCGCCGTCTTTGCGCGTGAGATCGTAAGTAAt                                            >  1:77490/1‑69 (MQ=255)
                    gCCGCGATCGTACACTTCACGTTTATGGCGGATAACTTCGCCGTCTTTGCGCGTGAGATCGTAAGTAAt                                            >  1:784841/1‑69 (MQ=255)
                    gCCGCGATCGTACACTTCACGTTTATGGCGGATAACTTCGCCGTCTTTGCGCGTGAGATCGTAAGTAAt                                            >  1:837733/1‑69 (MQ=255)
                    gCCGCGATCGTACACTTCACGTTTATGGCGGATAACTTCGCCGTCTTTGCGCGTGAGATCGTAAGTAAt                                            >  1:1783821/1‑69 (MQ=255)
                    gCCGCGATCGTACACTTCACGTTTATGGCGGATAACTTCGCCGTCTTTGCGCGTGAGATCGTAAGTAAt                                            >  1:844593/1‑69 (MQ=255)
                    gCCGCGATCGTACACTTCACGTTTATGGCGGATAACTTCGCCGTCTTTGCGCGTGAGATCGTAAGTAAt                                            >  1:1202021/1‑69 (MQ=255)
                    gCCGCGATCGTACACTTCACGTTTATGGCGGATAACTTCGCCGTCTTTGCGCGTGAGATCGTAAGTAAt                                            >  1:1795857/1‑69 (MQ=255)
                    gCCGCGATCGTACACTTCACGTTTATGGCGGATAACTTCGCCGTCTTTGCGCGTGAGATCGTAAGTAAt                                            >  1:2098592/1‑69 (MQ=255)
                          aTCGTACACTTCACGTTTATGGCGGATAACTTCGCCGTCTTTGCGCGTGAGATCGTAAGTAATGTtgtg                                      <  1:3028407/69‑1 (MQ=255)
                          aTCGTACACTTCACGTTTATGGCGGATAACTTCGCCGTCTTTGCGCGTGAGATCGTAAGTAATGTtgtg                                      <  1:2419507/69‑1 (MQ=255)
                            cGTACACTTCACGTTTATGGCGGATAACTTCGCCGTCTTTGCGCGTGAGATCGTAAGTAATGTTGTGCa                                    >  1:2524592/1‑69 (MQ=255)
                            cGTACACTTCACGTTTATGGCGGATAACTTCGCCGTCTTTGCGCGTGAGATCGTAAGTAATGTTGTGCa                                    >  1:883082/1‑69 (MQ=255)
                            cGTACACTTCACGTTTATGGCGGATAACTTCGCCGTCTTTGCGCGTGAGATCGTAAGTAATGTTGTGCa                                    <  1:1820057/69‑1 (MQ=255)
                            cGTACACTTCACGTTTATGGCGGATAACTTCGCCGTCTTTGCGCGTGAGATCGTAAGTAATGTTGTGCa                                    >  1:2889211/1‑69 (MQ=255)
                            cGTACACTTCACGTTTATGGCGGATAACTTCGCCGTCTTTGCGCGTGAGATCGTAAGTAATGTTGTGCa                                    >  1:2815092/1‑69 (MQ=255)
                            cGTACACTTCACGTTTATGGCGGATAACTTCGCCGTCTTTGCGCGTGAGATCGTAAGTAATGTTGTGCa                                    >  1:1693794/1‑69 (MQ=255)
                            cGTACACTTCACGTTTATGGCGGATAACTTCGCCGTCTTTGCGCGTGAGATCGTAAGTAATGTTGTGCa                                    >  1:2493097/1‑69 (MQ=255)
                            cGTACACTTCACGTTTATGGCGGATAACTTCGCCGTCTTTGCGCGTGAGATCGTAAGTAATGTTGTGCa                                    >  1:752173/1‑69 (MQ=255)
                            cGTACACTTCACGTTTATGGCGGATAACTTCGCCGTCTTTGCGCGTGAGATCGTAAGTAATGTTGTGCa                                    >  1:2298566/1‑69 (MQ=255)
                            cGTACACTTCACGTTTATGGCGGATAACTTCGCCGTCTTTGCGCGTGAGATCGTAAGTAATGTTGTGCa                                    >  1:789382/1‑69 (MQ=255)
                                    tCACGTTTATGGCGGATAACTTCGCCGTCTTTGCGCGTGAGATCGTAAGTAATGTTGTGCAGGGTGaaa                            <  1:145474/69‑1 (MQ=255)
                                    tCACGTTTATGGCGGATAACTTCGCCGTCTTTGCGCGTGAGATCGTAAGTAATGTTGTGCAGGGTGa                              >  1:522723/1‑67 (MQ=255)
                                    tCACGTTTATGGCGGATAACTTCGCCGTCTTTGCGCGTGAGATCGTAAGTAATGTTGTGCAGGGTGa                              >  1:3334516/1‑67 (MQ=255)
                                    tCACGTTTATGGCGGATAACTTCGCCGTCTTTGCGCGTGAGATCGTAAGTAATGTTGTGCAGGGTGa                              >  1:3024440/1‑67 (MQ=255)
                                    tCACGTTTATGGCGGATAACTTCGCCGTCTTTGCGCGTGAGATCGTAAGTAATGTTGTGCAGGGTGa                              >  1:678055/1‑67 (MQ=255)
                                    tCACGTTTATGGCGGATAACTTCGCCGTCTTTGCGCGTGAGATCGTAAGTAATGTTGTGCAGGGTGa                              >  1:1725797/1‑67 (MQ=255)
                                               gCGGATAACTTCGCCATCTTTGCGCGTGAGATCTTAAGTAATGTTGTGCAGGGTGAAATAGTTATCGga                 <  1:735871/69‑1 (MQ=255)
                                               gCGGATAACTTCGCCATCTTTGCGCGTGAGATCGTAAGTAATGTTGTGCAGGGTGAAATAGTTATCGga                 <  1:1332160/69‑1 (MQ=255)
                                               gCGGATAACTTCGCCATCTTTGCGCGTGAGATCGTAAGTAATGTTGTGCAGGGTGAAATAGTTATCGga                 <  1:1833273/69‑1 (MQ=255)
                                               gCGGATAACTTCGCCATCTTTGCGCGTGAGATCGTAAGTAATGTTGTGCAGGGTGAAATAGTTATCGga                 <  1:1859765/69‑1 (MQ=255)
                                               gCGGATAACTTCGCCATCTTTGCGCGTGAGATCGTAAGTAATGTTGTGCAGGGTGAAATAGTTATCGga                 <  1:1962132/69‑1 (MQ=255)
                                               gCGGATAACTTCGCCATCTTTGCGCGTGAGATCGTAAGTAATGTTGTGCAGGGTGAAATAGTTATCGga                 <  1:2264912/69‑1 (MQ=255)
                                               gCGGATAACTTCGCCATCTTTGCGCGTGAGATCGTAAGTAATGTTGTGCAGGGTGAAATAGTTATCGga                 <  1:3005378/69‑1 (MQ=255)
                                               gCGGATAACTTCGCCATCTTTGCGCGTGAGATCGTAAGTAATGTTGTGCAGGGTGAAATAGTTATCGga                 <  1:2427575/69‑1 (MQ=255)
                                               gCGGATAACTTCGCCATCTTTGCGCGTGAGATCGTAAGTAATGTTGTGCAGGGTGAAATAGTTATCGga                 <  1:1607231/69‑1 (MQ=255)
                                               gCGGATAACTTCGCCATCTTTGCGCGTGAGATCGTAAGTAATGTTGTGCAGGGTGAAATAGTTATCGga                 <  1:427527/69‑1 (MQ=255)
                                               gCGGATAACTTCGCCATCTTTGCGCGTGAGATCGTAAGTAATGTTGTGCAGGGTGAAATAGTTATCGga                 <  1:2925908/69‑1 (MQ=255)
                                                cgtaTAACTTCGCCATCTTTGCGCGTGAGATCGTAAGTAATGTTGTGCAGGGTGAAATAGTTATCGga                 <  1:1509225/65‑1 (MQ=255)
                                                cGGATAACTTCGCCATCTTTGCGCGTGAGATCGTAAGTAATGTTGTGCAGGGTGAAATAGTTATCGga                 <  1:431105/68‑1 (MQ=255)
                                                cGGATAACTTCGCCATCTTTGCGCGTGAGATCGTAAGTAATGTTGTGCAGGGTGAAATAGTTATCGga                 <  1:2392106/68‑1 (MQ=255)
                                                cGGATAACTTCGCCATCTTTGCGCGTGAGATCGTAAGTAATGTTGTGCAGGGTGAAATAGTTATCGga                 <  1:419473/68‑1 (MQ=255)
                                                cGGATAACTTCGCCATCTTTGCGCGTGAGATCGTAAGTAATGTTGTGCAGGGTGAAATAGTTATCGga                 <  1:137839/68‑1 (MQ=255)
                                                cGGATAACTTCGCCATCTTTGCGCGTGAGATCGTAAGTAATGTTGTGCAGGGTGAAATAGTTATCGga                 <  1:210099/68‑1 (MQ=255)
                                                    tAACTTCGCCGTCTTTGCGCGTGAGATCGTAAGTAATGTTGTGCAGGGTGAAATAGTTATCGGAGAGaa            >  1:485348/1‑69 (MQ=255)
                                                    tAACTTCGCCGTCTTTGCGCGTGAGATCGTAAGTAATGTTGTGCAGGGTGAAATAGTTATCGGAGAGaa            >  1:2074395/1‑69 (MQ=255)
                                                    tAACTTCGCCGTCTTTGCGCGTGAGATCGTAAGTAATGTTGTGCAGGGTGAAATAGTTATCGGAGAGaa            >  1:3201313/1‑69 (MQ=255)
                                                              gTCTTTGCGCGTGAGATCGTAAGTAATGTTGTGCAGGGTGAAATAGTTATCGGAGAGAATTTTGTCttt  <  1:124733/69‑1 (MQ=255)
                                                              gTCTTTGCGCGTGAGATCGTAAGTAATGTTGTGCAGGGTGAAATAGTTATCGGAGAGAATTTTGTCttt  >  1:2188347/1‑69 (MQ=255)
                                                              gTCTTTGCGCGTGAGATCGTAAGTAATGTTGTGCAGGGTGAAATAGTTATCGGAGAGAATTTTGTCttt  <  1:836022/69‑1 (MQ=255)
                                                              gTCTTTGCGCGTGAGATCGTAAGTAATGTTGTGCAGGGTGAAATAGTTATCGGAGAGAATTTTGTCttt  <  1:2871249/69‑1 (MQ=255)
                                                              |                                                                    
AGGAGGATCGTCGCACCATTGCCGCGATCGTACACTTCACGTTTATGGCGGATAACTTCGCCGTCTTTGCGCGTGAGATCGTAAGTAATGTTGTGCAGGGTGAAATAGTTATCGGAGAGAATTTTGTCTTT  >  minE/1629592‑1629722

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: