Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,058,356 T→C 56.9% I91I (ATT→ATC ttdB → L‑tartrate dehydratase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,058,3560TC56.9% 13.1 / 36.2 51I91I (ATT→ATCttdBL‑tartrate dehydratase, beta subunit
Reads supporting (aligned to +/- strand):  ref base T (10/12);  new base C (20/9);  total (30/21)
Fisher's exact test for biased strand distribution p-value = 1.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.83e-01

TGGGAGATGGTCTCCGTCGGCCCGACAACCAGTATGCGTATGGAAAGTTTTGAACGTGAATTTATTGAGCAGACCGGCGTGAAACTGGTGGTTGGCAAAGGTGGTATGGGGCCG  >  minE/2058291‑2058404
                                                                 |                                                
tGGGAGATGGTCTCCGTCGGCCCGACAACCAGTATGCGTATGGAAAGTTTTGAACGTGAATTTATTGAg                                               <  1:2427866/69‑1 (MQ=255)
          tgtcCGTCGGCCCGACAACCAGTATGGGTATGGAAAGTTTTGAACGTGAATTTATCGAGCAGACCGGc                                      >  1:2314170/3‑68 (MQ=255)
          tgtcCGTCGGCCCGACAACCAGTATGCGTATGGAAAGTTTTGAACGTGAATTTATCGAGCAGACCgg                                       >  1:2423756/3‑67 (MQ=255)
          tgtcCGTCGGCCCGACAACCAGTATGCGTATGGAAAGTTTTGAACGTGAATTTATCGAGCAGACCGGc                                      >  1:1960843/3‑68 (MQ=255)
          tgtcCGTCGGCCCGACAACCAGTATGCGTATGGAAAGTTTTGAACGTGAATTTATCGAGCAGACCGGc                                      >  1:2515730/3‑68 (MQ=255)
          tgtcCGTCGGCCCGACAACCAGTATGCGTATGGAAAGTTTTGAACGTGAATTTATCGAGCAGACCGGc                                      >  1:2438618/3‑68 (MQ=255)
          tgtcCGTCGGCCCGACAACCAGTATGCGTATGGAAAGTTTTGAACGTGAATTTATCGAGCAGACCGGc                                      >  1:2401962/3‑68 (MQ=255)
          tgtcCGTCGGCCCGACAACCAGTATGCGTATGGAAAGTTTTGAACGTGAATTTATCGAGCAGACCGGc                                      >  1:1037922/3‑68 (MQ=255)
          tgtcCGTCGGCCCGACAACCAGTATGCGTATGGAAAGTTTTGAACGTGAATTTATCGAGCAGACCGGc                                      >  1:2052860/3‑68 (MQ=255)
          tgtcCGTCGGCCCGACAACCAGTATGCGTATGGAAAGTTTTGAACGTGAATTTATCGAGCAGACCGGc                                      >  1:2979171/3‑68 (MQ=255)
          tgtcCGTCGGCCCGACAACCAGTATGCGTATGGAAAGTTTTGAACGTGAATTTATCGAGCAGACCGGc                                      >  1:302303/3‑68 (MQ=255)
          tgtcCGTCGGCCCGACAACCAGTATGCGTATGGAAAGTTTTGAACGTGAATTTATCGAGCAGACCGGc                                      >  1:3052377/3‑68 (MQ=255)
          tgtcCGTCGGCCCGACAACCAGTATGCGTATGGAAAGTTTTGAACGTGAATTTATCGAGCAGACCGGc                                      >  1:1839672/3‑68 (MQ=255)
          tgtcCGTCGGCCCGACAACCAGTATGCGTATGGAAAGTTTTGAACGTGAATTTATCGAGCAGACCGGc                                      >  1:1149435/3‑68 (MQ=255)
          tgtcCGTCGGCCCGACAACCAGTATGCGTATGGAAAGTTTTGAACGTGAATTTATCGAGCAGACCGGc                                      >  1:412068/3‑68 (MQ=255)
          tgtcCGTCGGCCCGACAACCAGTATGCGTATGGAAAGTTTTGAACGTGAATTTATCGAGCAGACCGGc                                      >  1:1171428/3‑68 (MQ=255)
          tgtcCGTCGGCCCGACAACCAGTATGCGTATGGAAAGTTTTGAACGTGAATTTATCGAGCAGACCGGc                                      >  1:1753332/3‑68 (MQ=255)
          tgtcCGTCGGCCCGACAACCAGTATGCGTATGGAAAGTTTTGAACGTGAATTTATCGAGCAGACCGGc                                      >  1:1688499/3‑68 (MQ=255)
          tgtcCGTCGGCCCGACAACCAGTATGCGTATGGAAAGTTTTGAACGTGAATTTATCGAGCAGACCGGc                                      >  1:1580866/3‑68 (MQ=255)
          tgtcCGTCGGCCCGACAACCAGTATGCGTATGGAAAGTTTTGAACGTGAATTTATCGAGCAGACCGGc                                      >  1:3176810/3‑68 (MQ=255)
          tgtcCGTCGGCCCGACAACCAGTATGCGTATGGAAAGTTTTGAACGTGAATTTATCGAGCAGACCGGc                                      >  1:473946/3‑68 (MQ=255)
                tCGGCCCGACAACCAGTATGCGTATGGAAAGTTTTGAACGTGAATTTATTGAGCAGACCGGCGTGa                                  <  1:2756728/66‑1 (MQ=255)
                           aCCAGTATGCGTATGGAAAGTTTTGAACGTGAATTTATTGAGCAGACCGGCGTGa                                  <  1:496492/55‑1 (MQ=255)
                             cAGTATGCGTATGGAAAGTTTTGAACGTGAATTTATTGAGCAGACTGGCGTGAAACTGGTGGTTGGCaa                  >  1:2347759/1‑69 (MQ=255)
                             cAGTATGCGTATGGAAAGTTTTGAACGTGAATTTATTGAGCAGACCGGCGTGAAACTGGTGGTTGGCaa                  >  1:1676470/1‑69 (MQ=255)
                             cAGTATGCGTATGGAAAGTTTTGAACGTGAATTTATTGAGCAGACCGGCGTGAAACTGGTGGTTGGCaa                  >  1:2440257/1‑69 (MQ=255)
                             cAGTATGCGTATGGAAAGTTTTGAACGTGAATTTATTGAGCAGACCGGCGTGAAACTGGTGGTTGGCaa                  >  1:2198845/1‑69 (MQ=255)
                             cAGTATGCGTATGGAAAGTTTTGAACGTGAATTTATTGAGCAGACCGGCGTGAAACTGGTGGTTGGCaa                  >  1:2804500/1‑69 (MQ=255)
                             cAGTATGCGTATGGAAAGTTTTGAACGTGAATTTATTGAGCAGACCGGCGTGAAACTGGTGGTTGGCaa                  >  1:1990395/1‑69 (MQ=255)
                             cAGTATGCGTATGGAAAGTTTTGAACGTGAATTTATTGAGCAGACCGGCGTGAAACTGGTGGTTGGCaa                  >  1:3316311/1‑69 (MQ=255)
                                  tGCGTATGGAAAGTTTTGAACGTGAATTTATTGAGCAGACCGGCGTGa                                  >  1:550015/1‑48 (MQ=255)
                                  tGCGTATGGAAAGTTTTGAACGTGAATTTATTGAGCAGACCGGCGTGAAACTGGTGGTTGGCAAAggt              >  1:2376512/1‑68 (MQ=255)
                                  tGCGTATGGAAAGTTTTGAACGTGAATTTATTGAGCAGACCGGCGTGAAACTGGTGGTTGGCAAAggt              >  1:1918359/1‑68 (MQ=255)
                                    cGTATGGAAAGTTTTGAACGTGAATTTATTGAGCAGACCGGCGTGAAACTGGTGGTTGGCAAAggtggt           <  1:698276/69‑1 (MQ=255)
                                    cGTATGGAAAGTTTTGAACGTGAATTTATTGAGCAGACCGGCGTGAAACTGGTGGTTGGCAAAggtggt           <  1:1810567/69‑1 (MQ=255)
                                    cGTATGGAAAGTTTTGAACGTGAATTTATTGAGCAGACCGGCGTGAAACTGGTGGTTGGCAAAggtggt           <  1:1800980/69‑1 (MQ=255)
                                    cGTATGGAAAGTTTTGAACGTGAATTTATTGAGCAGACCGGCGTGAAACTGGTGGTTGGCAAAggtggt           <  1:1796937/69‑1 (MQ=255)
                                    cGTATGGAAAGTTTTGAACGTGAATTTATTGAGCAGACCGGCGTGAAACTGGTGGTTGGCAAAggtgg            <  1:2004873/68‑1 (MQ=255)
                                             aGTTTTGAACGTGAATTTATTGAGCAGACCGGCGTGAAACTGGTGGTTGGCAAAggtgg            <  1:2369633/59‑1 (MQ=255)
                                             aGTTTTGAACGTGAATTTATTGAGCAGACCGGCGTGAAACTGGTGGTTGGCAAAggtgg            <  1:1860590/59‑1 (MQ=255)
                                             aGTTTTGAACGTGAATTTATTGAGCAGACCGGCGTGAAACTGGTGGTTGGCAAAggtgg            <  1:319431/59‑1 (MQ=255)
                                             aGTTTTGAACGTGAATTTATTGAGCAGACCGGCGTGAAACTGGTGGTTGGCAAAggtgg            <  1:145846/59‑1 (MQ=255)
                                             aGTTTTGAACGTGAATTTATCGAGCAGACCGGCGTGAAACTGGTGGTCGGCAAAGGCGGAATGGGGCCg  <  1:899520/69‑1 (MQ=255)
                                             aGTTTTGAACGTGAATTTATCGAGCAGACCGGCGTGAAACTGGTGGTCGGCAAAGGCGGAATGGGGCCg  <  1:2512113/69‑1 (MQ=255)
                                             aGTTTTGAACGTGAATTTATCGAGCAGACCGGCGTGAAACTGGTGGTCGGCAAAGGCGGAATGGGGCCg  <  1:784600/69‑1 (MQ=255)
                                             aGTTTTGAACGTGAATTTATCGAGCAGACCGGCGTGAAACTGGTGGTCGGCAAAGGCGGAATGGGGCCg  <  1:13947/69‑1 (MQ=255)
                                             aGTTTTGAACGTGAATTTATCGAGCAGACCGGCGTGAAACTGGTGGTCGGCAAAGGCGGAATGGGGCCg  <  1:3341077/69‑1 (MQ=255)
                                             aGTTTTGAACGTGAATTTATCGAGCAGACCGGCGTGAAACTGGTGGTCGGCAAAGGCGGAATGGGGCCg  <  1:3116664/69‑1 (MQ=255)
                                             aGTTTTGAACGTGAATTTATCGAGCAGACCGGCGTGAAACTGGTGGTCGGCAAAGGCGGAATGGGGCCg  <  1:2135314/69‑1 (MQ=255)
                                              gTTTTGAACGTGAATTTATCGAGCAGACCGGCGTGAAACTGGTGGTCGGCAAAGGCGGAATGGGGCCg  <  1:1652453/68‑1 (MQ=255)
                                              gTTTTGAACGTGAATTTATCGAGCAGACCGGCGTGAAACTGGTGGTCGGCAAAGGCGGAATGGGGCCg  <  1:2590862/68‑1 (MQ=255)
                                                                 |                                                
TGGGAGATGGTCTCCGTCGGCCCGACAACCAGTATGCGTATGGAAAGTTTTGAACGTGAATTTATTGAGCAGACCGGCGTGAAACTGGTGGTTGGCAAAGGTGGTATGGGGCCG  >  minE/2058291‑2058404

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: