Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 220,157 A→G 41.8% D657D (GAT→GAC fadE ← acyl coenzyme A dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE220,1570AG41.8% 19.9 / 48.2 67D657D (GAT→GACfadEacyl coenzyme A dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (24/15);  new base G (23/5);  total (47/20)
Fisher's exact test for biased strand distribution p-value = 1.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCGGGAAGATCACCACATTCAGCAGCCCGGCAACCACGCGGTTCGGGAAGTTTTGCAGTAAATCATCCATCGCCTGTTCAGCCTGATACAGCGCATCTTGTACGCCCCAGTGCACCAGCGGCAGGTCGG  >  minE/220092‑220221
                                                                 |                                                                
gTCGGGAAGATCACCACATTCAGCAGCCCGGCAACCACGCGGTTCGGGAAGTTTTGCAGTAAATcatcc                                                               <  1:3312239/69‑1 (MQ=255)
   gggAAGATCACCACATTCAGCAGCCCGGCAACCACGCGGTTCGGGAAGTTTTGCAGTAAATcatccatc                                                            >  1:3256516/1‑69 (MQ=255)
     gAAGATCACCACATTCAGCAGCCCGGCAACCACGCGGTTCGGGAAGTTTTGCAGTAAATCGTCCATCg                                                           <  1:1077559/68‑1 (MQ=255)
     gAAGATCACCACATTCAGCAGCCCGGCAACCACGCGGTTCGGGAAGTTTTGCAGTAAATCGTCCATCg                                                           <  1:1228942/68‑1 (MQ=255)
     gAAGATCACCACATTCAGCAGCCCGGCAACCACGCGGTTCGGGAAGTTTTGCAGTAAATCGTCCATCg                                                           <  1:139056/68‑1 (MQ=255)
     gAAGATCACCACATTCAGCAGCCCGGCAACCACGCGGTTCGGGAAGTTTTGCAGTAAATCGTCCATCg                                                           <  1:1417619/68‑1 (MQ=255)
     gAAGATCACCACATTCAGCAGCCCGGCAACCACGCGGTTCGGGAAGTTTTGCAGTAAATCGTCCATCg                                                           <  1:2374486/68‑1 (MQ=255)
       aGATCACCACATTCAGCAGCCCGGCAACCCCGCGGTTCGGGAAGTTTTGCAGTAAATCATCCATCGCCt                                                        >  1:763275/1‑69 (MQ=255)
       aGATCACCACATTCAGCAGCCCGGCAACCACGCGGTTCGGGAAGTTTTGCAGTAAATCATCCATCGCCt                                                        >  1:2130932/1‑69 (MQ=255)
       aGATCACCACATTCAGCAGCCCGGCAACCACGCGGTTCGGGAAGTTTTGCAGTAAATCATCCATCGCCt                                                        >  1:3308136/1‑69 (MQ=255)
       aGATCACCACATTCAGCAGCCCGGCAACCACGCGGTTCGGGAAGTTTTGCAGTAAATCATCCATCGCCt                                                        >  1:2315279/1‑69 (MQ=255)
       aGATCACCACATTCAGCAGCCCGGCAACCACGCGGTTCGGGAAGTTTTGCAGTAAATCATCCATCGCCt                                                        >  1:1675948/1‑69 (MQ=255)
       aGATCACCACATTCAGCAGCCCGGCAACCACGCGGTTCGGGAAGTTTTGCAGTAAATCATCCATCGCCt                                                        >  1:1489706/1‑69 (MQ=255)
       aGATCACCACATTCAGCAGCCCGGCAACCACGCGGTTCGGGAAGTTTTGCAGTAAATCATCCATCGCCt                                                        >  1:1299067/1‑69 (MQ=255)
         aTCACCACATTCAGCAGCCCGGCAACCACGCGGTTCGGGAAGTTTTGCAGTAAATCATCCATCGCCTg                                                       >  1:982653/1‑68 (MQ=255)
         aTCACCACATTCAGCAGCCCGGCAACCACGCGGTTCGGGAAGTTTTGCAGTAAATCATCCATCGCCTg                                                       >  1:19525/1‑68 (MQ=255)
         aTCACCACATTCAGCAGCCCGGCAACCACGCGGTTCGGGAAGTTTTGCAGTAAATCATCCATCGCCTg                                                       >  1:2223199/1‑68 (MQ=255)
           caccacATTCAGCAGCCCGGCAACCACGCGGTTCGGGAAGTTTTGCAGTAGATCATCCATCGCCTGTTc                                                    >  1:1448361/1‑69 (MQ=255)
           caccacATTCAGCAGCCCGGCAACCACGCGGTTCGGGAAGTTTTGCAGTAAATCATCCATCGCCTGTTc                                                    >  1:882173/1‑69 (MQ=255)
           caccacATTCAGCAGCCCGGCAACCACGCGGTTCGGGAAGTTTTGCAGTAAATCATCCATCGCCTGTTc                                                    >  1:2114466/1‑69 (MQ=255)
           caccacATTCAGCAGCCCGGCAACCACGCGGTTCGGGAAGTTTTGCAGTAAATCATCCATCGCCTGTTc                                                    >  1:1225854/1‑69 (MQ=255)
           caccacATTCAGCAGCCCGGCAACCACGCGGTTCGGGAAGTTTTGCAGTAAATCATCCATCGCCTGTTc                                                    >  1:3403579/1‑69 (MQ=255)
           caccacATTCAGCAGCCCGGCAACCACGCGGTTCGGGAAGTTTTGCAGTAAATCATCCATCGCCTGTTc                                                    >  1:3054701/1‑69 (MQ=255)
           caccacATTCAGCAGCCCGGCAACCACGCGGTTCGGGAAGTTTTGCAGTAAATCATCCATCGCCTGTTc                                                    >  1:884952/1‑69 (MQ=255)
           caccacATTCAGCAGCCCGGCAACCACGCGGTTCGGGAAGTTTTGCAGTAAATCATCCATCGCCTGTTc                                                    >  1:188054/1‑69 (MQ=255)
           caccacATTCAGCAGCCCGGCAACCACGCGGTTCGGGAAGTTTTGCAGTAAATCATCCATCGCCTGTTc                                                    >  1:1679423/1‑69 (MQ=255)
           caccacATTCAGCAGCCCGGCAACCACGCGGTTCGGGAAGTTTTGCAGTAAATCATCCATCGCCTGTTc                                                    >  1:1677222/1‑69 (MQ=255)
                caTTCAGCAGCCCGGCAACCACGCGGTTCGGGAAGTTTTGCAGTAAATCATCCATCGCCTGTTCAGCCt                                               <  1:1463371/69‑1 (MQ=255)
                caTTCAGCAGCCCGGCAACCACGCGGTTCGGGAAGTTTTGCAGTAAATCATCCATCGCCTGTTCAGCCt                                               <  1:1774700/69‑1 (MQ=255)
                caTTCAGCAGCCCGGCAACCACGCGGTTCGGGAAGTTTTGCAGTAAATCATCCATCGCCTGTTCAGCCt                                               <  1:2556118/69‑1 (MQ=255)
                caTTCAGCAGCCCGGCAACCACGCGGTTCGGGAAGTTTTGCAGTAAATCATCCATCGCCTGTTCAGCCt                                               <  1:1484780/69‑1 (MQ=255)
                  ttCAGCAGCCCGGCAACCACGCGGTTCGGGAAGTTTTGCAGTAAATCATCCATCGCCTGTTCAGCCTGa                                             <  1:362587/69‑1 (MQ=255)
                  ttCAGCAGCCCGGCAACCACGCGGTTCGGGAAGTTTTGCAGTAAATCATCCATCGCCTGTTCAGCCTGa                                             <  1:2673162/69‑1 (MQ=255)
                  ttCAGCAGCCCGGCAACCACGCGGTTCGGGAAGTTTTGCAGTAAATCATCCATCGCCTGTTCAGCCTGa                                             <  1:2109125/69‑1 (MQ=255)
                  ttCAGCAGCCCGGCAACCACGCGGTTCGGGAAGTTTTGCAGTAAATCATCCATCGCCTGGTCAGCCTGa                                             <  1:2407435/69‑1 (MQ=255)
                           ccGGCAACCACGCGGTTCGGGAAGTTTTGCAGTAAATCATCCATCGCCTGTTCAGCCTGATACAGCGCa                                    <  1:1557888/69‑1 (MQ=255)
                                aaCCACGCGGTTCGGGAAGTTTTGCAGTAAATCATCCATCGCCTGTTCAGCCTGATACAGCGCATCTTg                               <  1:1979500/69‑1 (MQ=255)
                                aaCCACGCGGTTCGGGAAGTTTTGCAGTAAATCATCCATCGCCTGTTCAGCCTGATACAGCGCATCTTg                               <  1:502896/69‑1 (MQ=255)
                                aaCCACGCGGTTCGGGAAGTTTTGCAGTAAATCATCCATCGCCTGTTCAGCCTGATACAGCGCATCTTg                               <  1:627586/69‑1 (MQ=255)
                                aaCCACGCGGTTCGGGAAGTTTTGCAGTAAATCATCCATCGCCTGTTCAGCCTGATACAGCGCATCTTg                               <  1:985903/69‑1 (MQ=255)
                                 aCCACGCGGTTCGGGAAGTTTTGCAGTAAATCATCCATCGCCTGTTCAGCCTGATACAGCGCATCTTg                               <  1:1891999/68‑1 (MQ=255)
                                     cgcgGTTCGGGAAGTTTTGCAGTAAATCATCCATCGCCTGTTCAGCCTGATACAGCGCATCTTGTACGc                          <  1:1541038/69‑1 (MQ=255)
                                       cgGTTCGGGAAGTTTTGCAGTAAATCGTCCATCGCCTGTTCTGCCTGATa                                           >  1:1891564/1‑50 (MQ=255)
                                       cgGTTCGGGAAGTTTTGCAGTAAATCGTCCATCGCCTGTTCTGCCTGATACAGCGCATCCTGTAcac                          >  1:1132225/1‑65 (MQ=255)
                                       cgGTTCGGGAAGTTTTGCAGTAAATCGTCCATCGCCTGTTCTGCCTGATACAGCGCATCCTGTAcac                          >  1:896303/1‑65 (MQ=255)
                                       cgGTTCGGGAAGTTTTGCAGTAAATCGTCCATCGCCTGTTCTGCCTGATACAGCGCATCCTGTAcac                          >  1:1880553/1‑65 (MQ=255)
                                       cgGTTCGGGAAGTTTTGCAGTAAATCGTCCATCGCCTGTTCTGCCTGATACAGCGCATCCTGTAcac                          >  1:1887662/1‑65 (MQ=255)
                                       cgGTTCGGGAAGTTTTGCAGTAAATCGTCCATCGCCTGTTCTGCCTGATACAGCGCATCCTGTAcac                          >  1:865182/1‑65 (MQ=255)
                                       cgGTTCGGGAAGTTTTGCAGTAAATCGTCCATCGCCTGTTCTGCCTGATACAGCGCATCCTGTAcac                          >  1:1111516/1‑65 (MQ=255)
                                       cgGTTCGGGAAGTTTTGCAGTAAATCGTCCATCGCCTGTTCTGCCTGATACAGCGCATCCTGTAcac                          >  1:907031/1‑65 (MQ=255)
                                       cgGTTCGGGAAGTTTTGCAGTAAATCGTCCATCGCCTGTTCTGCCTGATACAGCGCATCCTGTAcac                          >  1:593447/1‑65 (MQ=255)
                                       cgGTTCGGGAAGTTTTGCAGTAAATCGTCCATCGCCTGTTCTGCCTGATACAGCGCATCCTGTAcac                          >  1:191151/1‑65 (MQ=255)
                                       cgGTTCGGGAAGTTTTGCAGTAAATCGTCCATCGCCTGTTCTGCCTGATACAGCGCATCCTGTAcac                          >  1:3432016/1‑65 (MQ=255)
                                       cgGTTCGGGAAGTTTTGCAGTAAATCGTCCATCGCCTGTTCTGCCTGATACAGCGCATCCTGTAcac                          >  1:2009566/1‑65 (MQ=255)
                                       cgGTTCGGGAAGTTTTGCAGTAAATCGTCCATCGCCTGTTCTGCCTGATACAGCGCATCCTGTAcac                          >  1:2021883/1‑65 (MQ=255)
                                       cgGTTCGGGAAGTTTTGCAGTAAATCGTCCATCGCCTGTTCTGCCTGATACAGCGCATCCTGTAcac                          >  1:1059105/1‑65 (MQ=255)
                                       cgGTTCGGGAAGTTTTGCAGTAAATCGTCCATCGCCTGTTCTGCCTGATACAGCGCATCCTGTAcac                          >  1:3201898/1‑65 (MQ=255)
                                       cgGTTCGGGAAGTTTTGCAGTAAATCGTCCATCGCCTGTTCTGCCTGATACAGCGCATCCTGTAcac                          >  1:3055078/1‑65 (MQ=255)
                                       cgGTTCGGGAAGTTTTGCAGTAAATCGTCCATCGCCTGTTCTGCCTGATACAGCGCATCCTGTAcac                          >  1:2178427/1‑65 (MQ=255)
                                       cgGTTCGGGAAGTTTTGCAGTAAATCGTCCATCGCCTGTTCTGCCTGATACAGCGCATCCTGTAcac                          >  1:1611609/1‑65 (MQ=255)
                                       cgGTTCGGGAAGTTTTGCAGTAAATCGTCCATCGCCTGTTCTGCCTGATACAGCGCATCCTGTAcac                          >  1:2884486/1‑65 (MQ=255)
                                       cgGTTCGGGAAGTTTTGCAGTAAATCGTCCATCGCCTGTTCTGCCTGATACAGCGCATCCTGTAcac                          >  1:2881968/1‑65 (MQ=255)
                                       cgGTTCGGGAAGTTTTGCAGTAAATCGTCCATCGCCTGTTCTGCCTGATACAGCGCATCCTGTAcac                          >  1:2336832/1‑65 (MQ=255)
                                       cgGTTCGGGAAGTTTTGCAGTAAATCGTCCATCGCCTGTTCTGCCTGATACAGCGCATCCTGTAca                           >  1:2830073/1‑65 (MQ=255)
                                       cgGTTCGGGAAGTTTTGCAGTAAATCGTCCATCGCCTGTTCTCCCTGATACAGCGCATCCTGTAcac                          >  1:996610/1‑65 (MQ=255)
                                                             aaTCATCCATCGCCTGTTCAGCCTGATACAGCGCATCTTGTACGCCCCAGTGCACCAGCGGCAGGTCgg  >  1:1890003/1‑69 (MQ=255)
                                                             aaTCATCCATCGCCTGTTCAGCCTGATACAGCGCATCTTGTACGCCCCAGTGCACCAGCGGCAGGTCgg  >  1:1474317/1‑69 (MQ=255)
                                                             aaTCATCCATCGCCTGTTCAGCCTGATACAGCGCATCTTGTACGCCCCAGTGCACCAGCGGCAGGTCgg  >  1:3030353/1‑69 (MQ=255)
                                                             aaTCATCCATCGCCTGTTCAGCCTGATACAGCGCATCTTGTACGCCCCAGTGCACCAGCGGCAGGTCgg  >  1:2199534/1‑69 (MQ=255)
                                                                 |                                                                
GTCGGGAAGATCACCACATTCAGCAGCCCGGCAACCACGCGGTTCGGGAAGTTTTGCAGTAAATCATCCATCGCCTGTTCAGCCTGATACAGCGCATCTTGTACGCCCCAGTGCACCAGCGGCAGGTCGG  >  minE/220092‑220221

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: