Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,015,709 A→G 21.2% R120R (CGT→CGC fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported1,015,7090AG21.2% 124.3 / 28.7 67R120R (CGT→CGCfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base A (26/26);  new base G (2/12);  total (29/38)
Fisher's exact test for biased strand distribution p-value = 3.05e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.43e-01

TATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCC  >  Exported/1015642‑1015778
                                                                   |                                                                     
tATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACgcgc                                                                    >  1:2326450/1‑71 (MQ=255)
  tCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTg                                                                  >  1:3654795/1‑71 (MQ=255)
       gATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCtt                                                             <  1:1704682/71‑1 (MQ=255)
        atatAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCttt                                                            <  1:1357801/71‑1 (MQ=255)
        atatAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCttt                                                            <  1:2118962/71‑1 (MQ=255)
        atatAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCttt                                                            <  1:1647556/71‑1 (MQ=255)
         tataGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCttt                                                            <  1:1981250/70‑1 (MQ=255)
          ataAGTGTTATAGACACCGCGCGCCAGCGCCGCTCCATCACCACCACCGGTCCATACACGCTGCCCttttt                                                          >  1:4085079/1‑71 (MQ=255)
             aGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTAc                                                       >  1:2649153/1‑71 (MQ=255)
              gtgtTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTAcc                                                      <  1:5075048/71‑1 (MQ=255)
                   ataGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACcaaca                                                  >  1:5127351/1‑70 (MQ=255)
                    taGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACcaacaa                                                 >  1:3402511/1‑70 (MQ=255)
                    taGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACcaacaa                                                 >  1:2528372/1‑70 (MQ=255)
                    taGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACaat                                                <  1:4727397/71‑1 (MQ=255)
                      gACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACaataa                                              >  1:3811583/1‑71 (MQ=255)
                       acacCGCGCGCCAGCGCCGCTTCATCACCACCCCCGGTCCATACCCGCTGCCCTTTTTTACCAACaataat                                             >  1:520643/1‑71 (MQ=255)
                        cacCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATc                                            <  1:3447993/71‑1 (MQ=255)
                         acCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATcg                                           <  1:2695996/71‑1 (MQ=255)
                          ccGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATcgc                                          >  1:3982361/1‑71 (MQ=255)
                           cgcgcgCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATcg                                           <  1:3658047/69‑1 (MQ=255)
                            gcgcgcCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCgg                                        >  1:232508/1‑71 (MQ=255)
                             cgcgcCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCCttttt                                                          >  1:5052715/1‑52 (MQ=255)
                                gcCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGccc                                                               <  1:3002665/44‑1 (MQ=255)
                                 ccAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGccc                                                               <  1:4759638/43‑1 (MQ=255)
                                   aGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCgg                                  >  1:1551041/1‑70 (MQ=255)
                                     cgccgcTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTAt                               >  1:4722473/1‑71 (MQ=255)
                                      gccgcTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACaat                                                >  1:1775721/1‑53 (MQ=255)
                                      gccgcTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACaat                                                >  1:5430502/1‑53 (MQ=255)
                                      gccgcTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATc                                            >  1:4629508/1‑57 (MQ=255)
                                       ccgcTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATc                              <  1:1063355/70‑1 (MQ=255)
                                        cgcTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCt                            <  1:3515801/71‑1 (MQ=255)
                                        cgcTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCt                            <  1:3896823/71‑1 (MQ=255)
                                         tcTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACcgac                                                   >  1:1429097/2‑47 (MQ=16)
                                          cTTCATCACCACCACCGGTCCATACACGCTGCCCtt                                                             <  1:623383/36‑1 (MQ=255)
                                          cTTCATCACCACCACCGGTCCATACACGCTGCCCtt                                                             <  1:5502752/36‑1 (MQ=255)
                                          cTTCATCACCACCACCGGTCCATACACGCTGCCCtt                                                             <  1:5361848/36‑1 (MQ=255)
                                           ttCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTg                           <  1:1524142/69‑1 (MQ=255)
                                            tcatcaCCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCgg                                  <  1:237579/61‑1 (MQ=255)
                                             catcaCCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTggc                         <  1:2661484/69‑3 (MQ=255)
                                             catcaCCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGAGGTGCCGGTATCCTgtc                         <  1:2363479/69‑3 (MQ=255)
                                              atcaCCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCgg                                        >  1:2489002/1‑53 (MQ=255)
                                              atcaCCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGa                          >  1:3467340/1‑67 (MQ=255)
                                              atcaCCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGa                          >  1:5142457/1‑67 (MQ=255)
                                              atcaCCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAgg                      <  1:770359/71‑1 (MQ=255)
                                               tcaCCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGt                     >  1:4035952/1‑71 (MQ=255)
                                               tcaCCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGt                     <  1:243615/71‑1 (MQ=255)
                                               tcaCCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGa                          <  1:646429/66‑1 (MQ=255)
                                                 accaccaccGGTCCATACACGCTGCCCTTTTTTACCAACaataat                                             >  1:4356031/1‑45 (MQ=255)
                                                 accaccaccGGTCCATACACGCTGCCCTTTTTTACCAACAATAATcg                                           >  1:1572158/1‑47 (MQ=255)
                                                  ccaccaccGGTCCATACACGCTGCCCTTTTTTACCAACaata                                               <  1:431388/42‑1 (MQ=255)
                                                  ccaccaccGGTCCATACACGCTGCCCTTTTTTACCAACAATAATc                                            <  1:1340906/45‑1 (MQ=255)
                                                     ccgccGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAg               <  1:2558444/68‑1 (MQ=255)
                                                     ccaccGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAg               <  1:1060718/71‑1 (MQ=255)
                                                      caccGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGa              <  1:3155127/71‑1 (MQ=255)
                                                      caccGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGa              <  1:304152/71‑1 (MQ=255)
                                                        ccGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCgg                                  >  1:3911609/1‑49 (MQ=255)
                                                        ccGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGaa             >  1:1963509/1‑70 (MQ=255)
                                                          ggTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGaa             >  1:3296863/1‑68 (MQ=255)
                                                          ggTCCACACGCGCTGGCCTTTTTTACCGACGATGATcgc                                          <  1:5383109/39‑1 (MQ=9)
                                                           gTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCgg                                  <  1:5573713/46‑1 (MQ=255)
                                                           gTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCgg                                  <  1:2627401/46‑1 (MQ=255)
                                                           gTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGcc         >  1:4816617/1‑71 (MQ=255)
                                                           gTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGc          <  1:37240/70‑1 (MQ=255)
                                                            tCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCgg                                  <  1:4286053/45‑1 (MQ=255)
                                                             ccATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTa                                >  1:3679966/1‑46 (MQ=255)
                                                             ccATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGa                          >  1:3388531/1‑52 (MQ=255)
                                                               aTACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAg               <  1:1308118/61‑1 (MQ=255)
                                                                tACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGc          <  1:3271662/65‑1 (MQ=255)
                                                                 acacGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGaa             >  1:3667765/1‑61 (MQ=255)
                                                                  cacGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCg                                   >  1:3166586/1‑38 (MQ=255)
                                                                  cacGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTcgcc  >  1:5557374/1‑71 (MQ=255)
                                                                   |                                                                     
TATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCC  >  Exported/1015642‑1015778

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: