Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 354,489 C→T 40.7% S268S (TCC→TCT cysS → cysteinyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE354,4890CT40.7% 33.1 / 51.7 76S268S (TCC→TCTcysScysteinyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base C (12/33);  new base T (19/12);  total (31/45)
Fisher's exact test for biased strand distribution p-value = 4.11e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.24e-01

ATGGTCAGTATGTGAACTACTGGATGCACTCGGGGATGGTGATGGTTGACCGCGAGAAGATGTCCAAATCGCTGGGTAACTTCTTTACCGTGCGCGATGTGCTGAAATACTACGACGCGGAAACCGTGCGTT  >  minE/354425‑354556
                                                                |                                                                   
atgGTCAGTATGTGAACTACTGGATGCACTCGGGGATGGTGATGGTTGACCGCGAGAAGATGTCCAAAt                                                                 >  1:19083/1‑69 (MQ=255)
atgGTCAGTATGTGAACTACTGGATGCACTCGGGGATGGTGATGGTTGACCGCGAGAAGATGTCCAAAt                                                                 >  1:1638398/1‑69 (MQ=255)
atgGTCAGTATGTGAACTACTGGATGCACTCGGGGATGGTGATGGTTGACCGCGAGAAGATGTCCAAAt                                                                 >  1:3321352/1‑69 (MQ=255)
        tATGTGAACTACTGGATGCACTCGGGGATGGTGATGGTTGACCGCGAGAAGATGTCCAAATCGCTGGGt                                                         >  1:3020852/1‑69 (MQ=255)
         aTGTGAACTACTGGATGCACTCGGGGATGGTGATGGTTGACCGCGAGAAGATGTCCAAATCGCTGGGTa                                                        <  1:2547101/69‑1 (MQ=255)
          tgtgAACTACTGGATGCACTCGGGGATGGTGATGGTTGACCGCGAGAAGATGTCCAAATCGCTGGGTaa                                                       >  1:985796/1‑69 (MQ=255)
                          cACTCGGGGATGGTGATGGTTGACCGCGAGAAGATGTCCAAATCGCTGGGTAACTTCTTTACCGTgcgc                                       <  1:1831887/69‑1 (MQ=255)
                          cACTCGGGGATGGTGATGGTTGACCGCGAGAAGATGTCCAAATCGCTGGGTAACTTCTTTACCGTgcgc                                       >  1:2261986/1‑69 (MQ=255)
                          cACTCGGGGATGGTGATGGTTGACCGCGAGAAGATGTCCAAATCGCTGGGTAACTTCTTTACCGTgcgc                                       >  1:2440654/1‑69 (MQ=255)
                          cACTCGGGGATGGTGATGGTTGACCGCGAGAAGATGTCCAAATCGCTGGGTAACTTCTTTACCGTgcgc                                       <  1:2687502/69‑1 (MQ=255)
                          cACTCGGGGATGGTGATGGTTGACCGCGAGAAGATGTCCAAATCGCTGGGTAACTTCTTTACCGTgcgc                                       >  1:883165/1‑69 (MQ=255)
                                  gatggtgatggtTGACCGCGAGAAGATGTCCAAATCGCTGGGTAACTTCTTTACCGTGCGCGATGTGc                                >  1:3434034/1‑68 (MQ=255)
                                      gtgatggtTGACCGCGAGAAGATGTCCAAATCGCTTGGTAACTTCTTTACCGTGCGCGATGTGCTGaaa                           <  1:2064157/69‑1 (MQ=255)
                                      gtgatggtTGACCGCGAGAAGATGTCCAAATCGCTGGGTAACTTCTTTACCGTGCGCGATGTGGTGaaa                           <  1:2625940/69‑1 (MQ=255)
                                      gtgatggtTGACCGCGAGAAGATGTCCAAATCGCTGGGTAACTTCTTTACCGTGCGCGATGTGCTGaaa                           <  1:3434243/69‑1 (MQ=255)
                                      gtgatggtTGACCGCGAGAAGATGTCCAAATCGCTGGGTAACTTCTTTACCGTGCGCGATGTGCTGaaa                           <  1:2206775/69‑1 (MQ=255)
                                      gtgatggtTGACCGCGAGAAGATGTCCAAATCGCTGGGTAACTTCTTTACCGTGCGCGATGTGCTGaaa                           <  1:2186396/69‑1 (MQ=255)
                                      gtgatggtTGACCGCGAGAAGATGTCCAAATCGCTGGGTAACTTCTTTACCGTGCGCGATGTGCTGaaa                           <  1:371494/69‑1 (MQ=255)
                                      gtgatggtTGACCGCGAGAAGATGTCCAAATCGCTGGGTAACTTCTTTACCGTGCGCGATGTGCTGaaa                           <  1:998890/69‑1 (MQ=255)
                                      gtgatggtTGACCGCGAGAAGATGTCCAAATCGCTGGGTAACTTCTTTACCGTGCGCGATGTGCTGaaa                           <  1:1428794/69‑1 (MQ=255)
                                      gtgatggtTGACCGCGAGAAGATGTCCAAATCGCTGGGTAACTTCTTTACCGTGCGCGATGTGCTGaaa                           <  1:1208413/69‑1 (MQ=255)
                                      gtgatggtTGACCGCGAGAAGATGTCCAAATCGCTGGGTAACTTCTTTACCGTGCGCGATGTGCTGaaa                           <  1:1639643/69‑1 (MQ=255)
                                        gatggtTGACCGCGAGAAGATGTCCAAATCGCTGGGTAActtc                                                   <  1:979254/43‑1 (MQ=255)
                                        gatggtTGACCGCGAGAAGATGTCCAAATCGCTGGGTAActtc                                                   <  1:2879063/43‑1 (MQ=255)
                                        gatggtTGACCGCGAGAAGATGTCCAAATCGCTGGGTAActtc                                                   <  1:2820974/43‑1 (MQ=255)
                                        gatggtTGACCGCGAGAAGATGTCCAAATCGCTGGGTAActtc                                                   <  1:1489747/43‑1 (MQ=255)
                                        gatggtTGACCGCGAGAAGATGTCCAAATCGCTGGGTAActtc                                                   <  1:428049/43‑1 (MQ=255)
                                        gatggtTGACCGCGAGAAGATGTCCAAATCGCTGGGTAActtc                                                   <  1:2203867/43‑1 (MQ=255)
                                        gatggtTGACCGCGAGAAGATGTCCAAATCGCTGGGTAActtc                                                   <  1:698668/43‑1 (MQ=255)
                                        gatggtTGACCGCGAGAAGATGTCCAAATCGCTGGGTAActtc                                                   <  1:718221/43‑1 (MQ=255)
                                        gatggtTGACCGCGAGAAGATGTCCAAATCGCTGGGTAActtc                                                   <  1:1830949/43‑1 (MQ=255)
                                        gatggtTGACCGCGAGAAGATGTCCAAATCGCTGGGTAActtc                                                   <  1:1714773/43‑1 (MQ=255)
                                        gatggtTGACCGCGAGAAGATGTCCAAATCGCTGGGTAActtc                                                   <  1:2042608/43‑1 (MQ=255)
                                        gatggtTGACCGCGAGAAGATGTCCAAATCGCTGGGTAACTTCTTTACCGTGCGCGATGTGCTGaaa                           <  1:3065381/67‑1 (MQ=255)
                                        gatggtTGACCGCGAGAAGATGGCCAAATCGCTGGGTAActtc                                                   <  1:1436642/43‑1 (MQ=255)
                                         attgtTGACCGCGAGAAGATGTCCAAATCGCTGGGTAACTTCTTTACCGTGCGCGATGTGCTGAAAtac                        <  1:2101674/66‑1 (MQ=255)
                                         atggtTGACCGCGAGAAGATGTCCAAATCGCTGGGTAACTTCTTTACCGTGCGCGATGTGCTGAAAtac                        <  1:825544/69‑1 (MQ=255)
                                         atggtTGACCGCGAGAAGATGTCCAAATCGCTGGGTAACTTCTTTACCGTGCGCGATGTGCTGAAAtac                        <  1:1285522/69‑1 (MQ=255)
                                         atggtTGACCGCGAGAAGATGTCCAAATCGCTGGGTAACTTCTTTACCGTGCGCGATGTGCTGAAAtac                        <  1:3269909/69‑1 (MQ=255)
                                         atggtTGACCGCGAGAAGATGTCCAAATCGCTGGGTAACTTCTTTACCGTGCGCGATGTGCTGAAAtac                        <  1:2818854/69‑1 (MQ=255)
                                         atggtTGACCGCGAGAAGATGTCCAAATCGCTGGGTAACTTCTTTACCGTGCGCGATGTGCTGAAAtac                        <  1:3154716/69‑1 (MQ=255)
                                          tggtTGACCGCGAGAAGATGTCTAAATCGCTGGGTAACTTCTTCACCGTGCGTGATGTGCTGAAAtac                        >  1:433740/1‑68 (MQ=255)
                                          tggtTGACCGCGAGAAGATGTCTAAATCGCTGGGTAACTTCTTCACCGTGCGTGATGTGCTGAAAtac                        >  1:346634/1‑68 (MQ=255)
                                          tggtTGACCGCGAGAAGATGTCTAAATCGCTGGGTAACTTCTTCACCGTGCGTGATGTGCTGAAAtac                        >  1:975862/1‑68 (MQ=255)
                                          tggtTGACCGCGAGAAGATGTCTAAATCGCTGGGTAACTTCTTCACCGTGCGTGATGTGCTGAAAtac                        >  1:450180/1‑68 (MQ=255)
                                          tggtTGACCGCGAGAAGATGTCTAAATCGCTGGGTAACTTCTTCACCGTGCGTGATGTGCTGAAAtac                        >  1:3393803/1‑68 (MQ=255)
                                          tggtTGACCGCGAGAAGATGTCTAAATCGCTGGGTAACTTCTTCACCGTGCGTGATGTGCTGAAAtac                        >  1:3327219/1‑68 (MQ=255)
                                          tggtTGACCGCGAGAAGATGTCTAAATCGCTGGGTAACTTCTTCACCGTGCGTGATGTGCTGAAAtac                        >  1:1791371/1‑68 (MQ=255)
                                          tggtTGACCGCGAGAAGATGTCTAAATCGCTGGGTAACTTCTTCACCGTGCGTGATGTGCTGAAAtac                        >  1:3084287/1‑68 (MQ=255)
                                          tggtTGACCGCGAGAAGATGTCTAAATCGCTGGGTAACTTCTTCACCGTGCGTGATGTGCTGAAAtac                        >  1:3013652/1‑68 (MQ=255)
                                          tggtTGACCGCGAGAAGATGTCTAAATCGCTGGGTAACTTCTTCACCGTGCGTGATGTGCTGAAAtac                        >  1:1724244/1‑68 (MQ=255)
                                          tggtTGACCGCGAGAAGATGTCTAAATCGCTGGGTAACTTCTTCACCGTGCGTGATGTGCTGAAAtac                        >  1:1704700/1‑68 (MQ=255)
                                          tggtTGACCGCGAGAAGATGTCTAAATCGCTGGGTAACTTCTTCACCGTGCGTGATGTGCTGAAAtac                        >  1:2064074/1‑68 (MQ=255)
                                          tggtTGACCGCGAGAAGATGTCTAAATCGCTGGGTAACTTCTTCACCGTGCGTGATGTGCTGAAAtac                        >  1:1623802/1‑68 (MQ=255)
                                          tggtTGACCGCGAGAAGATGTCTAAATCGCTGGGTAACTTCTTCACCGTGCGTGATGTGCTGAAAtac                        >  1:1560094/1‑68 (MQ=255)
                                          tggtTGACCGCGAGAAGATGTCTAAATCGCTGGGTAACTTCTTCACCGTGCGTGATGTGCTGAAAtac                        >  1:2289937/1‑68 (MQ=255)
                                          tggtTGACCGCGAGAAGATGTCTAAATCGCTGGGTAACTTCTTCACCGTGCGTGATGTGCTGAAAtac                        >  1:2769218/1‑68 (MQ=255)
                                          tggtTGACCGCGAGAAGATGTCTAAATCGCTGGGTAACTTCTTCACCGTGCGTGATGTGCTGAAAtac                        >  1:2782569/1‑68 (MQ=255)
                                          tggtTGACCGCGAGAAGATGTCTAAATCGCTGGGTAACTTCTTCACCGTGCGTGATGTGCTGAAAtac                        >  1:1372647/1‑68 (MQ=255)
                                          tggtTGACCGCGAGAAGATGTCTAAATCGCTGGGTAACTTCTTCACCGTGCGTGATGTGCTGAAAtac                        >  1:2802933/1‑68 (MQ=255)
                                            gtTGACCGCGAGAAGATGTCCAAATCGCTGGGTAACTTCTTTACCGTGCGCGATGTGCTGAAAtactac                     >  1:1271981/1‑69 (MQ=255)
                                              tgcccGCGAGAAGATGTCCAAATCGCTGGGTAACTTCTTTACCGTGCGCGATGTGCTGaaa                           <  1:12639/58‑1 (MQ=255)
                                              tGACCGCGAGAAGATGTCTAAATCGCTGGGTAACTTCTTCACCGTGCGTGATGTGCTGAAAta                         <  1:630835/63‑1 (MQ=255)
                                              tGACCGCGAGAAGATGTCTAAATCGCTGGGTAACTTCTTCACCGTGCGTGATGTGCTGAAAta                         <  1:1145980/63‑1 (MQ=255)
                                              tGACCGCGAGAAGATGTCTAAATCGCTGGGTAACTTCTTCACCGTGCGTGATGTGCTGAAAta                         <  1:1191972/63‑1 (MQ=255)
                                              tGACCGCGAGAAGATGTCTAAATCGCTGGGTAACTTCTTCACCGTGCGTGATGTGCTGAAAta                         <  1:3219873/63‑1 (MQ=255)
                                              tGACCGCGAGAAGATGTCTAAATCGCTGGGTAACTTCTTCACCGTGCGTGATGTGCTGAAAta                         <  1:798050/63‑1 (MQ=255)
                                              tGACCGCGAGAAGATGTCTAAATCGCTGGGTAACTTCTTCACCGTGCGTGATGTGCTGAAAta                         <  1:2139808/63‑1 (MQ=255)
                                              tGACCGCGAGAAGATGTCTAAATCGCTGGGTAACTTCTTCACCGTGCGTGATGTGCTGAAAta                         <  1:3369772/63‑1 (MQ=255)
                                              tGACCGCGAGAAGATGTCTAAATCGCTGGGTAACTTCTTCACCGTGCGTGATGTGCTGAAAta                         <  1:3304961/63‑1 (MQ=255)
                                              tGACCGCGAGAAGATGTCTAAATCGCTGGGTAACTTCTTCACCGTGCGTGATGTGCTGAAAta                         <  1:439861/63‑1 (MQ=255)
                                              tGACCGCGAGAAGATGTCTAAATCGCTGGGTAACTTCTTCACCGTGCGTGATGTGCTGAAAta                         <  1:1438285/63‑1 (MQ=255)
                                              tGACCGCGAGAAGATGTCTAAATCGCTGGGTAACTTCTTCACCGTGCGTGATGTGCTGAAAta                         <  1:3426158/63‑1 (MQ=255)
                                               gACCGCGAGAAGATGTCTAAATCGCTGGGTAACTTCTTCACCGTGCGTGATGTGCTGAAAta                         <  1:1936163/62‑1 (MQ=255)
                                                           aTGTCCAAATCGCTGGGTAACTTCTTTACCGTGCGCGATGTGCTGAAATACTACGACGCGGAAACCGTg      >  1:1580406/1‑69 (MQ=255)
                                                           aTGTCCAAATCGCTGGGTAACTTCTTTACCGTGCGCGATGTGCTGAAATACTACGACGCGGAAACCGTg      >  1:1303297/1‑69 (MQ=255)
                                                                cAAATCGCTGGGTAACTTCTTTACCGTGCGCGATGTGCTGAAATACTACGACGCGGAAACCGTGCGtt  <  1:2073716/68‑1 (MQ=255)
                                                                |                                                                   
ATGGTCAGTATGTGAACTACTGGATGCACTCGGGGATGGTGATGGTTGACCGCGAGAAGATGTCCAAATCGCTGGGTAACTTCTTTACCGTGCGCGATGTGCTGAAATACTACGACGCGGAAACCGTGCGTT  >  minE/354425‑354556

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: