Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,037,015 A→C 29.3% V61V (GTA→GTC rsxD → predicted inner membrane oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,037,0150AC29.3% 57.2 / 36.7 58V61V (GTA→GTCrsxDpredicted inner membrane oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (24/17);  new base C (11/6);  total (35/23)
Fisher's exact test for biased strand distribution p-value = 7.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

CTCGTTCAGATCCTGTTGGCATCGGTTAGTGCTCTGTTAGCCGAAGCTCTCGTACTCAAACTACGCAAGCAGTCGGTAGCCGCAACGTTGAAAGATAACTCAGCATTGCTGACAGGCT  >  minE/1036962‑1037079
                                                     |                                                                
ctcGTTCAGATCCTGTTGGCATCGGTTAGTGCTCTGTTAGCCGAAGCTCTCGTACTCAAACTACGCa                                                     >  1:1756757/1‑67 (MQ=255)
ctcGTTCAGATCCTGTTGGCATCGGTTAGTGCTCTGTTAGCCGAAGCTCTCGTACTCAAACTACGCa                                                     >  1:1258474/1‑67 (MQ=255)
ctcGTTCAGATCCTGTTGGCATCGGTTAGTGCTCTGTTAGCCGAAGCTCTCGTACTCAAACTACGCa                                                     >  1:1696699/1‑67 (MQ=255)
ctcGTTCAGATCCTGTTGGCATCGGTTAGTGCTCTGTTAGCCGAAGCTCTCGTACTCAAACTACGCa                                                     >  1:1370038/1‑67 (MQ=255)
ctcGTTCAGATCCTGTTGGCATCGGTTAGTGCTCTGTTAGCCGAAGCTCTCGTACTCAAACTACGCa                                                     >  1:474290/1‑67 (MQ=255)
ctcGTTCAGATCCTGTTGGCATCGGTTAGTGCTCTGTTAGCCGAAGCTCTCGTACTCAAACTACGCa                                                     >  1:1461171/1‑67 (MQ=255)
ctcGTTCAGATCCTGTTGGCATCGGTTAGTGCTCTGTTAGCCGAAGCTCTCGTACTCAAACTACGCa                                                     >  1:386943/1‑67 (MQ=255)
ctcGTTCAGATCCTGTTGGCATCGGTTAGTGCTCTGTTAGCCGAAGCTCTCGTACTCAAACTACGCa                                                     >  1:3400495/1‑67 (MQ=255)
ctcGTTCAGATCCTGTTGGCATCGGTTAGTGCTCTGTTAGCCGAAGCTCTCGTACTCAAACTACGCa                                                     >  1:3398719/1‑67 (MQ=255)
ctcGTTCAGATCCTGTTGGCATCGGTTAGTGCTCTGTTAGCCGAAGCTCTCGTACTCAAACTACGCa                                                     >  1:2853124/1‑67 (MQ=255)
ctcGTTCAGATCCTGTTGGCATCGGTTAGTGCTCTGTTAGCCGAAGCTCTAGTACTCAAACTACGc                                                      >  1:2939594/1‑66 (MQ=255)
     tCAGATCCTGTTGGCATCGGTTAGTGCTCTGTTAGCCGAAGCTCTCGTACTCAAACTACGCAAGCAGTc                                              <  1:2317680/69‑1 (MQ=255)
          tCCTGTTGGCATCGGTTAGTGCTCTGTTAGCCGAAGCTCTCGTACTCAAACTACGCAAGCa                                                 >  1:2206695/1‑61 (MQ=255)
          tCCTGTTGGCATCGGTTAGTGCTCTGTTAGCCGAAGCTCTCGTACTCAAACTACGCAAGCa                                                 >  1:1047117/1‑61 (MQ=255)
          tCCTGTTGGCATCGGTTAGTGCTCTGTTAGCCGAAGCTCTCGTACTCAAACTACGCAAGCa                                                 >  1:25905/1‑61 (MQ=255)
          tCCTGTTGGCATCGGTTAGTGCTCTGTTAGCCGAAGCTCTCGTACTCAAACTACGCAAGCa                                                 >  1:660476/1‑61 (MQ=255)
          tCCTGTTGGCATCGGTTAGTGCTCTGTTAGCCGAAGCTCTCGTACTCAAACTACGCAAGCa                                                 >  1:807694/1‑61 (MQ=255)
          tCCTGTTGGCATCGGTTAGTGCTCTGTTAGCCGAAGCTCTCGTACTCAAACTACGCAAGCa                                                 >  1:1249717/1‑61 (MQ=255)
          tCCTGTTGGCATCGGTTAGTGCTCTGTTAGCCGAACCTCTCGTACTCAAACTACGCAAGCa                                                 >  1:1587423/1‑61 (MQ=255)
          tCCTGTTGGCATCGTTAGTGCTCTGTTAGCCGAAGCTCTCGTACTCAAACTACGCAAGCa                                                 >  1:2317609/1‑60 (MQ=255)
                 ggCATCGGTTAGTGCTCTGTTAGCCGAAGCTCTCGTACTCAAACTACGCAAGCAGTCGGTAGCCGCAAc                                  >  1:65755/1‑69 (MQ=255)
                      cGGTCAGTGCTCTTCTGGCCGAAGCTCTCGTCCTCAAACTACGCAAGCAGTCGGTAGCTGCAACGTTGa                             <  1:3010080/69‑1 (MQ=255)
                      cGGTCAGTGCTCTTCTGGCCGAAGCTCTCGTCCTCAAACTACGCAAGCAGTCGGTAGCTGCAACGTTGa                             <  1:99564/69‑1 (MQ=255)
                      cGGTCAGTGCTCTTCTGGCCGAAGCTCTCGTCCTCAAACTACGCAAGCAGTCGGTAGCTGCAACGTTGa                             <  1:1350729/69‑1 (MQ=255)
                      cGGTCAGTGCTCTTCTGGCCGAAGCTCTCGTCCTCAAACTACGCAAGCAGTCGGTAGCTGCAACGTTGa                             <  1:571809/69‑1 (MQ=255)
                      cGGTCAGTGCTCTTCTGGCCGAAGCTCTCGTCCTCAAACTACGCAAGCAGTCGGTAGCTGCAACGTTGa                             <  1:2772357/69‑1 (MQ=255)
                       ggTCAGTGCTCTTCTGGCCGAAGCTCTCGTCCTCAAACTACGCAAGCAGTCGGTAGCTGCAACGTTGa                             <  1:3377708/68‑1 (MQ=255)
                         ttAGTGCTCTGTTAGCCGAAGCTCTCGTACTCAAACTACGCAAGCAGTCGGTAGCCGCAACGTTGaaa                           <  1:659864/68‑1 (MQ=255)
                         ttAGTGCTATGTTAGCCGAAGCTCTCGTACTCAAACTACGCAAGCAGTCGGTAGCCGCAACGTTGaaa                           <  1:2297133/68‑1 (MQ=255)
                                      aGCCGAAGCTCTCGTACTCAAACTACGCAAGCAGTCGGTAGCCGCAACGTTGAAAGATAACTCAGc                <  1:1808305/66‑1 (MQ=255)
                                      aGCCGAAGCTCTCGTACTCAAACTACGCAAGCAGTCGGTAGCCGCAACGTTGAAAGATAACTCAGCAtt             <  1:2877315/69‑1 (MQ=255)
                                      aGCCGAAGCTCTCGTACTCAAACTACGCAAGCAGTCGGTAGCCGCAACGTTGAAAGATAACTCAGCAtt             <  1:3004360/69‑1 (MQ=255)
                                         cGAAGCTCTCGTACTCAAACTACGCAAGCAGTCTGTAGCCGCAACGTTGAAAGATAACTCAGCATTGCt          <  1:221015/69‑1 (MQ=255)
                                         cGAAGCTCTCGTACTCAAACTACGCAAGCAGTCGGTAGCCGCAACGTTGAAAGATAACTCAGCATTGCt          <  1:359722/69‑1 (MQ=255)
                                         cGAAGCTCTCGTACTCAAACTACGCAAGCAGTCGGTAGCCGCAACGTTGAAAGATAACTCAGCATTGCt          <  1:15712/69‑1 (MQ=255)
                                         cGAAGCTCTCGTACTCAAACTACGCAAGCAGTCGGTAGCCGCAACGTTGAAAGATAACTCAGCATTGCt          <  1:728430/69‑1 (MQ=255)
                                         cGAAGCTCTCGTACTCAAACTACGCAAGCAGTCGGTAGCCGCAACGTTGAAAGATAACTCAGCATTGCt          <  1:3027985/69‑1 (MQ=255)
                                         cGAAGCTCTCGTACTCAAACTACGCAAGCAGTCGGTAGCCGCAACGTTGAAAGATAACTCAGCATTGCt          <  1:1709890/69‑1 (MQ=255)
                                          gAAGCTCTCGTCCTCAAACTACGCAAGCAGTCGGTAGCTGCaa                                   >  1:1840134/1‑43 (MQ=38)
                                          gAAGCTCTCGTCCTCAAACTACGCAAGCAGTCGGTAGCTGCAACGTTGAAAGATAACTCAGCATTGCt          >  1:3409183/1‑68 (MQ=255)
                                          gAAGCTCTCGTCCTCAAACTACGCAAGCAGTCGGTAGCTGCAACGTTGAAAGATAACTCAGCATTGCt          >  1:2006612/1‑68 (MQ=255)
                                          gAAGCTCTCGTCCTCAAACTACGCAAGCAGTCGGTAGCTGCAACGTTGAAAGATAACTCAGCATTGCt          >  1:147662/1‑68 (MQ=255)
                                          gAAGCTCTCGTCCTCAAACTACGCAAGCAGTCGGTAGCTGCAACGTTGAAAGATAACTCAGCATTGCt          >  1:1469886/1‑68 (MQ=255)
                                          gAAGCTCTCGTCCTCAAACTACGCAAGCAGTCGGTAGCTGCAACGTTGAAAGATAACTCAGCATTGCt          >  1:1397988/1‑68 (MQ=255)
                                          gAAGCTCTCGTCCTCAAACTACGCAAGCAGTCGGTAGCTGCAACGTTGAAAGATAACTCAGCATTGCt          >  1:3139442/1‑68 (MQ=255)
                                          gAAGCTCTCGTCCTCAAACTACGCAAGCAGTCGGTAGCTGCAACGTTGAAAGATAACTCAGCATTGCt          >  1:1251899/1‑68 (MQ=255)
                                          gAAGCTCTCGTCCTCAAACTACGCAAGCAGTCGGTAGCTGCAACGTTGAAAGATAACTCAGCATTGCt          >  1:813178/1‑68 (MQ=255)
                                          gAAGCTCTCGTCCTCAAACTACGCAAGCAGTCGGTAGCTGCAACGTTGAAAGATAACTCAGCATTGCt          >  1:2331652/1‑68 (MQ=255)
                                          gAAGCTCTCGTCCTCAAACTACGCAAGCAGTCGGTAGCTGCAACGTTGAAAGATAACTCAGCATTGCt          >  1:1651220/1‑68 (MQ=255)
                                          gAAGCTCTCGTACTCAAACTACGCAAGCAGTCGGTAGCCGCAACGTTGAAAGATAACTCAGCATTGCt          <  1:1886223/68‑1 (MQ=255)
                                          gAAGCTCTCGTACTCAAACTACGCAAGCAGTCGGTAGCCGCAACGTTGAAAGATAACTCAGCATTGCt          >  1:3343237/1‑68 (MQ=255)
                                          gAAGCTCTCGTACTCAAACTACGCAAGCAGTCGGTAGCCGCAACGTTGAAAGATAACTCAGCATTGCt          >  1:2083345/1‑68 (MQ=255)
                                          gAAGCTCTCGTACTCAAACTACGCAAGCAGTCGGTAGCCGCAACGTTGAAAGATAACTCAGCATTGCt          >  1:1533595/1‑68 (MQ=255)
                                          gAAGCTCTCGTACTCAAACTACGCAAGCAGTCGGTAGCCGCAACGTTGAAAGATAACTCAGCATTGCt          >  1:1600081/1‑68 (MQ=255)
                                                   gTACTCAAACTACGCAAGCAGTCGGTAGCCGCAACGTTGAAAGATAACTCAGCATTGCTGACAGGCt  <  1:1985353/67‑1 (MQ=255)
                                                   gTACTCAAACTACGCAAGCAGTCGGTAGCCGCAACGTTGAAAGATAACTCAGCATTGCTGACAGGCt  <  1:2972689/67‑1 (MQ=255)
                                                   gTACTCAAACTACGCAAGCAGTCGGTAGCCGCAACGTTGAAAGATAACTCAGCATTGCTGACAGGCt  <  1:1547181/67‑1 (MQ=255)
                                                   gTACTCAAACTACGCAAGCAGTCGGTAGCCGCAACGTTGAAAGATAACTCAGCATTGCTGACAGGCt  <  1:2294686/67‑1 (MQ=255)
                                                     |                                                                
CTCGTTCAGATCCTGTTGGCATCGGTTAGTGCTCTGTTAGCCGAAGCTCTCGTACTCAAACTACGCAAGCAGTCGGTAGCCGCAACGTTGAAAGATAACTCAGCATTGCTGACAGGCT  >  minE/1036962‑1037079

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: