Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,384,433 G→A 62.4% A376A (GCC→GCT metE ← 5‑methyltetrahydropteroyltriglutamate‑ homocysteine S‑methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,384,4330GA62.4% 58.2 / 71.2 109A376A (GCC→GCTmetE5‑methyltetrahydropteroyltriglutamate‑ homocysteine S‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (7/34);  new base A (45/23);  total (52/57)
Fisher's exact test for biased strand distribution p-value = 5.32e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.12e-01

Only 100 of 109 total aligned reads displayed.

CAGACGCTTTTCTACCGCCGGATTATGTACGCGGGTAGAGTGACGACGTGCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTT  >  minE/2384373‑2384497
                                                            |                                                                
cAGACGCTTTTCTACCGCCGGATTATGTACGCGGGTAGAGTGACGACGTGCCTGAATCGGGGCGCTCCa                                                          <  1:746934/69‑1 (MQ=255)
 aGACGCTTTTCTACCGCCGGATTATGTACGCGGGTAGAGTGACGACGTGCCTGAATCGGGGCGCTCCa                                                          <  1:509640/68‑1 (MQ=255)
      cTTTTCTACCGCCGGATTATGTACTCGGGTAGAGTGACGACGTGCCTGAATCGGGGCGCTCCACTCTGc                                                    <  1:2449721/69‑1 (MQ=255)
                      ttATGTACGCGGGTAGAGTGACGACGTGCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTc                                     >  1:2688981/1‑68 (MQ=255)
                       tATGTACGCGGGTAGAGTGACGACGTGCCTGAATCGGGGCGCTCCACTCTTCCAGAGCTGCCGTGTcac                                   <  1:889212/69‑1 (MQ=255)
                       tATGTACGCGGGTAGAGTGACGACGTGCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTcac                                   <  1:2114753/69‑1 (MQ=255)
                             cgcgGGTAGAGTGACGACGTGCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTg                              >  1:1100307/1‑68 (MQ=255)
                             cgcgGGTAGAGTGACGACGTGCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTg                              >  1:537693/1‑68 (MQ=255)
                             cgcgGGTAGAGTGACGACGTGCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTg                              >  1:684275/1‑68 (MQ=255)
                                gagTAGAGTGACGACGCGCCTGCATCGGAGCGCTCCActc                                                       >  1:3309604/3‑40 (MQ=21)
                                gagTAGAGTGACGACGCGCCTGAATCGGAGCGCTCCACTCTg                                                     >  1:677050/3‑42 (MQ=37)
                                gagTAGAGTGACGACGCGCCTGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTc                           >  1:1225691/3‑68 (MQ=255)
                                gagTAGAGTGACGACGCGCCTGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTc                           >  1:2827367/3‑68 (MQ=255)
                                gagTAGAGTGACGACGCGCCTGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTc                           >  1:192924/3‑68 (MQ=255)
                                gagTAGAGTGACGACGCGCCTGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTc                           >  1:2458030/3‑68 (MQ=255)
                                gagTAGAGTGACGACGCGCCTGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTc                           >  1:1675073/3‑68 (MQ=255)
                                gagTAGAGTGACGACGCGCCTGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCa                          >  1:2258163/3‑69 (MQ=255)
                                gagTAGAGTGACGACGCGCCTGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCa                          >  1:2075871/3‑69 (MQ=255)
                                gagTAGAGTGACGACGCGCCTGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCa                          >  1:2313119/3‑69 (MQ=255)
                                gagTAGAGTGACGACGCGCCTGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCa                          >  1:2350687/3‑69 (MQ=255)
                                gagTAGAGTGACGACGCGCCTGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCa                          >  1:2525020/3‑69 (MQ=255)
                                gagTAGAGTGACGACGCGCCTGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCa                          >  1:2485390/3‑69 (MQ=255)
                                gagTAGAGTGACGACGCGCCTGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCa                          >  1:623790/3‑69 (MQ=255)
                                gagTAGAGTGACGACGCGCCTGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCa                          >  1:2569832/3‑69 (MQ=255)
                                gagTAGAGTGACGACGCGCCTGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCa                          >  1:2676352/3‑69 (MQ=255)
                                gagTAGAGTGACGACGCGCCTGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCa                          >  1:818742/3‑69 (MQ=255)
                                gagTAGAGTGACGACGCGCCTGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCa                          >  1:3171284/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:2543751/69‑1 (MQ=255)
                                          acgacgTGCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATcgcg                <  1:3189019/69‑1 (MQ=255)
                                          acgacgTGCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATcgcg                <  1:316537/69‑1 (MQ=255)
                                           cgacgTGCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATcgcg                <  1:3134861/68‑1 (MQ=255)
                                              cgTGCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCgcagc            <  1:2750741/69‑1 (MQ=255)
                                              cgTGCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCgcagc            <  1:3173894/69‑1 (MQ=255)
                                                 gCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGGCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGt         <  1:1043159/69‑1 (MQ=255)
                                                 gCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCTGTGTCACCACTGTTCAGCGCATCGCGCAGCAGt         <  1:1806855/69‑1 (MQ=255)
                                                 gCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGt         <  1:1564577/69‑1 (MQ=255)
                                                 gCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGt         <  1:1474263/69‑1 (MQ=255)
                                                 gCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGt         <  1:3180275/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:2832076/69‑1 (MQ=255)
                                                    tGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATCGCGCAGCAATGcc      <  1:1323286/69‑1 (MQ=255)
                                                    tGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATCGCGCAGCAATGcc      <  1:2670592/69‑1 (MQ=255)
                                                    tGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATCGCGCAGCAATGcc      <  1:3173104/69‑1 (MQ=255)
                                                    tGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATCGCGCAGCAATGcc      <  1:1544387/69‑1 (MQ=255)
                                                    tGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATCGCGCAGCAATGcc      <  1:2225330/69‑1 (MQ=255)
                                                    tGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATCGCGCAGCAATGcc      <  1:675024/69‑1 (MQ=255)
                                                    tGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATCGCGCAGCAATGcc      <  1:1873636/69‑1 (MQ=255)
                                                    tGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATCGCGCAGCAATGcc      <  1:1931304/69‑1 (MQ=255)
                                                    tGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATCGCGCAGCAATGcc      <  1:1992821/69‑1 (MQ=255)
                                                    tGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATCGCGCAGCAATGcc      <  1:2002631/69‑1 (MQ=255)
                                                    tGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATCGCGCAGCAATGcc      <  1:537219/69‑1 (MQ=255)
                                                    tGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATCGCGCAGCAATGcc      <  1:371147/69‑1 (MQ=255)
                                                    tGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATCGCGCAGCAATGcc      <  1:375238/69‑1 (MQ=255)
                                                    tGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATCGCGCAGCAATGcc      <  1:2080915/69‑1 (MQ=255)
                                                    tGAATCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATCGCGCAGCAATGcc      <  1:2254476/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:428748/1‑54 (MQ=255)
                                                       aTCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATcg                  >  1:75447/1‑54 (MQ=255)
                                                       aTCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATcg                  >  1:58626/1‑54 (MQ=255)
                                                       aTCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATcg                  >  1:2275473/1‑54 (MQ=255)
                                                       aTCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATcg                  >  1:1547434/1‑54 (MQ=255)
                                                       aTCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATcg                  >  1:1641733/1‑54 (MQ=255)
                                                       aTCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATcg                  >  1:1650004/1‑54 (MQ=255)
                                                       aTCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATcg                  >  1:1677807/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:2332204/1‑54 (MQ=255)
                                                       aTCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATcg                  >  1:1046215/1‑54 (MQ=255)
                                                       aTCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATcg                  >  1:3045893/1‑54 (MQ=255)
                                                       aTCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATcg                  >  1:2561458/1‑54 (MQ=255)
                                                       aTCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATcg                  >  1:2633215/1‑54 (MQ=255)
                                                       aTCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATcg                  >  1:2851264/1‑54 (MQ=255)
                                                       aTCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGTTCAACGCATcg                  >  1:2885997/1‑54 (MQ=255)
                                                       aTCGGAGCGCTCCACTCTGCCAGAGCTGCCGTATCACCACTGATCAACGCATcg                  >  1:2801065/1‑54 (MQ=255)
                                                        tCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGTAGTGCCAGtt  <  1:743517/69‑1 (MQ=255)
                                                        tCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGtt  <  1:1060171/69‑1 (MQ=255)
                                                        tCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGtt  <  1:1042581/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:2129274/69‑1 (MQ=255)
                                                        tCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGtt  <  1:2236404/69‑1 (MQ=255)
                                                        tCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGtt  <  1:365368/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)
                                                            |                                                                
CAGACGCTTTTCTACCGCCGGATTATGTACGCGGGTAGAGTGACGACGTGCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTT  >  minE/2384373‑2384497

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: