Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,385,126 G→C 30.9% D145E (GAC→GAG metE ← 5‑methyltetrahydropteroyltriglutamate‑ homocysteine S‑methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,385,1260GC30.9% 78.2 / 44.9 68D145E (GAC→GAGmetE5‑methyltetrahydropteroyltriglutamate‑ homocysteine S‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (20/27);  new base C (7/14);  total (27/41)
Fisher's exact test for biased strand distribution p-value = 5.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

CAGGTAACCGGCCCCAGCAGCACAGGTTTCACCTTGTGGCCCAGCGCCAGCGCCTCGTCCACTTCGTCCAGCAGCTGCGTCCAGGTCAGTTTGAACTGTTGGCCTTTAACGAACTCCGGCACCATGTA  >  minE/2385061‑2385188
                                                                 |                                                              
caGGTAACCGGCCCCAGCAGCACAGGTTTCACCTTGTGGCCCAGCGCCAGCGCCTCGTCCACTTCGTcc                                                             >  1:1025031/1‑69 (MQ=255)
caGGTAACCGGCCCCAGCAGCACAGGTTTCACCTTGTGGCCCAGCGCCAGCGCCTCGTCCACTTCGTcc                                                             >  1:1180535/1‑69 (MQ=255)
caGGTAACCGGCCCCAGCAGCACAGGTTTCACCTTGTGGCCCAGCGCCAGCGCCTCGTCCACTTCGTcc                                                             >  1:2209657/1‑69 (MQ=255)
caGGTAACCGGCCCCAGCAGCACAGGTTTCACCTTGTGGCCCAGCGCCAGCGCCTCGTCCACTTCGTcc                                                             >  1:1306977/1‑69 (MQ=255)
                        ggTTTCACCTTGTGGCCCAGCGCCAGCGCCTCGTCCACTTCGTCCAGCAGCTGCGTCCAGGTCAGTTTg                                     <  1:2169414/69‑1 (MQ=255)
                        ggTTTCACCTTGTGGCCCAGCGCCAGCGCCTCGTCCACTTCGTCCAGCAGCTGCGTCCAGGTCAGTTTg                                     <  1:3422807/69‑1 (MQ=255)
                        ggTTTCACCTTGTGGCCCAGCGCCAGCGCCTCGTCCACTTCGTCCAGCAGCTGCGTCCAGGTCAGTTTg                                     <  1:1291833/69‑1 (MQ=255)
                        ggTTTCACCTTGTGGCCCAGCGCCAGCGCCTCGTCCACTTCGTCCAGCAGCTGCGTCCAGGTCAGTTTg                                     <  1:1947363/69‑1 (MQ=255)
                        ggTTTCACCTTGTGGCCCAGCGCCAGCGCCTCGTCCACTTCGTCCAGCAGCTGCGTCCAGGTCAGTTTg                                     <  1:3155107/69‑1 (MQ=255)
                        ggTTTCACCTTGTGGCCCAGCGCCAGCGCCTCGTCCACTTCGTCCAGCAGCTGCGTCCAGGTCAGTTTg                                     <  1:3332092/69‑1 (MQ=255)
                          tttCACCTTGTGGCCCAGCGCCAGCGCCTCGTCCACTTCGTCCAGCAGCTGCGTCCAGGTCAGTTTGaa                                   <  1:597868/69‑1 (MQ=255)
                            tCACCTTGTGGCCCAGCGCCAGCGCCTCGTCCACTTCCTCCAGCAGCTGCGTCCAGGTCAGTTTGAACt                                 <  1:2182407/69‑1 (MQ=255)
                            tCACCTTGTGGCCCAGCGCCAGCGCCTCGTCCACTTCCTCCAGCAGCTGCGTCCAGGTCAGTTTGAACt                                 <  1:2331271/69‑1 (MQ=255)
                            tCACCTTGTGGCCCAGCGCCAGCGCCTCGTCCACTTCCTCCAGCAGCTGCGTCCAGGTCAGTTTGAACt                                 <  1:3400244/69‑1 (MQ=255)
                            tCACCTTGTGGCCCAGCGCCAGCGCCTCGTCCACTTCCTCCAGCAGCTGCGTCCAGGTCAGTTTGAACt                                 <  1:3385107/69‑1 (MQ=255)
                            tCACCTTGTGGCCCAGCGCCAGCGCCTCGTCCACTTCCTCCAGCAGCTGCGTCCAGGTCAGTTTGAACt                                 <  1:1686033/69‑1 (MQ=255)
                            tCACCTTGTGGCCCAGCGCCAGCGCCTCGTCCACTTCCTCCAGCAGCTGCGTCCAGGTCAGTTTGAACt                                 <  1:2756724/69‑1 (MQ=255)
                            tCACCTTGTGGCCCAGCGCCAGCGCCTCGTCCACTTCCTCCAGCAGCTGCGTCCAGGTCAGTTTGAACt                                 <  1:845692/69‑1 (MQ=255)
                            tCACCTTGTGGCCCAGCGCCAGCGCCTCGTCCACTTCCTCCAGCAGCTGCGTCCAGGTCAGTTTGAACt                                 <  1:1431298/69‑1 (MQ=255)
                            tCACCTTGTGGCCCAGCGCCAGCGCCTCGTCCACTTCCTCCAGCAGCTGCGTCCAGGTCAGTTTGAACt                                 <  1:1343535/69‑1 (MQ=255)
                            tCACCTTGTGGCCCAGCGCCAGCGCCTCGTCCACTTCCTCCAGCAGCTGCGTCCAGGTCAGTTTGAACt                                 <  1:952757/69‑1 (MQ=255)
                            tCACCTTGTGGCCCAGCGCCAGCGCCTCGTCCACTTCCTCCAGCAGCTGCGTCCAGGTCAGTTTGAACt                                 <  1:1212536/69‑1 (MQ=255)
                            tCACCTTGTGGCCCAGCGCCAGCGCCTCGTCCACTTCCTCCAGCAGCTGCGTCCAGGTCAGTTTGAACt                                 <  1:995536/69‑1 (MQ=255)
                            tCACCTTGTGGCCCAGCGCCAGCGCCTCGTCCACTTCCTCCAGCAGCTGCGTCCAGGTCAGTTTGAACt                                 <  1:1093534/69‑1 (MQ=255)
                                  tgtgGCCCAGCGCCAGCGCCTCGTCCACTTCGTCCAGCAGCTGCGTCCAGGTCAGTTTGAACTGTTGGc                           <  1:945201/69‑1 (MQ=255)
                                  tgtgGCCCAGCGCCAGCGCCTCGTCCACTTCGTCCAGCAGCTGCGTCCAGGTCAGTTTGAACTGTTGGc                           <  1:377160/69‑1 (MQ=255)
                                  tgtgGCCCAGCGCCAGCGCCTCGTCCACTTCGTCCAGCAGCTGCGTCCAGGTCAGTTTGAACTGTTGGc                           <  1:1182318/69‑1 (MQ=255)
                                  tgtgGCCCAGCGCCAGCGCCTCGTCCACATCGTCCAGCAGCTGCGTCCAGGTCAGTTTGAACTGTTGGc                           <  1:1904262/69‑1 (MQ=255)
                                      gCCCAGCGCCAGCGCCTCGTCCACTTCGTCCAGCAGCTGCGTCCAGGTCAGTTTGAACTGTTGGCCttt                       <  1:2844722/69‑1 (MQ=255)
                                      gCCCAGCGCCAGCGCCTCGTCCACTTCGTCCAGCAGCTGCGTCCAGGTCAGTTTGAACTGTTGGCCttt                       <  1:505130/69‑1 (MQ=255)
                                      gCCCAGCGCCAGCGCCTCGTCCACTTCGTCCAGCAGCTGCGTCCAGGTCAGTTTGAACTGTTGGCCttt                       <  1:791418/69‑1 (MQ=255)
                                      gCCCAGCGCCAGCGCCTCGTCCACTTCGTCCAGCAGCTGCGTCCAGGTCAGTTTGAACTGTTGGCCttt                       <  1:1670752/69‑1 (MQ=255)
                                        ccagcgccagcgCCTCGTCCACTTCGTCCAGCAGCTGCGTCCAGGTCAGTTTGAACTGTTGGCCTTTa                      >  1:2151397/1‑68 (MQ=255)
                                        ccagcgccagcgCCTCGTCCACTTCGTCCAGCAGCTGCGTCCAGGTCAGTTTGAACTGTTGGCCTTTa                      >  1:1950263/1‑68 (MQ=255)
                                        ccagcgccagcgCCTCGTCCACTTCGTCCAGCAGCTGCGTCCAGGTCAGTTTGAACTGTTGGCCTTTa                      >  1:1841008/1‑68 (MQ=255)
                                        ccagcgccagcgCCTCGTCCACTTCGTCCAGCAGCTGCGTCCAGGTCAGTTTGAACTGTTGGCCTTTa                      >  1:85687/1‑68 (MQ=255)
                                        ccagcgccagcgCCTCGTCCACTTCGTCCAGCAGCTGCGTCCAGGTCAGTTTGAACTGTTGGCCGTTa                      >  1:2569150/1‑68 (MQ=255)
                                          agcgccagcgccTCGTCCACTTCGTCCAGCAGCTGCGTCCAGGTCAGTTTGaa                                   <  1:1083156/53‑1 (MQ=255)
                                          agcgccagcgccTCGTCCACTTCGTCCAGCAGCTGCGTCCAGGTCAGTTTGAACTGTTGGCCTTTAACg                   <  1:3043347/69‑1 (MQ=255)
                                           gcgccagcgccTCGTCCACTTCCTCCAGCAGCTGCGTCCAGGTCAg                                         >  1:2876273/1‑46 (MQ=255)
                                           gcgccagcgccTCGTCCACTTCCTCCAGCAGCTGCGTCCAGGTCAg                                         >  1:67248/1‑46 (MQ=255)
                                           gcgccagcgccTCGTCCACTTCCTCCAGCAGCTGCGTCCAGGTCAg                                         >  1:540693/1‑46 (MQ=255)
                                           gcgccagcgccTCGTCCACTTCCTCCAGCAGCTGCGTCCAGGTCAg                                         >  1:2781088/1‑46 (MQ=255)
                                           gcgccagcgccTCGTCCACTTCCTCCAGCAGCTGCGTCCAGGTCAg                                         >  1:866166/1‑46 (MQ=255)
                                           gcgccagcgccTCGTCCACTTCCTCCAGCAGCTGCGTCCAGGTCAg                                         >  1:726766/1‑46 (MQ=255)
                                           gcgccagcgccTCGTCCACTTCCTCCAGCAGCTGCGTCCAGGTCAg                                         >  1:370445/1‑46 (MQ=255)
                                             gccagcgccTCGTCCACTTCGTCCAGCAGCTGCGTCCAGGTCAGTTTGAACTGTTGGCCTTTAACGAAc                >  1:1552509/1‑69 (MQ=255)
                                             gccagcgccTCGTCCACTTCGTCCAGCAGCTGCGTCCAGGTCAGTTTGAACTGTTGGCCTTTAACGAAc                >  1:2468982/1‑69 (MQ=255)
                                                   gccTCGTCCACTTCGTCCAGCAGCTGCGTCCAGGTCAg                                         >  1:1892628/1‑38 (MQ=255)
                                                   gccTCGTCCACTTCGTCCAGCAGCTGCGTCCAGGTCAGTTTGAACTGTTGGCCTTTAACGAACTCCGGc          >  1:1768429/1‑69 (MQ=255)
                                                   gccTCGTCCACTTCGTCCAGCAGCTGCGTCCAGGTCAGTTTGAACTGTTGGCCTTTAACGAACTCCGGc          >  1:1887391/1‑69 (MQ=255)
                                                   gccTCGTCCACTTCGTCCAGCAGCTGCGTCCAGGTCAGTTTGAACTGTTGGCCTTTAACGAACTCCGGc          >  1:3148931/1‑69 (MQ=255)
                                                   gccTCGTCCACTTCGTCCAGCAGCTGCGTCCAGGTCAGTTTGAACTGTTGGCCTTTAACGAACTCCGGc          >  1:1828495/1‑69 (MQ=255)
                                                      tCGTCCACTTCGTCCAGCAGCTGCGTCCAGGTCAGTTTGAACTGTTGGCCttt                       <  1:1196508/53‑1 (MQ=255)
                                                      tCGTCCACTTCGTCCAGCAGCTGCGTCCAGGTCAGTTTGAACTGTTGGCCttt                       <  1:1488366/53‑1 (MQ=255)
                                                      tCGTCCACTTCGTCCAGCAGCTGCGTCCAGGTCAGTTTGAACTGTTGGCCttt                       <  1:1323765/53‑1 (MQ=255)
                                                      tCGTCCACTTCGTCCAGCAGCTGCGTCCAGGTCAGTTTGAACTGTTGGCCttt                       <  1:659541/53‑1 (MQ=255)
                                                      tCGTCCACTTCGTCCAGCAGCTGCGTCCAGGTCAGTTTGAACTGTTGGCCttt                       <  1:643897/53‑1 (MQ=255)
                                                      tCGTCCACTTCGTCCAGCAGCTGCGTCCAGGTCAGTTTGAACTGTTGGCCttt                       <  1:501739/53‑1 (MQ=255)
                                                      tCGTCCACTTCGTCCAGCAGCTGCGTCCAGGTCAGTTTGAACTGTTGGCCttt                       <  1:2154755/53‑1 (MQ=255)
                                                      tCGTCCACTTCGTCCAGCAGCTGCGTCCAGGTCAGTTTGAACTGTTGGCCttt                       <  1:2549422/53‑1 (MQ=255)
                                                      tCGTCCACTTCGTCCAGCAGCTGCGTCCAGGTCAGTTTGAACTGTTGGCCttt                       <  1:2974912/53‑1 (MQ=255)
                                                        gTCCACTTCCTCCAGCAGCTTCGTCCAGGTCAGTTTGAACt                                 <  1:900458/41‑1 (MQ=38)
                                                           cACTTCGTCCAGCAGCTGCGTCCAGGTCAGTTTGAACTGTTGGCCTTTAACGAACTCCGGCACCATGTa  >  1:657054/1‑69 (MQ=255)
                                                           cACTTCGTCCAGCAGCTGCGTCCAGGTCAGTTTGAACTGTTGGCCTTTAACGAACTCCGGCACCATGTa  >  1:499332/1‑69 (MQ=255)
                                                           cACTTCGTCCAGCAGCTGCGTCCAGGTCAGTTTGAACTGTTGGCCTTTAACGAACTCCGGCACCATGTa  >  1:1409826/1‑69 (MQ=255)
                                                           cACTTCGTCCAGCAGCTGCGTCCAGGTCAGTTTGAACTGTTGGCCTTTAACGAACTCCGGCACCATGTa  >  1:3047572/1‑69 (MQ=255)
                                                             cTTCGTCCAGCAGCTGCGTCCAGGTCAGTTTGAACTGTTGGc                           <  1:2610550/42‑1 (MQ=255)
                                                                 |                                                              
CAGGTAACCGGCCCCAGCAGCACAGGTTTCACCTTGTGGCCCAGCGCCAGCGCCTCGTCCACTTCGTCCAGCAGCTGCGTCCAGGTCAGTTTGAACTGTTGGCCTTTAACGAACTCCGGCACCATGTA  >  minE/2385061‑2385188

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: