Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 49,231 A→G 44.8% A18A (GCT→GCC apaG ← protein associated with Co2+ and Mg2+ efflux

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE49,2310AG44.8% 21.5 / 96.0 116A18A (GCT→GCCapaGprotein associated with Co2+ and Mg2+ efflux
Reads supporting (aligned to +/- strand):  ref base A (33/31);  new base G (34/18);  total (67/49)
Fisher's exact test for biased strand distribution p-value = 1.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

Only 100 of 116 total aligned reads displayed.

TCGCCCCAGATTGCGTATGGTTACGGTATAAGCAAAAACGTAACGTTCATTATCAGGTGAAGATTGAGCCTCAATGTAGACGCTTTGAACCTGAATACACACTCGGGGCGAATTGATCATCGTTAACTCTCC  >  minE/49165‑49296
                                                                  |                                                                 
tCGCCCCAGATTGCGTATGGTTACGGTATAAGCAAAAACGTAACGTTCATTATCAGGTGAAGATTGAGc                                                                 <  1:513979/69‑1 (MQ=255)
tCGCCCCAGATTGCGTATGGTTACGGTATAAGCAAAAACGTAACGTTCATTATCAGGTGAAGATTGAGc                                                                 <  1:2121036/69‑1 (MQ=255)
tCGCCCCAGATTGCGTATGGTTACGGTATAAGCAAAAACGTAACGTTCATTATCAGGTGAAGATTGAGc                                                                 <  1:2143061/69‑1 (MQ=255)
tCGCCCCAGATTGCGTATGGTTACGGTATAAGCAAAAACGTAACGTTCATTATCAGGTGAAGATTGAGc                                                                 <  1:2471013/69‑1 (MQ=255)
tCGCCCCAGATTGCGTATGGTTACGGTATAAGCAAAAACGTAACGTTCATTATCAGGTGAAGATTGAGc                                                                 <  1:911730/69‑1 (MQ=255)
tCGCCCCAGATTGCGTATGGTTACGGTATAAGCAAAAACGTAACGTTCATTATCAGGTGAAGATTGAGc                                                                 <  1:764677/69‑1 (MQ=255)
tCGCCCCAGATTGCGTATGGTTACGGTATAAGCAAAAACGTAACGTTCATTATCAGGTGAAGATTGAGc                                                                 <  1:690041/69‑1 (MQ=255)
tCGCCCCAGATTGCGTATGGTTACGGTATAAGCAAAAACGTAACGTTCATTATCAGGTGAAGATTGAGc                                                                 <  1:3005129/69‑1 (MQ=255)
tCGCCCCAGATTGCGTATGGTTACGGTATAAGCAAAAACGTAACGTTCATTATCAGGTGAAGATTGAGc                                                                 <  1:686174/69‑1 (MQ=255)
tCGCCCCAGATTGCGTATGGTTACGGTATAAGCAAAAACGTAACCTTCATTATCAGGTGAAGATTGAGc                                                                 <  1:150806/69‑1 (MQ=255)
       tgATTGCGTATGGTTACGGTATAAGCAAAAACGTAACGTTCATTATCAGGTGAAGATTGGGCCTCAATg                                                          <  1:168256/68‑1 (MQ=255)
       aGATTGCTTATGGTTACGGTATAAGCAAAAACGTAACGTTCATTATCAGGTGAAGATTGGGCCTCAATg                                                          <  1:893285/69‑1 (MQ=255)
       aGATTGCGTATGGTTACGGTATAAGCAAAAACGTAACGTTCATTATCAGGTGAAGATTGGGCCTCAATg                                                          <  1:1394783/69‑1 (MQ=255)
       aGATTGCGTATGGTTACGGTATAAGCAAAAACGTAACGTTCATTATCAGGTGAAGATTGGGCCTCAATg                                                          <  1:872140/69‑1 (MQ=255)
       aGATTGCGTATGGTTACGGTATAAGCAAAAACGTAACGTTCATTATCAGGTGAAGATTGGGCCTCAATg                                                          <  1:454463/69‑1 (MQ=255)
       aGATTGCGTATGGTTACGGTATAAGCAAAAACGTAACGTTCATTATCAGGTGAAGATTGGGCCTCAATg                                                          <  1:2547613/69‑1 (MQ=255)
       aGATTGCGTATGGTTACGGTATAAGCAAAAACGTAACGTTCATTATCAGGTGAAGATTGGGCCTCAATg                                                          <  1:2557333/69‑1 (MQ=255)
       aGATTGCGTATGGTTACGGTATAAGCAAAAACGTAACGTTCATTATCAGGTGAAGATTGGGCCTCAATg                                                          <  1:1647673/69‑1 (MQ=255)
       aGATTGCGTATGGTTACGGTATAAGCAAAAACGTAACGTTCATTATCAGGTGAAGATTGGGCCTCAATg                                                          <  1:384603/69‑1 (MQ=255)
       aGATTGCGTATGGTTACGGTATAAGCAAAAACGTAACGTTCATTATCAGGTGAAGATTGGGCCTCAATg                                                          <  1:1575909/69‑1 (MQ=255)
       aGATTGCGTATGGTTACGGTATAAGCAAAAACGTAACGTTCATTATCAGGTGAAGATTGGGCCTCAATg                                                          <  1:1565703/69‑1 (MQ=255)
       aGATTGCGTATGGTTACGGTATAAGCAAAAACGTAACGTTCATTATCAGGTGAAGATTGGGCCTCAATg                                                          <  1:1241987/69‑1 (MQ=255)
       aGATTGCGTATGGTTACGGTATAAGCAAAAACGTAACGTTCATTATCAGGTGAAGATTGGGCCTCAATg                                                          <  1:3325692/69‑1 (MQ=255)
       aGATTGCGTATGGTTACGGTATAAGCAAAAACGTAACGTTCATTATCAGGTGAAGATTGGGCCTCAATg                                                          <  1:3188908/69‑1 (MQ=255)
       aGATTGCGTATGGTTACGGTATAAGCAAAAACGTAACGTTCATTATCAGGTGAAGATTGGGCCTCAATg                                                          <  1:1369600/69‑1 (MQ=255)
       aGATTGCGTATGGTTACGGTATAAGCAAAAACGTAACGTTCATTATCAGGTGAAGATTGAGCCTCAATg                                                          <  1:257111/69‑1 (MQ=255)
       aGATTGCGTATGGTTACGGTATAAGCAAAAACGTAACGTTCATTATCAGGTGAAGATTGAGCCTCAATg                                                          <  1:292229/69‑1 (MQ=255)
       aGATTGCGTATGGTTACGGTATAAGCAAAAACGTAACGTTCATTATCAGGTGAAGATTGAGCCTCAATg                                                          <  1:506520/69‑1 (MQ=255)
        gATTGCGTATGGTTACGGTATAAGCAAAAACGTAACGTTCATTATCAGGTGAAGATTGGGCCTCAATg                                                          <  1:3345646/68‑1 (MQ=255)
                 tGGTTACGGTATAAGCAAAAACGTAACGTTCATTATCAGGTGAAGATTGAGCCTCAATGTAg                                                       >  1:1041603/1‑62 (MQ=255)
                   gTTACGGTATAAGCAAAAACGTAACGTTCATTATCAGGTGAAGATTGAGCCTCAATGTAGACGCTTTGa                                              >  1:1138505/1‑69 (MQ=255)
                   gTTACGGTATAAGCAAAAACGTAACGTTCATTATCAGGTGAAGATTGAGCCTCAATGTAGACGCTTTGa                                              >  1:1157645/1‑69 (MQ=255)
                         gTATAAGCAAAAACGTAACGTTCATTATCAGGTGAAGATTGAGCCTCAATGTAGACGCTTTGAACCTg                                         <  1:1818908/68‑1 (MQ=255)
                         gTATAAGCAAAAACGTAACGTTCATTATCAGGTGAAGATTGAGCCTCAATGTAGACGCTTTGAACCTg                                         <  1:1253518/68‑1 (MQ=255)
                         gTATAAGCAAAAACGTAACGTTCATTATCAGGTGAAGATTGAGCCTCAATGTAGACGCTTTGAACCTg                                         <  1:2263278/68‑1 (MQ=255)
                          tataAGCAAAAACGTAACGTTCATTATCAGGTGAAGATTGAGCCTCAATGTAGACGCTTTGAACCTg                                         <  1:2294124/67‑1 (MQ=255)
                          tataAGCAAAAACGTAACGTTCATTATCAGGTGAAGATTGAGCCTCAATGTAGACGCTTTGAACCTGaa                                       >  1:1699779/1‑69 (MQ=255)
                          tataAGCAAAAACGTAACGTTCATTATCAGGTGAAGATTGAGCCTCAATGTAGACGCTTTGAACCTGaa                                       >  1:1744598/1‑69 (MQ=255)
                          tataAGCAAAAACGTAACGTTCATTATCAGGTGAAGATTGAGCCTCAATGTAGACGCTTTGAACCTGaa                                       >  1:610865/1‑69 (MQ=255)
                          tataAGCAAAAACGTAACGTTCATTATCAGGTGAAGATTGAGCCTCAATGTAGACGCTTTGAACCTGaa                                       >  1:1306673/1‑69 (MQ=255)
                              aGCAAAAACGTAACTTTCATTATCAGGTGAAGATTGAGCCCCAATGTAGACGc                                                   <  1:2408512/53‑1 (MQ=255)
                                 aaaaaCGTAACGTTCATTATCAGGTGAAGATTGAGCCTCAATGTAGACGCTTTGAACCTGAATaca                                   <  1:1947110/66‑1 (MQ=255)
                                         aacgtTCATTATCAGGTGAAGATTGGGCCTCAATGTa                                                        >  1:2628616/1‑37 (MQ=255)
                                         aacgtTCATTATCAGGTGAAGATTGGGCCTCAATGTAGACGCTTTGaa                                             >  1:66399/1‑48 (MQ=255)
                                         aacgtTCATTATCAGGTGAAGATTGGGCCTCAATGTAGACGCTTTGAACCTGAATACACACTCGGGGc                         >  1:917987/1‑68 (MQ=255)
                                         aacgtTCATTATCAGGTGAAGATTGGGCCTCAATGTAGACGCTTTGAACCTGAATACACACTCGGGGCg                        >  1:518099/1‑69 (MQ=255)
                                         aacgtTCATTATCAGGTGAAGATTGGGCCTCAATGTAGACGCTTTGAACCTGAATACACACTCGGGGCg                        >  1:517766/1‑69 (MQ=255)
                                         aacgtTCATTATCAGGTGAAGATTGGGCCTCAATGTAGACGCTTTGAACCTGAATACACACTCGGGGCg                        >  1:424045/1‑69 (MQ=255)
                                         aacgtTCATTATCAGGTGAAGATTGGGCCTCAATGTAGACGCTTTGAACCTGAATACACACTCGGGGCg                        >  1:495934/1‑69 (MQ=255)
                                         aacgtTCATTATCAGGTGAAGATTGGGCCTCAATGTAGACGCTTTGAACCTGAATACACACTCGGGGCg                        >  1:2770725/1‑69 (MQ=255)
                                         aacgtTCATTATCAGGTGAAGATTGGGCCTCAATGTAGACGCTTTGAACCTGAATACACACTCGGGGCg                        >  1:1465627/1‑69 (MQ=255)
                                         aacgtTCATTATCAGGTGAAGATTGGGCCTCAATGTAGACGCTTTGAACCTGAATACACACTCGGGGCg                        >  1:1475562/1‑69 (MQ=255)
                                         aacgtTCATTATCAGGTGAAGATTGGGCCTCAATGTAGACGCTTTGAACCTGAATACACACTCGGGGCg                        >  1:718389/1‑69 (MQ=255)
                                         aacgtTCATTATCAGGTGAAGATTGGGCCTCAATGTAGACGCTTTGAACCTGAATACACACTCGGGGCg                        >  1:1523718/1‑69 (MQ=255)
                                         aacgtTCATTATCAGGTGAAGATTGGGCCTCAATGTAGACGCTTTGAACCTGAATACACACTCGGGGCg                        >  1:2677068/1‑69 (MQ=255)
                                         aacgtTCATTATCAGGTGAAGATTGGGCCTCAATGTAGACGCTTTGAACCTGAATACACACTCGGGGCg                        >  1:20180/1‑69 (MQ=255)
                                         aacgtTCATTATCAGGTGAAGATTGGGCCTCAATGTAGACGCTTTGAACCTGAATACACACTCGGGGCg                        >  1:2199275/1‑69 (MQ=255)
                                         aacgtTCATTATCAGGTGAAGATTGGGCCTCAATGTAGACGCTTTGAACCTGAATACACACTCGGGGCg                        >  1:770128/1‑69 (MQ=255)
                                         aacgtTCATTATCAGGTGAAGATTGGGCCTCAATGTAGACGCTTTGAACCTGAATACACACTCGGGGCg                        >  1:2271903/1‑69 (MQ=255)
                                         aacgtTCATTATCAGGTGAAGATTGGGCCTCAATGTAGACGCTTTGAACCTGAATACACACTCGGGGCg                        >  1:2394620/1‑69 (MQ=255)
                                         aacgtTCATTATCAGGTGAAGATTGGGCCTCAATGTAGACGCTTTGAACCTGAATACACACTCGGGGCg                        >  1:2061938/1‑69 (MQ=255)
                                         aacgtTCATTATAAGGTGAAGATTGGGCCTCAATg                                                          >  1:2735726/1‑35 (MQ=38)
                                               cATTATCAGGTGAAGATTGAGCCTCAATGTAGACGCTTTGAACCTGAATACACTCTCGGGGCGAATTg                   >  1:3140677/1‑68 (MQ=255)
                                               cATTATCAGGTGAAGATTGAGCCTCAATGTAGACGCTTTGAACCTGAATACACACTCGGGGCGAATTg                   >  1:1427999/1‑68 (MQ=255)
                                               cATTATCAGGTGAAGATTGAGCCTCAATGTAGACGCTTTGAACCTGAATACACACTCGGGGCGAATTg                   >  1:1883919/1‑68 (MQ=255)
                                               cATTATCAGGTGAAGATTGAGCCTCAATGTAGACGCTTTGAACCTGAATACACACTCGGGGCGAATTg                   >  1:785111/1‑68 (MQ=255)
                                               cATTATCAGGTGAAGATTGAGCCTCAATGTAGACGCTTTGAACCTGAATACACACTCGGGGCGAATTg                   >  1:1665109/1‑68 (MQ=255)
                                               cATTATCAGGTGAAGATTGAGCCTCAATGTAGACGCTTTGAACCTGAATACACACTCGGGGCGAATTg                   >  1:2904286/1‑68 (MQ=255)
                                               cATTATCAGGTGAAGATTGAGCCTCAATGTAGACGCTTTGAACCTGAATACACACTCGGGGCGAATTg                   >  1:2852713/1‑68 (MQ=255)
                                                          gAAGATTGGGCCTCAATGTAGACGCTTTGAACCTGAATACACACTCg                             >  1:1909173/1‑47 (MQ=255)
                                                          gAAGATTGGGCCTCAATGTAGACGCTTTGAACCTGAATACACACTCg                             >  1:186530/1‑47 (MQ=255)
                                                          gAAGATTGGGCCTCAATGTAGACGCTTTGAACCTGAATACACACTCg                             >  1:919053/1‑47 (MQ=255)
                                                          gAAGATTGGGCCTCAATGTAGACGCTTTGAACCTGAATACACACTCg                             >  1:2207375/1‑47 (MQ=255)
                                                          gAAGATTGGGCCTCAATGTAGACGCTTTGAACCTGAATACACACTCg                             >  1:1194406/1‑47 (MQ=255)
                                                          gAAGATTGGGCCTCAATGTAGACGCTTTGAACCTGAATACACACTCg                             >  1:1857673/1‑47 (MQ=255)
                                                          gAAGATTGGGCCTCAATGTAGACGCTTTGAACCTGAATACACACTCg                             >  1:2498461/1‑47 (MQ=255)
                                                          gAAGATTGGGCCTCAATGTAGACGCTTTGAACCTGAATACACACTCg                             >  1:810886/1‑47 (MQ=255)
                                                          gAAGATTGGGCCTCAATGTAGACGCTTTGAACCTGAATACACACTCg                             >  1:2834349/1‑47 (MQ=255)
                                                          gAAGATTGAGCCTCAATGTAGACGCTTTGAACCTGAATACACACTCGGGGCg                        >  1:3300906/1‑52 (MQ=255)
                                                          gAAGATTGAGCCTCAATGTAGACGCTTTGAACCTGAATACACACTCGGGGCGAATTg                   <  1:1108331/57‑1 (MQ=255)
                                                          gAAGATTGAGCCTCAATGTAGACGCTTTGAACCTGAATACACACTCGGGGCGAATTg                   <  1:1409349/57‑1 (MQ=255)
                                                             gATTGAGCCTCAATGTAGACGCTTTGAACCTGAATACACACTCGGGGCGAATTGATCATCGTTAActc     <  1:622953/68‑1 (MQ=255)
                                                             gATTGAGCCTCAATGTAGACGCTTTGAACCTGAATACACACTCGGGGCGAATTGATCATCGTTAActc     <  1:3157568/68‑1 (MQ=255)
                                                             gATTGAGCCTCAATGTAGACGCTTTGAACCTGAATACACACTCGGGGCGAATTGATCATCGTTAActc     <  1:2923046/68‑1 (MQ=255)
                                                             gATTGAGCCTCAATGTAGACGCTTTGAACCTGAATACACACTCGGGGCGAATTGATCATCGTTAActc     <  1:2248390/68‑1 (MQ=255)
                                                              aTTGAGCCTCAATGTAGACGCTTTGCACCTGAATACACACTCGGGGCGAATTGATCATCGTTAAct      >  1:927619/1‑66 (MQ=255)
                                                              aTTGAGCCTCAATGTAGACGCTTTGAACCTGAATACACACTCGGGGCGAATTGATCATCGTTAAct      >  1:2801670/1‑66 (MQ=255)
                                                              aTTGAGCCTCAATGTAGACGCTTTGAACCTGAATACACACTCGGGGCGAATTGATCATCGTTAAct      >  1:1190254/1‑66 (MQ=255)
                                                              aTTGAGCCTCAATGTAGACGCTTTGAACCTGAATACACACTCGGGGCGAATTGATCATCGTTAAct      >  1:1646637/1‑66 (MQ=255)
                                                              aTTGAGCCTCAATGTAGACGCTTTGAACCTGAATACACACTCGGGGCGAATTGATCATCGTTAAct      >  1:1865274/1‑66 (MQ=255)
                                                              aTTGAGCCTCAATGTAGACGCTTTGAACCTGAATACACACTCGGGGCGAATTGATCATCGTTAAct      >  1:2552861/1‑66 (MQ=255)
                                                              aTTGAGCCTCAATGTAGACGCTTTGAACCTGAATACACACTCGGGGCGAATTGATCATCGTTAAct      >  1:2559291/1‑66 (MQ=255)
                                                              aTTGAGCCTCAATGTAGACGCTTTGAACCTGAATACACACTCGGGGCGAATTGATCATCGTTAAct      >  1:2627808/1‑66 (MQ=255)
                                                              aTTGAGCCTCAATGTAGACGCTTTGAACCTGAATACACACTCGGGGCGAATTGATCATCGTTAAct      >  1:388747/1‑66 (MQ=255)
                                                              aTTGAGCCTCAATGTAGACGCTTTGAACCTGAATACACACTCGGGGCGAATTGATCATCGTTAAct      >  1:734655/1‑66 (MQ=255)
                                                              aTTGAGCCTCAATGTAGACGCTTTGAACCTGAATACACACTCGGGGCGAATTGATCATCGTTAAct      >  1:3029672/1‑66 (MQ=255)
                                                              aTTGAGCCTCAATGTAGACGCTTTGAACCTGAATACACACTCGGGGCGAATTGATCATCGTTAAct      >  1:651232/1‑66 (MQ=255)
                                                              aTTGAGCCTCAATGTAGACGCTTTGAACCTGAATACACACTCGGGGCGAATTGATCATCGTTAAct      >  1:316189/1‑66 (MQ=255)
                                                              aTTGAGCCTCAATGTAGACGCTTTGAACCTGAATACACACTCGGGGCGAATTGATCATCGTTAAct      >  1:3349854/1‑66 (MQ=255)
                                                                tGAGCCTCAATGTAGACGCTTTGAACCTGAATACACACTCGGGGCGAATTGATCATCGTTAACTCTcc  >  1:2856152/1‑68 (MQ=255)
                                                                  |                                                                 
TCGCCCCAGATTGCGTATGGTTACGGTATAAGCAAAAACGTAACGTTCATTATCAGGTGAAGATTGAGCCTCAATGTAGACGCTTTGAACCTGAATACACACTCGGGGCGAATTGATCATCGTTAACTCTCC  >  minE/49165‑49296

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: