Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 597,053 T→C 46.6% T309T (ACT→ACC 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,0530TC46.6% 4.2 / 62.0 73T309T (ACT→ACCserSseryl‑tRNA synthetase, also charges selenocysteinyl‑tRNA with serine
Reads supporting (aligned to +/- strand):  ref base T (33/6);  new base C (15/19);  total (48/25)
Fisher's exact test for biased strand distribution p-value = 4.38e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAGTTGAAATGGTGCAGATCGTGCGCCCAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGCACTGGCGA  >  minE/596991‑597121
                                                              |                                                                    
aaaGTTGAAATGGTGCAGATCGTGCGCCCAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCa                                                                 >  1:2649528/1‑68 (MQ=255)
aaaGTTGAAATGGTGCAGATCGTGCGCCCAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCa                                                                 >  1:2594684/1‑68 (MQ=255)
aaaGTTGAAATGGTGCAGATCGTGCGCCCAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCa                                                                 >  1:1478627/1‑68 (MQ=255)
   gTTGAAATGGTGCAGATCGTGCGCCCAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGc                                                              >  1:3408077/1‑68 (MQ=255)
        aaTGGTGCAGATCGTGCGCCCGGAAGACTCAATGGCGGCGCTGGAAGAGATGACCGGTCATGCGGaaa                                                         <  1:2961236/68‑1 (MQ=255)
        aaTGGTGCAGATCGTGCGCCCGGAAGACTCAATGGCGGCGCTGGAAGAGATGACCGGTCATGCGGaaa                                                         <  1:2057391/68‑1 (MQ=255)
        aaTGGTGCAGATCGTGCGCCCGGAAGACTCAATGGCGGCGCTGGAAGAGATGACCGGTCATGCGGaaa                                                         <  1:1953708/68‑1 (MQ=255)
        aaTGGTGCAGATCGTGCGCCCGGAAGACTCAATGGCGGCGCTGGAAGAGATGACCGGTCATGCGGaaa                                                         <  1:2073164/68‑1 (MQ=255)
        aaTGGTGCAGATCGTGCGCCCGGAAGACTCAATGGCGGCGCTGGAAGAGATGACCGGTCATGCGGaaa                                                         <  1:2098/68‑1 (MQ=255)
        aaTGGTGCAGATCGTGCGCCCGGAAGACTCAATGGCGGCGCTGGAAGAGATGACCGGTCATGCGGaaa                                                         <  1:3305346/68‑1 (MQ=255)
        aaTGGTGCAGATCGTGCGCCCGGAAGACTCAATGGCGGCGCTGGAAGAGATGACCGGTCATGCGGaaa                                                         <  1:336563/68‑1 (MQ=255)
        aaTGGTGCAGATCGTGCGCCCGGAAGACTCAATGGCGGCGCTGGAAGAGATGACCGGTCATGCGGaaa                                                         <  1:1642137/68‑1 (MQ=255)
        aaTGGTGCAGATCGTGCGCCCGGAAGACTCAATGGCGGCGCTGGAAGAGATGACCGGTCATGCGGaaa                                                         <  1:2926376/68‑1 (MQ=255)
        aaTGGTGCAGATCGTGCGCCCGGAAGACTCAATGGCGGCGCTGGAAGAGATGACCGGTCATGCGGaaa                                                         <  1:1373341/68‑1 (MQ=255)
        aaTGGTGCAGATCGTGCGCCCGGAAGACTCAATGGCGGCGCTGGAAGAGATGACCGGTCATGCGGaaa                                                         <  1:645218/68‑1 (MQ=255)
        aaTGGTGCAGATCGTGCGCCCGGAAGACTCAATGGCGGCGCTGGAAGAGATGACCGGTCATGCGGaaa                                                         <  1:964825/68‑1 (MQ=255)
            gTGCAGATCGTGCGCCCAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTc                                                    >  1:2156999/1‑69 (MQ=255)
                    cGTGCGCCCAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGc                                                >  1:2774315/1‑65 (MQ=255)
                    cGTGCGCCCAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGc                                                >  1:3153444/1‑65 (MQ=255)
                    cGTGCGCCCAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGc                                                >  1:1046781/1‑65 (MQ=255)
                    cGTGCGCCCAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGc                                                >  1:2693629/1‑65 (MQ=255)
                    cGTGCGCCCAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGc                                                >  1:2836124/1‑65 (MQ=255)
                    cGTGCGCCCAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGc                                                >  1:1183030/1‑65 (MQ=255)
                    cGTGCGCCCAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGc                                                >  1:3219299/1‑65 (MQ=255)
                    cGTGCGCCCAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGc                                                >  1:1419537/1‑65 (MQ=255)
                    cGTGCGCCCAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGc                                                >  1:3396136/1‑65 (MQ=255)
                    cGTGCGCCCAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGc                                                >  1:163066/1‑65 (MQ=255)
                    cGTGCGCCCAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGtt                                            >  1:2671979/1‑69 (MQ=255)
                    cGTGCGCCCAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGtt                                            >  1:538749/1‑69 (MQ=255)
                    cGTGCGCCCAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGtt                                            >  1:452045/1‑69 (MQ=255)
                    cGTGCGCCCAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGtt                                            >  1:2935544/1‑69 (MQ=255)
                    cGTGCGCCCAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGtt                                            >  1:245689/1‑69 (MQ=255)
                          cccAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTggg                                      <  1:2038615/69‑1 (MQ=255)
                          cccAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTggg                                      <  1:2606129/69‑1 (MQ=255)
                          cccAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTggg                                      <  1:2490958/69‑1 (MQ=255)
                          cccAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTggg                                      <  1:2365301/69‑1 (MQ=255)
                          cccAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTggg                                      <  1:3179590/69‑1 (MQ=255)
                          cccAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGTCATGCAGAAAAAGTCCTGCAGTTGCTggg                                      <  1:2033566/68‑1 (MQ=255)
                               aagaCTCAATGGCGGCGCTGGAAGAGATGACCGGTCATGCGGAAAAAGTCCTGCAGCTGCTGGGCCTg                                  >  1:2188414/1‑68 (MQ=255)
                               aagaCTCAATGGCGGCGCTGGAAGAGATGACCGGTCATGCGGAAAAAGTCCTGCAGCTGCTGGGCCTGc                                 >  1:3428343/1‑69 (MQ=255)
                               aagaCTCAATGGCGGCGCTGGAAGAGATGACCGGTCATGCGGAAAAAGTCCTGCAGCTGCTGGGCCTGc                                 >  1:3299549/1‑69 (MQ=255)
                               aagaCTCAATGGCGGCGCTGGAAGAGATGACCGGTCATGCGGAAAAAGTCCTGCAGCTGCTGGGCCTGc                                 >  1:3189026/1‑69 (MQ=255)
                               aagaCTCAATGGCGGCGCTGGAAGAGATGACCGGTCATGCGGAAAAAGTCCTGCAGCTGCTGGGCCTGc                                 >  1:3184778/1‑69 (MQ=255)
                               aagaCTCAATGGCGGCGCTGGAAGAGATGACCGGTCATGCGGAAAAAGTCCTGCAGCTGCTGGGCCTGc                                 >  1:674390/1‑69 (MQ=255)
                               aagaCTCAATGGCGGCGCTGGAAGAGATGACCGGTCATGCGGAAAAAGTCCTGCAGCTGCTGGGCCTGc                                 >  1:902301/1‑69 (MQ=255)
                               aagaCTCAATGGCGGCGCTGGAAGAGATGACCGGTCATGCGGAAAAAGTCCTGCAGCTGCTGGGCCTGc                                 >  1:2530660/1‑69 (MQ=255)
                               aagaCTCAATGGCGGCGCTGGAAGAGATGACCGGTCATGCGGAAAAAGTCCTGCAGCTGCTGGGCCTGc                                 >  1:2543834/1‑69 (MQ=255)
                               aagaCTCAATGGCGGCGCTGGAAGAGATGACCGGTCATGCGGAAAAAGTCCTGCAGCTGCTGGGCCTGc                                 >  1:1606220/1‑69 (MQ=255)
                               aagaCTCAATGGCGGCGCTGGAAGAGATGACCGGTCATGCGGAAAAAGTCCTGCAGCTGCTGGGCCTGc                                 >  1:1765082/1‑69 (MQ=255)
                               aagaCTCAATGGCGGCGCTGGAAGAGATGACCGGTCATGCGGAAAAAGTCCTGCAGCTGCTGGGCCTGc                                 >  1:1786972/1‑69 (MQ=255)
                               aagaCTCAATGGCGGCGCTGGAAGAGATGACCGGTCATGCGGAAAAAGTCCTGCAGCTGCTGGGCCTGc                                 >  1:1182704/1‑69 (MQ=255)
                               aagaCTCAATGGCGGCGCTGGAAGAGATGACCGGTCATGCGGAAAAAGTCCTGCAGCTGCTGGGCCTGc                                 >  1:2492564/1‑69 (MQ=255)
                               aagaCTCAATGGCGGCGCTGGAAGAGATGACCGGTCATGCGGAAAAAGTCCTGCAGCTGCTAGGCCTGc                                 >  1:3361426/1‑69 (MQ=255)
                                    tCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGccgtac                            >  1:624058/1‑69 (MQ=255)
                                    tCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGccgtac                            >  1:1758087/1‑69 (MQ=255)
                                    tCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGccgtac                            >  1:381335/1‑69 (MQ=255)
                                         ggcggcGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTaa                       >  1:1996801/1‑69 (MQ=255)
                                          gcggcgCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTaaa                      >  1:2054611/1‑69 (MQ=255)
                                                      gagaTGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGc          >  1:661577/1‑69 (MQ=255)
                                                      gagaTGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGc          >  1:1169369/1‑69 (MQ=255)
                                                      gagaTGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGc          >  1:586642/1‑69 (MQ=255)
                                                      gagaTGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGc          >  1:1435293/1‑69 (MQ=255)
                                                      gagaTGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGc          >  1:1485938/1‑69 (MQ=255)
                                                      gagaTGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGc          >  1:1695906/1‑69 (MQ=255)
                                                      gagaTGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGc          >  1:1486179/1‑69 (MQ=255)
                                                      gagaTGACTGGTCATGCAGAAAAAGTCATGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGc          >  1:2213891/1‑69 (MQ=255)
                                                      gagaTGACCGGTCATGCGGAAAAAGTCCTGCAGCTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTgt          <  1:1660813/69‑2 (MQ=255)
                                                      gagaTGACCGGTCATGCGGAAAAAGTCCTGCAGCTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTgt          <  1:2726769/69‑2 (MQ=255)
                                                      gagaTGACCGGTCATGCGGAAAAAGTCCTGCAGCTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTgt          <  1:3266827/69‑2 (MQ=255)
                                                      gagaTGACCGGTCATGCGGAAAAAGTCCTGCAGCTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTgt          <  1:2224994/69‑2 (MQ=255)
                                                      gagaTGACCGGTCATGCGGAAAAAGTCCTGCAGCTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTgt          <  1:3060016/69‑2 (MQ=255)
                                                      gagaTGACCGGTCATGCGGAAAAAGTCCTGCAGCTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTgt          <  1:2373031/69‑2 (MQ=255)
                                                      gagaTGACCGGTCATGCGGAAAAAGTCCTGCAGCTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTgt          <  1:1115666/69‑2 (MQ=255)
                                                              tGGTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGCACTGGCGa  <  1:1228971/69‑1 (MQ=255)
                                                              |                                                                    
AAAGTTGAAATGGTGCAGATCGTGCGCCCAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGCACTGGCGA  >  minE/596991‑597121

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: