Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,111,527 G→A 21.6% C287Y (TGT→TAT)  ydiD → short chain acyl‑CoA synthetase, anaerobic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,111,5270GA21.6% 107.9 / 25.1 51C287Y (TGT→TAT) ydiDshort chain acyl‑CoA synthetase, anaerobic
Reads supporting (aligned to +/- strand):  ref base G (19/21);  new base A (5/6);  total (24/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

TAACGGCACCATTCTTAATTGGCGCTCGCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTGTCTCGCGCTGCTTGAGCAGCAGCGTTGCACCTGTATGCTCGGCGCAACGCCGTTTGTCTATGATCT  >  minE/1111464‑1111593
                                                               |                                                                  
tAACGGCACCATTCTTAATTGGCGCTCGCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTGtctc                                                                <  1:1027878/68‑1 (MQ=255)
tAACGGCACCATTCTTAATTGGCGCTCGCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTGtctc                                                                <  1:82632/68‑1 (MQ=255)
tAACGGCACCATTCTTAATTGGCGCTCGCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTGtctc                                                                <  1:146123/68‑1 (MQ=255)
tAACGGCACCATTCTTAATTGGCGCTCGCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTGtctc                                                                <  1:363530/68‑1 (MQ=255)
tAACGGCACCATTCTTAATTGGCGCTCGCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTGtctc                                                                <  1:2868158/68‑1 (MQ=255)
tAACGGCACCATTCTTAATTGGCGCTCGCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTGtctc                                                                <  1:2026555/68‑1 (MQ=255)
        ccATTCTTAATTGGCGCTCGCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTGTCTcgcg                                                             <  1:2463168/63‑1 (MQ=255)
        ccATTCTTAATTGGCGCTCGCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTATCTCGCGCTGCtt                                                       <  1:1536126/69‑1 (MQ=255)
        ccATTCTTAATTGGCGCTCGCAGCGTGTTGTTAGATATTTTCACTCCAGATGCGTGTCTcgcg                                                             <  1:1899772/63‑1 (MQ=255)
          aTTCTTAATTGGCGCTCGCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTATCTCGCGCTGCTTGa                                                     >  1:3604630/1‑69 (MQ=255)
            tCTTAATTGGCGCTCGCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTGTCTCGCGCTGCTTGagcag                                                 >  1:710766/1‑71 (MQ=255)
            tCTTAATTGGCGCTCGCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTGTCTCGCGCTGCTTGagc                                                   >  1:97594/1‑69 (MQ=255)
            tCTTAATTGGCGCTCGCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTGTCTCGCGCTGCTTGagc                                                   >  1:477396/1‑69 (MQ=255)
            tCTTAATTGGCGCTCGCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTGTCTCGCGCTGCTTGagc                                                   >  1:2942069/1‑69 (MQ=255)
            tCTTAATTGGCGCTCGCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTGTCTCGCGCTGCTTGagc                                                   >  1:2268834/1‑69 (MQ=255)
             cTTAATTGGCGCTCGCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTGTCTCGCGCTGCTTGagcag                                                 >  1:3652417/1‑70 (MQ=255)
             cTTAATTGGCGCTCGCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTATCTCGCGCTGCTTGagcag                                                 <  1:3525118/70‑1 (MQ=255)
             cTTAATTGGCGCTCGCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTATCTCGCGCTGCTTGagcag                                                 <  1:1261346/70‑1 (MQ=255)
             cTTAATTGGCGCTCGCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTATCTCGCGCTGCTTGagcag                                                 <  1:1519331/70‑1 (MQ=255)
             cTTAATTGGCGCCCGCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTATCTCGCGCTGCTTGagcag                                                 <  1:3704128/70‑1 (MQ=255)
                aaTTGGCGCTCGCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTGTCTCGCGCTGCTTGagcagcag                                              <  1:2377803/70‑1 (MQ=255)
                aaTTGGCGCTCGCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTGTCTCGCGCTGCTTGagcagcag                                              <  1:2544078/70‑1 (MQ=255)
                aaTTGGCGCTCGCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTGTCTCGCGCTGCTTGagcagc                                                <  1:3273906/68‑1 (MQ=255)
                           gCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTGTCTCGCGCTGCTTGAGCAGCAGCGTTGCACCTg                                   <  1:3181339/70‑1 (MQ=255)
                           gCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTGTCTCGCGCTGCTTGAGCAGCAGCGTTGCACCTg                                   <  1:620573/70‑1 (MQ=255)
                           gCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTGTCTCGCGCTGCTTAGCAGCAGCGTTGCACCTg                                   <  1:1080983/69‑1 (MQ=255)
                                    tgttAGATATTTTCACTCCTGATGCGTGTCTCGCGCTGCTTGAGCAGCAGCGTTGCACCTGTATGCTCGgc                         >  1:2193438/1‑71 (MQ=255)
                                     gttAGATATTTTCACTCCTGATGCGTGTCTCGCGCTGCTTGAGCAGCAGCGTTGCACCTGTATGCTCGgcg                        >  1:792876/1‑71 (MQ=255)
                                            atTTTCACTCCTGATGCGTGTCTCGCGCTGCTTGagcag                                                 >  1:4057097/1‑39 (MQ=255)
                                            atTTTCACTCCTGATGCGTGTCTCGCGCTGCTTGagcag                                                 >  1:805828/1‑39 (MQ=255)
                                            atTTTCACTCCTGATGCGTGTCTCGCGCTGCTTGagcag                                                 >  1:3450550/1‑39 (MQ=255)
                                            atTTTCACTCCTGATGCGTGTCTCGCGCTGCTTGagcag                                                 >  1:2180993/1‑39 (MQ=255)
                                            atTTTCACTCCTGATGCGTGTCTCGCGCTGCTTGagcag                                                 >  1:1999345/1‑39 (MQ=255)
                                            atTTTCACTCCTGATGCGTGTCTCGCGCTGCTTGAGCAGCAGCGTTGCACCTGTATGCTCGGCGCAACGcc                 >  1:3320408/1‑71 (MQ=255)
                                            atTTTCACTCCTGATGCGTGTCTCGCGCTGCTTGAGCAGCAGCGTTGCACCTGTATGCTCGGCGCAACGcc                 >  1:1557539/1‑71 (MQ=255)
                                              tttCACTCCTGATGCGTGTCTCGCGCTGCTTGAGCAGCAGCGTTGCACCTGTATGCTCGGCGCAACGCCGt               >  1:2253229/1‑71 (MQ=255)
                                              tttCACTCCTGATGCGTATCTCGCGCTGCTTGAGCAGCAGCGTTGCACCTGTATGCTCGGCGCAACGCCGt               >  1:324690/1‑71 (MQ=255)
                                              tttCACTCCTGATGCGTATCTCGCGCTGCTTGAGCAGCAGCGTTGCACCTGTATGCTCGGCGCAACGCCGt               >  1:2230031/1‑71 (MQ=255)
                                              tttCACTCCTGATGCGTATCTCGCGCTGCTTGAGCAGCAGCGTTGCACCTGTATGCTCGGCGCAACGCCGt               >  1:3090007/1‑71 (MQ=255)
                                              tttCACTCCTGATGCGTATCTCGCGCTGCTTGAGCAGCAGCGTTGCACCTGTATGCTCGGCGCAACGCCGt               >  1:4109996/1‑71 (MQ=255)
                                               ttCACTCCTGATGCGTGTCTCGCGCTGCTTGAGCAGCAGCGTTGCACCTGTATGCTCGGCGCAACGCCGtt              >  1:20340/1‑71 (MQ=255)
                                               ttCACTCCTGATGCGTATCTCGCGCTGCTTGAGCAGCAGCGTTGCACCTGTATGCTCGGCGCAACGCCGt               <  1:2346892/70‑1 (MQ=255)
                                                     ccTGATGCGTGTCTCGCGCTGCTTGAGCAGCAGCGTTGCAc                                      >  1:3291423/1‑41 (MQ=255)
                                                     ccTGATGCGTGTCTCGCGCTGCTTGAGCAGCAGCGTTGCAc                                      >  1:665974/1‑41 (MQ=255)
                                                         aTGCGTGTCTCGCGCTGCTTGAGCAGCAGCGTTGCACCTGTATGCTCGGCGCAACGGCGTTTGTCTATg      <  1:719771/69‑1 (MQ=255)
                                                         aTGCGTGTCTCGCGCTGCTTGAGCAGCAGCGTTGCACCTGTATGCTCGGCGCAACGCCGTTTGTCTa        <  1:2184815/67‑1 (MQ=255)
                                                         aTGCGTGTCTCGCGCTGCTTGAGCAGCAGCGTTGCACCTGTATGCTCGGCGCAACGCCGTTTGTCTATg      <  1:1790030/69‑1 (MQ=255)
                                                         aTGCGTGTCTCGCGCTGCTTGAGCAGCAGCGTTGCACCTGTATGCTCGGCGCAACGCCGTTTGTCTATg      <  1:1752845/69‑1 (MQ=255)
                                                         aTGCGTGTCTCGCGCTGCTTGAGCAGCAGCGTTGCACCTGTATGCTCGGCGCAACGCCGTTTGTCTATg      <  1:1625018/69‑1 (MQ=255)
                                                            cGTGTCTCGCGCTGCTTGAGCAGCAGCGTTGCACCTGTATGCTCGGCGCAACGCCGTTTGTCTATGATCt  <  1:553392/70‑1 (MQ=255)
                                                            cGTGTCTCGCGCTGCTTGAGCAGCAGCGTTGCACCTGTATGCTCGGCGCAACGCCGTTTGTCTATGATCt  <  1:1446339/70‑1 (MQ=255)
                                                               |                                                                  
TAACGGCACCATTCTTAATTGGCGCTCGCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTGTCTCGCGCTGCTTGAGCAGCAGCGTTGCACCTGTATGCTCGGCGCAACGCCGTTTGTCTATGATCT  >  minE/1111464‑1111593

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: