Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 948,935 C→A 35.3% V634V (GTG→GTT dos ← cAMP phosphodiesterase, heme‑regulated

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE948,9350CA35.3% 80.3 / 59.4 85V634V (GTG→GTTdoscAMP phosphodiesterase, heme‑regulated
Reads supporting (aligned to +/- strand):  ref base C (42/13);  new base A (27/3);  total (69/16)
Fisher's exact test for biased strand distribution p-value = 1.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.02e-01

TGCATTGCATCAGACACCTGATTAGGCAGTTGATTACTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGCCGGGATATGAATATTCTGGCTACGCCATTCTGCTAACTGACGGCAAGCTTCCGCGAT  >  minE/948873‑949003
                                                              |                                                                    
tgcattgcatCAGACACCTGATTAGGCAGTTGATTACTGCGAAAGTGCAGCGCCGACAAGTTCACGGAt                                                                <  1:3143618/69‑1 (MQ=255)
tgcattgcatCAGACACCTGATTAGGCAGTTGATTACTGCGAAAGTGCAGCGCCGACAAGTTCACGGAt                                                                <  1:2827979/69‑1 (MQ=255)
      gcatCAGACACCTGATTAGGCAGTTGATTACTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGcc                                                          >  1:89570/1‑69 (MQ=255)
      gcatCAGACACCTGATTAGGCAGTTGATTACTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGcc                                                          >  1:3298363/1‑69 (MQ=255)
      gcatCAGACACCTGATTAGGCAGTTGATTACTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGcc                                                          >  1:323680/1‑69 (MQ=255)
      gcatCAGACACCTGATTAGGCAGTTGATTACTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGcc                                                          >  1:3176572/1‑69 (MQ=255)
      gcatCAGACACCTGATTAGGCAGTTGATTACTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGcc                                                          >  1:3149166/1‑69 (MQ=255)
      gcatCAGACACCTGATTAGGCAGTTGATTACTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGcc                                                          >  1:3145819/1‑69 (MQ=255)
      gcatCAGACACCTGATTAGGCAGTTGATTACTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGcc                                                          >  1:285776/1‑69 (MQ=255)
      gcatCAGACACCTGATTAGGCAGTTGATTACTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGcc                                                          >  1:2684881/1‑69 (MQ=255)
      gcatCAGACACCTGATTAGGCAGTTGATTACTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGcc                                                          >  1:3416460/1‑69 (MQ=255)
      gcatCAGACACCTGATTAGGCAGTTGATTACTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGcc                                                          >  1:2314089/1‑69 (MQ=255)
      gcatCAGACACCTGATTAGGCAGTTGATTACTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGcc                                                          >  1:2165060/1‑69 (MQ=255)
      gcatCAGACACCTGATTAGGCAGTTGATTACTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGcc                                                          >  1:2010599/1‑69 (MQ=255)
      gcatCAGACACCTGATTAGGCAGTTGATTACTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGcc                                                          >  1:1961920/1‑69 (MQ=255)
      gcatCAGACACCTGATTAGGCAGTTGATTACTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGcc                                                          >  1:1407670/1‑69 (MQ=255)
      gcatCAGACACCTGATTAGGCAGTTGATTACTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGcc                                                          >  1:1426157/1‑69 (MQ=255)
      gcatCAGACACCTGATTAGGCAGTTGATTACTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGcc                                                          >  1:1102671/1‑69 (MQ=255)
      gcatCAGACACCTGATTAGGCAGTTGATTACTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGcc                                                          >  1:1236023/1‑69 (MQ=255)
      gcatCAGACACCTGATTAGGCAGTTGATTACTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGcc                                                          >  1:1180442/1‑69 (MQ=255)
      gcatCAGACACCTGATTAGGCAGTTGATTACTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGcc                                                          >  1:508804/1‑69 (MQ=255)
          cAGATACCTGATTAGGCAGCTGATTACTGCGAAAGTGCAGCGCAGACAAGTTAACGGATAACGCCGGGa                                                      >  1:2261067/1‑69 (MQ=255)
          cAGATACCTGATTAGGCAGCTGATTACTGCGAAAGTGCAGCGCAGACAAGTTAACGGATAACGCCGGGa                                                      >  1:1017650/1‑69 (MQ=255)
          cAGATACCTGATTAGGCAGCTGATTACTGCGAAAGTGCAGCGCAGACAAGTTAACGGATAACGCCGGGa                                                      >  1:2020147/1‑69 (MQ=255)
          cAGATACCTGATTAGGCAGCTGATTACTGCGAAAGTGCAGCGCAGACAAGTTAACGGATAACGCCGGGa                                                      >  1:2835972/1‑69 (MQ=255)
          cAGATACCTGATTAGGCAGCTGATTACTGCGAAAGTGCAGCGCAGACAAGTTAACGGATAACGCCGGGa                                                      >  1:1222052/1‑69 (MQ=255)
          cAGATACCTGATTAGGCAGCTGATTACTGCGAAAGTGCAGCGCAGACAAGTTAACGGATAACGCCGGGa                                                      >  1:3182548/1‑69 (MQ=255)
          cAGATACCTGATTAGGCAGCTGATTACTGCGAAAGTGCAGCGCAGACAAGTTAACGGATAACGCCGGGa                                                      >  1:1380376/1‑69 (MQ=255)
          cAGATACCTGATTAGGCAGCTGATTACTGCGAAAGTGCAGCGCAGACAAGTTAACGGATAACGCCGGGa                                                      >  1:1749792/1‑69 (MQ=255)
                  tGATTAGGCAGTTGATTACTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGTCGGGATATGa                                                 >  1:1806132/1‑66 (MQ=255)
                  tGATTAGGCAGTTGATTACTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGCCGGGATATGaa                                                >  1:634536/1‑67 (MQ=255)
                  tGATTAGGCAGTTGATTACTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGCCGGGATATGaa                                                >  1:1862288/1‑67 (MQ=255)
                  tGATTAGGCAGTTGATTACTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGCCGGGATATGaa                                                >  1:141927/1‑67 (MQ=255)
                  tGATTAGGCAGTTGATTACTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGCCGGGATATGaa                                                >  1:1770848/1‑67 (MQ=255)
                  tGATTAGGCAGTTGATTACTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGCCGGGATATGAata                                              >  1:2843639/1‑69 (MQ=255)
                  tGATTAGGCAGTTGATTACTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGCCGGGATATGAata                                              >  1:2256763/1‑69 (MQ=255)
                  tGATTAGGCAGTTGATTACTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGCCGGGATATGAata                                              >  1:2738563/1‑69 (MQ=255)
                  tGATTAGGCAGTTGATTACTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGCCGGGATATGAata                                              >  1:1786563/1‑69 (MQ=255)
                  tGATTAGGCAGTTGATTACTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGCCGGGATATGAat                                               >  1:2256365/1‑68 (MQ=255)
                   gATTAGGCAGCTGATTACTGCGAAAGTGCAGCGCAGACAAGTTAACgg                                                                  >  1:1725317/1‑48 (MQ=255)
                   gATTAGGCAGCTGATTACTGCGAAAGTGCAGCGCAGACAAGTTAACgg                                                                  >  1:551769/1‑48 (MQ=255)
                   gATTAGGCAGCTGATTACTGCGAAAGTGCAGCGCAGACAAGTTAACgg                                                                  >  1:125136/1‑48 (MQ=255)
                   gATTAGGCAGCTGATTACTGCGAAAGTGCAGCGCAGACAAGTTAACgg                                                                  >  1:2892755/1‑48 (MQ=255)
                   gATTAGGCAGCTGATTACTGCGAAAGTGCAGCGCAGACAAGTTAACgg                                                                  >  1:330669/1‑48 (MQ=255)
                   gATTAGGCAGCTGATTACTGCGAAAGTGCAGCGCAGACAAGTTAACgg                                                                  >  1:296130/1‑48 (MQ=255)
                   gATTAGGCAGCTGATTACTGCGAAAGTGCAGCGCAGACAAGTTAACgg                                                                  >  1:3143944/1‑48 (MQ=255)
                   gATTAGGCAGCTGATTACTGCGAAAGTGCAGCGCAGACAAGTTAACgg                                                                  >  1:3301473/1‑48 (MQ=255)
                   gATTAGGCAGCTGATTACTGCGAAAGTGCAGCGCAGACAAGTTAACgg                                                                  >  1:1620361/1‑48 (MQ=255)
                   gATTAGGCAGCTGATTACTGCGAAAGTGCAGCGCAGACAAGTTAACgg                                                                  >  1:1780168/1‑48 (MQ=255)
                   gATTAGGCAGCTGATTACTGCGAAAGTGCAGCGCAGACAAGTTAACgg                                                                  >  1:2576892/1‑48 (MQ=255)
                   gATTAGGCAGCTGATTACTGCGAAAGTGCAGCGCAGACAAGTTAACgg                                                                  >  1:2413977/1‑48 (MQ=255)
                   gATTAGGCAGCTGATTACTGCGAAAGTGCAGCGCAGACAAGTTAACgg                                                                  >  1:2318926/1‑48 (MQ=255)
                   gATTAGGCAGCTGATTACTGCGAAAGTGCAGCGCAGACAAGTTAACgg                                                                  >  1:1790447/1‑48 (MQ=255)
                   gATTAGGCAGCTGATTACTGCGAAAGTGCAGCGCAGACAAGTTAACgg                                                                  >  1:2019015/1‑48 (MQ=255)
                   gATTAGGCAGCTGATTACTGCGAAAGTGCAGCGCAGACAAGTTAACgg                                                                  >  1:1907217/1‑48 (MQ=255)
                   gATTAGGCAGCTGATTACTGCGAAAGTGCAGCGCAGACAAGTTAACgg                                                                  >  1:547797/1‑48 (MQ=255)
                   gATTAGGCAGCTGATTACTGCGAAAGTGCAGCGCAGACAAGTTAACgg                                                                  >  1:197229/1‑48 (MQ=255)
                   gATTAGGCAGCTGATTACTGCGAAAGTGCAGCGCAGACAAGTTAACgg                                                                  >  1:468293/1‑48 (MQ=255)
                              tGATTACTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGCCGGGATATGAatat                                             >  1:3404000/1‑58 (MQ=255)
                              tGATTACTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGCCGGGATATGAATATTCTGGCTACGc                                  <  1:3032982/69‑1 (MQ=255)
                              tGATTACTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGCCGGGATATGAATATTCTGGCTACGc                                  <  1:441573/69‑1 (MQ=255)
                              tGATTACTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGCCGGGATATGAATATTCTGGCTACGc                                  >  1:501837/1‑69 (MQ=255)
                              tGATTACTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGCCGGGATATGAATATTCTGGCTACGc                                  <  1:3382778/69‑1 (MQ=255)
                              tGATTACTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGCCGGGATATGAATATTCTGGCTACGc                                  <  1:2826341/69‑1 (MQ=255)
                               gATTACTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGCCGGGATATGAATATTCTGGCTACGc                                  <  1:2698493/68‑1 (MQ=255)
                                   aCTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGCCGGGATATGAATATTCTGGCTACGCCAtt                              >  1:1085672/1‑68 (MQ=255)
                                   aCTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGCCGGGATATGAATATTCTGGCTACGCCAtt                              >  1:2748241/1‑68 (MQ=255)
                                   aCTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGCCGGGATATGAATATTCTGGCTACGCCAt                               >  1:3260162/1‑67 (MQ=255)
                                   aCTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGCCGGGATATGAATATTCTGGCTACGCCATTc                             >  1:2345520/1‑69 (MQ=255)
                                   aCTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGCCGGGATATGAATATTCTGGCTACGCCATTc                             >  1:2123533/1‑69 (MQ=255)
                                   aCTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGCCGGGATATGAATATTCTGGCTACGCCATTc                             >  1:1910858/1‑69 (MQ=255)
                                   aCTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGCCGGGATATGAATATTCTGGCTACGCCATTc                             >  1:1789878/1‑69 (MQ=255)
                                      gCGAAAGTGCAGCGCAGACAAGTTAACGGATAACGCCGGGATATGAATATTCTGGCTACGCCATTCTg                           <  1:3287630/68‑1 (MQ=255)
                                      gCGAAAGTGCAGCGCAGACAAGTTAACGGATAACGCCGGGATATGAATATTCTGGCTACGCCATTCTg                           <  1:1304765/68‑1 (MQ=255)
                                        gAAAGTGCAGCTCCGACAAGTTCACGGATAACGCCGGGATATGAATATTCTGGCTACGCCATTCTGCt                         <  1:844297/68‑1 (MQ=255)
                                         aaaGTGCAGCGCCGACAAGTTCACGGATAACGCCGGGATATGAATATTCTGGCTACGCCATTCTGCt                         <  1:1577354/67‑1 (MQ=255)
                                         aaaGTGCAGCGCCGACAAGTTCACGGATAACGCCGGGATATGAATATTCTGGCTACGCCATTCTGCt                         <  1:1532897/67‑1 (MQ=255)
                                         aaaGTGCAGCGCCGACAAGTTCACGGATAACGCCGGGATATGAATATTCTGGCTACGCCATTCTGCt                         <  1:41496/67‑1 (MQ=255)
                                         aaaGTGCAGCGCCGACAAGTTCACGGATAACGCCGGGATATGAATATTCTGGCTACGCCATTCTGCTaa                       >  1:1183675/1‑69 (MQ=255)
                                         aaaGTGCAGCGCCGACAAGTTCACGGATAACGCCGGGATATGAATATTCTGGCTACGCCATTCTGCTaa                       >  1:705303/1‑69 (MQ=255)
                                         aaaGTGCAGCGCCGACAAGTTCACGGATAACGCCGGGATATGAATATTCTGGCTACGCCATTCTGCTaa                       >  1:729855/1‑69 (MQ=255)
                                         aaaGTGCAGCGCCGACAAGTTCACGGATAACGCCGGGATATGAATATTCTGGCTACGCCATTCTGCTaa                       >  1:2899987/1‑69 (MQ=255)
                                          aaGTGCAGCGCCGACAAGTTCACGGATAACGCCGGGATATGAATATTCTGGCTACGCCATTCTGCt                         <  1:855126/66‑1 (MQ=255)
                                               cAGCGCAGACAAGTTAACGGATAACGCCGGGATATGAATATTCTGGCTACGCCATTCTg                           <  1:392012/59‑1 (MQ=255)
                                                  cgcCGACAAGTTCACGGATAACGCCGGGATATGAATATTCTGGCTACGCCATTCTGCTAACTGACGGCa              <  1:284204/69‑1 (MQ=255)
                                                              cACGGATAACGCCGGGATATGAATATTCTGGCTACGCCATTCTGCTAACTGACGGCAAGCTTCCGCGAt  >  1:501102/1‑69 (MQ=255)
                                                              |                                                                    
TGCATTGCATCAGACACCTGATTAGGCAGTTGATTACTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGCCGGGATATGAATATTCTGGCTACGCCATTCTGCTAACTGACGGCAAGCTTCCGCGAT  >  minE/948873‑949003

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: