Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,014,521 A→G 38.3% I516I (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,5210AG38.3% 60.7 / 80.4 81I516I (ATT→ATCfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base A (29/21);  new base G (16/15);  total (45/36)
Fisher's exact test for biased strand distribution p-value = 6.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.60e-01

TGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCC  >  minE/1014453‑1014586
                                                                    |                                                                 
tGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAAt                                                                 >  1:2839048/1‑71 (MQ=255)
 gTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAAt                                                                 >  1:565999/1‑70 (MQ=255)
   aTCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAAtttt                                                              >  1:535536/1‑71 (MQ=255)
        ctgGAAGAAGTCGTTACCTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAg                                                          <  1:3389495/70‑1 (MQ=255)
           gaagaaGTCGTTACCTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGa                                                         >  1:4708865/1‑68 (MQ=255)
           gaagaaGTCGTTACCTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGAtag                                                      >  1:4840571/1‑69 (MQ=255)
           gaagaaGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATg                                                       <  1:3396548/70‑1 (MQ=255)
           gaagaaGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATg                                                       >  1:1281765/1‑70 (MQ=255)
           caagaagTCGTTACCTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGAtag                                                      >  1:4840576/2‑69 (MQ=255)
               aagTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTc                                                  >  1:3585077/1‑71 (MQ=255)
                agTCGTTACCTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGCTTc                                                  <  1:1851522/70‑1 (MQ=255)
                  tCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCAt                                               >  1:1554285/1‑71 (MQ=255)
                   cATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCAtt                                              >  1:3199116/1‑71 (MQ=255)
                    aTTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAAtt                                                                <  1:3480619/52‑1 (MQ=255)
                    aTTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGAt                                                        >  1:5655099/1‑60 (MQ=255)
                     tttCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTcc                                            >  1:489388/1‑71 (MQ=255)
                      taccTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTc                                                             <  1:384221/51‑1 (MQ=25)
                          tttGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGCTTCCAt                                               <  1:5383765/63‑1 (MQ=255)
                           ttgtcatcGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGCTTCCATACCCAGCTc                                      >  1:1254032/4‑71 (MQ=255)
                            tgtcatcGACCAGGATAAACGCCGGGAAATCTTCTACTTCGAtttt                                                              >  1:2888147/3‑46 (MQ=21)
                            tgtcatcGACCAGGATAAACGCCGGGAAATCTTCTACTTCGAtttt                                                              >  1:1219716/3‑46 (MQ=21)
                                atcGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGCTTCCATACCCAGCTCCGGa                                  >  1:2639127/1‑70 (MQ=255)
                                atcGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGCTTCCATACCCAGCTCCGGa                                  >  1:5086069/1‑70 (MQ=255)
                                 tcCACAAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTcc                                                            <  1:3869038/43‑1 (MQ=255)
                                 tcCACAAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTcc                                                            <  1:2663123/43‑1 (MQ=255)
                                  ccACAAGGATAAACGCCGGGAAATCTTCCACTTCaa                                                                  <  1:995606/36‑1 (MQ=255)
                                  ccACAAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTccc                                           >  1:2747386/1‑59 (MQ=255)
                                   cacaAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGatat                               >  1:1079912/1‑70 (MQ=255)
                                   cacaAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGatat                               >  1:1175426/1‑70 (MQ=255)
                                   cacaAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATAtt                              >  1:453674/1‑71 (MQ=255)
                                    acaAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTcc                                            <  1:2981343/56‑1 (MQ=255)
                                     caAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGa                                                         <  1:1626221/42‑1 (MQ=255)
                                      aaGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATATTc                             >  1:5501114/1‑69 (MQ=255)
                                       aGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGa                                                         <  1:3495772/40‑1 (MQ=255)
                                       aGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCa                                                >  1:1172567/1‑49 (MQ=255)
                                           tAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGCtt                                                   <  1:4758803/42‑1 (MQ=255)
                                           tAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGCTTCCATACCCa                                          >  1:119644/1‑51 (MQ=255)
                                           tAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGCTTCCATACCCa                                          >  1:126505/1‑51 (MQ=255)
                                           tAAACGCCGGGAAATCTTCTACTTCGATTTTCAAGATAGCTTCCATACCCAGCTCCGGa                                  >  1:275311/1‑59 (MQ=255)
                                           tAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCg                                    <  1:1877283/57‑1 (MQ=255)
                                           tAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATATTCAACACAt                      >  1:2076490/1‑71 (MQ=255)
                                             aaCGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATATTCaa                           <  1:10312/64‑1 (MQ=255)
                                              aCGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGCTTCCa                                                >  1:1908448/1‑42 (MQ=255)
                                              aCGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGCTTCCATAcc                                            >  1:4521493/1‑46 (MQ=255)
                                               cGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATATTCAACACATTCCa                  <  1:4026277/71‑1 (MQ=255)
                                               cGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATATTCAACACATTCCa                  <  1:4032262/71‑1 (MQ=255)
                                               cGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATATTCAACACATTCCa                  <  1:358922/71‑1 (MQ=255)
                                                gCCGGGAAATCTTCTACTTCGATTTTCCAGATAGCtt                                                   <  1:4636530/37‑1 (MQ=37)
                                                gCCGGGAAATCTTCTACTTCGATTTTCCAGATAGCTTCCATACCCAg                                         >  1:4561960/1‑47 (MQ=255)
                                                  cGGGAAATCTTCTACTTCGATTTTCCAGATAGCTTCCATACCCAGCTCCGGATAAGCGACGCACTCCAgat               <  1:1093234/71‑3 (MQ=255)
                                                  cGGGAAATCTTCTACTTCGATTTTCCAGATAGATTcc                                                 >  1:1111159/1‑37 (MQ=25)
                                                   gggAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATATTCAACACATTCCAGGct              <  1:2595090/71‑1 (MQ=255)
                                                    ggAAATCTTCTACTTCGATTTTCCAGATAGCTTCCATACCCAgc                                        <  1:878322/44‑2 (MQ=255)
                                                    ggAAATCTTCTACTTCGATTTTCCAGATAGCTTCCATACCCAgc                                        <  1:450747/44‑2 (MQ=255)
                                                    ggAAATCTTCTACTTCGATTTTCCAGATAGCTTCCATACCCAgc                                        <  1:3154114/44‑2 (MQ=255)
                                                    ggAAATCTTCTACTTCGATTTTCCAGATAGCTTCCATACCCAGCTCCg                                    <  1:3488393/48‑1 (MQ=255)
                                                    ggAAATCTTCTACTTCGATTTTCCAGATAGATTCCATACCCAGCTCCGGATAAGCGACGCACTCCAgatg              >  1:5477308/1‑67 (MQ=255)
                                                      aaaTCTTCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTc                                      <  1:4764458/44‑1 (MQ=255)
                                                      aaaTCTTCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTc                                      <  1:1354722/44‑1 (MQ=255)
                                                      aaaTCTTCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATATTCAACACATTCCAGGCTCtt           <  1:1585056/71‑1 (MQ=255)
                                                          cTTCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATATTCAACACATTCCAGGCTCTTAATa       >  1:5071271/1‑71 (MQ=255)
                                                           ttCTACTTCGATTTTCCAGATAGCTTCCATACCCAGCTcc                                     <  1:4732056/40‑1 (MQ=255)
                                                           ttCTACTTCGATTTTCCAGATAGCTTCCATACCCAGCTcc                                     <  1:2219373/40‑1 (MQ=255)
                                                           ttCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATATTc                             <  1:1886663/48‑1 (MQ=255)
                                                           ttCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATATTc                             <  1:351022/48‑1 (MQ=255)
                                                           ttCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATATTCAACACATTCCAGGCTCTTAATa       >  1:3467154/1‑70 (MQ=255)
                                                              cACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATATTCAACACATTCCAgg                >  1:5163477/1‑58 (MQ=255)
                                                              cACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATATTCAACACATTCCAGGCTCTTAATa       >  1:3412052/1‑67 (MQ=255)
                                                              cACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATATTCAACACATTCCAGGCTCTTAATa       >  1:388147/1‑67 (MQ=255)
                                                              cACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACtt    <  1:4335428/70‑1 (MQ=255)
                                                              cACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATAATCAACACATTCCAGGCTCTTAATACtt    <  1:2827364/70‑1 (MQ=255)
                                                              cACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCAGGATATTCAACACATTCCAgg                >  1:3877106/1‑58 (MQ=255)
                                                              cACTTCAATTTTCCAGATGGCTTCCACTCCCAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTc   >  1:5353610/1‑71 (MQ=255)
                                                               aCTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATATTGAACACATTCCAGGCTCTTAATACTTcc  >  1:5670052/1‑71 (MQ=255)
                                                               aCTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTcc  >  1:3431314/1‑71 (MQ=255)
                                                                cTTCGATTTTCCAGATAGCTTCCATACCCAGCTCCGGat                                 <  1:5003720/39‑1 (MQ=25)
                                                                cTTCGATTTTCCAGATAGCTTCCATACCCAGCTCCGGat                                 <  1:3316729/39‑1 (MQ=25)
                                                                cTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATATTCAACACATTCCAGGc               <  1:499907/57‑1 (MQ=255)
                                                                 ttCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATATTCAACACATTCCAgg                >  1:1280102/1‑55 (MQ=255)
                                                                   cAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATATTCaa                           >  1:3713759/1‑42 (MQ=255)
                                                                   cAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATATTCAACACATTCCAgg                >  1:1895583/1‑53 (MQ=255)
                                                                   cAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATATTCAACACATTCCAg                 >  1:4543578/1‑52 (MQ=255)
                                                                    |                                                                 
TGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCC  >  minE/1014453‑1014586

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: