Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,031,586 G→A 43.0% S232S (TCC→TCT cpdA ← cyclic 3',5'‑adenosine monophosphate phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,031,5860GA43.0% 24.8 / 58.2 79S232S (TCC→TCTcpdAcyclic 3',5'‑adenosine monophosphate phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base G (15/30);  new base A (14/20);  total (29/50)
Fisher's exact test for biased strand distribution p-value = 4.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.96e-01

CATGTAACTCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCA  >  minE/2031534‑2031648
                                                    |                                                              
cATGTAACTCTAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTa                                                 <  1:131271/68‑1 (MQ=255)
cATGTAACTCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTa                                                 <  1:1042988/68‑1 (MQ=255)
cATGTAACTCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTa                                                 <  1:2861100/68‑1 (MQ=255)
cATGTAACTCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTa                                                 <  1:2533831/68‑1 (MQ=255)
cATGTAACTCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTa                                                 <  1:2387141/68‑1 (MQ=255)
cATGTAACTCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTa                                                 <  1:2052716/68‑1 (MQ=255)
cATGTAACTCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTa                                                 <  1:2050287/68‑1 (MQ=255)
cATGTAACTCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTa                                                 <  1:1066567/68‑1 (MQ=255)
 aTGTAACTCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTa                                                 <  1:2334101/67‑1 (MQ=255)
  tGTAACTCGAGAGTACGCCAGCCGGGCTCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAAc                                              <  1:1391765/69‑1 (MQ=255)
  tGTAACTCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAAc                                              <  1:1325663/69‑1 (MQ=255)
  tGTAACTCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAAc                                              <  1:1493668/69‑1 (MQ=255)
  tGTAACTCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAAc                                              <  1:391528/69‑1 (MQ=255)
  tGTAACTCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAAc                                              <  1:854737/69‑1 (MQ=255)
  tGTAACTCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAAc                                              <  1:1921702/69‑1 (MQ=255)
    tAACTCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTg                                            >  1:2556940/1‑69 (MQ=255)
    tAACTCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTg                                            >  1:241112/1‑69 (MQ=255)
    tAACTCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTg                                            >  1:2325333/1‑69 (MQ=255)
    tAACTCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTg                                            >  1:2632713/1‑69 (MQ=255)
    tAACTCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTg                                            >  1:3376192/1‑69 (MQ=255)
    tAACTCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTg                                            >  1:707393/1‑69 (MQ=255)
    tAACTCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTg                                            >  1:2049450/1‑69 (MQ=255)
    tAACTCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTg                                            >  1:1849149/1‑69 (MQ=255)
    tAACTCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTg                                            >  1:1017379/1‑69 (MQ=255)
    tAACTCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTg                                            >  1:3155201/1‑69 (MQ=255)
    tAACTCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTg                                            >  1:1642368/1‑69 (MQ=255)
    tAACTCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTg                                            >  1:2775761/1‑69 (MQ=255)
    tAACTCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTg                                            >  1:3105886/1‑69 (MQ=255)
    tAACTCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTg                                            >  1:1542318/1‑69 (MQ=255)
                 cGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGAc                               >  1:2626624/1‑69 (MQ=255)
                   ccAGCCGGGCGCGATGGTGTCCAGCGTAAAGTTAGAACAGTGCGGCTTAAACTGCACACAGGTCGACg                              <  1:2655059/68‑1 (MQ=255)
                   ccAGCCGGGCGCGATGGTGTCCAGCGTAAAGTTAGAACAGTGCGGCTTAAACTGCACACAGGTCGACg                              <  1:941262/68‑1 (MQ=255)
                   ccAGCCGGGCGCGATGGTGTCCAGCGTAAAGTTAGAACAGTGCGGCTTAAACTGCACACAGGTCGACg                              <  1:3186463/68‑1 (MQ=255)
                   ccAGCCGGGCGCGATGGTGTCCAGCGTAAAGTTAGAACAGTGCGGCTTAAACTGCACACAGGTCGACg                              <  1:2568227/68‑1 (MQ=255)
                   ccAGCCGGGCGCGATGGTGTCCAGCGTAAAGTTAGAACAGTGCGGCTTAAACTGCACACAGGTCGACg                              <  1:305473/68‑1 (MQ=255)
                   ccAGCCGGGCGCGATGGTGTCCAGCGTAAAGTTAGAACAGTGCGGCTTAAACTGCACACAGGTCGACg                              <  1:2643034/68‑1 (MQ=255)
                   ccAGCCGGGCGCGATGGTGTCCAGCGTAAAGTTAGAACAGTGCGGCTTAAACTGCACACAGGTCGACg                              <  1:2717331/68‑1 (MQ=255)
                   ccAGCCGGGCGCGATGGTGTCCAGCGTAAAGTTAGAACAGTGCGGCTTAAACTGCACACAGGTCGACg                              <  1:2738161/68‑1 (MQ=255)
                   ccAGCCGGGCGCGATGGTGTCCAGCGTAAAGTTAGAACAGTGCGGCTTAAACTGCACACAGGTCGACg                              <  1:2911755/68‑1 (MQ=255)
                   ccAGCCGGGCGCGATGGTGTCCAGCGTAAAGTTAGAACAGTGCGGCTTAAACTGCACACAGGTCGACg                              <  1:3328728/68‑1 (MQ=255)
                   ccAGCCGGGCGCGATGGTGTCCAGCGTAAAGTTAGAACAGTGCGGCTTAAACTGCACACAGGTCGACg                              <  1:3349640/68‑1 (MQ=255)
                   ccAGCCGGGCGCGATGGTGTCCAGCGTAAAGTTAGAACAGTGCGGCTTAAACTGCACACAGGTCGACg                              <  1:1169891/68‑1 (MQ=255)
                   ccAGCCGGGCGCGATGGTGTCCAGCGTAAAGTTAGAACAGTGCGGCTTAAACTGCACACAGGTCGACg                              <  1:2226240/68‑1 (MQ=255)
                   ccAGCCGGGCGCGATGGTGTCCAGCGTAAAGTTAGAACAGTGCGGCTTAAACTGCACACAGGTCGACg                              <  1:349343/68‑1 (MQ=255)
                   ccAGCCGGGCGCGATGGTGTCCAGCGTAAAGTTAGAACAGTGCGGCTTAAACTGCACACAGGTCGACg                              <  1:2210564/68‑1 (MQ=255)
                   ccAGCCGGGCGCGATGGTGTCCAGCGTAAAGTTAGAACAGTGCGGCTTAAACTGCACACAGGTCGACg                              <  1:2173316/68‑1 (MQ=255)
                   ccAGCCGGGCGCGATGGTGTCCAGCGTAAAGTTAGAACAGTGCGGCTTAAACTGCACACAGGTCGACg                              <  1:2055198/68‑1 (MQ=255)
                   ccAGCCGGGCGCGATGGTGTCCAGCGTAAAGTTAGAACAGTGCGGCTTAAACTGCACACAGGTCGACg                              <  1:1889881/68‑1 (MQ=255)
                   ccAGCCGGGCGCGATGGTGTCCAGCGTAAAGTTAGAACAGTGCGGCTTAAACTGCACACAGGTCGACg                              <  1:1490288/68‑1 (MQ=255)
                                   gtgtCCAGCGTAAAGTTAGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAagcg              >  1:2348941/4‑68 (MQ=255)
                                   gtgtCCAGCGTAAAGTTAGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAagcg              >  1:3187425/4‑68 (MQ=255)
                                   gtgtCCAGCGTAAAGTTAGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAagcg              >  1:1350919/4‑68 (MQ=255)
                                   gtgtCCAGCGTAAAGTTAGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAagcg              >  1:1595195/4‑68 (MQ=255)
                                   gtgtCCAGCGTAAAGTTAGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAagcg              >  1:3029936/4‑68 (MQ=255)
                                   gtgtCCAGCGTAAAGTTAGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAagcg              >  1:170311/4‑68 (MQ=255)
                                   gtgtCCAGCGTAAAGTTAGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAagcg              >  1:2868406/4‑68 (MQ=255)
                                   gtgtCCAGCGTAAAGTTAGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAagcg              >  1:869512/4‑68 (MQ=255)
                                   gtgtCCAGCGTAAAGTTAGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAagcg              >  1:2855675/4‑68 (MQ=255)
                                   gtgtCCAGCGTAAAGTTAGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAagcg              >  1:28197/4‑68 (MQ=255)
                                   gtgtCCAGCGTAAAGTTAGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAagcg              >  1:1817335/4‑68 (MQ=255)
                                   gtgtCCAGCGTAAAGTTAGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAagcg              >  1:1856776/4‑68 (MQ=255)
                                   gtgtCCAGCGTAAAGTTAGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAagcg              >  1:273016/4‑68 (MQ=255)
                                   gtgtCCAGCGTAAAGTTAGAACAGTGCGGCTTAAACTACACACAGGTCGACGGCGTTGCCAGCAagcg              >  1:286766/4‑68 (MQ=255)
                                     aTCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAggcggc            <  1:32399/68‑1 (MQ=255)
                                     aTCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAggcggc            <  1:477445/68‑1 (MQ=255)
                                     aTCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAggcggc            <  1:98181/68‑1 (MQ=255)
                                     aTCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAggcggc            <  1:2830018/68‑1 (MQ=255)
                                     aTCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAggcggc            <  1:2773055/68‑1 (MQ=255)
                                     aTCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAggcggc            <  1:2768861/68‑1 (MQ=255)
                                     aTCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAggcggc            <  1:1171559/68‑1 (MQ=255)
                                     aTCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAggcggc            <  1:1777807/68‑1 (MQ=255)
                                        cAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGAc         <  1:1186151/68‑1 (MQ=255)
                                        cAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGAc         <  1:1757410/68‑1 (MQ=255)
                                        cAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGAc         <  1:1194504/68‑1 (MQ=255)
                                         aGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGAc         <  1:2210872/67‑1 (MQ=255)
                                           cGTAAAGTTAGAACAGTGCGGCTTAAACTGCACACAGGTCGACg                              <  1:281481/44‑1 (MQ=255)
                                              aaaGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCa  <  1:3160199/69‑1 (MQ=255)
                                              aaaGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCa  <  1:2052669/69‑1 (MQ=255)
                                               aaGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCa  <  1:1207586/68‑1 (MQ=255)
                                                    |                                                              
CATGTAACTCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCA  >  minE/2031534‑2031648

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: