Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,014,533 G→A 33.5% A512A (GCC→GCT fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,014,5330GA33.5% 80.0 / 55.9 72A512A (GCC→GCTfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base G (28/19);  new base A (11/13);  total (39/33)
Fisher's exact test for biased strand distribution p-value = 3.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATA  >  minE/1014464‑1014597
                                                                     |                                                                
gaagaaGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATg                                                                  <  1:3396548/70‑1 (MQ=255)
gaagaaGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATg                                                                  >  1:1281765/1‑70 (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)
          tttCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTcc                                                       >  1:489388/1‑71 (MQ=255)
               tttGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGCTTCCAt                                                          <  1:5383765/63‑1 (MQ=255)
                ttgtcatcGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGCTTCCATACCCAGCTc                                                 >  1:1254032/4‑71 (MQ=255)
                     atcGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGCTTCCATACCCAGCTCCGGa                                             >  1:2639127/1‑70 (MQ=255)
                     atcGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGCTTCCATACCCAGCTCCGGa                                             >  1:5086069/1‑70 (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)
                           aaGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATATTc                                        >  1:5501114/1‑69 (MQ=255)
                            aGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCa                                                           >  1:1172567/1‑49 (MQ=255)
                                tAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGCtt                                                              <  1:4758803/42‑1 (MQ=255)
                                tAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGCTTCCATACCCa                                                     >  1:126505/1‑51 (MQ=255)
                                tAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGCTTCCATACCCa                                                     >  1:119644/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:4032262/71‑1 (MQ=255)
                                    cGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATATTCAACACATTCCa                             <  1:4026277/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:450747/44‑2 (MQ=255)
                                         ggAAATCTTCTACTTCGATTTTCCAGATAGCTTCCATACCCAgc                                                   <  1:878322/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:1354722/44‑1 (MQ=255)
                                           aaaTCTTCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTc                                                 <  1:4764458/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:351022/48‑1 (MQ=255)
                                                ttCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATATTc                                        <  1:1886663/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)
                                                           ttttCCAGATGGCTTCCATTCCCAGTTCCGGATATTCAACACATTCCa                             <  1:2060825/48‑1 (MQ=255)
                                                           ttttCCAGATGGCTTCCATTCCCAGTTCCGGATATTCAACACATTCCa                             <  1:3468184/48‑1 (MQ=255)
                                                           ttttCCAGATGGCTTCCATTCCCAGTTCCGGATATTCAACACATTCCAGGCTCTTAATa                  >  1:4637325/1‑59 (MQ=255)
                                                           ttttCCAGATGGCTTCCATTCCCAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTgcgc       >  1:4608584/1‑70 (MQ=255)
                                                            tttCCAGATGGCTTCCATTCCCAGTTCCGGATATTc                                        <  1:5380084/36‑1 (MQ=255)
                                                             ttCCAGATGGCTTCCATTCCCAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACtt               >  1:5137394/1‑60 (MQ=255)
                                                             ttCCAGATGCTTCCATTCCCAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCaa    <  1:839897/70‑1 (MQ=255)
                                                                cAGATGGCTTCCATTCCCAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATa  <  1:4875235/70‑1 (MQ=255)
                                                                   aTGGCTTCCATTCCCAGTTCCGGATATTCAACACATTCCAGGCTCTTaa                    >  1:1433669/1‑49 (MQ=255)
                                                                     |                                                                
GAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATA  >  minE/1014464‑1014597

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: