Predicted mutation
evidence seq id position mutation annotation gene description
RA minE 2,367,537 C→T intergenic (‑85/‑105) pepQ ← / → fadB proline dipeptidase/fused 3‑hydroxybutyryl‑CoA epimerase/delta(3)‑cis‑delta(2)‑trans‑enoyl‑CoA isomerase/enoyl‑CoA hydratase and 3‑hydroxyacyl‑CoA dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,367,5370CT100.0% 87.8 / NA 27intergenic (‑85/‑105)pepQ/fadBproline dipeptidase/fused 3‑hydroxybutyryl‑CoA epimerase/delta(3)‑cis‑delta(2)‑trans‑enoyl‑CoA isomerase/enoyl‑CoA hydratase and 3‑hydroxyacyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (16/11);  total (16/11)

GACCTTAAGTCTCCGCATCTTAGCACATCGTTCATCCAGAGCGTGATTTCTGCCGAGCGTGATCAGATCGGCATTTCTTTAATCTTTTGTTTGCATATTTTTAACACAAAATACACACTTCGACTCATC  >  minE/2367474‑2367602
                                                               |                                                                 
gACCTTAAGTCTCCGCATCTTAGCACATCGTTCATCCAGAGCGTGATTTCTGCCGAGCGTGATTAGATCgg                                                            >  1:979797/1‑71 (MQ=255)
  ccTTAAGTCTCCGCATCTTAGCACATCGTTCATCCAGAGCGTGATTTCTGCCGAGCGTGATTAGATCGGCa                                                          >  1:2095058/1‑71 (MQ=255)
   cTTAAGTCTCCGCATCTTAGCACATCGTTCATCCAGAGCGTGATTTCTGCCGAGCGTGATTAGATCGGCAt                                                         <  1:1019620/71‑1 (MQ=255)
   cTTAAGTCTCCGCATCTTAGCACATCGTTCATCCAGAGCGTGATTTCTGCCGAGCGTGATTAGATCGGCAt                                                         <  1:237801/71‑1 (MQ=255)
     tAAGTCTCCGCATCTTAGCACATCGTTCATCCAGAGCGTGATTTCTGCCGAGCGTGATTAGATCGGCAttt                                                       >  1:2946625/1‑71 (MQ=255)
           tcCGCATCTTAGCACATCGTTCATCCAGAGCGTGATTTCTGCCGAGCGTGATTAGATCGGCATTTCTTTa                                                  >  1:2243556/1‑70 (MQ=255)
              gCATCTTAGCACATCGTTCATCCAGAGCGTGATTTCTGCCGAGCGTGATTAGATCGGCATTTCtt                                                    >  1:876275/1‑65 (MQ=255)
              gCATCTTAGCACATCGTTCATCCAGAGCGTGATTTCTGCCGAGCGTGATTAGATCGGCATTTCtt                                                    >  1:677922/1‑65 (MQ=255)
                  cTTAGCACATCGTTCATACAGAGCGTGATTTCTGCCGAGCGTGATTAGATCGGCa                                                          >  1:2125818/1‑55 (MQ=255)
                    tAGCACATCGTTCATCCAGAGCGTGATTTCTGCCGAGCGTGATTAGATCGGCATTTCTTTAATCTtttgtt                                        >  1:1850958/1‑71 (MQ=255)
                     aGCACATCGTTCATCCAGAGCGTGATTTCTGCCGAGCGTGATTAGATCGGCATTTCTTTAATCtttt                                           <  1:1020452/67‑1 (MQ=255)
                     aGCACATCGTTCATCCAGAGCGTGATTTCTGCCGAGCGTGATTAGATCGGCATTTCTTTAATCTtttgt                                         <  1:3177712/69‑1 (MQ=255)
                     aGCACATCGTTCATCCAGAGCGTGATTTCTGCCGAGCGTGATTAGATCGGCATTTCTTTAATCTtttgt                                         <  1:1141509/69‑1 (MQ=255)
                     aGCACATCGTTCATCCAGAGCGTGATTTCTGCCGAGCGTGATTAGATCGGCATTTCTTTAATCTtttg                                          <  1:2688625/68‑1 (MQ=255)
                           tCGTTCATCCAGAGCGTGATTTCTGCCGAGCGTGATTAGATCGGCATTTCTTTAATCTTTTGTTTGCatat                                 <  1:3046236/71‑1 (MQ=255)
                           tCGTTCATCCAGAGCGTGATTTCTGCCGAGCGTGATTAGATCGGCATTTCTTTAATCTTTTGTTTGCata                                  >  1:1875903/1‑70 (MQ=255)
                                  tCCAGAGCGTGATTTCTGCCGAGCGTGATTAGATCGGCATTTCTTTAATCTtttg                                          <  1:1537677/55‑1 (MQ=255)
                                     agagCGTGATTTCTGCCGAGCGTGATTAGATCGGCATTTCTTTAATCTTTTGTTTGCATATTTTTAacac                        >  1:2387317/1‑70 (MQ=255)
                                        gCGTGATTTCTGCCGAGCGTGATTAGATCGGCATTTCTTTAATCTTTTGTTTGCATATTTTTAACACaaaa                    <  1:831303/71‑1 (MQ=255)
                                              tttCTGCCGAGCGTGATTAGATCGGCATTTCTTTAATCTTTTGTTTGCATATTTTTAACACAAAATacaca              >  1:1683810/1‑71 (MQ=255)
                                               ttCTGCCGAGCGTGATTAGATCGGCATTTCTTTAATCTTTTGTTTGCATATTTTTAAcaca                       >  1:2452543/1‑61 (MQ=255)
                                                tCTGCCGAGCGTGATTAGATCGGCATTTCTTTAATCTTTTGTTTGCATATTTTTAACACAAAATACACACt            >  1:391621/1‑71 (MQ=255)
                                                tCTGCCGAGCGTGATTAGATCGGCATTTCTTTAATCTTTTGTTTGCATATTTTTAACACAAAATACACACt            >  1:1529794/1‑71 (MQ=255)
                                                tCTGCCGAGCGTGATTAGATCGGCATTTCTTTAATCTTTTGTTTGCATATTTTTAACACAAAATACACACt            >  1:1052307/1‑71 (MQ=255)
                                                      gAGCGTGATTAGATCGGCATTTCTTTAATCTTTTGTTTGCATATTTTTAACACAAAATACACACTTCGACt      >  1:1847716/1‑71 (MQ=255)
                                                       aGCGTGATTAGATCGGCATTTCTTTAATCTTTTGTTTGCATATTTTTAACACAAAATACACACTTCGACTc     <  1:2314643/71‑1 (MQ=255)
                                                          gTGATTAGATCGGCATTTCTTTAATCTTTTGTTTGCATATTTTTAACACAAAATACACACTTCGACTCATc  <  1:778741/71‑1 (MQ=255)
                                                               |                                                                 
GACCTTAAGTCTCCGCATCTTAGCACATCGTTCATCCAGAGCGTGATTTCTGCCGAGCGTGATCAGATCGGCATTTCTTTAATCTTTTGTTTGCATATTTTTAACACAAAATACACACTTCGACTCATC  >  minE/2367474‑2367602

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: