Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 219,996 C→T 32.8% G711D (GGT→GAT)  fadE ← acyl coenzyme A dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE219,9960CT32.8% 71.5 / 49.8 61G711D (GGT→GAT) fadEacyl coenzyme A dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (19/22);  new base T (10/10);  total (29/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.74e-01

GGCAATCACATCCACCAGCGCCTCTTCCAGCAAGCCAACCGGATTATGCTCGCTCGGCGTCAGGTACTGACCGCGACCAATGCGGGAACGGGTGGCGTTCGGCACTTGTAAAATCTTCGCCACTTTATGATCCAGC  >  minE/219926‑220061
                                                                      |                                                                 
ggcAATCACATCCACCAGCGCCTCTTCCAGCAAGCCAACCGGATTATGCTCGCTCGGCGTCAGGTACTGAc                                                                   >  1:1995698/1‑71 (MQ=255)
   aaTCACATCCACCAGCGCCTCTTCCAGCAAGCCAACCGGATTATGCTCGCTCGGCGTCAGGTACTGACcgc                                                                <  1:644941/71‑1 (MQ=255)
   aaTCACATCCACCAGCGCCTCTTCCAGCAAGCCAACCGGATTATGCTCGCTCGGCGTCAGGTACTGACcg                                                                 >  1:552899/1‑70 (MQ=255)
   aaTCACATCCACCAGCGCCTCTTCCAGCAAGCCAACCGGATTATGCGCGCGCGGCGTCAGGTACTGACcgc                                                                <  1:2722277/71‑1 (MQ=255)
    aTCACATCCACCAGCGCCTCTTCCAGCAAGCCAACCGGATTATGCTCGCTCGGCGTCAGGTACTGACcgc                                                                >  1:2547847/1‑70 (MQ=255)
    aTCACATCCACCAGCGCCTCTTCCAGCAAGCCAACCGGATTATGCTCGCTCGGCGTCAGGTACTGACcgc                                                                >  1:3456952/1‑70 (MQ=255)
    aTCACATCCACCAGCGCCTCTTCCAGAAAGCCAACCGGATTATGCTCGCTCGGCGTCAGGTACTGACcgc                                                                >  1:1097003/1‑70 (MQ=255)
     tCACATCCACCAGCGCCTCTTCCAGCAAGCCAACCGGATTATGCTCGCTCGGCGTCAGGTACTGACCGCGa                                                              >  1:722938/1‑71 (MQ=255)
       acaTCCACCAGCGCCTCTTCCAGCAAGCCAACCGGATTATGCTCGCTCGGCGTCAGGTACTGACCGCGAc                                                             >  1:2979205/1‑70 (MQ=255)
           ccaccaGCGCCTCTTCCAGCAAGCCAACCGGATTATGCTCGCTCGGCGTCAGGTACTGACCGCGACCAATg                                                        >  1:717154/1‑71 (MQ=255)
           ccaccaGCGCCTCTTCCAGCAAGCCAACCGGATTATGCTCGCTCGGCGTCAGGTACTGACCGCGACCAATg                                                        >  1:1609063/1‑71 (MQ=255)
              ccaGCGCCTCTTCCAGCAAGCCAACCGGATTATGCTCGCTCGGCGTCAGGTACTGATCGCGACCAAt                                                         <  1:39509/67‑1 (MQ=255)
              ccaGCGCCTCTTCCAGCAAGCCAACCGGATTATGCTCGCTCGGCGTCAGGTACTGATCGCGACCAATGCgg                                                     <  1:1160176/71‑1 (MQ=255)
                  cgcCTCTTCCAGCAAGCCAACCGGATTATGCTCGCTCGGCGTCAGGTACTGACCGCGACCAATGCGGGAAc                                                 <  1:226224/71‑1 (MQ=255)
                   gcCTCTTCCAGCAAGCCAACCGGATTATGCTCGCTCGGCGTCAGGTACTGACCGCGACCAATGCGGGAACg                                                <  1:2502034/71‑1 (MQ=255)
                     ctctTCCAGCAAGCCAACCGGATTATGCTCGCTCGGCGTCAGGTACTGACCGCGACCAATGCGGGAACggg                                              <  1:1814882/71‑1 (MQ=255)
                     ctctTCCAGCAAGCCAACCGGATTATGCTCGCTCGGCGTCAGGTACTGACCGCGACCAATGCGGGAACgg                                               <  1:144061/70‑1 (MQ=255)
                      tctTCCAGCAAGCCAACCGGATTATGCTCGCTCGGCGTCAGGTACTGACCGCGACCAATGCGGGAACGGGt                                             >  1:2511753/1‑71 (MQ=255)
                      tctTCCAGCAAGCCAACCGGATTATGCTCGCTCGGCGTCAGGTACTGACCGCGACCAATGCGGGAACGGGt                                             >  1:2137723/1‑71 (MQ=255)
                      tctTCCAGCAAGCCAACCGGATTATGCTCGCTCGGCGTCAGGTACTGACCGCGACCAATGCGGGAACGGGt                                             >  1:2732576/1‑71 (MQ=255)
                      tctTCCAGCAAGCCAACCGGATTATGCTCGCTCGGCGTCAGGTACTGACCGCGACCAATGCGGGAACGGGt                                             >  1:1357341/1‑71 (MQ=255)
                       ctTCCAGCAAGCCAACCGGATTATGCTCGCTCGGCGTCAGGTACTGACCGCGACCAATGCGGGAACGGGTg                                            >  1:3868489/1‑71 (MQ=255)
                        ttCCAGCAAGCCAACCGGATTATGCTCGCTCGGCGTCAGGTACTGATCGCGACCAATGCGGGa                                                   >  1:1266367/1‑63 (MQ=255)
                         tCCAGCAAGCCAACCGGATTATGCTCGCTCGGCGTCAGGTACTGAcc                                                                  >  1:983170/1‑47 (MQ=255)
                          ccAGCAAGCCAACCGGATTATGCTCGCTCGGCGTCAGGTACTGATCGCGACCa                                                           >  1:2178921/1‑53 (MQ=255)
                          ccAGCAAGCCAACCGGATTATGCTCGCTCGGCGTCAGGTACTGATCGCGACCa                                                           >  1:3862446/1‑53 (MQ=255)
                          ccAGCAAGCCAACCGGATTATGCTCGCTCGGCGTCAGGTACTGATCGCGACCAATGCGGGAACGGGTGGCg                                         >  1:457572/1‑71 (MQ=255)
                           cAGCAAGCCAACCGGATTATGCTCGCTCGGCGTCAGGTACTGATCGCGACCAATGCGGGAACGGGTGGCg                                         <  1:204308/70‑1 (MQ=255)
                             gCAAGCCAACCGGATTATGCTCGCTCGGCGTCAGGTACTGATCGCGACCAATGc                                                       <  1:1124176/54‑1 (MQ=255)
                             gCAAGCCAACCGGATTATGCTCGCTCGGCGTCAGGTACTGATCGCGACCAATGCGGGAACGGGTg                                            <  1:3073086/65‑1 (MQ=255)
                              cAAGCCAACCGGATTATGCTCGCTCGGCGTCAGGTACTGACcgc                                                                <  1:919426/44‑1 (MQ=255)
                              cAAGCCAACCGGATTATGCTCGCTCGGCGTCAGGTACTGACCGCGACCAATGCGGGAACGGGTGGCGTTCg                                     >  1:3319098/1‑71 (MQ=255)
                               aaGCCAACCGGATTATGCTCGCTCGGCGTCAGGTACTGATCGCGACCAATGCGGGAACGGGTGGCGTTCgg                                    >  1:2152902/1‑71 (MQ=255)
                               aaGCCAACCGGATTATGCTCGCTCGGCGTCAGGTACTGACCGCGACCAATGCGGGAACGGGTGGCGTTCgg                                    <  1:392242/71‑1 (MQ=255)
                                aGCCAACCGGATTATGCTCGCTCGGCGTCAGGTACTGACCGCGACCAATGCggg                                                    <  1:912427/54‑1 (MQ=255)
                                aGCCAACCGGATTATGCTCGCTCGGCGTCAGGTACTGACCGCGACCAATGCggg                                                    <  1:527136/54‑1 (MQ=255)
                                aGCCAACCGGATTATGCTCGCTCGGCGTCAGGTACTGACCGCGACCAATGCggg                                                    <  1:3381616/54‑1 (MQ=255)
                                    aaCCGGATTATGCTCGCTCGGCGTCAGGTACTGACCGCGACCAATGCGGGAACGGGTGGCGTTCGGCACtt                               >  1:2864650/1‑71 (MQ=255)
                                    aaCCGGATTATGCTCGCTCGGCGTCAGGTACTGACCGCGACCAATGCGGGAACGGGTGGCGTTCGGCACtt                               <  1:3677117/71‑1 (MQ=255)
                                    aaCCGGATTATGCTCGCTCGGCGTCAGGTACTGACCGCGACCAATGCGGGAACGGGTGGCGTTCGGCACt                                >  1:2323210/1‑70 (MQ=255)
                                     aCCGGATTATGCTCGCTCGGCGTCAGGTACTGATCGCGACCAATGCGGGAACGGGTGGCGTTCGGCACTTg                              >  1:2667353/1‑71 (MQ=255)
                                       cGGATTATGCTCGCTCGGCGTCAGGTACTGATCGCGACCAATGAGGGAACGGGTGGCGTTCGGCACTTGTa                            <  1:778230/71‑1 (MQ=255)
                                       cGGATTATGCTCGCTCGGCGTCAGGTACTGACCGCGAcc                                                            <  1:2020041/39‑1 (MQ=255)
                                       cGGATTATGCTCGCTCGGCGTCAGGTACTGACCGCGACCAATGCGGGAACGGGTg                                            >  1:2492891/1‑55 (MQ=255)
                                        ggATTATGCTCGCTCGGCGTCAGGTACTGACCGCGACCAATGCGGGAACGGGt                                             <  1:310770/53‑1 (MQ=255)
                                              tGCTCGCTCGGCGTCAGGTACTGACCGCGACCAATGCGGGAACGGGTGGCGTTCGGCACTTGTAAAATCt                      <  1:1843732/70‑1 (MQ=255)
                                               gctcgctcGGCGTCAGGTACTGATCGCGACCAATGCGGGAACGGGTGGCGTTCGGCACTTGTAAAATCTTc                    >  1:2059319/1‑71 (MQ=255)
                                                 tcgctcgGCGTCAGGTACTGACCGCGACCAATGCGGGAACGGGTGGCGTTCGGCACTTGTAAAATCTTc                    <  1:2216023/69‑1 (MQ=255)
                                                    ctcgGCGTCAGGTACTGATCGCGACCAATGCGGGAACGGGTGGCGTTCGGCACTTGTaaa                          <  1:2010088/60‑1 (MQ=255)
                                                    ctcgGCGTCAGGTACTGATCGCGACCAATGCGGGAACGGGTGGCGTTCGGCACTTGTAAAATCTTCGCCAc               >  1:3731307/1‑71 (MQ=255)
                                                       ggCGTCAGGTACTGACCGCTACCAATGCGGGAACGGGTGGCGTTCGGCACTTGTAAAATCTTCGCCACt              <  1:643879/69‑1 (MQ=255)
                                                       ggCGTCAGGTACTGACCGCGACCAATGCGGGAACGGGTGGCGTTCGGCACTTGTAAAATCTTCGCCACt              <  1:987381/69‑1 (MQ=255)
                                                        gagTCAGGTACTGACCGCGACCAATGCGGGAACGGGTGGCGTTCGGCACTTGTAAAATCTTCGCCACttt            <  1:3267438/68‑1 (MQ=255)
                                                           tCAGGTACTGATCGCGACCAATGCGGGAATGGGTGGCGTTCGGCACTTGTAAAATCTTCGCCACTTTATGa        >  1:669994/1‑71 (MQ=255)
                                                           tCAGGTACTGATCGCGACCAATGCGGGAACGGGTGGCGTTCGGCACTTGTAAAATCTTCGCCACTTTATGa        <  1:687767/71‑1 (MQ=255)
                                                             aGGTACTGATCGCGACCAATGCGGGAACGGGTGGCGTTCGGCACTTGTAAAATCTTCGCCACTTTATGATc      <  1:3938695/71‑1 (MQ=255)
                                                             aGGTACTGACCGCGACCAATGCGGGAACGGGTGGCGTTCGGCACTTGTAAAATCTTCGCCACTTTATGAt       <  1:2960908/70‑1 (MQ=255)
                                                              ggTACTGACCGCGACCAATGCGGGAACGGGTGGCGTTCGGCACTTGTAAAATCTTCGCCACTTTATGATcc     <  1:1787010/71‑1 (MQ=255)
                                                               gTACTGACCGCGACCAATGCGGGAACGGGTGGCGTTCGGCACTTGTAAAATCTTCGCCACTTTATGATc      >  1:2708969/1‑69 (MQ=255)
                                                                 aCTGATAGCGACCAATGCGGGAACGGGTGGCGTTCGGCACTTGTAAAATCTTCGCCACtt             >  1:510199/1‑60 (MQ=255)
                                                                 aCTGACCGCGACCAATGCGGGAACGGGTGGCGTTCGGCACTTGTAAAATCTTCGCCACTTTATg         >  1:1746318/1‑64 (MQ=255)
                                                                  cTGATCGCGACCAATGCGGGAACGGGTGGCGTTCGGCACTTGTAAAATCTTCGCCACTTTATGATCCAGc  <  1:308402/70‑1 (MQ=255)
                                                                   tGACCGCGACCAATGCGGGAACGGGTGGCGTTCGGCACTTGTAAAATCt                      <  1:3893192/49‑1 (MQ=255)
                                                                      ccGCGACCAATGCGGGAACGGGTGGCGTTCGGCACTTGTAAAAt                        >  1:3002189/1‑44 (MQ=255)
                                                                      |                                                                 
GGCAATCACATCCACCAGCGCCTCTTCCAGCAAGCCAACCGGATTATGCTCGCTCGGCGTCAGGTACTGACCGCGACCAATGCGGGAACGGGTGGCGTTCGGCACTTGTAAAATCTTCGCCACTTTATGATCCAGC  >  minE/219926‑220061

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: