Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 948,926 C→A 30.6% S637S (TCG→TCT dos ← cAMP phosphodiesterase, heme‑regulated

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE948,9260CA30.6% 128.4 / 57.1 98S637S (TCG→TCTdoscAMP phosphodiesterase, heme‑regulated
Reads supporting (aligned to +/- strand):  ref base C (55/13);  new base A (27/3);  total (82/16)
Fisher's exact test for biased strand distribution p-value = 3.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.66e-01

AATACCCCAGGCGTGCATTGCATCAGACACCTGATTAGGCAGTTGATTACTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGCCGGGATATGAATATTCTGGCTACGCCATTCTGCTAACTGACGGCA  >  minE/948860‑948991
                                                                  |                                                                 
aaTACCCCAGGCGTGCATTGCATCAGACACCTGATTAGGCAGTTGATTACTGCGAAAGTGCAGCGCCg                                                                  >  1:965629/1‑68 (MQ=255)
aaTACCCCAGGCGTGCATTGCATCAGACACCTGATTAGGCAGTTGATTACTGCGAAAGTGCAGCGCCg                                                                  >  1:2620692/1‑68 (MQ=255)
aaTACCCCAGGCGTGCATTGCATCAGACACCTGATTAGGCAGTTGATTACTGCGAAAGTGCAGCGCCg                                                                  >  1:2867113/1‑68 (MQ=255)
aaTACCCCAGGCGTGCATTGCATCAGACACCTGATTAGGCAGTTGATTACTGCGAAAGTGCAGCGCCg                                                                  >  1:2174366/1‑68 (MQ=255)
aaTACCCCAGGCGTGCATTGCATCAGACACCTGATTAGGCAGTTGATTACTGCGAAAGTGCAGCGCCg                                                                  >  1:1231562/1‑68 (MQ=255)
aaTACCCCAGGCGTGCATTGCATCAGACACCTGATTAGGCAGTTGATTACTGCGAAAGTGCAGCGCCg                                                                  >  1:2115824/1‑68 (MQ=255)
aaTACCCCAGGCGTGCATTGCATCAGACACCTGATTAGGCAGTTGATTACTGCGAAAGTGCAGCGCCg                                                                  >  1:2972145/1‑68 (MQ=255)
aaTACCCCAGGCGTGCATTGCATCAGACACCTGATTAGGCAGTTGATTACTGCGAAAGTGCAGCGCCg                                                                  >  1:3008092/1‑68 (MQ=255)
aaTACCCCAGGCGTGCATTGCATCAGACACCTGATTAGGCAGTTGATTACTGCGAAAGTGCAGCGCCg                                                                  >  1:3090874/1‑68 (MQ=255)
aaTACCCCAGGCGTGCATTGCATCAGACACCTGATTAGGCAGTTGATTACTGCGAAAGTGCAGCGCCg                                                                  >  1:1770637/1‑68 (MQ=255)
aaTACCCCAGGCGTGCATTGCATCAGACACCTGATTAGGCAGTTGATTACTGCGAAAGTGCAGCGCCg                                                                  >  1:1734732/1‑68 (MQ=255)
aaTACCCCAGGCGTGCATTGCATCAGACACCTGATTAGGCAGTTGATTACTGCGAAAGTGCAGCGCCg                                                                  >  1:3169044/1‑68 (MQ=255)
aaTACCACAGGCGTGCATTGCATCAGACACCTGATTAGGCAGTTGATTACTGCGAAAGTGCAGCGCCg                                                                  >  1:4705/1‑68 (MQ=255)
             tgcattgcatCAGACACCTGATTAGGCAGTTGATTACTGCGAAAGTGCAGCGCCGACAAGTTCACGGAt                                                    <  1:2827979/69‑1 (MQ=255)
             tgcattgcatCAGACACCTGATTAGGCAGTTGATTACTGCGAAAGTGCAGCGCCGACAAGTTCACGGAt                                                    <  1:3143618/69‑1 (MQ=255)
                   gcatCAGACACCTGATTAGGCAGTTGATTACTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGcc                                              >  1:3176572/1‑69 (MQ=255)
                   gcatCAGACACCTGATTAGGCAGTTGATTACTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGcc                                              >  1:2684881/1‑69 (MQ=255)
                   gcatCAGACACCTGATTAGGCAGTTGATTACTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGcc                                              >  1:1102671/1‑69 (MQ=255)
                   gcatCAGACACCTGATTAGGCAGTTGATTACTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGcc                                              >  1:3416460/1‑69 (MQ=255)
                   gcatCAGACACCTGATTAGGCAGTTGATTACTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGcc                                              >  1:285776/1‑69 (MQ=255)
                   gcatCAGACACCTGATTAGGCAGTTGATTACTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGcc                                              >  1:2314089/1‑69 (MQ=255)
                   gcatCAGACACCTGATTAGGCAGTTGATTACTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGcc                                              >  1:1180442/1‑69 (MQ=255)
                   gcatCAGACACCTGATTAGGCAGTTGATTACTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGcc                                              >  1:89570/1‑69 (MQ=255)
                   gcatCAGACACCTGATTAGGCAGTTGATTACTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGcc                                              >  1:2165060/1‑69 (MQ=255)
                   gcatCAGACACCTGATTAGGCAGTTGATTACTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGcc                                              >  1:1236023/1‑69 (MQ=255)
                   gcatCAGACACCTGATTAGGCAGTTGATTACTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGcc                                              >  1:3298363/1‑69 (MQ=255)
                   gcatCAGACACCTGATTAGGCAGTTGATTACTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGcc                                              >  1:2010599/1‑69 (MQ=255)
                   gcatCAGACACCTGATTAGGCAGTTGATTACTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGcc                                              >  1:1961920/1‑69 (MQ=255)
                   gcatCAGACACCTGATTAGGCAGTTGATTACTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGcc                                              >  1:323680/1‑69 (MQ=255)
                   gcatCAGACACCTGATTAGGCAGTTGATTACTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGcc                                              >  1:1426157/1‑69 (MQ=255)
                   gcatCAGACACCTGATTAGGCAGTTGATTACTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGcc                                              >  1:3145819/1‑69 (MQ=255)
                   gcatCAGACACCTGATTAGGCAGTTGATTACTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGcc                                              >  1:508804/1‑69 (MQ=255)
                   gcatCAGACACCTGATTAGGCAGTTGATTACTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGcc                                              >  1:1407670/1‑69 (MQ=255)
                   gcatCAGACACCTGATTAGGCAGTTGATTACTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGcc                                              >  1:3149166/1‑69 (MQ=255)
                       cAGATACCTGATTAGGCAGCTGATTACTGCGAAAGTGCAGCGCAGACAAGTTAACGGATAACGCCGGGa                                          >  1:3182548/1‑69 (MQ=255)
                       cAGATACCTGATTAGGCAGCTGATTACTGCGAAAGTGCAGCGCAGACAAGTTAACGGATAACGCCGGGa                                          >  1:2835972/1‑69 (MQ=255)
                       cAGATACCTGATTAGGCAGCTGATTACTGCGAAAGTGCAGCGCAGACAAGTTAACGGATAACGCCGGGa                                          >  1:1017650/1‑69 (MQ=255)
                       cAGATACCTGATTAGGCAGCTGATTACTGCGAAAGTGCAGCGCAGACAAGTTAACGGATAACGCCGGGa                                          >  1:1749792/1‑69 (MQ=255)
                       cAGATACCTGATTAGGCAGCTGATTACTGCGAAAGTGCAGCGCAGACAAGTTAACGGATAACGCCGGGa                                          >  1:2261067/1‑69 (MQ=255)
                       cAGATACCTGATTAGGCAGCTGATTACTGCGAAAGTGCAGCGCAGACAAGTTAACGGATAACGCCGGGa                                          >  1:1222052/1‑69 (MQ=255)
                       cAGATACCTGATTAGGCAGCTGATTACTGCGAAAGTGCAGCGCAGACAAGTTAACGGATAACGCCGGGa                                          >  1:2020147/1‑69 (MQ=255)
                       cAGATACCTGATTAGGCAGCTGATTACTGCGAAAGTGCAGCGCAGACAAGTTAACGGATAACGCCGGGa                                          >  1:1380376/1‑69 (MQ=255)
                               tGATTAGGCAGTTGATTACTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGTCGGGATATGa                                     >  1:1806132/1‑66 (MQ=255)
                               tGATTAGGCAGTTGATTACTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGCCGGGATATGaa                                    >  1:1770848/1‑67 (MQ=255)
                               tGATTAGGCAGTTGATTACTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGCCGGGATATGaa                                    >  1:141927/1‑67 (MQ=255)
                               tGATTAGGCAGTTGATTACTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGCCGGGATATGaa                                    >  1:634536/1‑67 (MQ=255)
                               tGATTAGGCAGTTGATTACTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGCCGGGATATGaa                                    >  1:1862288/1‑67 (MQ=255)
                               tGATTAGGCAGTTGATTACTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGCCGGGATATGAata                                  >  1:2738563/1‑69 (MQ=255)
                               tGATTAGGCAGTTGATTACTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGCCGGGATATGAata                                  >  1:2256763/1‑69 (MQ=255)
                               tGATTAGGCAGTTGATTACTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGCCGGGATATGAata                                  >  1:2843639/1‑69 (MQ=255)
                               tGATTAGGCAGTTGATTACTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGCCGGGATATGAata                                  >  1:1786563/1‑69 (MQ=255)
                               tGATTAGGCAGTTGATTACTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGCCGGGATATGAat                                   >  1:2256365/1‑68 (MQ=255)
                                gATTAGGCAGCTGATTACTGCGAAAGTGCAGCGCAGACAAGTTAACgg                                                      >  1:197229/1‑48 (MQ=255)
                                gATTAGGCAGCTGATTACTGCGAAAGTGCAGCGCAGACAAGTTAACgg                                                      >  1:551769/1‑48 (MQ=255)
                                gATTAGGCAGCTGATTACTGCGAAAGTGCAGCGCAGACAAGTTAACgg                                                      >  1:547797/1‑48 (MQ=255)
                                gATTAGGCAGCTGATTACTGCGAAAGTGCAGCGCAGACAAGTTAACgg                                                      >  1:1620361/1‑48 (MQ=255)
                                gATTAGGCAGCTGATTACTGCGAAAGTGCAGCGCAGACAAGTTAACgg                                                      >  1:1725317/1‑48 (MQ=255)
                                gATTAGGCAGCTGATTACTGCGAAAGTGCAGCGCAGACAAGTTAACgg                                                      >  1:1780168/1‑48 (MQ=255)
                                gATTAGGCAGCTGATTACTGCGAAAGTGCAGCGCAGACAAGTTAACgg                                                      >  1:1790447/1‑48 (MQ=255)
                                gATTAGGCAGCTGATTACTGCGAAAGTGCAGCGCAGACAAGTTAACgg                                                      >  1:3143944/1‑48 (MQ=255)
                                gATTAGGCAGCTGATTACTGCGAAAGTGCAGCGCAGACAAGTTAACgg                                                      >  1:1907217/1‑48 (MQ=255)
                                gATTAGGCAGCTGATTACTGCGAAAGTGCAGCGCAGACAAGTTAACgg                                                      >  1:330669/1‑48 (MQ=255)
                                gATTAGGCAGCTGATTACTGCGAAAGTGCAGCGCAGACAAGTTAACgg                                                      >  1:2019015/1‑48 (MQ=255)
                                gATTAGGCAGCTGATTACTGCGAAAGTGCAGCGCAGACAAGTTAACgg                                                      >  1:296130/1‑48 (MQ=255)
                                gATTAGGCAGCTGATTACTGCGAAAGTGCAGCGCAGACAAGTTAACgg                                                      >  1:2892755/1‑48 (MQ=255)
                                gATTAGGCAGCTGATTACTGCGAAAGTGCAGCGCAGACAAGTTAACgg                                                      >  1:125136/1‑48 (MQ=255)
                                gATTAGGCAGCTGATTACTGCGAAAGTGCAGCGCAGACAAGTTAACgg                                                      >  1:2318926/1‑48 (MQ=255)
                                gATTAGGCAGCTGATTACTGCGAAAGTGCAGCGCAGACAAGTTAACgg                                                      >  1:468293/1‑48 (MQ=255)
                                gATTAGGCAGCTGATTACTGCGAAAGTGCAGCGCAGACAAGTTAACgg                                                      >  1:2413977/1‑48 (MQ=255)
                                gATTAGGCAGCTGATTACTGCGAAAGTGCAGCGCAGACAAGTTAACgg                                                      >  1:2576892/1‑48 (MQ=255)
                                gATTAGGCAGCTGATTACTGCGAAAGTGCAGCGCAGACAAGTTAACgg                                                      >  1:3301473/1‑48 (MQ=255)
                                           tGATTACTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGCCGGGATATGAatat                                 >  1:3404000/1‑58 (MQ=255)
                                           tGATTACTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGCCGGGATATGAATATTCTGGCTACGc                      >  1:501837/1‑69 (MQ=255)
                                           tGATTACTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGCCGGGATATGAATATTCTGGCTACGc                      <  1:3382778/69‑1 (MQ=255)
                                           tGATTACTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGCCGGGATATGAATATTCTGGCTACGc                      <  1:2826341/69‑1 (MQ=255)
                                           tGATTACTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGCCGGGATATGAATATTCTGGCTACGc                      <  1:3032982/69‑1 (MQ=255)
                                           tGATTACTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGCCGGGATATGAATATTCTGGCTACGc                      <  1:441573/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:2899987/1‑69 (MQ=255)
                                                      aaaGTGCAGCGCCGACAAGTTCACGGATAACGCCGGGATATGAATATTCTGGCTACGCCATTCTGCTaa           >  1:705303/1‑69 (MQ=255)
                                                      aaaGTGCAGCGCCGACAAGTTCACGGATAACGCCGGGATATGAATATTCTGGCTACGCCATTCTGCTaa           >  1:729855/1‑69 (MQ=255)
                                                      aaaGTGCAGCGCCGACAAGTTCACGGATAACGCCGGGATATGAATATTCTGGCTACGCCATTCTGCTaa           >  1:1183675/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)
                                                                  |                                                                 
AATACCCCAGGCGTGCATTGCATCAGACACCTGATTAGGCAGTTGATTACTGCGAAAGTGCAGCGCCGACAAGTTCACGGATAACGCCGGGATATGAATATTCTGGCTACGCCATTCTGCTAACTGACGGCA  >  minE/948860‑948991

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: