Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 597,078 T→C 33.7% L318L (TTG→CTG)  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,0780TC33.7% 48.5 / 35.7 68L318L (TTG→CTG) serSseryl‑tRNA synthetase, also charges selenocysteinyl‑tRNA with serine
Reads supporting (aligned to +/- strand):  ref base T (21/24);  new base C (15/8);  total (36/32)
Fisher's exact test for biased strand distribution p-value = 2.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTGCGCCCAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGCACTGGCGACATGGGCTTTGGCGCTTGCAAA  >  minE/597011‑597143
                                                                   |                                                                 
cGTGCGCCCAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGtt                                                                  >  1:2671979/1‑69 (MQ=255)
cGTGCGCCCAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGtt                                                                  >  1:2935544/1‑69 (MQ=255)
cGTGCGCCCAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGtt                                                                  >  1:452045/1‑69 (MQ=255)
cGTGCGCCCAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGtt                                                                  >  1:538749/1‑69 (MQ=255)
cGTGCGCCCAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGtt                                                                  >  1:245689/1‑69 (MQ=255)
      cccAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTggg                                                            <  1:3179590/69‑1 (MQ=255)
      cccAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTggg                                                            <  1:2038615/69‑1 (MQ=255)
      cccAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTggg                                                            <  1:2365301/69‑1 (MQ=255)
      cccAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTggg                                                            <  1:2490958/69‑1 (MQ=255)
      cccAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTggg                                                            <  1:2606129/69‑1 (MQ=255)
      cccAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGTCATGCAGAAAAAGTCCTGCAGTTGCTggg                                                            <  1:2033566/68‑1 (MQ=255)
           aagaCTCAATGGCGGCGCTGGAAGAGATGACCGGTCATGCGGAAAAAGTCCTGCAGCTGCTGGGCCTg                                                        >  1:2188414/1‑68 (MQ=255)
           aagaCTCAATGGCGGCGCTGGAAGAGATGACCGGTCATGCGGAAAAAGTCCTGCAGCTGCTGGGCCTGc                                                       >  1:902301/1‑69 (MQ=255)
           aagaCTCAATGGCGGCGCTGGAAGAGATGACCGGTCATGCGGAAAAAGTCCTGCAGCTGCTGGGCCTGc                                                       >  1:2492564/1‑69 (MQ=255)
           aagaCTCAATGGCGGCGCTGGAAGAGATGACCGGTCATGCGGAAAAAGTCCTGCAGCTGCTGGGCCTGc                                                       >  1:3184778/1‑69 (MQ=255)
           aagaCTCAATGGCGGCGCTGGAAGAGATGACCGGTCATGCGGAAAAAGTCCTGCAGCTGCTGGGCCTGc                                                       >  1:2530660/1‑69 (MQ=255)
           aagaCTCAATGGCGGCGCTGGAAGAGATGACCGGTCATGCGGAAAAAGTCCTGCAGCTGCTGGGCCTGc                                                       >  1:2543834/1‑69 (MQ=255)
           aagaCTCAATGGCGGCGCTGGAAGAGATGACCGGTCATGCGGAAAAAGTCCTGCAGCTGCTGGGCCTGc                                                       >  1:3189026/1‑69 (MQ=255)
           aagaCTCAATGGCGGCGCTGGAAGAGATGACCGGTCATGCGGAAAAAGTCCTGCAGCTGCTGGGCCTGc                                                       >  1:3299549/1‑69 (MQ=255)
           aagaCTCAATGGCGGCGCTGGAAGAGATGACCGGTCATGCGGAAAAAGTCCTGCAGCTGCTGGGCCTGc                                                       >  1:3428343/1‑69 (MQ=255)
           aagaCTCAATGGCGGCGCTGGAAGAGATGACCGGTCATGCGGAAAAAGTCCTGCAGCTGCTGGGCCTGc                                                       >  1:1786972/1‑69 (MQ=255)
           aagaCTCAATGGCGGCGCTGGAAGAGATGACCGGTCATGCGGAAAAAGTCCTGCAGCTGCTGGGCCTGc                                                       >  1:1182704/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)
           aagaCTCAATGGCGGCGCTGGAAGAGATGACCGGTCATGCGGAAAAAGTCCTGCAGCTGCTAGGCCTGc                                                       >  1:3361426/1‑69 (MQ=255)
                tCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGccgtac                                                  >  1:381335/1‑69 (MQ=255)
                tCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGccgtac                                                  >  1:1758087/1‑69 (MQ=255)
                tCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGccgtac                                                  >  1:624058/1‑69 (MQ=255)
                     ggcggcGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTaa                                             >  1:1996801/1‑69 (MQ=255)
                      gcggcgCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTaaa                                            >  1:2054611/1‑69 (MQ=255)
                                  gagaTGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGc                                >  1:1485938/1‑69 (MQ=255)
                                  gagaTGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGc                                >  1:1169369/1‑69 (MQ=255)
                                  gagaTGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGc                                >  1:661577/1‑69 (MQ=255)
                                  gagaTGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGc                                >  1:1435293/1‑69 (MQ=255)
                                  gagaTGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGc                                >  1:1486179/1‑69 (MQ=255)
                                  gagaTGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGc                                >  1:586642/1‑69 (MQ=255)
                                  gagaTGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGc                                >  1:1695906/1‑69 (MQ=255)
                                  gagaTGACTGGTCATGCAGAAAAAGTCATGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGc                                >  1:2213891/1‑69 (MQ=255)
                                  gagaTGACCGGTCATGCGGAAAAAGTCCTGCAGCTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTgt                                <  1:3266827/69‑2 (MQ=255)
                                  gagaTGACCGGTCATGCGGAAAAAGTCCTGCAGCTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTgt                                <  1:3060016/69‑2 (MQ=255)
                                  gagaTGACCGGTCATGCGGAAAAAGTCCTGCAGCTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTgt                                <  1:1115666/69‑2 (MQ=255)
                                  gagaTGACCGGTCATGCGGAAAAAGTCCTGCAGCTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTgt                                <  1:2726769/69‑2 (MQ=255)
                                  gagaTGACCGGTCATGCGGAAAAAGTCCTGCAGCTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTgt                                <  1:2373031/69‑2 (MQ=255)
                                  gagaTGACCGGTCATGCGGAAAAAGTCCTGCAGCTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTgt                                <  1:1660813/69‑2 (MQ=255)
                                  gagaTGACCGGTCATGCGGAAAAAGTCCTGCAGCTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTgt                                <  1:2224994/69‑2 (MQ=255)
                                          tGGTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGCACTGGCGa                        <  1:1228971/69‑1 (MQ=255)
                                          cggTCATGCGGAAAAAGTCCTGCAGCTGCTGGGCCTGCCGTACCGTAAAATCATGCTTTgt                                <  1:372160/60‑2 (MQ=255)
                                           ggTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGCACTGGCGAc                       <  1:2674514/69‑1 (MQ=255)
                                           ggTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGCACTGGCGAc                       <  1:997984/69‑1 (MQ=255)
                                           ggTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGCACTGGCGAc                       <  1:789595/69‑1 (MQ=255)
                                           ggTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGCACTGGCGAc                       <  1:713607/69‑1 (MQ=255)
                                           ggTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGCACTGGCGAc                       <  1:136985/69‑1 (MQ=255)
                                           ggTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGCACTGGCGAc                       <  1:2430499/69‑1 (MQ=255)
                                           ggTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGCACTGGCGAc                       <  1:2305991/69‑1 (MQ=255)
                                            gTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAatca                                        <  1:945173/51‑1 (MQ=255)
                                            gTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAatca                                        <  1:1066129/51‑1 (MQ=255)
                                            gTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAatca                                        <  1:1806073/51‑1 (MQ=255)
                                            gTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAatca                                        <  1:3260277/51‑1 (MQ=255)
                                            gTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAatca                                        <  1:1728658/51‑1 (MQ=255)
                                            gTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAatca                                        <  1:696814/51‑1 (MQ=255)
                                            gTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAatca                                        <  1:470000/51‑1 (MQ=255)
                                            gTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAatca                                        <  1:3067264/51‑1 (MQ=255)
                                            gTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGCACCGTAAAatca                                        <  1:1412614/51‑1 (MQ=255)
                                                  cAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTaaa                                            >  1:2688980/1‑41 (MQ=255)
                                                  cAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGCACTGGCGACATGGGCt                >  1:1673981/1‑69 (MQ=255)
                                                  cAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGCACTGGCGACATGGGCt                >  1:528350/1‑69 (MQ=255)
                                                  cAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGCACTGGCGACATGGGCt                >  1:1685174/1‑69 (MQ=255)
                                                  cAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGCACTGGCGACATGGGCt                >  1:29032/1‑69 (MQ=255)
                                                  cAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGCACTGGCGACATGGGCt                >  1:182658/1‑69 (MQ=255)
                                                  cAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGCACTGGCGACATGGGCt                >  1:1949306/1‑69 (MQ=255)
                                                  cAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGCACTGGCGACATGGGCt                >  1:2163121/1‑69 (MQ=255)
                                                                  gTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGCACTGGCGACATGGGCTTTGGCGCTTGCaaa  <  1:2590995/67‑1 (MQ=255)
                                                                   |                                                                 
CGTGCGCCCAGAAGACTCAATGGCGGCGCTGGAAGAGATGACTGGTCATGCAGAAAAAGTCCTGCAGTTGCTGGGCCTGCCGTACCGTAAAATCATCCTTTGCACTGGCGACATGGGCTTTGGCGCTTGCAAA  >  minE/597011‑597143

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: