Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,014,542 T→A 25.3% G509G (GGA→GGT fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,014,5420TA25.3% 132.3 / 43.2 71G509G (GGA→GGTfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base T (30/23);  new base A (9/9);  total (39/32)
Fisher's exact test for biased strand distribution p-value = 7.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGA  >  minE/1014472‑1014608
                                                                      |                                                                  
cATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCAtt                                                                    >  1:3199116/1‑71 (MQ=255)
  tttCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTcc                                                                  >  1:489388/1‑71 (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)
                        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)
                           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)
                             gCCGGGAAATCTTCTACTTCGATTTTCCAGATAGCTTCCATACCCAg                                                               >  1:4561960/1‑47 (MQ=255)
                               cGGGAAATCTTCTACTTCGATTTTCCAGATAGCTTCCATACCCAGCTCCGGATAAGCGACGCACTCCAgat                                     <  1:1093234/71‑3 (MQ=255)
                                gggAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATATTCAACACATTCCAGGct                                    <  1:2595090/71‑1 (MQ=255)
                                 ggAAATCTTCTACTTCGATTTTCCAGATAGCTTCCATACCCAgc                                                              <  1:450747/44‑2 (MQ=255)
                                 ggAAATCTTCTACTTCGATTTTCCAGATAGCTTCCATACCCAgc                                                              <  1:3154114/44‑2 (MQ=255)
                                 ggAAATCTTCTACTTCGATTTTCCAGATAGCTTCCATACCCAgc                                                              <  1:878322/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:3316729/39‑1 (MQ=25)
                                             cTTCGATTTTCCAGATAGCTTCCATACCCAGCTCCGGat                                                       <  1:5003720/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)
                                                              gCTTCCATTCCCAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGc      >  1:838174/1‑71 (MQ=255)
                                                              gCTTCCATTCCCAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGc      >  1:2413168/1‑71 (MQ=255)
                                                              gCTTCCATTCCCAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGc      <  1:3358408/71‑1 (MQ=255)
                                                               cTTCCATTCCCAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGcc     <  1:1383428/71‑1 (MQ=255)
                                                               cTTCCATTCCCAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGcc     <  1:1576348/71‑1 (MQ=255)
                                                               cTTCCATTCCCAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGcc     >  1:836335/1‑71 (MQ=255)
                                                                 tCCATTCCCAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGc         >  1:4193030/1‑65 (MQ=255)
                                                                  ccATTCCCAGTTCCGGATATTCAACACATTCCAgg                                      >  1:652606/1‑35 (MQ=255)
                                                                  ccATTCCCAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGa  >  1:3342986/1‑71 (MQ=255)
                                                                     ttCCCAGTTCCGGATATTCAACACATTCCAGGCTCtt                                 <  1:381481/37‑1 (MQ=255)
                                                                     ttCCCAGTTCCGGATATTCAACACATTCCAGGCTCtt                                 <  1:716008/37‑1 (MQ=255)
                                                                     ttCCCAGTTCCGGATATTCAACACATTCCAGGCTCtt                                 <  1:1178778/37‑1 (MQ=255)
                                                                      |                                                                  
CATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGA  >  minE/1014472‑1014608

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: