Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,384,489 G→A 57.3% L358L (CTG→TTG)  metE ← 5‑methyltetrahydropteroyltriglutamate‑ homocysteine S‑methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,384,4890GA57.3% 18.5 / 57.3 75L358L (CTG→TTG) metE5‑methyltetrahydropteroyltriglutamate‑ homocysteine S‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (7/25);  new base A (20/23);  total (27/48)
Fisher's exact test for biased strand distribution p-value = 3.19e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCA  >  minE/2384422‑2384544
                                                                   |                                                       
gCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGGCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGt                                                        <  1:1043159/69‑1 (MQ=255)
gCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCTGTGTCACCACTGTTCAGCGCATCGCGCAGCAGt                                                        <  1:1806855/69‑1 (MQ=255)
gCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGt                                                        <  1:3180275/69‑1 (MQ=255)
gCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGt                                                        <  1:1564577/69‑1 (MQ=255)
gCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGt                                                        <  1:1474263/69‑1 (MQ=255)
 ccTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGt                                                        <  1:992183/68‑1 (MQ=255)
 ccTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGt                                                        >  1:460003/1‑68 (MQ=255)
 ccTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGt                                                        >  1:2570648/1‑68 (MQ=255)
 ccTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGt                                                        >  1:2025025/1‑68 (MQ=255)
   tGAATCGTAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATCGCGCAGCAATGcc                                                     <  1:2456861/69‑1 (MQ=255)
   tGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACTCATCGCGCAGCAATGcc                                                     <  1:2691258/69‑1 (MQ=255)
   tGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATCGCGCAGCAATGcc                                                     <  1:2002631/69‑1 (MQ=255)
   tGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATCGCGCAGCAATGcc                                                     <  1:2670592/69‑1 (MQ=255)
   tGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATCGCGCAGCAATGcc                                                     <  1:2832076/69‑1 (MQ=255)
   tGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATCGCGCAGCAATGcc                                                     <  1:2977538/69‑1 (MQ=255)
   tGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATCGCGCAGCAATGcc                                                     <  1:3173104/69‑1 (MQ=255)
   tGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATCGCGCAGCAATGcc                                                     <  1:2254476/69‑1 (MQ=255)
   tGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATCGCGCAGCAATGcc                                                     <  1:2225330/69‑1 (MQ=255)
   tGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATCGCGCAGCAATGcc                                                     <  1:2121573/69‑1 (MQ=255)
   tGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATCGCGCAGCAATGcc                                                     <  1:2080915/69‑1 (MQ=255)
   tGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATCGCGCAGCAATGcc                                                     <  1:371147/69‑1 (MQ=255)
   tGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATCGCGCAGCAATGcc                                                     <  1:1992821/69‑1 (MQ=255)
   tGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATCGCGCAGCAATGcc                                                     <  1:1931304/69‑1 (MQ=255)
   tGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATCGCGCAGCAATGcc                                                     <  1:375238/69‑1 (MQ=255)
   tGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATCGCGCAGCAATGcc                                                     <  1:1544387/69‑1 (MQ=255)
   tGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATCGCGCAGCAATGcc                                                     <  1:1873636/69‑1 (MQ=255)
   tGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATCGCGCAGCAATGcc                                                     <  1:537219/69‑1 (MQ=255)
   tGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATCGCGCAGCAATGcc                                                     <  1:1323286/69‑1 (MQ=255)
   tGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATCGCGCAGCAATGcc                                                     <  1:603508/69‑1 (MQ=255)
   tGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATCGCGCAGCAATGcc                                                     <  1:675024/69‑1 (MQ=255)
   tGAATCGGAGCGCTACACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATCGCGCAGCAATGcc                                                     <  1:898860/69‑1 (MQ=255)
    gAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATCGCGCAGCAATGcc                                                     <  1:1388945/68‑1 (MQ=255)
       tCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGTAGTGCCAGtt                                                 <  1:743517/69‑1 (MQ=255)
       tCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGtt                                                 <  1:3268265/69‑1 (MQ=255)
       tCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGtt                                                 <  1:1060171/69‑1 (MQ=255)
       tCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGtt                                                 <  1:3389613/69‑1 (MQ=255)
       tCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGtt                                                 <  1:1042581/69‑1 (MQ=255)
       tCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGtt                                                 <  1:365368/69‑1 (MQ=255)
       tCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGtt                                                 <  1:1054153/69‑1 (MQ=255)
       tCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGtt                                                 <  1:1325541/69‑1 (MQ=255)
       tCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGtt                                                 <  1:1805521/69‑1 (MQ=255)
       tCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGtt                                                 <  1:1812269/69‑1 (MQ=255)
       tCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGtt                                                 <  1:2129274/69‑1 (MQ=255)
       tCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGtt                                                 <  1:2236404/69‑1 (MQ=255)
        cGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGtt                                                 <  1:1946987/68‑1 (MQ=255)
        cGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGtt                                                 <  1:2515399/68‑1 (MQ=255)
                                 gCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCa                       <  1:265513/69‑1 (MQ=255)
                                 gCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCa                       <  1:914856/69‑1 (MQ=255)
                                 gCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCa                       <  1:723924/69‑1 (MQ=255)
                                 gCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCa                       <  1:65973/69‑1 (MQ=255)
                                         accacTGTTCAACGCATCGCGCCGCAATGCCAGTTCATGACATTTTTGCAGGGCGAAGGCAAACCAGCt               >  1:541092/1‑69 (MQ=255)
                                         accacTGTTCAACGCATCGCGCAGCAATGCCAGTTCATGACATTTTTGCAGGGCGAAGGCAAACCAGCt               >  1:1866971/1‑69 (MQ=255)
                                         accacTGTTCAACGCATCGCGCAGCAATGCCAGTTCATGACATTTTTGCAGGGCGAAGGCAAACCAGCt               >  1:1926499/1‑69 (MQ=255)
                                         accacTGTTCAACGCATCGCGCAGCAATGCCAGTTCATGACATTTTTGCAGGGCGAAGGCAAACCAGCt               >  1:621657/1‑69 (MQ=255)
                                         accacTGTTCAACGCATCGCGCAGCAATGCCAGTTCATGACATTTTTGCAGGGCGAAGGCAAACCAGCt               >  1:652246/1‑69 (MQ=255)
                                         accacTGTTCAACGCATCGCGCAGCAATGCCAGTTCATGACATTTTTGCAGGGCGAAGGCAAACCAGCt               >  1:698449/1‑69 (MQ=255)
                                         accacTGTTCAACGCATCGCGCAGCAATGCCAGTTCATGACATTTTTGCAGGGCGAAGGCAAACCAGCt               >  1:1778753/1‑69 (MQ=255)
                                         accacTGTTCAACGCATCGCGCAGCAATGCCAGTTCATGACATTTTTGCAGGGCGAAGGCAAACCAGCt               >  1:95710/1‑69 (MQ=255)
                                         accacTGTTCAACGCATCGCGCAGCAATGCCAGTTCATGACATTTTTGCAGGGCGAAGGCAAACCAGCt               >  1:989854/1‑69 (MQ=255)
                                         accacTGTTCAACGCATCGCGCAGCAATGCCAGTTCATGACATTTTTGCAGGGCGAAGGCAAACCAGCt               >  1:1606479/1‑69 (MQ=255)
                                         accacTGTTCAACGCATCGCGCAGCAATGCCAGTTCATGACATTTTTGCAGGGCGAAGGCAAACCAGCt               >  1:511112/1‑69 (MQ=255)
                                         accacTGTTCAACGCATCGCGCAGCAATGCCAGTTCATGACATTTTTGCAGGGCGAAGGCAAACCAGCt               >  1:4610/1‑69 (MQ=255)
                                         accacTGTTCAACGCATCGCGCAGCAATGCCAGTTCATGACATTTTTGCAGGGCGAAGGCAAACCAGCt               >  1:1098774/1‑69 (MQ=255)
                                         accacTGTTCAACGCATCGCGCAGCAATGCCAGTTCATGACATTTTTGCAGGGCGAAGGCAAACCAGCt               >  1:2235025/1‑69 (MQ=255)
                                         accacTGTTCAACGCATCGCGCAGCAATGCCAGTTCATGACATTTTTGCAGGGCGAAGGCAAACCAGCt               >  1:3163558/1‑69 (MQ=255)
                                         accacTGTTCAACGCATCGCGCAGCAATGCCAGTTCATGACATTTTTGCAGGGCGAAGGCAAACCAGCt               >  1:2864477/1‑69 (MQ=255)
                                         accacTGTTCAACGCATCGCGCAGCAATGCCAGTTCATGACATTTTTGCAGGGCGAAGGCAAACCAGCt               >  1:2846979/1‑69 (MQ=255)
                                         accacTGTTCAACGCATCGCGCAGCAATGCCAGTTCATGACATTTTTGCAGGGCGAAGGCAAACCAGCt               >  1:1072375/1‑69 (MQ=255)
                                         accacTGTTCAACGCATCGCGCAGCAATGCCAGTTCATGACATTTTTGCAGGGCGAAGGCAAACCAGCt               >  1:2587076/1‑69 (MQ=255)
                                         accacTGTTCAACGCATCGCGCAGCAATGCCAGTTCATGACATTATTGCAGGGCGAAGGCAAACCAGCt               >  1:1183079/1‑69 (MQ=255)
                                            acTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCtttt            >  1:2691100/1‑69 (MQ=255)
                                            acTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCtttt            >  1:2687263/1‑69 (MQ=255)
                                               gTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTc           <  1:2385407/67‑1 (MQ=255)
                                                         tCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCa  >  1:2338028/1‑66 (MQ=255)
                                                         tCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCa  >  1:2672639/1‑66 (MQ=255)
                                                                   |                                                       
GCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCA  >  minE/2384422‑2384544

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: