Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,014,818 A→G 6.9% H417H (CAT→CAC fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,014,8180AG 6.9% 164.8 / 3.8 58H417H (CAT→CACfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base A (27/27);  new base G (2/2);  total (29/29)
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.38e-01

GTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATAT  >  minE/1014753‑1014888
                                                                 |                                                                      
gTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATc                                                                   >  1:1537813/1‑71 (MQ=255)
 ttGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCt                                                                  <  1:2487878/71‑1 (MQ=255)
  tGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCtt                                                                 >  1:2524295/1‑71 (MQ=255)
   gggCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCtt                                                                 <  1:4091952/70‑1 (MQ=255)
    ggCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTTATCttt                                                                <  1:8670/70‑1 (MQ=255)
    ggCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCttt                                                                <  1:498954/70‑1 (MQ=255)
     gCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGa                                                              <  1:1929100/71‑1 (MQ=255)
       cAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATg                                                            <  1:1344157/71‑1 (MQ=255)
          gagaACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGga                                                                        >  1:944244/1‑56 (MQ=255)
            gaACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTg                                                       >  1:833732/1‑71 (MQ=255)
             aaCCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCttt                                                                >  1:551054/1‑61 (MQ=255)
             aaCCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCttt                                                                >  1:4093745/1‑61 (MQ=255)
             aaCCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGa                                                              >  1:2627726/1‑63 (MQ=255)
                     gCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGa                                                              <  1:1054043/55‑1 (MQ=255)
                     gCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGa                                                              <  1:270241/55‑1 (MQ=255)
                       aTAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGa                                                              >  1:2907868/1‑53 (MQ=255)
                        tAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCt                                                                  <  1:2888106/48‑1 (MQ=255)
                        tAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCtt                                            >  1:423521/1‑70 (MQ=255)
                         aaCCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCa                                          >  1:3085683/1‑71 (MQ=255)
                          aCCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCt                                                                  <  1:504573/46‑1 (MQ=255)
                          aCCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCAc                                         >  1:2407409/1‑71 (MQ=255)
                          aCCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCAc                                         >  1:916124/1‑71 (MQ=255)
                            cTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCttt                                                                >  1:3807066/1‑46 (MQ=255)
                             ttCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCt                                                                  <  1:2291594/43‑1 (MQ=255)
                             ttCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCt                                                                  <  1:2107492/43‑1 (MQ=255)
                             ttCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATg                                                            <  1:1870067/49‑1 (MQ=255)
                               ccGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTa                                    >  1:1097324/1‑71 (MQ=255)
                                cGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTat                                   <  1:3468916/71‑1 (MQ=255)
                                cGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTa                                    <  1:556199/70‑1 (MQ=255)
                                  gCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCtt                                                                 <  1:2954730/39‑1 (MQ=255)
                                       ttGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGa                                                              >  1:2898046/1‑37 (MQ=255)
                                       ttGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGa                                                              >  1:1232856/1‑37 (MQ=255)
                                       ttGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCg                            >  1:2590988/1‑71 (MQ=255)
                                       ttGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCg                            >  1:4077071/1‑71 (MQ=255)
                                       ttGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCg                            >  1:2827994/1‑71 (MQ=255)
                                       ttGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCg                            >  1:3619533/1‑71 (MQ=255)
                                        tGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGa                                                              >  1:2293454/1‑36 (MQ=255)
                                        tGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGa                                                              >  1:3995564/1‑36 (MQ=255)
                                        tGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGt                                      <  1:1855589/60‑1 (MQ=255)
                                        tGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGt                                      <  1:3240787/60‑1 (MQ=255)
                                            cGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCtt                       >  1:607638/1‑71 (MQ=255)
                                            cGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCtt                       >  1:2561392/1‑71 (MQ=255)
                                             ggACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCttt                      <  1:1227072/71‑1 (MQ=255)
                                               aCCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGt                                      <  1:3873275/53‑1 (MQ=255)
                                                   gcgTAGTAGATCGGGTGATCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGctt                >  1:189733/1‑71 (MQ=255)
                                                       agtaAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCttt                      <  1:2832601/61‑1 (MQ=255)
                                                         taAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTatcc                                 <  1:4058249/48‑1 (MQ=255)
                                                           gaTCGGGTGATCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTgtg        >  1:4086970/2‑69 (MQ=255)
                                                           aaTCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGtt                 <  1:1008676/62‑1 (MQ=255)
                                                           aaTCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAg        >  1:1604695/1‑71 (MQ=255)
                                                           aaTCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAg        >  1:609762/1‑71 (MQ=255)
                                                           aaTCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAg        >  1:3266590/1‑71 (MQ=255)
                                                            aTCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGc       <  1:518570/71‑1 (MQ=255)
                                                             tCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGa         <  1:2043162/68‑1 (MQ=255)
                                                               ggATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCtt                       <  1:1294073/52‑1 (MQ=255)
                                                               ggATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGttt                <  1:760663/59‑1 (MQ=255)
                                                               ggATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCaa     >  1:3321036/1‑70 (MQ=255)
                                                               ggATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAat    <  1:137192/71‑1 (MQ=255)
                                                               ggATGATCTTTGATGTACTGCGGCAGCCCTTCACAGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAat    <  1:2608323/71‑1 (MQ=255)
                                                                gatgatCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCtt                       <  1:992132/51‑1 (MQ=255)
                                                                gatgatCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAata   <  1:52035/71‑1 (MQ=255)
                                                                 atgatCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAatat  <  1:3800675/71‑1 (MQ=255)
                                                                 |                                                                      
GTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATAT  >  minE/1014753‑1014888

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