Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,014,812 A→G 6.1% I419I (ATT→ATC fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,014,8120AG 6.1% 189.1 / 3.5 65I419I (ATT→ATCfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base A (26/35);  new base G (2/2);  total (28/37)
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.51e-01

CGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGT  >  minE/1014745‑1014879
                                                                   |                                                                   
cggcggTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGGTTTGGCCGGACCCGCGTAGTAAATc                                                                  >  1:721406/1‑71 (MQ=255)
         ttGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTaa                                                                     <  1:1460445/59‑1 (MQ=255)
                aaGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACGCGCGTAGTAGATCGGGTGATCTTTGATGt                                                  <  1:2619708/71‑1 (MQ=255)
                aaGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGATg                                                   >  1:101962/1‑70 (MQ=255)
                aaGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGATg                                                   >  1:3567164/1‑70 (MQ=255)
                aaGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGt                                                  <  1:113606/71‑1 (MQ=255)
                aaGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGt                                                  <  1:1732437/71‑1 (MQ=255)
                 agagAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGt                                                  <  1:3653243/70‑1 (MQ=255)
                    gaACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTg                                              >  1:1948386/1‑71 (MQ=255)
                    gaACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTg                                              <  1:611043/71‑1 (MQ=255)
                     aaCCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGc                                             >  1:2491710/1‑71 (MQ=255)
                       ccGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCg                                                                 <  1:3653738/49‑1 (MQ=255)
                       ccGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGt                                                  <  1:736646/64‑1 (MQ=255)
                          gAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTa                                                 >  1:584844/1‑62 (MQ=255)
                          gAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCgg                                           >  1:1275875/1‑68 (MQ=255)
                          gAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCgg                                           >  1:2820544/1‑68 (MQ=255)
                            ggCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGcc                                      >  1:1218710/1‑71 (MQ=255)
                             gCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCttt                                                       <  1:1371049/53‑1 (MQ=255)
                             gCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCttt                                                       <  1:1933503/53‑1 (MQ=255)
                             gCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGccc                                     >  1:323446/1‑71 (MQ=255)
                               aTAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCtt                                   >  1:1414733/1‑71 (MQ=255)
                                tAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATg                                                   <  1:812547/54‑1 (MQ=255)
                                 aaCCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCt                                    >  1:551355/1‑68 (MQ=255)
                                  aCCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGgatgat                                                           >  1:329019/1‑44 (MQ=255)
                                  aCCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTa                                                 >  1:993509/1‑54 (MQ=255)
                                  aCCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCtt                                   >  1:1103119/1‑68 (MQ=255)
                                     ttCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCtt                                                        <  1:2461649/44‑1 (MQ=255)
                                       ccGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGtt                            <  1:2773770/70‑1 (MQ=255)
                                       ccGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGtt                            <  1:2375409/70‑1 (MQ=255)
                                          gCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGc                                          <  1:776656/53‑1 (MQ=255)
                                          gCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGc                                          <  1:1796592/53‑1 (MQ=255)
                                          gCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTat                          >  1:2638307/1‑69 (MQ=255)
                                           cGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTg                                                      <  1:1316758/40‑1 (MQ=255)
                                           cGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTg                                                      <  1:1724404/40‑1 (MQ=255)
                                           cGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTg                                                      <  1:541624/40‑1 (MQ=255)
                                           cGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTg                                                      <  1:2583768/40‑1 (MQ=255)
                                               ttGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGa                                                     >  1:367144/1‑37 (MQ=255)
                                               ttGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCg                   >  1:2329281/1‑71 (MQ=255)
                                                tGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCg                   <  1:3097906/70‑1 (MQ=255)
                                                tGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCg                   <  1:2192979/70‑1 (MQ=255)
                                                  gCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGctc                <  1:1342874/71‑1 (MQ=255)
                                                  gCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGctc                <  1:1733072/71‑1 (MQ=255)
                                                    cGGACCCGCGTAGTAGATCGGGTGATCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCtt              <  1:2672719/71‑1 (MQ=255)
                                                    cGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTTCAGCAGCCCTTCACCGTTATCCATCCGCTCtt              >  1:2900480/1‑71 (MQ=255)
                                                     ggACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCttt             >  1:139355/1‑71 (MQ=255)
                                                     ggACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCttt             >  1:1813195/1‑71 (MQ=255)
                                                     ggACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCttt             >  1:2808968/1‑71 (MQ=255)
                                                     ggACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCttt             >  1:1026715/1‑71 (MQ=255)
                                                     ggACCCGCGTAGTAAATAGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCttt             >  1:1368526/1‑71 (MQ=255)
                                                      gACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCa                                         >  1:2155113/1‑42 (MQ=255)
                                                        cccGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCa           >  1:1869584/1‑70 (MQ=255)
                                                         ccGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGt         <  1:2727957/71‑1 (MQ=255)
                                                         ccGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGt         <  1:130702/71‑1 (MQ=255)
                                                          cgcgTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGtt        <  1:1907587/71‑1 (MQ=255)
                                                           gcgTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGttt       >  1:2412296/1‑71 (MQ=255)
                                                           gcgTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGttt       <  1:2126898/71‑1 (MQ=255)
                                                            cgTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTg      <  1:3120526/71‑1 (MQ=255)
                                                             gtagtaAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTc            <  1:2240471/64‑1 (MQ=255)
                                                                gtaAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTGACCGtt                            <  1:475581/45‑1 (MQ=255)
                                                                gtaAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGtt                            <  1:1022622/45‑1 (MQ=255)
                                                                gtaAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGtt                            <  1:2298116/45‑1 (MQ=255)
                                                                gtaAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGt  >  1:1711950/1‑71 (MQ=255)
                                                                gtaAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGt  <  1:3463488/71‑1 (MQ=255)
                                                                gtaAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGt  <  1:2272110/71‑1 (MQ=255)
                                                                gtaAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGt  <  1:1039109/71‑1 (MQ=255)
                                                                gtaAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGt  <  1:2491895/71‑1 (MQ=255)
                                                                   aaTCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCg                              >  1:1540795/1‑40 (MQ=255)
                                                                   |                                                                   
CGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGT  >  minE/1014745‑1014879

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: