Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,015,631 G→A 15.6% Y146Y (TAC→TAT fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,015,6310GA15.6% 163.3 / 19.5 71Y146Y (TAC→TATfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base G (32/27);  new base A (6/5);  total (39/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

ATCGATCTGCGCTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGG  >  minE/1015562‑1015701
                                                                     |                                                                      
atcgatcTGCGCTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGt                                                                       <  1:2303055/71‑1 (MQ=255)
         cgcTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGa                                                               >  1:4642723/1‑70 (MQ=255)
         cgcTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGa                                                               <  1:3757661/70‑1 (MQ=255)
         cgcTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGAt                                                              >  1:4709590/1‑71 (MQ=255)
          gcTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGAtt                                                             <  1:4079060/71‑1 (MQ=255)
          gcTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGAtt                                                             <  1:3908554/71‑1 (MQ=255)
          gcTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGAtt                                                             <  1:465080/71‑1 (MQ=255)
              gCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCt                                                         <  1:1982053/71‑1 (MQ=255)
                aGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAg                                                                <  1:1692001/62‑1 (MQ=255)
                   ttGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGat                                                    <  1:2594335/71‑1 (MQ=255)
                     agTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGata                                                   <  1:3778340/69‑1 (MQ=255)
                      gTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGatat                                                  <  1:4143419/70‑1 (MQ=255)
                      gTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGatat                                                  <  1:3461270/70‑1 (MQ=255)
                      gTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGata                                                   >  1:5581302/1‑69 (MQ=255)
                      gTGCCGGTGGTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGatat                                                  <  1:3454570/70‑1 (MQ=255)
                      gTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGatat                                                  <  1:3564154/70‑1 (MQ=255)
                      gTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGatat                                                  <  1:1775278/70‑1 (MQ=255)
                           ggTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGata                                                   >  1:1171213/1‑64 (MQ=255)
                           ggTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGata                                                   >  1:534164/1‑64 (MQ=255)
                           ggTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGgtgt                                            <  1:4440195/71‑1 (MQ=255)
                            gTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTc                                                       <  1:2761376/59‑1 (MQ=255)
                            gTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGtt                                           <  1:4724312/71‑1 (MQ=255)
                              aTTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTa                                                            >  1:171411/1‑52 (MQ=255)
                               ttCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTtata                                        <  1:5054075/71‑1 (MQ=255)
                                tCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGtt                                           >  1:294279/1‑67 (MQ=255)
                                 cACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTc                                                       >  1:5493322/1‑54 (MQ=255)
                                 cACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGa                                      >  1:2754096/1‑71 (MQ=255)
                                 cACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGa                                      >  1:609669/1‑71 (MQ=255)
                                 cACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGa                                      >  1:1504845/1‑71 (MQ=255)
                                   cttcttTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCtt                                                        >  1:2849303/1‑51 (MQ=255)
                                   cttcttTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGAtata                                                 <  1:5037207/58‑1 (MQ=255)
                                      cttTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGcct                                 >  1:3645695/1‑70 (MQ=255)
                                      cttTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGa                                                               >  1:5372689/1‑41 (MQ=255)
                                      cttTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGa                                                               >  1:5617610/1‑41 (MQ=255)
                                      cttTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGa                                                               >  1:5365871/1‑41 (MQ=255)
                                      cttTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGa                                                               >  1:5193590/1‑41 (MQ=255)
                                      cttTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGa                                                               >  1:5226723/1‑41 (MQ=255)
                                      cttTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTAt                                                           >  1:1919862/1‑45 (MQ=255)
                                      cttTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTAt                                                           >  1:1387125/1‑45 (MQ=255)
                                      cttTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTAt                                                           >  1:3926092/1‑45 (MQ=255)
                                         tataCATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACcgc                                <  1:4172446/69‑1 (MQ=255)
                                         tataCATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACcgc                                <  1:1219466/69‑1 (MQ=255)
                                          gtaCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGata                                                   >  1:99921/2‑49 (MQ=255)
                                          gtaCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGata                                                   >  1:3952652/2‑49 (MQ=255)
                                           taCATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCgcgcgc                            >  1:4243451/1‑71 (MQ=255)
                                            aCATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGAtata                                                 <  1:301726/49‑1 (MQ=255)
                                            aCATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGtt                                           <  1:931394/55‑1 (MQ=255)
                                            aCATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGtt                                           <  1:5163592/55‑1 (MQ=255)
                                            aCATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGtt                                           <  1:3652130/55‑1 (MQ=255)
                                            aCATATCCAGCGGCGCGTTTTGCGACTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCgcgcgc                            >  1:3526384/1‑70 (MQ=255)
                                             cATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCa                          <  1:4582949/71‑1 (MQ=255)
                                                atCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTc                                                       >  1:3748015/1‑39 (MQ=255)
                                                atCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGatat                                                  <  1:1678795/44‑1 (MQ=255)
                                                atCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGatat                                                  <  1:4993292/44‑1 (MQ=255)
                                                atCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAgcg                       >  1:3766262/1‑71 (MQ=255)
                                                 tCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATaa                                                >  1:3536829/1‑45 (MQ=255)
                                                 tCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGaca                                    >  1:3149096/1‑57 (MQ=255)
                                                  ccAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGcgcc                     >  1:1071204/1‑71 (MQ=255)
                                                  ccAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGcgcc                     >  1:1919440/1‑71 (MQ=255)
                                                  ccAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGcgcc                     >  1:5688093/1‑71 (MQ=255)
                                                   cAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGcgccg                    >  1:4531729/1‑71 (MQ=255)
                                                   cAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGcgccg                    >  1:967218/1‑71 (MQ=255)
                                                     gcggcgCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGac                                     <  1:760103/52‑1 (MQ=255)
                                                      cggcgCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCtt                 <  1:2243750/71‑1 (MQ=255)
                                                       ggcgCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGacac                                   >  1:389997/1‑52 (MQ=255)
                                                              tttGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATcacca         >  1:327459/1‑71 (MQ=255)
                                                              tttGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATcacca         >  1:163289/1‑71 (MQ=255)
                                                               ttGCGAGTAGCGCAGATTATCTTCGATATAAGTGTta                                          >  1:4225252/1‑37 (MQ=255)
                                                               ttGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATcaccac        >  1:4547975/1‑71 (MQ=255)
                                                                 gCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGa                                      <  1:1262096/39‑1 (MQ=255)
                                                                    aGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAgcgc                      <  1:3941957/52‑1 (MQ=255)
                                                                     gTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCgg  >  1:1569488/1‑71 (MQ=255)
                                                                     |                                                                      
ATCGATCTGCGCTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGG  >  minE/1015562‑1015701

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: