Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,384,460 G→A 61.1% D367D (GAC→GAT metE ← 5‑methyltetrahydropteroyltriglutamate‑ homocysteine S‑methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,384,4600GA61.1% 51.0 / 73.9 108D367D (GAC→GATmetE5‑methyltetrahydropteroyltriglutamate‑ homocysteine S‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (7/35);  new base A (43/23);  total (50/58)
Fisher's exact test for biased strand distribution p-value = 6.30e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

Only 100 of 108 total aligned reads displayed.

TTATGTACGCGGGTAGAGTGACGACGTGCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCA  >  minE/2384395‑2384523
                                                                 |                                                               
ttATGTACGCGGGTAGAGTGACGACGTGCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTc                                                               >  1:2688981/1‑68 (MQ=255)
 tATGTACGCGGGTAGAGTGACGACGTGCCTGAATCGGGGCGCTCCACTCTTCCAGAGCTGCCGTGTcac                                                             <  1:889212/69‑1 (MQ=255)
 tATGTACGCGGGTAGAGTGACGACGTGCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTcac                                                             <  1:2114753/69‑1 (MQ=255)
       cgcgGGTAGAGTGACGACGTGCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTg                                                        >  1:537693/1‑68 (MQ=255)
       cgcgGGTAGAGTGACGACGTGCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTg                                                        >  1:1100307/1‑68 (MQ=255)
       cgcgGGTAGAGTGACGACGTGCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTg                                                        >  1:684275/1‑68 (MQ=255)
          gagTAGAGTGACGACGCGCCTGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTc                                                     >  1:1675073/3‑68 (MQ=255)
          gagTAGAGTGACGACGCGCCTGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTc                                                     >  1:237381/3‑68 (MQ=255)
          gagTAGAGTGACGACGCGCCTGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTc                                                     >  1:2458030/3‑68 (MQ=255)
          gagTAGAGTGACGACGCGCCTGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTc                                                     >  1:1225691/3‑68 (MQ=255)
          gagTAGAGTGACGACGCGCCTGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTc                                                     >  1:2827367/3‑68 (MQ=255)
          gagTAGAGTGACGACGCGCCTGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCa                                                    >  1:3171284/3‑69 (MQ=255)
          gagTAGAGTGACGACGCGCCTGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCa                                                    >  1:2569832/3‑69 (MQ=255)
          gagTAGAGTGACGACGCGCCTGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCa                                                    >  1:2075871/3‑69 (MQ=255)
          gagTAGAGTGACGACGCGCCTGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCa                                                    >  1:2258163/3‑69 (MQ=255)
          gagTAGAGTGACGACGCGCCTGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCa                                                    >  1:2668622/3‑69 (MQ=255)
          gagTAGAGTGACGACGCGCCTGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCa                                                    >  1:2485390/3‑69 (MQ=255)
          gagTAGAGTGACGACGCGCCTGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCa                                                    >  1:345140/3‑69 (MQ=255)
          gagTAGAGTGACGACGCGCCTGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCa                                                    >  1:623790/3‑69 (MQ=255)
          gagTAGAGTGACGACGCGCCTGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCa                                                    >  1:2676352/3‑69 (MQ=255)
          gagTAGAGTGACGACGCGCCTGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCa                                                    >  1:1450387/3‑69 (MQ=255)
          gagTAGAGTGACGACGCGCCTGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCa                                                    >  1:2313119/3‑69 (MQ=255)
          gagTAGAGTGACGACGCGCCTGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCa                                                    >  1:818742/3‑69 (MQ=255)
          gagTAGAGTGACGACGCGCCTGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCa                                                    >  1:2350687/3‑69 (MQ=255)
           ggTAGAGTGACGACGTGCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAg                                                   <  1:353827/69‑1 (MQ=255)
            gTAGAGTGACGACGTGCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAg                                                   <  1:2748713/68‑1 (MQ=255)
                    acgacgTGCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATcgcg                                          <  1:58112/69‑1 (MQ=255)
                    acgacgTGCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATcgcg                                          <  1:316537/69‑1 (MQ=255)
                    acgacgTGCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATcgcg                                          <  1:2543751/69‑1 (MQ=255)
                     cgacgTGCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATcgcg                                          <  1:3134861/68‑1 (MQ=255)
                        cgTGCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCgcagc                                      <  1:3173894/69‑1 (MQ=255)
                        cgTGCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCgcagc                                      <  1:2750741/69‑1 (MQ=255)
                           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)
                            ccTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGt                                   >  1:2025025/1‑68 (MQ=255)
                            ccTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGt                                   >  1:2570648/1‑68 (MQ=255)
                            ccTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGt                                   >  1:460003/1‑68 (MQ=255)
                              tGAATCGTAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATCGCGCAGCAATGcc                                <  1:2456861/69‑1 (MQ=255)
                              tGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACTCATCGCGCAGCAATGcc                                <  1:2691258/69‑1 (MQ=255)
                              tGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATCGCGCAGCAATGcc                                <  1:2977538/69‑1 (MQ=255)
                              tGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATCGCGCAGCAATGcc                                <  1:2670592/69‑1 (MQ=255)
                              tGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATCGCGCAGCAATGcc                                <  1:2832076/69‑1 (MQ=255)
                              tGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATCGCGCAGCAATGcc                                <  1:1931304/69‑1 (MQ=255)
                              tGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATCGCGCAGCAATGcc                                <  1:1323286/69‑1 (MQ=255)
                              tGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATCGCGCAGCAATGcc                                <  1:1544387/69‑1 (MQ=255)
                              tGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATCGCGCAGCAATGcc                                <  1:675024/69‑1 (MQ=255)
                              tGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATCGCGCAGCAATGcc                                <  1:537219/69‑1 (MQ=255)
                              tGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATCGCGCAGCAATGcc                                <  1:1873636/69‑1 (MQ=255)
                              tGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATCGCGCAGCAATGcc                                <  1:3173104/69‑1 (MQ=255)
                              tGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATCGCGCAGCAATGcc                                <  1:1992821/69‑1 (MQ=255)
                              tGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATCGCGCAGCAATGcc                                <  1:2002631/69‑1 (MQ=255)
                              tGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATCGCGCAGCAATGcc                                <  1:375238/69‑1 (MQ=255)
                              tGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATCGCGCAGCAATGcc                                <  1:2080915/69‑1 (MQ=255)
                              tGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATCGCGCAGCAATGcc                                <  1:2121573/69‑1 (MQ=255)
                              tGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATCGCGCAGCAATGcc                                <  1:371147/69‑1 (MQ=255)
                              tGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATCGCGCAGCAATGcc                                <  1:2225330/69‑1 (MQ=255)
                              tGAATCGGAGCGCTACACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATCGCGCAGCAATGcc                                <  1:898860/69‑1 (MQ=255)
                               gAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATCGCGCAGCAATGcc                                <  1:1388945/68‑1 (MQ=255)
                                 aTCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATcg                                            >  1:915625/1‑54 (MQ=255)
                                 aTCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATcg                                            >  1:1030723/1‑54 (MQ=255)
                                 aTCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATcg                                            >  1:75447/1‑54 (MQ=255)
                                 aTCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATcg                                            >  1:58626/1‑54 (MQ=255)
                                 aTCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATcg                                            >  1:3045893/1‑54 (MQ=255)
                                 aTCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATcg                                            >  1:428748/1‑54 (MQ=255)
                                 aTCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATcg                                            >  1:2851264/1‑54 (MQ=255)
                                 aTCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATcg                                            >  1:1650004/1‑54 (MQ=255)
                                 aTCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATcg                                            >  1:1677807/1‑54 (MQ=255)
                                 aTCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATcg                                            >  1:1641733/1‑54 (MQ=255)
                                 aTCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATcg                                            >  1:1547434/1‑54 (MQ=255)
                                 aTCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATcg                                            >  1:1545984/1‑54 (MQ=255)
                                 aTCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATcg                                            >  1:2010083/1‑54 (MQ=255)
                                 aTCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATcg                                            >  1:1398881/1‑54 (MQ=255)
                                 aTCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATcg                                            >  1:2269511/1‑54 (MQ=255)
                                 aTCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATcg                                            >  1:2272616/1‑54 (MQ=255)
                                 aTCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATcg                                            >  1:2275473/1‑54 (MQ=255)
                                 aTCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATcg                                            >  1:2332204/1‑54 (MQ=255)
                                 aTCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATcg                                            >  1:1046215/1‑54 (MQ=255)
                                 aTCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATcg                                            >  1:2561458/1‑54 (MQ=255)
                                 aTCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATcg                                            >  1:2633215/1‑54 (MQ=255)
                                 aTCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATcg                                            >  1:2885997/1‑54 (MQ=255)
                                  tCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGTAGTGCCAGtt                            <  1:743517/69‑1 (MQ=255)
                                  tCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGtt                            <  1:1054153/69‑1 (MQ=255)
                                  tCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGtt                            <  1:1060171/69‑1 (MQ=255)
                                  tCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGtt                            <  1:1325541/69‑1 (MQ=255)
                                  tCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGtt                            <  1:1042581/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:365368/69‑1 (MQ=255)
                                  tCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGtt                            <  1:2236404/69‑1 (MQ=255)
                                  tCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGtt                            <  1:3389613/69‑1 (MQ=255)
                                  tCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGtt                            <  1:3268265/69‑1 (MQ=255)
                                   cGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGtt                            <  1:1946987/68‑1 (MQ=255)
                                   cGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGtt                            <  1:2515399/68‑1 (MQ=255)
                                                            gCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCa  <  1:65973/69‑1 (MQ=255)
                                                            gCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCa  <  1:723924/69‑1 (MQ=255)
                                                            gCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCa  <  1:914856/69‑1 (MQ=255)
                                                            gCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCa  <  1:265513/69‑1 (MQ=255)
                                                                 |                                                               
TTATGTACGCGGGTAGAGTGACGACGTGCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCA  >  minE/2384395‑2384523

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: