Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,336,886 C→T 43.0% L125L (CTG→TTG)  fdoI → formate dehydrogenase‑O, cytochrome b556 subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,336,8860CT43.0% 18.9 / 47.3 65L125L (CTG→TTG) fdoIformate dehydrogenase‑O, cytochrome b556 subunit
Reads supporting (aligned to +/- strand):  ref base C (12/25);  new base T (13/15);  total (25/40)
Fisher's exact test for biased strand distribution p-value = 3.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.79e-01

GACACCGGGCGTTATAACTTCGGTCAGAAATGCGTTTTCTGGGCGGCGATTATTTTCCTGGTTCTGCTGCTGGTGAGCGGCGTGATTATCTGGCGTCCTTATTTTGCGCCTGCTTTCTCAATCCCGGTGA  >  minE/2336820‑2336949
                                                                  |                                                               
gACACCGGGCGTTATAACTTCGGTCAGAAATGCGTTTTCTGGGCGGCGATTATTTTCCTGGTTgtgctg                                                               >  1:1932546/1‑69 (MQ=255)
gACACCGGGCGTTATAACTTCGGTCAGAAATGCGTTTTCTGGGCGGCGATTATTTTCCTGGTTctgctg                                                               >  1:1293303/1‑69 (MQ=255)
gACACCGGGCGTTATAACTTCGGTCAGAAATGCGTTTTCTGGGCGGCGATTATTTTCCTGGTTctgctg                                                               >  1:1907169/1‑69 (MQ=255)
             ataACTTCGGTCAGAAATGCGTTTTCTGGGCGGCGATTATTTTCCTGGTTCTGCTGCTGGTGAgcggc                                                   >  1:804889/1‑68 (MQ=255)
             ataACTTCGGTCAGAAATGCGTTTTCTGGGCGGCGATTATTTTCCTGGTTCTGCTGCTGGTGAgcggc                                                   >  1:1910252/1‑68 (MQ=255)
             ataACTTCGGTCAGAAATGCGTTTTCTGGGCGGCGATTATTTTCCTGGTTCTGCTGCTGGTGAgcggc                                                   >  1:1676408/1‑68 (MQ=255)
             ataACTTCGGTCAGAAATGAGTTTTCTGGGCGGCGATTATTTTCCTGGTTCTGCTGCTGGTGAgagga                                                   >  1:1933644/1‑67 (MQ=255)
                aCTTCGGTCAGAAATGCGTTTTCTGGGCGGCGATTATTTTCCTGGTCCTGTTGCTGGTGAGCGGCGTg                                                >  1:1727216/1‑68 (MQ=255)
                aCTTCGGTCAGAAATGCGTTTTCTGGGCGGCGATTATTTTCCTGGTCCTGTTGCTGGTGAGCGGCGTg                                                >  1:2151632/1‑68 (MQ=255)
                aCTTCGGTCAGAAATGCGTTTTCTGGGCGGCGATTATTTTCCTGGTCCTGTTGCTGGTGAGCGGCGTg                                                >  1:282419/1‑68 (MQ=255)
                aCTTCGGTCAGAAATGCGTTTTCTGGGCGGCGATTATTTTCCTGGTCCTGTTGCTGGTGAGCGGCGTg                                                >  1:2211449/1‑68 (MQ=255)
                aCTTCGGTCAGAAATGCGTTTTCTGGGCGGCGATTATTTTCCTGGTCCTGTTGCTGGTGAGCGGCGTg                                                >  1:3164352/1‑68 (MQ=255)
                aCTTCGGTCAGAAATGCGTTTTCTGGGCGGCGATTATTTTCCTGGTCCTGTTGCTGGTGAGCGGCGTg                                                >  1:3276896/1‑68 (MQ=255)
                aCTTCGGTCAGAAATGCGTTTTCTGGGCGGCGATTATTTTCCTGGTCCTGTTGCTGGTGAGCGGCGTg                                                >  1:3297644/1‑68 (MQ=255)
                aCTTCGGTCAGAAATGCGTTTTCTGGGCGGCGATTATTTTCCTGGTCCTGTTGCTGGTGAGCGGCGTg                                                >  1:2772792/1‑68 (MQ=255)
                aCTTCGGTCAGAAATGCGTTTTCTGGGCGGCGATTATTTTCCTGGTCCTGTTGCTGGTGAGCGGCGTg                                                >  1:512159/1‑68 (MQ=255)
                aCTTCGGTCAGAAATGCGTTTTCTGGGCGGCGATTATTTTCCTGGTCCTGTTGCTGGTGAGCGGCGTg                                                >  1:653483/1‑68 (MQ=255)
                aCTTCGGTCAGAAATGCGTTTTCTGGGCGGCGATTATTTTCCTGGTCCTGTTGCTGGTGAGCGGCGTg                                                >  1:701942/1‑68 (MQ=255)
                aCTTCGGTCAGAAATGCGTTTTCTGGGCGGCGATTATTTTCCTGGTCCTGTTGCTGGTGAGCGGCGTg                                                >  1:117559/1‑68 (MQ=255)
                aCTTCGGTCAGAAATGCGTTTTCTGGGCGGCGATTATTTTCCTGGTCCTGTTGCTGGTGAGCGGCGTg                                                >  1:2646071/1‑68 (MQ=255)
                aCTTCGGTCAGAAATGCGTTTTCTGGGCGGCGATTATTTTCCAGGTCctgttgctg                                                            >  1:1685492/1‑56 (MQ=255)
                  ttCGGTCAGAAATGCGTTTTCTGGGCGGCGATTATTTTCCTGGTTCTGCTGCTGGTGAGCGGCGTGAtt                                             <  1:2360844/69‑1 (MQ=255)
                  ttCGGTCAGAAATGCGTTTTCTGGGCGGCGATTATTTTCCTGGTTCTGCTGCTGGTGAGCGGCGTGAtt                                             <  1:2688084/69‑1 (MQ=255)
                  ttCGGTCAGAAATGCGTTTTCTGGGCGGCGATTATTTTCCTGGTTCTGCTGCTGGTGAGCGGCGTGAtt                                             <  1:2790516/69‑1 (MQ=255)
                  ttCGGTCAGAAATGCGTTTTCTGGGCGGCGATTATTTTCCTGGTTCTGCTGCTGGTGAGCGGCGTGAtt                                             <  1:2827806/69‑1 (MQ=255)
                  ttCGGTCAGAAATGCGTTTTCTGGGCGGCGATTATTTTCCTGGTTCTGCTGCTGGTGAGCGGCGTGAtt                                             <  1:3321765/69‑1 (MQ=255)
                  ttCGGTCAGAAATGCGTTTTCTGGGCGGCGATTATTTTCCTGGTTCTGCTGCTGGTGAGCGGCGTGAtt                                             <  1:538394/69‑1 (MQ=255)
                  ttCGGTCAGAAATGCGTTTTCTGGGCGGCGATTATTTTCCTGGTTCTGCTGCTGGTGAGCGGCGTGAtt                                             <  1:981962/69‑1 (MQ=255)
                             aTGCGTTTTCTGGGCGGCGATTATTTTCCTGGTTCTGCTGCTGGTGAGCGGCGTGATTATCTg                                        >  1:1157191/1‑63 (MQ=255)
                             aTGCGTTTTCTGGGCGGCGATTATTTTCCTGGTTCTGCTGCTGGTGAGCGGCGTGATTATCTGGCGTc                                   >  1:2237141/1‑68 (MQ=255)
                             aTGCGTTTTCTGGGCGGCGATTATTTTCCTGGTTCTGCTGCTGGTGAGCGGCGTGATTATCTGGCGTc                                   >  1:2977858/1‑68 (MQ=255)
                             aTGCGTTTTCTGGGCGGCGATTATTTTCCTGGTTCTGCTGCTGGTGAGCGGCGTGATTATCTGGCGTc                                   >  1:796991/1‑68 (MQ=255)
                             aTGCGTTTTCTGGGCGGCGATTATTTTCCTGGTTCTGCTGCTGGTGAGCGGCGTGATTATCTGGCGTc                                   >  1:2218338/1‑68 (MQ=255)
                             aTGCGTTTTCTGGGCGGCGATTATTTTCCTGGTTCTGCTGCTGGTGAGCGGCGTGATTATCTGGCGTc                                   >  1:3153827/1‑68 (MQ=255)
                             aTGCGTTTTCTGGGCGGCGATTATTTTCCTGGTTCTGCTGCTGGTGAGCGGCGTGATTATCTGGCGTc                                   >  1:3180581/1‑68 (MQ=255)
                             aTGCGTTTTCTGGGCGGCGATTATTTTACTGGTTCTGCTGCTGGTGAGCGGCGTGATTATCTGGCGTc                                   >  1:1346296/1‑68 (MQ=255)
                              tGCGTTTTCTGGGCGGCGATTATTTTCCTGGTTCTGCTGCTGGTGAGCGGCGTGATTATCTGGCGTCCt                                 <  1:225354/69‑1 (MQ=255)
                              tGCGTTTTCTGGGCGGCGATTATTTTCCTGGTTCTGCTGCTGGTGAGCGGCGTGATTATCTGGCGTCCt                                 <  1:392234/69‑1 (MQ=255)
                              tGCGTTTTCTGGGCGGCGATTATTTTCCTGGTTCTGCTGCTGGTGAGCGGCGTGATTATCTGGCGTCCt                                 <  1:1564487/69‑1 (MQ=255)
                              tGCGTTTTCTGGGCGGCGATTATTTTCCTGGTTCTGCTGCTGGTGAGCGGCGTGATTATCTGGCGTCCt                                 <  1:3080611/69‑1 (MQ=255)
                              tGCGTTTTCTGGGCGGCGATTATTTTCCTGGTTCTGCTGCTGGTGAGCGGCGTGATTATCTGGCGTCCt                                 <  1:1555000/69‑1 (MQ=255)
                               gCGTTTTCTGGGCGGCGATTATTTTCCTGGTTctgct                                                                >  1:3161674/1‑37 (MQ=255)
                               gCGTTTTCTGGGCGGCGATTATTTTCCTGGTTctgct                                                                >  1:3431555/1‑37 (MQ=255)
                               gCGTTTTCTGGGCGGCGATTATTTTCCTGGTTctgct                                                                >  1:1766628/1‑37 (MQ=255)
                               gCGTTTTCTGGGCGGCGATTATTTTCCTGGTTctgct                                                                >  1:173776/1‑37 (MQ=255)
                               gCGTTTTCTGGGCGGCGATTATTTTCCTGGTTctgct                                                                >  1:1728667/1‑37 (MQ=255)
                               gCGTTTTCTGGGCGGCGATTATTTTCCTGGTTctgct                                                                >  1:2530390/1‑37 (MQ=255)
                               gCGTTTTCTGGGCGGCGATTATTTTCCTGGTTctgct                                                                >  1:2350957/1‑37 (MQ=255)
                               gCGTTTTCTGGGCGGCGATTATTTTCCTGGTTctgct                                                                >  1:2264573/1‑37 (MQ=255)
                               gCGTTTTCTGGGCGGCGATTATTTTCCTGGTTCTGCTGCTGGTGAGCGGCGTGATTATCTGGCGTCCt                                 <  1:2608864/68‑1 (MQ=255)
                               gCGTTTTCTGGGCGGCGATTATTTTCCTGGTTCTGCTGCTGGGGAGCGGCGTGATTATCTGGCGTCCt                                 <  1:2006062/68‑1 (MQ=255)
                                cGTTTTCTGGGCGGCGATTATTTTCCTGGTTCTGCTGCTGGTGAGCGGCGTGATTATCTGGCGTCCTTa                               <  1:1893319/69‑1 (MQ=255)
                                cGTTTTCTGGGCGGCGATTATTTTCCTGGTTCTGCTGCTGGTGAGCGGCGTGATTATCTGGCGTCCTTa                               <  1:1619979/69‑1 (MQ=255)
                                cGTTTTCTGGGCGGCGATTATTTTCCTGGTTCTGCTGCTGGTGAGCGGCGTGATTATCTGGCGTCCTTa                               <  1:2262967/69‑1 (MQ=255)
                                        gggCGGCGATTATTTTCCTGGTTCTGCTGCTGGTGAGCGGCGTGATTATCTGGCGTCCTTATTTTgcg                        <  1:50634/68‑1 (MQ=255)
                                                  tattTTCCTGGTTCTGCTGCTGGTGAGCGGCGTGATTATCTGGCGTCCTTATTTTGCGCCTGCTTtctc             <  1:80329/69‑1 (MQ=255)
                                                  tattTTCCTGGTTCTGCTGCTGGTGAGCGGCGTGATTATCTGGCGTCCTTATTTTGCGCCTGCTTtctc             <  1:811948/69‑1 (MQ=255)
                                                  tattTTCCTGGTTCTGCTGCTGGTGAGCGGCGTGATTATCTGGCGTCCTTATTTTGCGCCTGCTTtctc             <  1:95222/69‑1 (MQ=255)
                                                  tattTTCCTGGTTCTGCTGCTGGTGAGCGGCGTGATTATCTGGCGTCCTTATTTTGCGCCTGCTTtctc             <  1:229472/69‑1 (MQ=255)
                                                  tattTTCCTGGTTCTGCTGCTGGTGAGCGGCGTGATTATCTGGCGTCCTTATTTTGCGCCTGCTTtctc             <  1:1894476/69‑1 (MQ=255)
                                                  tattTTCCTGGTTCTGCTGCTGGTGAGCGGCGTGATTATCTGGCGTCCTTATTTTGCGCCTGCTTtctc             <  1:1961974/69‑1 (MQ=255)
                                                  tattTTCCTGGTCCTGTTGCTGGTGAGCGGCGTGATTATCTGGCGTCCTTATTTTGCGCCTGCTTtctc             <  1:388940/69‑1 (MQ=255)
                                                  tattTTCCTGGTCCTGTTGCTGGTGAGCGGCGTGATTATCTGGCGTCCTTATTTTGCGCCTGCTTtctc             <  1:3174355/69‑1 (MQ=255)
                                                  tattTTCCTGGTCCTGTTGCTGGTGAGCGGCGTGATTATCTGGCGTCCTTATTTTGCGCCTGCTTtctc             <  1:3173907/69‑1 (MQ=255)
                                                  tattTTCCTGGTCCTGTTGCTGGTGAGCGGCGTGATTATCTGGCGTCCTTATTTTGCGCCTGCTTtctc             <  1:3116440/69‑1 (MQ=255)
                                                  tattTTCCTGGTCCTGTTGCTGGTGAGCGGCGTGATTATCTGGCGTCCTTATTTTGCGCCTGCTTtctc             <  1:615183/69‑1 (MQ=255)
                                                  tattTTCCTGGTCCTGTTGCTGGTGAGCGGCGTGATTATCTGGCGTCCTTATTTTGCGCCTGCTTtctc             <  1:657487/69‑1 (MQ=255)
                                                  tattTTCCTGGTCCTGTTGCTGGTGAGCGGCGTGATTATCTGGCGTCCTTATTTTGCGCCTGCTTtctc             <  1:710984/69‑1 (MQ=255)
                                                  tattTTCCTGGTCCTGTTGCTGGTGAGCGGCGTGATTATCTGGCGTCCTTATTTTGCGCCTGCTTtctc             <  1:2738229/69‑1 (MQ=255)
                                                  tattTTCCTGGTCCTGTTGCTGGTGAGCGGCGTGATTATCTGGCGTCCTTATTTTGCGCCTGCTTtctc             <  1:2724387/69‑1 (MQ=255)
                                                  tattTTCCTGGTCCTGTTGCTGGTGAGCGGCGTGATTATCTGGCGTCCTTATTTTGCGCCTGCTTtctc             <  1:264532/69‑1 (MQ=255)
                                                  tattTTCCTGGTCCTGTTGCTGGTGAGCGGCGTGATTATCTGGCGTCCTTATTTTGCGCCTGCTTtctc             <  1:851902/69‑1 (MQ=255)
                                                  tattTTCCTGGTCCTGTTGCTGGTGAGCGGCGTGATTATCTGGCGTCCTTATTTTGCGCCTGCTTtctc             <  1:904348/69‑1 (MQ=255)
                                                  tattTTCCTGGTCCTGTTGCTGGTGAGAGGCGTGATTATCTGGCGTCCTTATTTTGCGCCTGCTTtctc             <  1:1491424/69‑1 (MQ=37)
                                                  tattTTCCTGCTCCTGTTGCTGGTGAGCGGCGTGATTATCTGGCGTCCTTATTTTGCGCCTGCTTtctc             <  1:3192857/69‑1 (MQ=25)
                                                   attTTCCTGGTCCTGTTGCTGGTGAGCGGCGTGATTATCTGGCGTCCTTATTTTGCGCCTGCTTtctc             <  1:2715418/68‑1 (MQ=255)
                                                              tCTGCTGCTGGTGAGCGGCGTGATTATCTGGCGTCCTTATTTTGCGCCTGCTTTCTCAATCCCGGTGa  <  1:1125768/68‑1 (MQ=255)
                                                                  |                                                               
GACACCGGGCGTTATAACTTCGGTCAGAAATGCGTTTTCTGGGCGGCGATTATTTTCCTGGTTCTGCTGCTGGTGAGCGGCGTGATTATCTGGCGTCCTTATTTTGCGCCTGCTTTCTCAATCCCGGTGA  >  minE/2336820‑2336949

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: