Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 555,416 T→G 30.9% M211L (ATG→CTG)  moeB ← molybdopterin synthase sulfurylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE555,4160TG30.9% 119.7 / 87.6 123M211L (ATG→CTG) moeBmolybdopterin synthase sulfurylase
Reads supporting (aligned to +/- strand):  ref base T (26/59);  new base G (15/23);  total (41/82)
Fisher's exact test for biased strand distribution p-value = 4.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.34e-01

Only 100 of 123 total aligned reads displayed.

TCATCGCATCGTACATGACGATTTTCCCGCTGGCCGGTTTGCCATAACCTGCCAGCATTTTGATCGCTTCCATTGCCTGCAACGAACCAATTACGCCGATCAACGGTGCCATTACGCCTGCTTC  >  minE/555359‑555482
                                                         |                                                                  
tCATCGCATCGTACATGACGATTTTCCCGCTGGCCGGTTTGCCATAACCTGCCAGCATTTTGATCGCt                                                          >  1:70715/1‑68 (MQ=255)
tCATCGCATCGTACATGACGATTTTCCCGCTGGCCGGTTTGCCATAACCTGCCAGCATTTTGATCGCt                                                          >  1:86332/1‑68 (MQ=255)
tCATCGCATCGTACATGACGATTTTCCCGCTGGCCGGTTTGCCATAACCTGCCAGCATTTTGATCGCt                                                          >  1:781236/1‑68 (MQ=255)
tCATCGCATCGTACATGACGATTTTCCCGCTGGCCGGTTTGCCATAACCTGCCAGCATTTTGATCGCt                                                          >  1:787306/1‑68 (MQ=255)
 catcgcatcgTACATGACGATTTTCCCGCTGGCCGGTTTGCCATAACCTGCCAGCATTTTGATCGCTTc                                                        >  1:3112221/1‑69 (MQ=255)
 catcgcatcgTACATGACGATTTTCCCGCTGGCCGGTTTGCCATAACCTGCCAGCATTTTGATCGCTTc                                                        <  1:698586/69‑1 (MQ=255)
 catcgcatcgTACATGACGATTTTCCCGCTGGCCGGTTTGCCATAACCTGCCAGCATTTTGATCGCTTc                                                        >  1:1952244/1‑69 (MQ=255)
 catcgcatcgTACATGACGATTTTCCCGCTGGCCGGTTTGCCATAACCTGCCAGCATTTTGATCGCTTc                                                        >  1:1948148/1‑69 (MQ=255)
 catcgcatcgTACATGACGATTTTCCCGCTGGCCGGTTTGCCATAACCTGCCAGCATTTTGATCGCTTc                                                        >  1:545736/1‑69 (MQ=255)
 catcgcatcgTACATGACGATTTTCCCGCTGGCCGGTTTGCCATAACCTGCCAGCATTTTGATCGCTTc                                                        <  1:1863835/69‑1 (MQ=255)
 catcgcatcgTACATGACGATTTTCCCGCTGGCCGGTTTGCCATAACCTGCCAGCATTTTGATCGCTTc                                                        <  1:2490011/69‑1 (MQ=255)
 catcgcatcgTACATGACGATTTTCCCGCTGGCCGGTTTGCCATAACCTGCCAGCATTTTGATCGCTTc                                                        >  1:168957/1‑69 (MQ=255)
 catcgcatcgTACATGACGATTTTCCCGCTGGCCGGTTTGCCATAACCTGCCAGCATTTTGATCGCTTc                                                        >  1:793179/1‑69 (MQ=255)
 catcgcatcgTACATGACGATTTTCCCGCTGGCCGGTTTGCCATAACCTGCCAGCATTTTGATCGCTTc                                                        <  1:3163338/69‑1 (MQ=255)
 catcgcatcgTACATGACGATTTTCCCGCTGGCCGGTTTGCCATAACCTGCCAGCATTTTGATCGCTTc                                                        <  1:3002090/69‑1 (MQ=255)
 catcgcatcgTACATGACGATTTTCCCGCTGGCCGGTTTGCCATAACCTGCCAGCATTTTGATCGCTTc                                                        >  1:3404106/1‑69 (MQ=255)
 catcgcatcgTACATGACGATTTTCCCGCTGGCCGGTTTGCCATAACCTGCCAGCATTTTGATCGCTTc                                                        <  1:1356153/69‑1 (MQ=255)
 catcgcatcgTACATGACGATTTTCCCGCTGGCCGGTTTGCCATAACCTGCCAGCATTTTGATCGCTTc                                                        <  1:2899440/69‑1 (MQ=255)
    cgcatcgTACATGACGATTTTCCCGCTGGCCGGTTTGCCATAACCTGCCAGCATTTTGATCGCTTCCAt                                                     >  1:425577/1‑69 (MQ=255)
    cgcatcgTACATGACGATTTTCCCGCTGGCCGGTTTGCCATAACCTGCCAGCATTTTGATCGCTTCCAt                                                     >  1:2147294/1‑69 (MQ=255)
    cgcatcgTACATGACGATTTTCCCGCTGGCCGGTTTGCCATAACCTGCCAGCATTTTGATCGCTTCCAt                                                     >  1:1201135/1‑69 (MQ=255)
    cgcatcgTACATGACGATTTTCCCGCTGGCCGGTTTGCCATAACCTGCCAGCATTTTGATCGCTTCCAt                                                     >  1:1294622/1‑69 (MQ=255)
    cgcatcgTACATGACGATTTTCCCGCTGGCCGGTTTGCCATAACCTGCCAGCATTTTGATCGCTTCCAt                                                     >  1:2331384/1‑69 (MQ=255)
    cgcatcgTACATGACGATTTTCCCGCTGGCCGGTTTGCCATAACCTGCCAGCATTTTGATCGCTTCCAt                                                     >  1:1547680/1‑69 (MQ=255)
    cgcatcgTACATGACGATTTTCCCGCTGGCCGGTTTGCCATAACCTGCCAGCATTTTGATCGCTTCCAt                                                     >  1:1620866/1‑69 (MQ=255)
    cgcatcgTACATGACGATTTTCCCGCTGGCCGGTTTGCCATAACCTGCCAGCATTTTGATCGCTTCCAt                                                     >  1:1595102/1‑69 (MQ=255)
     gcatcgTACATGACGATTTTCCCGCTGGCCGGTTTGCCATAACCTGCCAGCATTTTGATCGCTTCCAtt                                                    >  1:2611240/1‑69 (MQ=255)
     gcatcgTACATGACGATTTTCCCGCTGGCCGGTTTGCCATAACCTGCCAGCATTTTGATCGCTTCCAtt                                                    >  1:19332/1‑69 (MQ=255)
     gcatcgTACATGACGATTTTCCCGCTGGCCGGTTTGCCATAACCTGCCAGCATTTTGATCGCTTCCAtt                                                    >  1:76987/1‑69 (MQ=255)
        tcgTACATGACGATTTTCCCGCTGGCCGGTTTGCCATAACCTGCCAGCATTTTGATCGCTTc                                                        <  1:1604117/62‑1 (MQ=255)
            aCATGACGATTTTCCCGCTGGCCGGTTTGCCATAACCTGCCAGCATTTTGATCGCTTCCATTGCCTGc                                              <  1:2004936/68‑1 (MQ=255)
                 aCGATTTTCCCGCTGGCCGGTTTGCCATAACCTGCCAGCATTTTg                                                                <  1:799511/45‑1 (MQ=255)
                 aCGATTTTCCCGCTGGCCGGTTTGCCATAACCTGCCAGCATTTTg                                                                <  1:1634442/45‑1 (MQ=255)
                 aCGATTTTCCCGCTGGCCGGTTTGCCATAACCTGCCAGCATTTTg                                                                <  1:747899/45‑1 (MQ=255)
                   gATTTTCCCGCTGGCCGGTTTGCCATAACCTGCCAGCATTTTGATCGCTTCCATTGCCTGCAACGAAcc                                      <  1:3192774/69‑1 (MQ=255)
                   gATTTTCCCGCTGGCCGGTTTGCCATAACCTGCCAGCATTTTGATCGCTTCCATTGCCTGCAACGAAcc                                      <  1:2929492/69‑1 (MQ=255)
                   gATTTTCCCGCTGGCCGGTTTGCCATAACCTGCCAGCATTTTGATCGCTTCCATTGCCTGCAACGAAcc                                      <  1:2769596/69‑1 (MQ=255)
                   gATTTTCCCGCTGGCCGGTTTGCCATAACCTGCCAGCATTTTGATCGCTTCCATTGCCTGCAACGAAcc                                      <  1:1729392/69‑1 (MQ=255)
                   gATTTTCCCGCTGGCCGGTTTGCCATAACCTGCCAGCATTTTGATCGCTTCCATTGCCTGCAACGAAcc                                      <  1:2200326/69‑1 (MQ=255)
                   gATTTTCCCGCTGGCCGGTTTGCCATAACCTGCCAGCATTTTGATCGCTTCCATTGCCTGCAACGAAcc                                      <  1:1968509/69‑1 (MQ=255)
                   gATTTTCCCGCTGGCCGGTTTGCCATAACCTGCCAGCATTATGATCGCTTCCATTGCCTGCAACGAAca                                      <  1:1805365/69‑2 (MQ=255)
                     ttttTCCGCTGGCCGGTTTGCCATAACCTGCCAGCAGTTTGATCGCTTCCATTGCCTGCAACGAACCa                                     >  1:2264259/1‑68 (MQ=255)
                     ttttTCCGCTGGCCGGTTTGCCATAACCTGCCAGCAGTTTGATCGCTTCCATTGCCTGCAACGAACCa                                     >  1:1875622/1‑68 (MQ=255)
                     ttttTCCGCTGGCCGGTTTGCCATAACCTGCCAGCAGTTTGATCGCTTCCATTGCCTGCAACGAACCa                                     >  1:111352/1‑68 (MQ=255)
                     ttttTCCGCTGGCCGGTTTGCCATAACCTGCCAGCAGTTTGATCGCTTCCATTGCCTGCAACGAACCa                                     >  1:2532871/1‑68 (MQ=255)
                     ttttTCCGCTGGCCGGTTTGCCATAACCTGCCAGCAGTTTGATCGCTTCCATTGCCTGCAACGAACCa                                     >  1:399132/1‑68 (MQ=255)
                     ttttTCCGCTGGCCGGTTTGCCATAACCTGCCAGCAGTTTGATCGCTTCCATTGCCTGCAACGAACCa                                     >  1:3432234/1‑68 (MQ=255)
                     ttttTCCGCTGGCCGGTTTGCCATAACCTGCCAGCAGTTTGATCGCTTCCATTGCCTGCAACGAACCa                                     >  1:897901/1‑68 (MQ=255)
                     ttttTCCGCTGGCCGGTTTGCCATAACCTGCCAGCAGTTTGATCGCTTCCATTGCCTGCAACGAACCa                                     >  1:3289827/1‑68 (MQ=255)
                     ttttTCCGCTGGCCGGTTTGCCATAACCTGCCAGCAGTTTGATCGCTTCCATTGCCTGCAACGAACCa                                     >  1:104171/1‑68 (MQ=255)
                     ttttTCCGCTGGCCGGTTTGCCATAACCTGCCAGCAGTTTGATCGCTTCCATTGCCTGCAACGAACCa                                     >  1:2877454/1‑68 (MQ=255)
                     ttttTCCGCTGGCCGGTTTGCCATAACCTGCCAGCAGTTTGATCGCTTCCATTGCCTGCAACGAACCa                                     >  1:844663/1‑68 (MQ=255)
                     ttttTCCGCTGGCCGGTTTGACATAACCTGCCAGCAGTTTGATCGCTTCCATTGCCTGCAACGAACCa                                     >  1:501749/1‑68 (MQ=255)
                     ttttTCCGCTGGCCGGTTTCCCATAACCTGCCAGCAGTTTGATCGCTTCCATTGCCTGCAACGAACCa                                     >  1:2139088/1‑68 (MQ=255)
                      tttCCCGCTGGCCGGTTTGCCATAACCTGCCAGCATTTTGATCGCTTCCATTGCCTGCAACGAACCAAt                                   <  1:2870954/69‑1 (MQ=255)
                      tttCCCGCTGGCCGGTTTGCCATAACCTGCCAGCATTTTGATCGCTTCCATTGCCTGCAACGAACCAAt                                   <  1:3172683/69‑1 (MQ=255)
                      tttCCCGCTGGCCGGTTTGCCATAACCTGCCAGCATTTTGATCGCTTCCATTGCCTGCAACGAACCAAt                                   <  1:987202/69‑1 (MQ=255)
                      tttCCCGCTGGCCGGTTTGCCATAACCTGCCAGCATTTTGATCGCTTCCATTGCCTGCAACGAACCAAt                                   <  1:3258464/69‑1 (MQ=255)
                      tttCCCGCTGGCCGGTTTGCCATAACCTGCCAGCATTTTGATCGCTTCCATTGCCTGCAACGAACCAAt                                   <  1:588616/69‑1 (MQ=255)
                      tttCCCGCTGGCCGGTTTGCCATAACCTGCCAGCATTTTGATCGCTTCCATTGCCTGCAACGAACCAAt                                   <  1:653174/69‑1 (MQ=255)
                      tttCCCGCTGGCCGGTTTGCCATAACCTGCCAGCATTTTGATCGCTTCCATTGCCTGCAACGAACCAAt                                   <  1:2845527/69‑1 (MQ=255)
                         cccGCTGGCCGGTTTGCCATAACCTGCCAGCATTTTGATCGCTACCATTGCCTGCAACGAACCAAt                                   <  1:1239965/66‑1 (MQ=255)
                                       tGCCATAACCTGCCAGCAGTTTGATCGCTTCCATTGCCTGCAACGAAcc                                      <  1:3159540/49‑1 (MQ=255)
                                       tGCCATAACCTGCCAGCAGTTTGATCGCTTCCATTGCCTGCAACGAAcc                                      <  1:903645/49‑1 (MQ=255)
                                       tGCCATAACCTGCCAGCAGTTTGATCGCTTCCATTGCCTGCAACGAAcc                                      <  1:1185151/49‑1 (MQ=255)
                                       tGCCATAACCTGCCAGCAGTTTGATCGCTTCCATTGCCTGCAACGAAcc                                      <  1:853574/49‑1 (MQ=255)
                                       tGCCATAACCTGCCAGCAGTTTGATCGCTTCCATTGCCTGCAACGAAcc                                      <  1:802676/49‑1 (MQ=255)
                                       tGCCATAACCTGCCAGCAGTTTGATCGCTTCCATTGCCTGCAACGAAcc                                      <  1:1463101/49‑1 (MQ=255)
                                       tGCCATAACCTGCCAGCAGTTTGATCGCTTCCATTGCCTGCAACGAAcc                                      <  1:1549776/49‑1 (MQ=255)
                                       tGCCATAACCTGCCAGCAGTTTGATCGCTTCCATTGCCTGCAACGAAcc                                      <  1:2046455/49‑1 (MQ=255)
                                       tGCCATAACCTGCCAGCAGTTTGATCGCTTCCATTGCCTGCAACGAAcc                                      <  1:2124257/49‑1 (MQ=255)
                                       tGCCATAACCTGCCAGCAGTTTGATCGCTTCCATTGCCTGCAACGAAcc                                      <  1:2207309/49‑1 (MQ=255)
                                       tGCCATAACCTGCCAGCAGTTTGATCGCTTCCATTGCCTGCAACGAAcc                                      <  1:2225620/49‑1 (MQ=255)
                                       tGCCATAACCTGCCAGCAGTTTGATCGCTTCCATTGCCTGCAACGAAcc                                      <  1:2686898/49‑1 (MQ=255)
                                       tGCCATAACCTGCCAGCAGTTTGATCGCTTCCATTGCCTGCAACGAAcc                                      <  1:2813956/49‑1 (MQ=255)
                                       tGCCATAACCTGCCAGCAGTTTGATCGCTTCCATTGCCTGCAACGAAcc                                      <  1:2819841/49‑1 (MQ=255)
                                       tGCCATAACCTGCCAGCAGTTTGATCGCTTCCATTGCCTGCAACGAAcc                                      <  1:336964/49‑1 (MQ=255)
                                       tGCCATAACCTGCCAGCAGTTTGATCGCTTCCATTGCCTGCAACGAAcc                                      <  1:3119993/49‑1 (MQ=255)
                                       tGCCATAACCTGCCAGCAGTTTGATCGCTTCCATTGCCTGCAACGAAcc                                      <  1:3159575/49‑1 (MQ=255)
                                       tGCCATAACCTGCCAGCAGTTTGATCGCTTCCATTGCCTGCAACGAAcc                                      <  1:3001369/49‑1 (MQ=255)
                                       tGCCATAACCTGCCAGCAGTTTGATCGCTTCCATTGCCTGCAACGAAcc                                      <  1:3089570/49‑1 (MQ=255)
                                       tGCCATAACCTGCCAGCAGTTTGATCGCTTCCATTGCCTGCAACGAAcc                                      <  1:300843/49‑1 (MQ=255)
                                           aTAACCTGCCAGCATTTTGATCGCTTCCATTGCCTGCAACGAACCAATTACGCCGATCAACGGTGCCAt              <  1:2442364/69‑1 (MQ=255)
                                           aTAACCTGCCAGCATTTTGATCGCTTCCATTGCCTGCAACGAACCAATTACGCCGATCAACGGTGCCAt              <  1:1916781/69‑1 (MQ=255)
                                           aTAACCTGCCAGCATTTTGATCGCTTCCATTGCCTGCAACGAACCAATTACGCCGATCAACGGTGCCAt              <  1:1886636/69‑1 (MQ=255)
                                           aTAACCTGCCAGCATTTTGATCGCTTCCATTGCCTGCAACGAACCAATTACGCCGATCAACGGTGCCAt              <  1:1499836/69‑1 (MQ=255)
                                           aTAACCTGCCAGCATTTTGATCGCTTCCATTCCTGCAACGAACCAATTACGCCGATCAACGGTGCCAt              <  1:865098/68‑1 (MQ=255)
                                            tAACCTGCCAGCATTTTGATCGCTTCCATTGCCTGCAACGAACCAATTACGCCGATCAACGGTGCCAt              <  1:1079631/68‑1 (MQ=255)
                                              aCCTGCCAGCATTTTGATCGCTTCCATTGCCTGCAACGAAcc                                      <  1:1577023/42‑1 (MQ=255)
                                              aCCTGCCAGCATTTTGATCGCTTCCATTGCCTGCAACGAACCAATTACGCCGATCAACGGTGCCAt              <  1:1605620/66‑1 (MQ=255)
                                              aCCTGCCAGCATTTTGATCGCTTCCATTGCCTGCAACGAACCAATTACGCCGATCAACGGTGCCAt              <  1:1333874/66‑1 (MQ=255)
                                              aCCTGCCAGCATTTTGATCGCTTCCATTGCCTGCAACGAACCAATTACGCCGATCAACGGTGCCAt              <  1:715420/66‑1 (MQ=255)
                                              aCCTGCCAGCATTTTGATCGCTTCCATTGCCTGCAACGAACCAATTACGCCGATCAACGGTGCCAt              <  1:3363587/66‑1 (MQ=255)
                                              aCCTGCCAGCATTTTGATCGCTTCCATTGCCTGCAACGAACCAATTACGCCGATCAACGGTGCCAt              <  1:639127/66‑1 (MQ=255)
                                              aCCTGCCAGCATTTTGATCGCTTCCATTGCCTGCAACGAACCAATTACGCCGATCAACGGTGCCAt              <  1:911460/66‑1 (MQ=255)
                                              aCCTGCCAGCATTTTGATCGCTTCCATTGCCTGCAACGAACCAATTACGCCGATCAACGGTGCCAt              <  1:2274988/66‑1 (MQ=255)
                                                       cATTTTGATCGCTTCCATTGCCTGCAACGAACCAATTACGCCGATCAACGGTGCCATTACGCCTGCTTc  <  1:3190920/69‑1 (MQ=255)
                                                       cATTTTGATCGCTTCCATTGCCTGCAACGAACCAATTACGCCGATCAACGGTGCCATTACGCCTGCTTc  <  1:3060611/69‑1 (MQ=255)
                                                       cATTTTGATCGCTTCCATTGCCTGCAACGAACCAATTACGCCGATCAACGGTGCCATTACGCCTGCTTc  <  1:3171108/69‑1 (MQ=255)
                                                       cATTTTGATCGCTTCCATTGCCTGCAACGAACCAATTACGCCGATCAACGGTGCCATTACGCCTGCTTc  <  1:674134/69‑1 (MQ=255)
                                                         |                                                                  
TCATCGCATCGTACATGACGATTTTCCCGCTGGCCGGTTTGCCATAACCTGCCAGCATTTTGATCGCTTCCATTGCCTGCAACGAACCAATTACGCCGATCAACGGTGCCATTACGCCTGCTTC  >  minE/555359‑555482

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: