Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,336,882 T→C 37.9% V123V (GTT→GTC 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,8820TC37.9% 39.6 / 50.6 79V123V (GTT→GTCfdoIformate dehydrogenase‑O, cytochrome b556 subunit
Reads supporting (aligned to +/- strand):  ref base T (25/24);  new base C (15/15);  total (40/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTGACACCGGGCGTTATAACTTCGGTCAGAAATGCGTTTTCTGGGCGGCGATTATTTTCCTGGTTCTGCTGCTGGTGAGCGGCGTGATTATCTGGCGTCCTTATTTTGCGCCTGCTTTCTCAATCCCGGTGA  >  minE/2336817‑2336949
                                                                 |                                                                   
ggTGACACCGGGCGTTATAACTTCGGTCAGAAATGCGTTTTCTGGGCGGCGATTATTTTCCTGGTTctg                                                                  >  1:220214/1‑69 (MQ=255)
ggTGACACCGGGCGTTATAACTTCGGTCAGAAATGCGTTTTCTGGGCGGCGATTATTTTCCTGGTTatg                                                                  >  1:2603886/1‑69 (MQ=255)
   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:1676408/1‑68 (MQ=255)
                ataACTTCGGTCAGAAATGCGTTTTCTGGGCGGCGATTATTTTCCTGGTTCTGCTGCTGGTGAgcggc                                                   >  1:1910252/1‑68 (MQ=255)
                ataACTTCGGTCAGAAATGAGTTTTCTGGGCGGCGATTATTTTCCTGGTTCTGCTGCTGGTGAgagga                                                   >  1:1933644/1‑67 (MQ=255)
                   aCTTCGGTCAGAAATGCGTTTTCTGGGCGGCGATTATTTTCCTGGTCctgtt                                                                >  1:581632/1‑50 (MQ=255)
                   aCTTCGGTCAGAAATGCGTTTTCTGGGCGGCGATTATTTTCCTGGTCCTGTTGCTGGTGAGCGGCGTg                                                >  1:2211449/1‑68 (MQ=255)
                   aCTTCGGTCAGAAATGCGTTTTCTGGGCGGCGATTATTTTCCTGGTCCTGTTGCTGGTGAGCGGCGTg                                                >  1:2151632/1‑68 (MQ=255)
                   aCTTCGGTCAGAAATGCGTTTTCTGGGCGGCGATTATTTTCCTGGTCCTGTTGCTGGTGAGCGGCGTg                                                >  1:282419/1‑68 (MQ=255)
                   aCTTCGGTCAGAAATGCGTTTTCTGGGCGGCGATTATTTTCCTGGTCCTGTTGCTGGTGAGCGGCGTg                                                >  1:2772792/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:512159/1‑68 (MQ=255)
                   aCTTCGGTCAGAAATGCGTTTTCTGGGCGGCGATTATTTTCCTGGTCCTGTTGCTGGTGAGCGGCGTg                                                >  1:1727216/1‑68 (MQ=255)
                   aCTTCGGTCAGAAATGCGTTTTCTGGGCGGCGATTATTTTCCTGGTCCTGTTGCTGGTGAGCGGCGTg                                                >  1:2646071/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)
                   aCTTCGGTCAGAAATGCGTTTTCTGGGCGGCGATTATTTTCCAGGTCctgttgctg                                                            >  1:1685492/1‑56 (MQ=255)
                     ttCGGTCAGAAATGCGTTTTCTGGGCGGCGATTATTTTCCTGGTTCTGCTGCTGGTGAGCGGCGTGAtt                                             <  1:2688084/69‑1 (MQ=255)
                     ttCGGTCAGAAATGCGTTTTCTGGGCGGCGATTATTTTCCTGGTTCTGCTGCTGGTGAGCGGCGTGAtt                                             <  1:2360844/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:796991/1‑68 (MQ=255)
                                aTGCGTTTTCTGGGCGGCGATTATTTTCCTGGTTCTGCTGCTGGTGAGCGGCGTGATTATCTGGCGTc                                   >  1:3153827/1‑68 (MQ=255)
                                aTGCGTTTTCTGGGCGGCGATTATTTTCCTGGTTCTGCTGCTGGTGAGCGGCGTGATTATCTGGCGTc                                   >  1:2977858/1‑68 (MQ=255)
                                aTGCGTTTTCTGGGCGGCGATTATTTTCCTGGTTCTGCTGCTGGTGAGCGGCGTGATTATCTGGCGTc                                   >  1:2218338/1‑68 (MQ=255)
                                aTGCGTTTTCTGGGCGGCGATTATTTTCCTGGTTCTGCTGCTGGTGAGCGGCGTGATTATCTGGCGTc                                   >  1:2237141/1‑68 (MQ=255)
                                aTGCGTTTTCTGGGCGGCGATTATTTTCCTGGTTCTGCTGCTGGTGAGCGGCGTGATTATCTGGCGTc                                   >  1:3180581/1‑68 (MQ=255)
                                aTGCGTTTTCTGGGCGGCGATTATTTTACTGGTTCTGCTGCTGGTGAGCGGCGTGATTATCTGGCGTc                                   >  1:1346296/1‑68 (MQ=255)
                                 tGCGTTTTCTGGGCGGCGATTATTTTCCTGGTTCTGCTGCTGGTGAGCGGCGTGATTATCTGGCGTCCt                                 <  1:392234/69‑1 (MQ=255)
                                 tGCGTTTTCTGGGCGGCGATTATTTTCCTGGTTCTGCTGCTGGTGAGCGGCGTGATTATCTGGCGTCCt                                 <  1:225354/69‑1 (MQ=255)
                                 tGCGTTTTCTGGGCGGCGATTATTTTCCTGGTTCTGCTGCTGGTGAGCGGCGTGATTATCTGGCGTCCt                                 <  1:1555000/69‑1 (MQ=255)
                                 tGCGTTTTCTGGGCGGCGATTATTTTCCTGGTTCTGCTGCTGGTGAGCGGCGTGATTATCTGGCGTCCt                                 <  1:1564487/69‑1 (MQ=255)
                                 tGCGTTTTCTGGGCGGCGATTATTTTCCTGGTTCTGCTGCTGGTGAGCGGCGTGATTATCTGGCGTCCt                                 <  1:3080611/69‑1 (MQ=255)
                                  gCGTTTTCTGGGCGGCGATTATTTTCCTGGTTctgct                                                                >  1:1728667/1‑37 (MQ=255)
                                  gCGTTTTCTGGGCGGCGATTATTTTCCTGGTTctgct                                                                >  1:3431555/1‑37 (MQ=255)
                                  gCGTTTTCTGGGCGGCGATTATTTTCCTGGTTctgct                                                                >  1:2350957/1‑37 (MQ=255)
                                  gCGTTTTCTGGGCGGCGATTATTTTCCTGGTTctgct                                                                >  1:3161674/1‑37 (MQ=255)
                                  gCGTTTTCTGGGCGGCGATTATTTTCCTGGTTctgct                                                                >  1:1766628/1‑37 (MQ=255)
                                  gCGTTTTCTGGGCGGCGATTATTTTCCTGGTTctgct                                                                >  1:173776/1‑37 (MQ=255)
                                  gCGTTTTCTGGGCGGCGATTATTTTCCTGGTTctgct                                                                >  1:2264573/1‑37 (MQ=255)
                                  gCGTTTTCTGGGCGGCGATTATTTTCCTGGTTctgct                                                                >  1:2530390/1‑37 (MQ=255)
                                  gCGTTTTCTGGGCGGCGATTATTTTCCTGGTTCTGCTGCTGGTGAGCGGCGTGATTATCTGGCGTCCt                                 <  1:2608864/68‑1 (MQ=255)
                                  gCGTTTTCTGGGCGGCGATTATTTTCCTGGTTCTGCTGCTGGGGAGCGGCGTGATTATCTGGCGTCCt                                 <  1:2006062/68‑1 (MQ=255)
                                   cGTTTTCTGGGCGGCGATTATTTTCCTGGTTCTGCTGCTGGTGAGCGGCGTGATTATCTGGCGTCCTTa                               <  1:1619979/69‑1 (MQ=255)
                                   cGTTTTCTGGGCGGCGATTATTTTCCTGGTTCTGCTGCTGGTGAGCGGCGTGATTATCTGGCGTCCTTa                               <  1:2262967/69‑1 (MQ=255)
                                   cGTTTTCTGGGCGGCGATTATTTTCCTGGTTCTGCTGCTGGTGAGCGGCGTGATTATCTGGCGTCCTTa                               <  1:1893319/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:229472/69‑1 (MQ=255)
                                                     tattTTCCTGGTTCTGCTGCTGGTGAGCGGCGTGATTATCTGGCGTCCTTATTTTGCGCCTGCTTtctc             <  1:1894476/69‑1 (MQ=255)
                                                     tattTTCCTGGTTCTGCTGCTGGTGAGCGGCGTGATTATCTGGCGTCCTTATTTTGCGCCTGCTTtctc             <  1:1961974/69‑1 (MQ=255)
                                                     tattTTCCTGGTTCTGCTGCTGGTGAGCGGCGTGATTATCTGGCGTCCTTATTTTGCGCCTGCTTtctc             <  1:95222/69‑1 (MQ=255)
                                                     tattTTCCTGGTCCTGTTGCTGGTGAGCGGCGTGATTATCTGGCGTCCTTATTTTGCGCCTGCTTtctc             <  1:657487/69‑1 (MQ=255)
                                                     tattTTCCTGGTCCTGTTGCTGGTGAGCGGCGTGATTATCTGGCGTCCTTATTTTGCGCCTGCTTtctc             <  1:904348/69‑1 (MQ=255)
                                                     tattTTCCTGGTCCTGTTGCTGGTGAGCGGCGTGATTATCTGGCGTCCTTATTTTGCGCCTGCTTtctc             <  1:851902/69‑1 (MQ=255)
                                                     tattTTCCTGGTCCTGTTGCTGGTGAGCGGCGTGATTATCTGGCGTCCTTATTTTGCGCCTGCTTtctc             <  1:264532/69‑1 (MQ=255)
                                                     tattTTCCTGGTCCTGTTGCTGGTGAGCGGCGTGATTATCTGGCGTCCTTATTTTGCGCCTGCTTtctc             <  1:2724387/69‑1 (MQ=255)
                                                     tattTTCCTGGTCCTGTTGCTGGTGAGCGGCGTGATTATCTGGCGTCCTTATTTTGCGCCTGCTTtctc             <  1:710984/69‑1 (MQ=255)
                                                     tattTTCCTGGTCCTGTTGCTGGTGAGCGGCGTGATTATCTGGCGTCCTTATTTTGCGCCTGCTTtctc             <  1:388940/69‑1 (MQ=255)
                                                     tattTTCCTGGTCCTGTTGCTGGTGAGCGGCGTGATTATCTGGCGTCCTTATTTTGCGCCTGCTTtctc             <  1:2738229/69‑1 (MQ=255)
                                                     tattTTCCTGGTCCTGTTGCTGGTGAGCGGCGTGATTATCTGGCGTCCTTATTTTGCGCCTGCTTtctc             <  1:615183/69‑1 (MQ=255)
                                                     tattTTCCTGGTCCTGTTGCTGGTGAGCGGCGTGATTATCTGGCGTCCTTATTTTGCGCCTGCTTtctc             <  1:3116440/69‑1 (MQ=255)
                                                     tattTTCCTGGTCCTGTTGCTGGTGAGCGGCGTGATTATCTGGCGTCCTTATTTTGCGCCTGCTTtctc             <  1:3173907/69‑1 (MQ=255)
                                                     tattTTCCTGGTCCTGTTGCTGGTGAGCGGCGTGATTATCTGGCGTCCTTATTTTGCGCCTGCTTtctc             <  1:3174355/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)
                                                                 |                                                                   
GGTGACACCGGGCGTTATAACTTCGGTCAGAAATGCGTTTTCTGGGCGGCGATTATTTTCCTGGTTCTGCTGCTGGTGAGCGGCGTGATTATCTGGCGTCCTTATTTTGCGCCTGCTTTCTCAATCCCGGTGA  >  minE/2336817‑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: