Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,014,488 C→G 25.8% V527V (GTG→GTC fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,014,4880CG25.8% 46.6 / 15.0 31V527V (GTG→GTCfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base C (9/14);  new base G (3/5);  total (12/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.14e-01

TATTTCACACAGCGGGTGCATTGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGAT  >  minE/1014423‑1014555
                                                                 |                                                                   
tATTTCACACAGCGGGTGCATTGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAg                                                                >  1:802756/1‑71 (MQ=255)
   ttCACACAGCGGGTGCATTGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACaa                                                                 <  1:574330/67‑1 (MQ=255)
     cacacaGCGGGTGCATTGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATa                                                            <  1:1409133/70‑1 (MQ=255)
       cacaGCGGGTGCATTGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAAc                                                         <  1:2397066/71‑1 (MQ=255)
        acaGCGGGTGCATTGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACg                                                        >  1:766524/1‑71 (MQ=255)
               gTGCATTGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATcc                                                                     <  1:2570528/51‑1 (MQ=255)
                   aTTGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCt                                             <  1:879644/71‑1 (MQ=255)
                        gtgAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTcc                                          <  1:3225121/69‑1 (MQ=255)
                          gAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACg                                                        >  1:520733/1‑53 (MQ=255)
                          gAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACg                                                        >  1:3318998/1‑53 (MQ=255)
                               gTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAAt                                  <  1:786651/70‑1 (MQ=255)
                                tATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAAttt                                <  1:3271043/71‑1 (MQ=255)
                                     gctgGAAGAAGTCGTTACCTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCa                            <  1:1875849/70‑1 (MQ=255)
                                     gctgGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGaa                                                 <  1:1684018/49‑1 (MQ=255)
                                      ctgGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACGCCCGGAAATCTTCCACTTCAATTTTCCAg                           <  1:1928436/70‑1 (MQ=255)
                                         gaagaaGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGAt                         <  1:1659224/69‑1 (MQ=255)
                                          aagaagTCGTTACCTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGc                      <  1:2995623/71‑1 (MQ=255)
                                           agaagTCGTTACCTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGCt                     <  1:3026090/71‑1 (MQ=255)
                                           agaagTCGTTACCTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGCt                     <  1:1022793/71‑1 (MQ=255)
                                             aagTCATTTCCTTTATCATCCACAAGGATAAACGc                                                       >  1:879320/1‑35 (MQ=255)
                                               gTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCt                                             <  1:2939095/43‑1 (MQ=255)
                                                tCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCa                 <  1:1330090/70‑1 (MQ=255)
                                                 cATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCAtt               >  1:22054/1‑71 (MQ=255)
                                                  aTTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTc              >  1:3510638/1‑71 (MQ=255)
                                                  aTTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTc              <  1:1984600/71‑1 (MQ=255)
                                                    ttCCTTTATCATCCACAAGGATAAACGCCGGGAGATCTTCCACTTCaa                                   >  1:2030932/1‑48 (MQ=255)
                                                     accTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGCTTCCATAcc             >  1:3306792/2‑69 (MQ=255)
                                                        tttGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGCTTCCATACCCAgct        >  1:199777/1‑69 (MQ=255)
                                                         ttgtcatcGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGCTTCCATACCCAGCTc       >  1:122209/4‑71 (MQ=255)
                                                              atcGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGCTTCCATACCCAGCTCCGGat  <  1:2037959/71‑1 (MQ=255)
                                                              atcCACAAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGa   <  1:3450025/70‑1 (MQ=255)
                                                              atcCACAAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGa   >  1:412057/1‑70 (MQ=255)
                                                                 |                                                                   
TATTTCACACAGCGGGTGCATTGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGAT  >  minE/1014423‑1014555

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: