Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,014,611 A→G 30.8% G486G (GGT→GGC fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,014,6110AG30.8% 60.5 / 37.4 52G486G (GGT→GGCfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base A (15/21);  new base G (8/8);  total (23/29)
Fisher's exact test for biased strand distribution p-value = 7.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.52e-01

CCCAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTT  >  minE/1014543‑1014676
                                                                    |                                                                 
cccAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGaccacc                                                                 >  1:5442996/1‑71 (MQ=255)
cccAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGaccacc                                                                 >  1:1637925/1‑71 (MQ=255)
 ccAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGaccacc                                                                 <  1:1525995/70‑1 (MQ=255)
     ttCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATAc                                                            <  1:607906/71‑1 (MQ=255)
         ggATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGcc                                                        <  1:3217333/71‑1 (MQ=255)
                 ctcacaTTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAg                                                 >  1:1360593/3‑70 (MQ=255)
                 aaCACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAg                                                 <  1:4134195/70‑1 (MQ=255)
                 aaCACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAg                                                 <  1:2641148/70‑1 (MQ=255)
                 aaCACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAg                                                 <  1:2520645/70‑1 (MQ=255)
                 aaCACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAg                                                 <  1:1980413/70‑1 (MQ=255)
                        tCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCaa                                                      >  1:5568517/1‑58 (MQ=255)
                        tCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAgccgcc                                         >  1:924183/1‑71 (MQ=255)
                           aGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCgtg                                      <  1:3329862/71‑1 (MQ=255)
                           aGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCgt                                       <  1:2919898/70‑1 (MQ=255)
                             gCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAgccgcc                                         <  1:1879519/66‑1 (MQ=255)
                               tctTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGa                                                               <  1:2680070/42‑1 (MQ=255)
                                ctTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGttttt                                 >  1:4765620/1‑71 (MQ=255)
                                 ttAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCa                                                       <  1:462459/48‑1 (MQ=255)
                                  aaaTACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATa                                                             >  1:4431457/2‑41 (MQ=255)
                                   aaTACTTCCCTGCGCCAATACAGCGGCCGGACCACCGAt                                                              >  1:5015040/1‑39 (MQ=255)
                                   aaTACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATa                                                             >  1:1266360/1‑40 (MQ=255)
                                   aaTACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATa                                                             >  1:4257964/1‑40 (MQ=255)
                                   aaTACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGa                                                >  1:1823041/1‑53 (MQ=255)
                                     tACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAgg                                                    <  1:3616730/47‑1 (MQ=255)
                                        ttCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAg                                                 <  1:162245/47‑1 (MQ=255)
                                          cccTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCa                         <  1:5435483/69‑1 (MQ=255)
                                          cccTGCGCCAATACAGCGGCCGGACCACCGATACAGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCa                         <  1:3040774/69‑1 (MQ=255)
                                            cTGCGCCAGTACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGcgt                        <  1:2813926/68‑3 (MQ=255)
                                            cTGCGCCAGTACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTcgg                     >  1:5648506/1‑69 (MQ=255)
                                            cTGCGCCAGTACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTcg                      >  1:5141058/1‑69 (MQ=255)
                                            cTGCGCCAGTACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTc                       <  1:227697/69‑1 (MQ=255)
                                            cTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCa                         >  1:5498553/1‑67 (MQ=255)
                                            cTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCg                     <  1:406601/71‑1 (MQ=255)
                                               cgcCAGTACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGtaa                  >  1:4676832/1‑69 (MQ=255)
                                               cgcCAGTACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGta                   <  1:2331422/70‑2 (MQ=255)
                                               cgcCAGTACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGta                   <  1:2271464/70‑2 (MQ=255)
                                               cgcCAGTACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGt                    <  1:3096114/69‑1 (MQ=255)
                                               cgcCAGTACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGt                    <  1:2188830/69‑1 (MQ=255)
                                               cgcCAGTACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGt                    <  1:4717182/69‑1 (MQ=255)
                                               cgcCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTc                   <  1:1776041/70‑1 (MQ=255)
                                                  cAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACct               >  1:5625313/1‑71 (MQ=255)
                                                  cAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACct               >  1:2937934/1‑71 (MQ=255)
                                                   agtACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACct               >  1:3344664/3‑70 (MQ=255)
                                                    gtACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACctgc             >  1:401206/2‑71 (MQ=255)
                                                    aTACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGccgccg                                        >  1:3263390/1‑44 (MQ=255)
                                                    aTACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCAt                        <  1:5582778/60‑1 (MQ=255)
                                                     tACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATcc                      <  1:2153796/61‑1 (MQ=255)
                                                     tACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACctgct            >  1:477448/1‑71 (MQ=255)
                                                      aCAGCGGCCGGACCACCGATACTGACAAGGTAGAAGCCGCCGTGTTTTTTAc                              >  1:2034702/1‑52 (MQ=255)
                                                         gccgCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACCTgctgact        >  1:2804282/4‑71 (MQ=255)
                                                          ccgCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACCTgctgactg       >  1:1923016/3‑71 (MQ=255)
                                                           ggCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTAc                              <  1:3324028/47‑1 (MQ=255)
                                                             ccGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACCTGCTGACTGCgg    >  1:1332908/1‑71 (MQ=255)
                                                               ggACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACCTGCTGACTGCGGtt  <  1:3406118/71‑1 (MQ=255)
                                                               ggACCACCGATACTGCCAAGGTAGAAGCCGCCGTGttttt                                 <  1:2221431/40‑1 (MQ=255)
                                                                   caccGATACTGCCAAGGTAGAAGCCGCCGTGttttt                                 <  1:5541609/36‑1 (MQ=255)
                                                                    accGATACTGCCAAGGTAGAAGCCTCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGtt  <  1:1485907/66‑1 (MQ=255)
                                                                    |                                                                 
CCCAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTT  >  minE/1014543‑1014676

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: