Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,542,575 T→G 46.0% I613L (ATC→CTC)  yfcX ← fused enoyl‑CoA hydratase, 3‑hydroxybutyryl‑CoA epimerase, dodecenoyl‑CoA D‑isomerase, and 3‑hydroxyacyl‑CoA dehydrogenase, subunit of anaerobic fatty acid oxidation complex

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,542,5750TG46.0% ‑6.2 / 41.0 37I613L (ATC→CTC) yfcXfused enoyl‑CoA hydratase, 3‑hydroxybutyryl‑CoA epimerase, dodecenoyl‑CoA D‑isomerase, and 3‑hydroxyacyl‑CoA dehydrogenase, subunit of anaerobic fatty acid oxidation complex
Reads supporting (aligned to +/- strand):  ref base T (9/11);  new base G (11/6);  total (20/17)
Fisher's exact test for biased strand distribution p-value = 3.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.61e-01

AACACAACGTACTGCTTCATTCAGCATCAACATCACACACCGTTCAGCAACCTGCGGTGCGGAGATTCGCCCCTGCCCTTGTGTGCCAATCAGCGGGTAAATGGCGGGATCGACCTGTTTTTTGCTTTTACGCC  >  minE/1542510‑1542643
                                                                 |                                                                    
aaCACAACGTACTGCTTCATTCAGCATCAACATCACACACCGTTCAGCAACCTGCGGTGCGGAGATTCg                                                                   >  1:1049020/1‑69 (MQ=255)
aaCACAACGTACTGCTTCATTCAGCATCAACATCACACACCGTTCAGCAACCTGCGGTGCGGAGATTCg                                                                   >  1:380510/1‑69 (MQ=255)
aaCACAACGTACTGCTTCATTCAGCATCAACATCACACACCGTTCAGCAACCTGCGGTGCGGAGATTCg                                                                   >  1:406281/1‑69 (MQ=255)
aaCACAACGTACTGCTTCATTCAGCATCAACATCACACACCGTTCAGCAACCTGCGGTGCGGAGATTCg                                                                   >  1:3296161/1‑69 (MQ=255)
aaCACAACGTACTGCTTCATTCAGCATCAACATCACACACCGTTCAGCAACCTGCGGTGCGGAGATTCg                                                                   >  1:2700751/1‑69 (MQ=255)
aaCACAACGTACTGCTTCATTCAGCATCAACATCACACACCGTTCAGCAACCTGCGGTGCGGAGATTCg                                                                   >  1:1905980/1‑69 (MQ=255)
aaCACAACGTACTGCTTCATTCAGCATCAACATCACACACCGTTCAGCAACCTGCGGTGCGGAGATTCg                                                                   >  1:929427/1‑69 (MQ=255)
aaCACAACGTACTGCTTCATTCAGCATCAACATCACACACCGTTCAGCAACCTGCGGTGCGGAGATTCg                                                                   >  1:1813086/1‑69 (MQ=255)
  cacaACGTACTGCTTCATTCAGCATCAACATCACACACCGTTCAGCAACCTGCGGTGCGGAGATTCGc                                                                  <  1:393783/68‑1 (MQ=255)
  cacaACGTACTGCTTCATTCAGCATCAACATCACACACCGTTCAGCAACCTGCGGTGCGGAGATTCGc                                                                  <  1:551019/68‑1 (MQ=255)
          aCTGCTTCATTCAGCATCAACATCACACACCGTTCAGCAACCTGCGGTGCGGAGAGTCGCCCCTGCCCt                                                         >  1:2183134/1‑69 (MQ=255)
          aCTGCTTCATTCAGCATCAACATCACACACCGTTCAGCAACCTGCGGTGCGGAGAGTCGCCCCTGCCCt                                                         >  1:2780468/1‑69 (MQ=255)
          aCTGCTTCATTCAGCATCAACATCACACACCGTTCAGCAACCTGCGGTGCGGAGAGTCGCCCCTGCCCt                                                         >  1:2827088/1‑69 (MQ=255)
          aCTGCTTCATTCAGCATCAACATCACACACCGTTCAGCAACCTGCGGTGCGGAGAGTCGCCCCTGCCCt                                                         >  1:150218/1‑69 (MQ=255)
          aCTGCTTCATTCAGCATCAACATCACACACCGTTCAGCAACCTGCGGTGCGGAGAGTCGCCCCTGCCCt                                                         >  1:3233005/1‑69 (MQ=255)
          aCTGCTTCATTCAGCATCAACATCACACACCGTTCAGCAACCTGCGGTGCGGAGAGTCGCCCCTGCCCt                                                         >  1:1887906/1‑69 (MQ=255)
              cTTCATTCAGCATCAACATCACACACCGTTCAGCAACCTGCGGTGCGGAGAGTCGCCCCTGCCCTtgtg                                                     >  1:1346963/1‑69 (MQ=255)
              cTTCATTCAGCATCAACATCACACACCGTTCAGCAACCTGCGGTGCGGAGAGTCGCCCCTGCCCTtgtg                                                     >  1:1530009/1‑69 (MQ=255)
              cTTCATTCAGCATCAACATCACACACCGTTCAGCAACCTGCGGTGCGGAGAGTCGCCCCTGCCCTtgtg                                                     >  1:2457504/1‑69 (MQ=255)
              cTTCATTCAGCATCAACATCACACACCGTTCAGCAACCTGCGGTGCGGAGAGTCGCCCCTGCCCTtgtg                                                     >  1:653415/1‑69 (MQ=255)
              cTTCATTCAGCATCAACATCACACACCGTTCAGCAACCTGCGGTGCGGAGAGTCGCCCCTGCCCTtgt                                                      >  1:1660329/1‑68 (MQ=255)
                            aaCATCACACACCGTTCAGCAACCTGCGGTGCGGAGATTCGCCCCTGCCCTTGTGTGCCAATCAGCgg                                        >  1:1002690/1‑68 (MQ=255)
                                       ccGTTCAGCAACCTGCGGTGCGGAGATTCGCCCCTGCCCTTTTGTGCCAATCAGCGGGTAAATGGCggg                            <  1:354106/69‑1 (MQ=255)
                                       ccGTTCAGCAACCTGCGGTGCGGAGATTCGCCCCTGCCCTTGTGTGCCAATCAGCGGGTAAATGGCggg                            <  1:1105553/69‑1 (MQ=255)
                                       ccGTTCAGCAACCTGCGGTGCGGAGATTCGCCCCTGCCCTTGTGTGCCAATCAGCGGGTAAATGGCggg                            <  1:543430/69‑1 (MQ=255)
                                       ccGTTCAGCAACCTGCGGTGCGGAGATTCGCCCCTGCCCTTGTGTGCCAATCAGCGGGTAAATGGCggg                            <  1:2940720/69‑1 (MQ=255)
                                       ccGTTCAGCAACCTGCGGTGCGGAGATTCGCCCCTGCCCTTGTGTGCCAATCAGCGGGTAAATGGCggg                            <  1:2817549/69‑1 (MQ=255)
                                       ccGTTCAGCAACCTGCGGTGCGGAGATTCGCCCCTGCCCTTGTGTGCCAATCAGCGGGTAAATGGCggg                            <  1:272613/69‑1 (MQ=255)
                                       ccGTTCAGCAACCTGCGGTGCGGAGATTCGCCCCTGCCCTTGTGTGCCAATCAGCGGGTAAATGGCggg                            <  1:936101/69‑1 (MQ=255)
                                       ccGTTCAGCAACCTGCGGTGCGGAGATTCGCCCCTGCCCTTGTGTGCCAATCAGCGGGTAAATGGCggg                            <  1:1457804/69‑1 (MQ=255)
                                       ccGTTCAGCAACCTGCGGTGCGGAGATTCGCCCCTGCCCTTGTGTGCCAATCAGCGGGTAAATGGCggg                            <  1:233247/69‑1 (MQ=255)
                                          ttCAGCAACCTGCGGTGCGGAGAGTCGCCCCTGCCCTTGTGCGCCAATGAGCGGGTAAATGGCGGGATc                         <  1:80198/69‑1 (MQ=255)
                                          ttCAGCAACCTGCGGTGCGGAGAGTCGCCCCTGCCCTTGTGCGCCAATGAGCGGGTAAATGGCGGGATc                         <  1:49859/69‑1 (MQ=255)
                                          ttCAGCAACCTGCGGTGCGGAGAGTCGCCCCTGCCCTTGTGCGCCAATGAGCGGGTAAATGGCGGGATc                         <  1:2029600/69‑1 (MQ=255)
                                          ttCAGCAACCTGCGGTGCGGAGAGTCGCCCCTGCCCTTGTGCGCCAATGAGCGGGTAAATGGCGGGATc                         <  1:84131/69‑1 (MQ=255)
                                          ttCAGCAACCTGCGGTGCGGAGAGTCGCCCCTGCCCTTGTGCGCCAATGAGCGGGTAAATGGCGGGATc                         <  1:1083789/69‑1 (MQ=255)
                                          ttCAGCAACCTGCGGTGCGGAGAGTCGCACCTGCCCTTGTGCGCCAATGAGCGGGTAAATGGCGGGATc                         <  1:3264353/69‑1 (MQ=255)
                                                                 ttCGCCCCTGCCCTTGTGTGCCAATCAGCGGGTAAATGGCGGGATCGACCTGTTTTTTGCTTTTACGcc  <  1:953813/69‑1 (MQ=255)
                                                                 ttCGCCCCTGCCCTTGTGTGCCAATCAGCGGGTAAATGGCGGGATCGACCTGTTTTTTGCTTTTACGcc  <  1:3393105/69‑1 (MQ=255)
                                                                 ttCGCCCCTGCCCTTGTGTGCCAATCAGCGGGTAAATGGCGGGATCGACCTGTTTTTTGCTTTTACGcc  <  1:3388186/69‑1 (MQ=255)
                                                                 ttCGCCCCTGCCCTTGTGTGCCAATCAGCGGGTAAATGGCGGGATCGACCTGTTTTTTGCTTTTACGcc  <  1:2709121/69‑1 (MQ=255)
                                                                 ttCGCCCCTGCCCTTGTGTGCCAATCAGCGGGTAAATGGCGGGATCGACCTGTTTTTTGCTTTTACGcc  <  1:25429/69‑1 (MQ=255)
                                                                 ttCGCCCCTGCCCTTGTGTGCCAATCAGCGGGTAAATGGCGGGATCGACCTGTTTTTTGCTTTTACGcc  <  1:2493824/69‑1 (MQ=255)
                                                                 ttCGCCCCTGCCCTTGTGTGCCAATCAGCGGGTAAATGGCGGGATCGACCTGTTTTTTGCTTTTACGcc  <  1:2062463/69‑1 (MQ=255)
                                                                 ttCGCCCCTGCCCTTGTGTGCCAATCAGCGGGTAAATGGCGGGATCGACCTGTTTTTTGCTTTTACGcc  <  1:2019823/69‑1 (MQ=255)
                                                                 ttCGCCCCTGCCCTTGTGTGCCAATCAGCGGGTAAATGGCGGGATCGACCTGTTTTTTGCTTTTACGcc  <  1:1847397/69‑1 (MQ=255)
                                                                 ttCGCCCCTGCCCTTGTGTGCCAATCAGCGGGTAAATGGCGGGATCGACCTGTTTTTTGCTTTTACGcc  <  1:1711490/69‑1 (MQ=255)
                                                                 ttCGCCCCTGCCCTTGTGTGCCAATCAGCGGGTAAATGGCGGGATCGACCTGTTTTTTGCTTTTACGcc  <  1:1660194/69‑1 (MQ=255)
                                                                 |                                                                    
AACACAACGTACTGCTTCATTCAGCATCAACATCACACACCGTTCAGCAACCTGCGGTGCGGAGATTCGCCCCTGCCCTTGTGTGCCAATCAGCGGGTAAATGGCGGGATCGACCTGTTTTTTGCTTTTACGCC  >  minE/1542510‑1542643

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: