Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,111,527 G→A 16.7% 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,5270GA16.7% 72.7 / 7.7 30C287Y (TGT→TAT) ydiDshort chain acyl‑CoA synthetase, anaerobic
Reads supporting (aligned to +/- strand):  ref base G (15/10);  new base A (3/2);  total (18/12)
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.48e-01

ACGGCACCATTCTTAATTGGCGCTCGCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTGTCTCGCGCTGCTTGAGCAGCAGCGTTGCACCTGTATGCTCGGCGCAACGCCGTTTGTCTATGATCT  >  minE/1111466‑1111593
                                                             |                                                                  
aCGGCACCATTCTTAATTGGCGCTCGCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTGTCTCGCgct                                                           >  1:1145079/1‑71 (MQ=255)
aCGGCACCATTCTTAATTGGCGCTCGCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTGTCTCGCgct                                                           >  1:471975/1‑71 (MQ=255)
aCGGCACCATTCTTAATTGGCGCTCGCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTGTCTCGCgct                                                           >  1:149413/1‑71 (MQ=255)
      ccATTCTTAATTGGCGCTCGCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTGTCTCGCGCTGCTTGa                                                     <  1:331381/71‑1 (MQ=255)
       cATTCTTAATTGGCGCTCGCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTGTCTCGCGCTGCTTGag                                                    >  1:1786288/1‑71 (MQ=255)
        aTTCTTAATTGGCGCTCGCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTGTCTCGCGCTGCTTGag                                                    >  1:1300060/1‑70 (MQ=255)
        aTTCTTAATTGGCGCTCGCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTGTCTCGCGCTGCTTGag                                                    >  1:1674798/1‑70 (MQ=255)
        aTTCTTAATTGGCGCTCGCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTATCTCGCGCTGCTTg                                                      <  1:503625/68‑1 (MQ=255)
           cTTAATTGGCGCTCGCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTGTCTCGCGCTGCTTGagcag                                                 >  1:167283/1‑70 (MQ=255)
           cTTAATTGGCGCTCGCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTGTCTCGCGCTGCTTGagca                                                  >  1:1438481/1‑69 (MQ=255)
                         gCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTGTCTCGCGCTGCTTGAGCAGCAGCGTTGCACCTg                                   >  1:2024784/1‑70 (MQ=255)
                           aGCGTGTTGTTAGATATTTTCACTCCTGATGCGTGTCTCGCGCTGCTTGAGCAGCAGCGTTGCACCTGTAt                                >  1:1650413/1‑71 (MQ=255)
                              gtgtTGTTAGATATTTTCACTCCTGATGCGTGTCTCGCGCTGCTTGAGCAGCAGCGTTGCACCTGTATGCt                             <  1:1046837/71‑1 (MQ=255)
                                    ttAGATATTTTCACTCCTGATGCGTGTCTCGCGCTGCtt                                                       <  1:1473771/39‑1 (MQ=255)
                                    ttAGATATTTTCACTCCTGATGCGTGTCTCGCGCTGCtt                                                       <  1:2129330/39‑1 (MQ=255)
                                          atTTTCACTCCTGATGCGTATCTCGCGCTGCTTGAGCAGCAGCGTTGCACCTGTATGCTCGGCGCAACGcc                 >  1:598313/1‑71 (MQ=255)
                                             ttCACTCCTGATGCGTGTCTCGCGCTGCTTGAGCAGCAGCGTTGCACCTGTATGCTCGGCGCAACGCCGtt              >  1:845917/1‑71 (MQ=255)
                                             ttCACTCCTGATGCGTGTCTCGCGCTGCTTGAGCAGCAGCGTTGCACCTGTATGCTCGGCGCAACGCCGtt              >  1:1593711/1‑71 (MQ=255)
                                             ttCACTCCTGATGCGTGTCTCGCGCTGCTTGAGCAGCAGCGTTGCACCTGTATGCTCGGCGCAACGCCGtt              >  1:109203/1‑71 (MQ=255)
                                              tCACTCCTGATGCGTATCTCGCGCTGCTTGAGCAGCAGCGTTGCACCTGTATGCTCGGCGCAACGCCGtt              >  1:2060092/1‑70 (MQ=255)
                                              tCACTCCTGATGCGTATCTCGCGCTGCTTGAGCAGCAGCGTTGCACCTGTATGCTCGGCGCAACGCCGtt              >  1:550381/1‑70 (MQ=255)
                                                aCTCCTGATGCGTGTCTCGCGCTGCTTGAGCAGCAGCGTTGCACCTGTATGCTCGTCGCAACGCCGTTTGt           <  1:570437/71‑1 (MQ=255)
                                                    cTGATGCGTGTCTCGCGCTGCTTGAGCAGCAGCGTTGCACCTGt                                  >  1:2174434/1‑44 (MQ=255)
                                                       aTGCGTGTCTCGCGCTGCTTGAGCAGCAGCGTTGCACCTGTATGCTCGGCGCAACGCCGTTTGTCTATGa     >  1:1400115/1‑70 (MQ=255)
                                                       aTGCGTATCTCGCGCTGCTTGAGCAGCAGCGTTGCACCTGTATGCTCGGCGCAACGCCGTTTGTCTATg      <  1:1863260/69‑1 (MQ=255)
                                                        tGCGTGTCTCGCGCTGCTTGAGCAGCAGCGTTGCACCTGTATGCTCGGCGCAACGCCGTTTGTCTATGATc   <  1:1914317/71‑1 (MQ=255)
                                                        tGCGTGTCTCGCGCTGCTTGAGCAGCAGCGTTGCACCTGTATGCTCGGCGCAACGCCGTTTGTCTATGATc   <  1:2026221/71‑1 (MQ=255)
                                                        tGCGTGTCTCGCGCTGCTTGAGCAGCAGCGTTGCACCTGTATGCTCGGCGCAACGCCGTTTGTCTATGATc   <  1:1256416/71‑1 (MQ=255)
                                                        tGCGTGTCTCGCGCTGCTTGAGCAGCAGCGTTGCACCTGTATGCTCGGCGCAACGCAGTTTGTCTATAATc   <  1:10814/71‑1 (MQ=255)
                                                          cGTGTCTCGCGCTGCTTGAGCAGCAGCGTTGCACCTGTATGCTCGGCGCAACGCCGTTTGTCTATGATCt  <  1:1730989/70‑1 (MQ=255)
                                                             |                                                                  
ACGGCACCATTCTTAATTGGCGCTCGCAGCGTGTTGTTAGATATTTTCACTCCTGATGCGTGTCTCGCGCTGCTTGAGCAGCAGCGTTGCACCTGTATGCTCGGCGCAACGCCGTTTGTCTATGATCT  >  minE/1111466‑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: