Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 354,510 T→C 39.3% F275F (TTT→TTC cysS → cysteinyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE354,5100TC39.3% 36.6 / 57.1 84F275F (TTT→TTCcysScysteinyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base T (11/40);  new base C (19/14);  total (30/54)
Fisher's exact test for biased strand distribution p-value = 1.09e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.15e-02

CACTCGGGGATGGTGATGGTTGACCGCGAGAAGATGTCCAAATCGCTGGGTAACTTCTTTACCGTGCGCGATGTGCTGAAATACTACGACGCGGAAACCGTGCGTTACTTCCTGATGTCGGGCC  >  minE/354451‑354574
                                                           |                                                                
cACTCGGGGATGGTGATGGTTGACCGCGAGAAGATGTCCAAATCGCTGGGTAACTTCTTTACCGTgcgc                                                         <  1:2687502/69‑1 (MQ=255)
cACTCGGGGATGGTGATGGTTGACCGCGAGAAGATGTCCAAATCGCTGGGTAACTTCTTTACCGTgcgc                                                         >  1:2440654/1‑69 (MQ=255)
cACTCGGGGATGGTGATGGTTGACCGCGAGAAGATGTCCAAATCGCTGGGTAACTTCTTTACCGTgcgc                                                         >  1:2261986/1‑69 (MQ=255)
cACTCGGGGATGGTGATGGTTGACCGCGAGAAGATGTCCAAATCGCTGGGTAACTTCTTTACCGTgcgc                                                         >  1:883165/1‑69 (MQ=255)
cACTCGGGGATGGTGATGGTTGACCGCGAGAAGATGTCCAAATCGCTGGGTAACTTCTTTACCGTgcgc                                                         <  1:1831887/69‑1 (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:998890/69‑1 (MQ=255)
            gtgatggtTGACCGCGAGAAGATGTCCAAATCGCTGGGTAACTTCTTTACCGTGCGCGATGTGCTGaaa                                             <  1:2206775/69‑1 (MQ=255)
            gtgatggtTGACCGCGAGAAGATGTCCAAATCGCTGGGTAACTTCTTTACCGTGCGCGATGTGCTGaaa                                             <  1:2186396/69‑1 (MQ=255)
            gtgatggtTGACCGCGAGAAGATGTCCAAATCGCTGGGTAACTTCTTTACCGTGCGCGATGTGCTGaaa                                             <  1:3434243/69‑1 (MQ=255)
            gtgatggtTGACCGCGAGAAGATGTCCAAATCGCTGGGTAACTTCTTTACCGTGCGCGATGTGCTGaaa                                             <  1:371494/69‑1 (MQ=255)
            gtgatggtTGACCGCGAGAAGATGTCCAAATCGCTGGGTAACTTCTTTACCGTGCGCGATGTGCTGaaa                                             <  1:1639643/69‑1 (MQ=255)
            gtgatggtTGACCGCGAGAAGATGTCCAAATCGCTGGGTAACTTCTTTACCGTGCGCGATGTGCTGaaa                                             <  1:1208413/69‑1 (MQ=255)
            gtgatggtTGACCGCGAGAAGATGTCCAAATCGCTGGGTAACTTCTTTACCGTGCGCGATGTGCTGaaa                                             <  1:1428794/69‑1 (MQ=255)
              gatggtTGACCGCGAGAAGATGTCCAAATCGCTGGGTAACTTCTTTACCGTGCGCGATGTGCTGaaa                                             <  1:3065381/67‑1 (MQ=255)
               attgtTGACCGCGAGAAGATGTCCAAATCGCTGGGTAACTTCTTTACCGTGCGCGATGTGCTGAAAtac                                          <  1:2101674/66‑1 (MQ=255)
               atggtTGACCGCGAGAAGATGTCCAAATCGCTGGGTAACTTCTTTACCGTGCGCGATGTGCTGAAAtac                                          <  1:2818854/69‑1 (MQ=255)
               atggtTGACCGCGAGAAGATGTCCAAATCGCTGGGTAACTTCTTTACCGTGCGCGATGTGCTGAAAtac                                          <  1:1285522/69‑1 (MQ=255)
               atggtTGACCGCGAGAAGATGTCCAAATCGCTGGGTAACTTCTTTACCGTGCGCGATGTGCTGAAAtac                                          <  1:3154716/69‑1 (MQ=255)
               atggtTGACCGCGAGAAGATGTCCAAATCGCTGGGTAACTTCTTTACCGTGCGCGATGTGCTGAAAtac                                          <  1:3269909/69‑1 (MQ=255)
               atggtTGACCGCGAGAAGATGTCCAAATCGCTGGGTAACTTCTTTACCGTGCGCGATGTGCTGAAAtac                                          <  1:825544/69‑1 (MQ=255)
                tggtTGACCGCGAGAAGATGTCTAAATCGCTGGGTAACTTCTTCACCGTGCGTGATGTGCTGAAAtac                                          >  1:1560094/1‑68 (MQ=255)
                tggtTGACCGCGAGAAGATGTCTAAATCGCTGGGTAACTTCTTCACCGTGCGTGATGTGCTGAAAtac                                          >  1:450180/1‑68 (MQ=255)
                tggtTGACCGCGAGAAGATGTCTAAATCGCTGGGTAACTTCTTCACCGTGCGTGATGTGCTGAAAtac                                          >  1:3084287/1‑68 (MQ=255)
                tggtTGACCGCGAGAAGATGTCTAAATCGCTGGGTAACTTCTTCACCGTGCGTGATGTGCTGAAAtac                                          >  1:2289937/1‑68 (MQ=255)
                tggtTGACCGCGAGAAGATGTCTAAATCGCTGGGTAACTTCTTCACCGTGCGTGATGTGCTGAAAtac                                          >  1:1704700/1‑68 (MQ=255)
                tggtTGACCGCGAGAAGATGTCTAAATCGCTGGGTAACTTCTTCACCGTGCGTGATGTGCTGAAAtac                                          >  1:975862/1‑68 (MQ=255)
                tggtTGACCGCGAGAAGATGTCTAAATCGCTGGGTAACTTCTTCACCGTGCGTGATGTGCTGAAAtac                                          >  1:2802933/1‑68 (MQ=255)
                tggtTGACCGCGAGAAGATGTCTAAATCGCTGGGTAACTTCTTCACCGTGCGTGATGTGCTGAAAtac                                          >  1:3013652/1‑68 (MQ=255)
                tggtTGACCGCGAGAAGATGTCTAAATCGCTGGGTAACTTCTTCACCGTGCGTGATGTGCTGAAAtac                                          >  1:1724244/1‑68 (MQ=255)
                tggtTGACCGCGAGAAGATGTCTAAATCGCTGGGTAACTTCTTCACCGTGCGTGATGTGCTGAAAtac                                          >  1:433740/1‑68 (MQ=255)
                tggtTGACCGCGAGAAGATGTCTAAATCGCTGGGTAACTTCTTCACCGTGCGTGATGTGCTGAAAtac                                          >  1:1623802/1‑68 (MQ=255)
                tggtTGACCGCGAGAAGATGTCTAAATCGCTGGGTAACTTCTTCACCGTGCGTGATGTGCTGAAAtac                                          >  1:1372647/1‑68 (MQ=255)
                tggtTGACCGCGAGAAGATGTCTAAATCGCTGGGTAACTTCTTCACCGTGCGTGATGTGCTGAAAtac                                          >  1:2064074/1‑68 (MQ=255)
                tggtTGACCGCGAGAAGATGTCTAAATCGCTGGGTAACTTCTTCACCGTGCGTGATGTGCTGAAAtac                                          >  1:3327219/1‑68 (MQ=255)
                tggtTGACCGCGAGAAGATGTCTAAATCGCTGGGTAACTTCTTCACCGTGCGTGATGTGCTGAAAtac                                          >  1:3393803/1‑68 (MQ=255)
                tggtTGACCGCGAGAAGATGTCTAAATCGCTGGGTAACTTCTTCACCGTGCGTGATGTGCTGAAAtac                                          >  1:2769218/1‑68 (MQ=255)
                tggtTGACCGCGAGAAGATGTCTAAATCGCTGGGTAACTTCTTCACCGTGCGTGATGTGCTGAAAtac                                          >  1:1791371/1‑68 (MQ=255)
                tggtTGACCGCGAGAAGATGTCTAAATCGCTGGGTAACTTCTTCACCGTGCGTGATGTGCTGAAAtac                                          >  1:2782569/1‑68 (MQ=255)
                tggtTGACCGCGAGAAGATGTCTAAATCGCTGGGTAACTTCTTCACCGTGCGTGATGTGCTGAAAtac                                          >  1:346634/1‑68 (MQ=255)
                  gtTGACCGCGAGAAGATGTCCAAATCGCTGGGTAACTTCTTTACCGTGCGCGATGTGCTGAAAtactac                                       >  1:1271981/1‑69 (MQ=255)
                    tgcccGCGAGAAGATGTCCAAATCGCTGGGTAACTTCTTTACCGTGCGCGATGTGCTGaaa                                             <  1:12639/58‑1 (MQ=255)
                    tGACCGCGAGAAGATGTCTAAATCGCTGGGTAACTTCTTCACCGTGCGTGATGTGCTGAAAta                                           <  1:1191972/63‑1 (MQ=255)
                    tGACCGCGAGAAGATGTCTAAATCGCTGGGTAACTTCTTCACCGTGCGTGATGTGCTGAAAta                                           <  1:3304961/63‑1 (MQ=255)
                    tGACCGCGAGAAGATGTCTAAATCGCTGGGTAACTTCTTCACCGTGCGTGATGTGCTGAAAta                                           <  1:3219873/63‑1 (MQ=255)
                    tGACCGCGAGAAGATGTCTAAATCGCTGGGTAACTTCTTCACCGTGCGTGATGTGCTGAAAta                                           <  1:3369772/63‑1 (MQ=255)
                    tGACCGCGAGAAGATGTCTAAATCGCTGGGTAACTTCTTCACCGTGCGTGATGTGCTGAAAta                                           <  1:3426158/63‑1 (MQ=255)
                    tGACCGCGAGAAGATGTCTAAATCGCTGGGTAACTTCTTCACCGTGCGTGATGTGCTGAAAta                                           <  1:798050/63‑1 (MQ=255)
                    tGACCGCGAGAAGATGTCTAAATCGCTGGGTAACTTCTTCACCGTGCGTGATGTGCTGAAAta                                           <  1:439861/63‑1 (MQ=255)
                    tGACCGCGAGAAGATGTCTAAATCGCTGGGTAACTTCTTCACCGTGCGTGATGTGCTGAAAta                                           <  1:630835/63‑1 (MQ=255)
                    tGACCGCGAGAAGATGTCTAAATCGCTGGGTAACTTCTTCACCGTGCGTGATGTGCTGAAAta                                           <  1:1145980/63‑1 (MQ=255)
                    tGACCGCGAGAAGATGTCTAAATCGCTGGGTAACTTCTTCACCGTGCGTGATGTGCTGAAAta                                           <  1:2139808/63‑1 (MQ=255)
                    tGACCGCGAGAAGATGTCTAAATCGCTGGGTAACTTCTTCACCGTGCGTGATGTGCTGAAAta                                           <  1:1438285/63‑1 (MQ=255)
                     gACCGCGAGAAGATGTCTAAATCGCTGGGTAACTTCTTCACCGTGCGTGATGTGCTGAAAta                                           <  1:1936163/62‑1 (MQ=255)
                                 aTGTCCAAATCGCTGGGTAACTTCTTTACCGTGCGCGATGTGCTGAAATACTACGACGCGGAAACCGTg                        >  1:1303297/1‑69 (MQ=255)
                                 aTGTCCAAATCGCTGGGTAACTTCTTTACCGTGCGCGATGTGCTGAAATACTACGACGCGGAAACCGTg                        >  1:1580406/1‑69 (MQ=255)
                                      cAAATCGCTGGGTAACTTCTTTACCGTGCGCGATGTGCTGAAATACTACGACGCGGAAACCGTGCGtt                    <  1:2073716/68‑1 (MQ=255)
                                            gCTGGGTAACTTCTTTACCGTGCGCGATGTGCTGAAATACTACGACGCGGAAACCGTGCGtt                    >  1:2377693/1‑62 (MQ=255)
                                            gCTGGGTAACTTCTTTACCGTGCGCGATGTGCTGAAATACTACGACGCGGAAACCGTGCGtt                    >  1:658294/1‑62 (MQ=255)
                                            gCTGGGTAACTTCTTTACCGTGCGCGATGTGCTGAAATACTACGACGCGGAAACCGTGCGtt                    >  1:1663877/1‑62 (MQ=255)
                                            gCTGGGTAACTTCTTTACCGTGCGCGATGTGCTGAAATACTACGACGCGGAAACCGTGCGtt                    >  1:1775102/1‑62 (MQ=255)
                                            gCTGGGTAACTTCTTTACCGTGCGCGATGTGCTGAAATACTACGACGCGGAAACCGTGCGTTACTTCCt             <  1:2653220/69‑1 (MQ=255)
                                            gCTGGGTAACTTCTTTACCGTGCGCGATGTGCTGAAATACTACGACGCGGAAACCGTGCGTTACTTCCt             <  1:1286440/69‑1 (MQ=255)
                                            gCTGGGTAACTTCTTTACCGTGCGCGATGTGCTGAAATACTACGACGCGGAAACCGTGCGTTACTTCCt             <  1:428542/69‑1 (MQ=255)
                                            gCTGGGTAACTTCTTTACCGTGCGCGATGTGCTGAAATACTACGACGCGGAAACCGTGCGTTACTTCCt             <  1:2139751/69‑1 (MQ=255)
                                            gCTGGGTAACTTCTTTACCGTGCGCGATGTGCTGAAATACTACGACGCGGAAACCGTGCGTTACTTCCt             <  1:2903534/69‑1 (MQ=255)
                                            gCTGGGTAACTTCTTTACCGTGCGCGATGTGCTGAAATACTACGACGCGGAAACCGTGCGTTACTTCCt             <  1:2303165/69‑1 (MQ=255)
                                            gCTGGGTAACTTCTTTACCGTGCGCGATGTGCTGAAATACTACGACGCGGAAACCGTGCGTTACTTCCt             <  1:3255532/69‑1 (MQ=255)
                                            gCTGGGTAACTTCTTTACCGTGCGCGATGTGCTGAAATACTACGACGCGGAAACCGTGCGTTACTTCCt             <  1:1838746/69‑1 (MQ=255)
                                            gCTGGGTAACTTCTTTACCGTGCGCGATGTGCTGAAATACTACGACGCGGAAACCGTGCGTTACTTCCt             <  1:329321/69‑1 (MQ=255)
                                            gCTGGGTAACTTCTTTACCGTGCGCGATGTGCTGAAATACTACGACGCGGAAACCGTGCGTTACTTCCt             <  1:2008539/69‑1 (MQ=255)
                                            gCTGGGTAACTTCTTTACCGTGCGCGATGTGCTGAAATACTACGACGCGGAAACCGTGCGTTACTTCCt             <  1:1971512/69‑1 (MQ=255)
                                             cTGGGTAACTTCTTCACCGTGCGTGATGTGCTGAAAta                                           <  1:1828409/38‑1 (MQ=38)
                                               gggTAACTTCTTCACCGTGCGTGATGTGCTGAAAta                                           <  1:2142185/36‑1 (MQ=37)
                                                 gTAACTTCTTTACCGTGCGCGATGTGCTGAAATACTACGACGCGGAAACCGTGCGTTACTTCCt             <  1:1937059/64‑1 (MQ=255)
                                                      ttcttTACCGTGCGCGATGTGCTGAAATACTACGACGCGGAAACCGTGCGTTACTTCCTGATGTCGGGc   <  1:677788/69‑1 (MQ=255)
                                                      ttcttTACCGTGCGCGATGTGCTGAAATACTACGACGCGGAAACCGTGCGTTACTTCCTGATGTCGGGc   <  1:1573101/69‑1 (MQ=255)
                                                      ttcttTACCGTGCGCGATGTGCTGAAATACTACGACGCGGAAACCGTGCGTTACTTCCTGATGTCGGGc   <  1:180755/69‑1 (MQ=255)
                                                      ttcttTACCGTGCGCGATGTGCTGAAATACTACGACGCGGAAACCGTGCGTTACTTCCTGATGTCGGGc   <  1:1838335/69‑1 (MQ=255)
                                                       tcttTACCGTGCGCGATGTGCTTAAATACTACGACGCGGAAACCGTGCGTTACTTCCTGATGTCGGGcc  <  1:745480/69‑1 (MQ=255)
                                                       tcttTACCGTGCGCGATGTGCTGAAATACTACGACGCGGAAACCGTGCGTTACTTCCTGATGTCGGGcc  <  1:1706025/69‑1 (MQ=255)
                                                       tcttTACCGTGCGCGATGTGCTGAAATACTACGACGCGGAAACCGTGCGTTACTTCCTGATGTCGGGcc  <  1:681498/69‑1 (MQ=255)
                                                           |                                                                
CACTCGGGGATGGTGATGGTTGACCGCGAGAAGATGTCCAAATCGCTGGGTAACTTCTTTACCGTGCGCGATGTGCTGAAATACTACGACGCGGAAACCGTGCGTTACTTCCTGATGTCGGGCC  >  minE/354451‑354574

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: