Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,015,742 G→A 12.1% G109G (GGC→GGT fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,015,7420GA12.1% 193.3 / 13.9 66G109G (GGC→GGTfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base G (31/27);  new base A (4/4);  total (35/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGC  >  minE/1015672‑1015805
                                                                      |                                                               
gcgcCAGCGCCGCTTCATAACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTg                                                                 <  1:3858535/71‑1 (MQ=255)
  gcCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGcc                                                               <  1:1555891/71‑1 (MQ=255)
  gcCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGcc                                                               <  1:1550816/71‑1 (MQ=255)
    cAGCGCCGCTTCATCTCCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCg                                                              <  1:4120867/70‑1 (MQ=255)
    cAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGTGTTACCgg                                                             >  1:3997553/1‑71 (MQ=255)
    cAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCgg                                                             >  1:3433433/1‑71 (MQ=255)
      gcgcCGCTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTa                                                           >  1:440759/1‑71 (MQ=255)
        gccgcTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATc                                                         >  1:1570852/1‑71 (MQ=255)
         ccgcTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATcc                                                        <  1:1267889/71‑1 (MQ=255)
           gcTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCt                                                       >  1:1993825/1‑70 (MQ=255)
             ttCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGa                                                     >  1:1292530/1‑70 (MQ=255)
              tcatcaCCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGa                                                     >  1:1416971/1‑69 (MQ=255)
              tcatcaCCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGa                                                     >  1:478473/1‑69 (MQ=255)
               catcaCCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGAc                                                    <  1:3461133/69‑1 (MQ=255)
                atcaCCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGcc                                                               >  1:904290/1‑57 (MQ=255)
                atcaCCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAg                                                  <  1:453733/70‑1 (MQ=255)
                 tcaCCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGt                                                <  1:2080477/71‑1 (MQ=255)
                       ccaccGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTgg                                             >  1:2748941/1‑68 (MQ=255)
                       ccaccGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAg                                          <  1:187228/71‑1 (MQ=255)
                       ccaccGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCCGACAGGTTGGCAg                                          <  1:916815/71‑1 (MQ=255)
                        caccGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTgg                                             <  1:2739178/67‑1 (MQ=255)
                         accGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAg                                                  <  1:288501/61‑1 (MQ=255)
                         accGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGaa                                        >  1:3190041/1‑71 (MQ=255)
                          ccGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGaa                                        >  1:883492/1‑70 (MQ=255)
                          ccGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGaa                                        >  1:3891143/1‑70 (MQ=255)
                           cGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAAc                                       <  1:3544685/70‑1 (MQ=255)
                            ggTCCATACACGCTGCCCTTTTTTACTAACAATAATCGCGGTGCCTGTATCCTGACAGGTTGGCAGAACGc                                     <  1:888784/71‑1 (MQ=255)
                            ggTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGc                                     >  1:3795709/1‑71 (MQ=255)
                            ggTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGc                                     >  1:1838300/1‑71 (MQ=255)
                                cATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCtt                                  <  1:434743/70‑1 (MQ=255)
                                  tACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGAc                                                    >  1:1635468/1‑50 (MQ=255)
                                    cacGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTcgcc                             >  1:1674049/1‑71 (MQ=255)
                                    cacGCTGCCCTTTTTTACCAACAATAATCGCGGTACCGGTa                                                           >  1:938722/1‑41 (MQ=255)
                                     acGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCg                                                              >  1:736277/1‑37 (MQ=255)
                                     acGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCg                                                              >  1:2487963/1‑37 (MQ=255)
                                     acGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCg                                                              >  1:852895/1‑37 (MQ=255)
                                      cGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTCGCAGAACGCCTTTcgccgc                           <  1:314583/71‑1 (MQ=255)
                                       gCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACg                                      >  1:2120855/1‑59 (MQ=255)
                                         tGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGa                                                     >  1:2795689/1‑42 (MQ=255)
                                           cccTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGAt                        <  1:4119048/69‑1 (MQ=255)
                                            ccTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAg                                          <  1:623639/50‑1 (MQ=255)
                                            ccTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTc                     <  1:302048/71‑1 (MQ=255)
                                             cTTTTTTACCAACAATAATCGCGGTACCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCg                    <  1:1655130/71‑1 (MQ=255)
                                               tttttACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCg                    >  1:1171732/1‑69 (MQ=255)
                                               tttttACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCg                    >  1:2380246/1‑69 (MQ=255)
                                                ttttACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCgg                   <  1:1737907/69‑1 (MQ=255)
                                                ttttACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAg                 >  1:1585987/1‑71 (MQ=255)
                                                 tttACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCgg                   <  1:3447688/68‑1 (MQ=255)
                                                 tttACCAACAATAATCGCGGTACCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGa                         >  1:3131853/1‑62 (MQ=255)
                                                 tttACCAACAATAATCGCGGTACCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGa                         >  1:1581633/1‑62 (MQ=255)
                                                    aCCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTAc             >  1:916538/1‑71 (MQ=255)
                                                     ccAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAg                                          <  1:595295/41‑1 (MQ=255)
                                                        acaaTAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAg         >  1:2982685/1‑71 (MQ=255)
                                                        acaaTAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAg         >  1:3062939/1‑71 (MQ=255)
                                                        acaaTAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAg         >  1:1587391/1‑71 (MQ=255)
                                                         caaTAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAgg        <  1:500676/71‑1 (MQ=255)
                                                         caaTAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAg         >  1:3006653/1‑70 (MQ=255)
                                                          aataatCGCTGTGCCGGTATCCTGACAGGTCGGCAGAACGCCTTTcgccgc                           <  1:3438468/51‑1 (MQ=39)
                                                          aataatCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAgg        >  1:2158346/1‑70 (MQ=255)
                                                            taatCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAAt     <  1:3171930/71‑1 (MQ=255)
                                                            taatCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAAt     <  1:1609567/71‑1 (MQ=255)
                                                              atCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTcgccgc                           <  1:1192979/47‑1 (MQ=255)
                                                              atCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTg   >  1:263159/1‑71 (MQ=255)
                                                              atCGCGGTACCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAAtt    <  1:208523/70‑1 (MQ=255)
                                                              atCGCGGTACCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTg   <  1:2743737/71‑1 (MQ=255)
                                                               tCGCGGTACCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGc  <  1:1907472/71‑1 (MQ=255)
                                                                     tGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGa                  <  1:142306/49‑1 (MQ=255)
                                                                      |                                                               
GCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGC  >  minE/1015672‑1015805

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: