Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 597,062 A→G 36.8% A312A (GCA→GCG serS → seryl‑tRNA synthetase, also charges selenocysteinyl‑tRNA with serine

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE597,0620AG36.8% 56.6 / 59.1 95A312A (GCA→GCGserSseryl‑tRNA synthetase, also charges selenocysteinyl‑tRNA with serine
Reads supporting (aligned to +/- strand):  ref base A (37/23);  new base G (15/20);  total (52/43)
Fisher's exact test for biased strand distribution p-value = 9.01e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.44e-01

AATGGTGCAGATCGTGCGCCCAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGCACTGGCGACATGGGCT  >  minE/596999‑597129
                                                               |                                                                   
aaTGGTGCAGATCGTGCGCCCGGAAGACTCAATGGCGGCGCTGGAAGAGATGACCGGTCATGCGGaaa                                                                 <  1:3305346/68‑1 (MQ=255)
aaTGGTGCAGATCGTGCGCCCGGAAGACTCAATGGCGGCGCTGGAAGAGATGACCGGTCATGCGGaaa                                                                 <  1:2098/68‑1 (MQ=255)
aaTGGTGCAGATCGTGCGCCCGGAAGACTCAATGGCGGCGCTGGAAGAGATGACCGGTCATGCGGaaa                                                                 <  1:2073164/68‑1 (MQ=255)
aaTGGTGCAGATCGTGCGCCCGGAAGACTCAATGGCGGCGCTGGAAGAGATGACCGGTCATGCGGaaa                                                                 <  1:336563/68‑1 (MQ=255)
aaTGGTGCAGATCGTGCGCCCGGAAGACTCAATGGCGGCGCTGGAAGAGATGACCGGTCATGCGGaaa                                                                 <  1:964825/68‑1 (MQ=255)
aaTGGTGCAGATCGTGCGCCCGGAAGACTCAATGGCGGCGCTGGAAGAGATGACCGGTCATGCGGaaa                                                                 <  1:1373341/68‑1 (MQ=255)
aaTGGTGCAGATCGTGCGCCCGGAAGACTCAATGGCGGCGCTGGAAGAGATGACCGGTCATGCGGaaa                                                                 <  1:2057391/68‑1 (MQ=255)
aaTGGTGCAGATCGTGCGCCCGGAAGACTCAATGGCGGCGCTGGAAGAGATGACCGGTCATGCGGaaa                                                                 <  1:1642137/68‑1 (MQ=255)
aaTGGTGCAGATCGTGCGCCCGGAAGACTCAATGGCGGCGCTGGAAGAGATGACCGGTCATGCGGaaa                                                                 <  1:2961236/68‑1 (MQ=255)
aaTGGTGCAGATCGTGCGCCCGGAAGACTCAATGGCGGCGCTGGAAGAGATGACCGGTCATGCGGaaa                                                                 <  1:1953708/68‑1 (MQ=255)
aaTGGTGCAGATCGTGCGCCCGGAAGACTCAATGGCGGCGCTGGAAGAGATGACCGGTCATGCGGaaa                                                                 <  1:645218/68‑1 (MQ=255)
aaTGGTGCAGATCGTGCGCCCGGAAGACTCAATGGCGGCGCTGGAAGAGATGACCGGTCATGCGGaaa                                                                 <  1:2926376/68‑1 (MQ=255)
    gTGCAGATCGTGCGCCCAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTc                                                            >  1:2156999/1‑69 (MQ=255)
            cGTGCGCCCAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGc                                                        >  1:3153444/1‑65 (MQ=255)
            cGTGCGCCCAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGc                                                        >  1:3396136/1‑65 (MQ=255)
            cGTGCGCCCAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGc                                                        >  1:2836124/1‑65 (MQ=255)
            cGTGCGCCCAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGc                                                        >  1:2774315/1‑65 (MQ=255)
            cGTGCGCCCAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGc                                                        >  1:2693629/1‑65 (MQ=255)
            cGTGCGCCCAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGc                                                        >  1:1046781/1‑65 (MQ=255)
            cGTGCGCCCAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGc                                                        >  1:163066/1‑65 (MQ=255)
            cGTGCGCCCAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGc                                                        >  1:1183030/1‑65 (MQ=255)
            cGTGCGCCCAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGc                                                        >  1:1419537/1‑65 (MQ=255)
            cGTGCGCCCAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGc                                                        >  1:3219299/1‑65 (MQ=255)
            cGTGCGCCCAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGtt                                                    >  1:452045/1‑69 (MQ=255)
            cGTGCGCCCAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGtt                                                    >  1:245689/1‑69 (MQ=255)
            cGTGCGCCCAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGtt                                                    >  1:2671979/1‑69 (MQ=255)
            cGTGCGCCCAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGtt                                                    >  1:2935544/1‑69 (MQ=255)
            cGTGCGCCCAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGtt                                                    >  1:538749/1‑69 (MQ=255)
                  cccAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTggg                                              <  1:2365301/69‑1 (MQ=255)
                  cccAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTggg                                              <  1:3179590/69‑1 (MQ=255)
                  cccAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTggg                                              <  1:2038615/69‑1 (MQ=255)
                  cccAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTggg                                              <  1:2606129/69‑1 (MQ=255)
                  cccAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTggg                                              <  1:2490958/69‑1 (MQ=255)
                  cccAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGTCATGCAGAAAAAGTCCTGCAGTTGCTggg                                              <  1:2033566/68‑1 (MQ=255)
                       aagaCTCAATGGCGGCGCTGGAAGAGATGACCGGTCATGCGGAAAAAGTCCTGCAGCTGCTGGGCCTg                                          >  1:2188414/1‑68 (MQ=255)
                       aagaCTCAATGGCGGCGCTGGAAGAGATGACCGGTCATGCGGAAAAAGTCCTGCAGCTGCTGGGCCTGc                                         >  1:2543834/1‑69 (MQ=255)
                       aagaCTCAATGGCGGCGCTGGAAGAGATGACCGGTCATGCGGAAAAAGTCCTGCAGCTGCTGGGCCTGc                                         >  1:2492564/1‑69 (MQ=255)
                       aagaCTCAATGGCGGCGCTGGAAGAGATGACCGGTCATGCGGAAAAAGTCCTGCAGCTGCTGGGCCTGc                                         >  1:2530660/1‑69 (MQ=255)
                       aagaCTCAATGGCGGCGCTGGAAGAGATGACCGGTCATGCGGAAAAAGTCCTGCAGCTGCTGGGCCTGc                                         >  1:3428343/1‑69 (MQ=255)
                       aagaCTCAATGGCGGCGCTGGAAGAGATGACCGGTCATGCGGAAAAAGTCCTGCAGCTGCTGGGCCTGc                                         >  1:1786972/1‑69 (MQ=255)
                       aagaCTCAATGGCGGCGCTGGAAGAGATGACCGGTCATGCGGAAAAAGTCCTGCAGCTGCTGGGCCTGc                                         >  1:1765082/1‑69 (MQ=255)
                       aagaCTCAATGGCGGCGCTGGAAGAGATGACCGGTCATGCGGAAAAAGTCCTGCAGCTGCTGGGCCTGc                                         >  1:674390/1‑69 (MQ=255)
                       aagaCTCAATGGCGGCGCTGGAAGAGATGACCGGTCATGCGGAAAAAGTCCTGCAGCTGCTGGGCCTGc                                         >  1:1606220/1‑69 (MQ=255)
                       aagaCTCAATGGCGGCGCTGGAAGAGATGACCGGTCATGCGGAAAAAGTCCTGCAGCTGCTGGGCCTGc                                         >  1:3184778/1‑69 (MQ=255)
                       aagaCTCAATGGCGGCGCTGGAAGAGATGACCGGTCATGCGGAAAAAGTCCTGCAGCTGCTGGGCCTGc                                         >  1:3189026/1‑69 (MQ=255)
                       aagaCTCAATGGCGGCGCTGGAAGAGATGACCGGTCATGCGGAAAAAGTCCTGCAGCTGCTGGGCCTGc                                         >  1:902301/1‑69 (MQ=255)
                       aagaCTCAATGGCGGCGCTGGAAGAGATGACCGGTCATGCGGAAAAAGTCCTGCAGCTGCTGGGCCTGc                                         >  1:1182704/1‑69 (MQ=255)
                       aagaCTCAATGGCGGCGCTGGAAGAGATGACCGGTCATGCGGAAAAAGTCCTGCAGCTGCTGGGCCTGc                                         >  1:3299549/1‑69 (MQ=255)
                       aagaCTCAATGGCGGCGCTGGAAGAGATGACCGGTCATGCGGAAAAAGTCCTGCAGCTGCTAGGCCTGc                                         >  1:3361426/1‑69 (MQ=255)
                            tCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGccgtac                                    >  1:381335/1‑69 (MQ=255)
                            tCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGccgtac                                    >  1:624058/1‑69 (MQ=255)
                            tCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGccgtac                                    >  1:1758087/1‑69 (MQ=255)
                                 ggcggcGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTaa                               >  1:1996801/1‑69 (MQ=255)
                                  gcggcgCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTaaa                              >  1:2054611/1‑69 (MQ=255)
                                              gagaTGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGc                  >  1:1169369/1‑69 (MQ=255)
                                              gagaTGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGc                  >  1:586642/1‑69 (MQ=255)
                                              gagaTGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGc                  >  1:661577/1‑69 (MQ=255)
                                              gagaTGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGc                  >  1:1695906/1‑69 (MQ=255)
                                              gagaTGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGc                  >  1:1435293/1‑69 (MQ=255)
                                              gagaTGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGc                  >  1:1485938/1‑69 (MQ=255)
                                              gagaTGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGc                  >  1:1486179/1‑69 (MQ=255)
                                              gagaTGACTGGTCATGCAGAAAAAGTCATGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGc                  >  1:2213891/1‑69 (MQ=255)
                                              gagaTGACCGGTCATGCGGAAAAAGTCCTGCAGCTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTgt                  <  1:2224994/69‑2 (MQ=255)
                                              gagaTGACCGGTCATGCGGAAAAAGTCCTGCAGCTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTgt                  <  1:1115666/69‑2 (MQ=255)
                                              gagaTGACCGGTCATGCGGAAAAAGTCCTGCAGCTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTgt                  <  1:3266827/69‑2 (MQ=255)
                                              gagaTGACCGGTCATGCGGAAAAAGTCCTGCAGCTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTgt                  <  1:1660813/69‑2 (MQ=255)
                                              gagaTGACCGGTCATGCGGAAAAAGTCCTGCAGCTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTgt                  <  1:3060016/69‑2 (MQ=255)
                                              gagaTGACCGGTCATGCGGAAAAAGTCCTGCAGCTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTgt                  <  1:2726769/69‑2 (MQ=255)
                                              gagaTGACCGGTCATGCGGAAAAAGTCCTGCAGCTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTgt                  <  1:2373031/69‑2 (MQ=255)
                                                      tGGTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGCACTGGCGa          <  1:1228971/69‑1 (MQ=255)
                                                      cggTCATGCGGAAAAAGTCCTGCAGCTGCTGGGCCTGCCGTACCGTAAAATCATGCTTTgt                  <  1:372160/60‑2 (MQ=255)
                                                       ggTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGCACTGGCGAc         <  1:789595/69‑1 (MQ=255)
                                                       ggTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGCACTGGCGAc         <  1:713607/69‑1 (MQ=255)
                                                       ggTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGCACTGGCGAc         <  1:997984/69‑1 (MQ=255)
                                                       ggTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGCACTGGCGAc         <  1:136985/69‑1 (MQ=255)
                                                       ggTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGCACTGGCGAc         <  1:2674514/69‑1 (MQ=255)
                                                       ggTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGCACTGGCGAc         <  1:2430499/69‑1 (MQ=255)
                                                       ggTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGCACTGGCGAc         <  1:2305991/69‑1 (MQ=255)
                                                        gTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAatca                          <  1:3067264/51‑1 (MQ=255)
                                                        gTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAatca                          <  1:1066129/51‑1 (MQ=255)
                                                        gTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAatca                          <  1:945173/51‑1 (MQ=255)
                                                        gTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAatca                          <  1:696814/51‑1 (MQ=255)
                                                        gTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAatca                          <  1:1728658/51‑1 (MQ=255)
                                                        gTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAatca                          <  1:1806073/51‑1 (MQ=255)
                                                        gTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAatca                          <  1:470000/51‑1 (MQ=255)
                                                        gTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAatca                          <  1:3260277/51‑1 (MQ=255)
                                                        gTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGCACCGTAAAatca                          <  1:1412614/51‑1 (MQ=255)
                                                              cAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTaaa                              >  1:2688980/1‑41 (MQ=255)
                                                              cAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGCACTGGCGACATGGGCt  >  1:1949306/1‑69 (MQ=255)
                                                              cAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGCACTGGCGACATGGGCt  >  1:182658/1‑69 (MQ=255)
                                                              cAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGCACTGGCGACATGGGCt  >  1:2163121/1‑69 (MQ=255)
                                                              cAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGCACTGGCGACATGGGCt  >  1:528350/1‑69 (MQ=255)
                                                              cAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGCACTGGCGACATGGGCt  >  1:1685174/1‑69 (MQ=255)
                                                              cAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGCACTGGCGACATGGGCt  >  1:1673981/1‑69 (MQ=255)
                                                              cAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGCACTGGCGACATGGGCt  >  1:29032/1‑69 (MQ=255)
                                                               |                                                                   
AATGGTGCAGATCGTGCGCCCAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGCACTGGCGACATGGGCT  >  minE/596999‑597129

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: