Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,814,214 A→G 46.5% G321G (GGT→GGC)  alaS ← alanyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,814,2140AG46.5% 9.2 / 86.1 99G321G (GGT→GGC) alaSalanyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base A (48/5);  new base G (26/20);  total (74/25)
Fisher's exact test for biased strand distribution p-value = 1.43e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.89e-01

Only 100 of 102 total aligned reads displayed.

ACGTCGATCAGCGGACCAACCAGTTTGTAGAAGAAGGTTTCTTTCGCGCCGAGCATATTACCGTGACGCACTGCGCGACGAATGATACGACGCAGTACATAACCACGGTTTTCATTGGACGGCATTAC  >  minE/1814155‑1814282
                                                           |                                                                    
aCGTCGATCAGCGGACCAACCAGTTTGTAGAAGAAGGTTTCTTTCGCGCCGAGCATATTACCGTGACg                                                              >  1:3088839/1‑68 (MQ=255)
aCGTCGATCAGCGGACCAACCAGTTTGTAGAAGAAGGTTTCTTTCGCGCCGAGCATATTACCGTGACGc                                                             >  1:2142227/1‑69 (MQ=255)
aCGTCGATCAGCGGACCAACCAGTTTGTAGAAGAAGGTTTCTTTCGCGCCGAGCATATTACCGTGACGc                                                             >  1:1844048/1‑69 (MQ=255)
aCGTCGATCAGCGGACCAACCAGTTTGTAGAAGAAGGTTTCTTTCGCGCCGAGCATATTACCGTGACGc                                                             >  1:320989/1‑69 (MQ=255)
aCGTCGATCAGCGGACCAACCAGTTTGTAGAAGAAGGTTTCTTTCGCGCCGAGCATATTACCGTGACGc                                                             >  1:1881581/1‑69 (MQ=255)
aCGTCGATCAGCGGACCAACCAGTTTGTAGAAGAAGGTTTCTTTCGCGCCGAGCATATTACCGTGACGc                                                             >  1:2749944/1‑69 (MQ=255)
aCGTCGATCAGCGGACCAACCAGTTTGTAGAAGAAGGTTTCTTTCGCGCCGAGCATATTACCGTGACGc                                                             >  1:3413795/1‑69 (MQ=255)
aCGTCGATCAGCGGACCAACCAGTTTGTAGAAGAAGGTTTCTTTCGCGCCGAGCATATTACCGTGACGc                                                             >  1:2612801/1‑69 (MQ=255)
aCGTCGATCAGCGGACCAACCAGTTTGTAGAAGAAGGTTTCTTTCGCGCCGAGCATATTACCGTGACGc                                                             >  1:2599672/1‑69 (MQ=255)
aCGTCGATCAGCGGACCAACCAGTTTGTAGAAGAAGGTTTCTTTCGCGCCGAGCATATTACCGTGACGc                                                             <  1:495156/69‑1 (MQ=255)
aCGTCGATCAGCGGACCAACCAGTTTGTAGAAGAAGGTTTCTTTCGCGCCGAGCATATTACCGTGACGc                                                             >  1:1451610/1‑69 (MQ=255)
aCGTCGATCAGCGGACCAACCAGTTTGTAGAAGAAGGTTTCTTTCGCGCCGAGCATATTACCGTGACGc                                                             >  1:1423741/1‑69 (MQ=255)
aCGTCGATCAGCGGACCAACCAGTTTGTAGAAGAAGGTTTCTTTCGCGCCGAGCATATTACCGTGACGc                                                             >  1:1079609/1‑69 (MQ=255)
aCGTCGATCAGCGGACCAACCAGTTTGTAGAAGAAGGTTTCTTTCGCGCCGAGCATATTACCGTGACGc                                                             >  1:6086/1‑69 (MQ=255)
aCGTCGATCAGCGGACCAACCAGTTTGTAGAAGAAGGTTTCTTTCGCGCCGAGCATATTACCGTGACGc                                                             >  1:1340310/1‑69 (MQ=255)
aCGTCGATCAGCGGACCAACCAGTTTGTAGAAGAAGGTTTCTTTCGCGCCGAGCATATTACCGTGACGc                                                             >  1:1268659/1‑69 (MQ=255)
aCGTCGATCAGCGGACCAACCAGTTTGTAGAAGAAGGTTTCTTTCGCGCCGAGCATATTACCGTGACGc                                                             >  1:2459319/1‑69 (MQ=255)
aCGTCGATCAGCGGACCAACCAGTTTGTAGAAGAAGGTTTCTTTCGCGCCGAGCATATTACCGTGACGc                                                             >  1:1100248/1‑69 (MQ=255)
aCGTCGAACAGCGGACCAACCAGTTTGTAGAAGAAGGTTTCTTTCGCGCCGAGCATATTACCGTGACGc                                                             >  1:2359762/1‑69 (MQ=255)
        cAGCGGACCAACCAGTTTGTAGAAGAAGGTTTCTTTCGCGCCGAGCATATTGCCGTGACGCACTgcgc                                                      >  1:1388039/1‑68 (MQ=255)
        cAGCGGACCAACCAGTTTGTAGAAGAAGGTTTCTTTCGCGCCGAGCATATTGCCGTGACGCACTgcgc                                                      >  1:2066394/1‑68 (MQ=255)
        cAGCGGACCAACCAGTTTGTAGAAGAAGGTTTCTTTCGCGCCGAGCATATTGCCGTGACGCACTGcgcg                                                     >  1:2313710/1‑69 (MQ=255)
        cAGCGGACCAACCAGTTTGTAGAAGAAGGTTTCTTTCGCGCCGAGCATATTGCCGTGACGCACTGcgcg                                                     >  1:2387306/1‑69 (MQ=255)
        cAGCGGACCAACCAGTTTGTAGAAGAAGGTTTCTTTCGCGCCGAGCATATTGCCGTGACGCACTGcgcg                                                     >  1:296374/1‑69 (MQ=255)
        cAGCGGACCAACCAGTTTGTAGAAGAAGGTTTCTTTCGCGCCGAGCATATTGCCGTGACGCACTGcgcg                                                     >  1:2801881/1‑69 (MQ=255)
        cAGCGGACCAACCAGTTTGTAGAAGAAGGTTTCTTTCGCGCCGAGCATATTGCCGTGACGCACTGcgcg                                                     >  1:1957440/1‑69 (MQ=255)
        cAGCGGACCAACCAGTTTGTAGAAGAAGGTTTCTTTCGCGCCGAGCATATTGCCGTGACGCACTGcgcg                                                     >  1:1961619/1‑69 (MQ=255)
        cAGCGGACCAACCAGTTTGTAGAAGAAGGTTTCTTTCGCGCCGAGCATATTGCCGTGACGCACTGcgcg                                                     >  1:217700/1‑69 (MQ=255)
        cAGCGGACCAACCAGTTTGTAGAAGAAGGTTTCTTTCGCGCCGAGCATATTGCCGTGACGCACTGcgcg                                                     >  1:3215099/1‑69 (MQ=255)
        cAGCGGACCAACCAGTTTGTAGAAGAAGGTTTCTTTCGCGCCGAGCATATTGCCGTGACGCACTGcgcg                                                     >  1:2528366/1‑69 (MQ=255)
        cAGCGGACCAACCAGTTTGTAGAAGAAGGTTTCTTTCGCGCCGAGCATATTGCCGTGACGCACTGcgcg                                                     >  1:1537572/1‑69 (MQ=255)
        cAGCGGACCAACCAGTTTGTAGAAGAAGGTTTCTTTCGCGCCGAGCATATTGCCGTGACGCACTGcgcg                                                     >  1:1023866/1‑69 (MQ=255)
        cAGCGGACCAACCAGTTTGTAGAAGAAGGTTTCTTTCGCGCCGAGCATATTGCCGTGACGCACTGcgcg                                                     >  1:831323/1‑69 (MQ=255)
        cAGCGGACCAACCAGTTTGTAGAAGAAGGTTTCTTTCGCGCCGAGCATATTGCCGTGACGCACTGcgcg                                                     >  1:787321/1‑69 (MQ=255)
        cAGCGGACCAACCAGTTTGTAGAAGAAGGTTTCTTTCGCGCCGAGCATATTGCCGTGACGCACTGcgcg                                                     >  1:125670/1‑69 (MQ=255)
        cAGCGGACCAACCAGTTTGTAGAAGAAGGTTTCTTTCGCGCCGAGCATATTGCCGTGACGCACTGcgcg                                                     >  1:1410346/1‑69 (MQ=255)
        cAGCGGACCAACCAGTTTGTAGAAGAAGGTTTCTTTCGCGCCGAGCATATTGCCGTGACGCACTGcgcg                                                     >  1:1410560/1‑69 (MQ=255)
        cAGCGGACCAACCAGTTTGTAGAAGAAGGTTTCTTTCGCGCCGAGCATATTGCCGTGACGCACTGcgcg                                                     >  1:1490723/1‑69 (MQ=255)
        cAGCGGACCAACCAGTTTGTAGAAGAAGGTTTCTTTCGCGCCGAGCATATTGCCGTGACGCACTGcgcg                                                     >  1:1700862/1‑69 (MQ=255)
        cAGCGGACCAACCAGTTTGTAGAAGAAGGTTTCTTTCGCGCCGAGCATATTGCCGTGACGCACTGcgcg                                                     >  1:1569392/1‑69 (MQ=255)
        cAGCGGACCAACCAGTTTGTAGAAGAAGGTTTCTTTCGCGCCGAGCATATTGCCGTGACGCACTGcgcg                                                     >  1:3394745/1‑69 (MQ=255)
        cAGCGGACCAACCAGTTTGTAGAAGAAGGTTTCTTTCGCGCCGAGCATATTGCCGTGACGCACTGcgcg                                                     >  1:1633766/1‑69 (MQ=255)
        cAGCGGACCAACCAGTTTGTAGAAGAAGGTTTCTTTCGCGCCGAGCATATTGCCGTGACGCACTGcgcg                                                     >  1:3289507/1‑69 (MQ=255)
        cAGCGGACCAACCAGTTTGTAGAAGAAGGTTTCTTTCGCGCCGAGCATATTGCCGTGACGCACTGcgcg                                                     >  1:1659455/1‑69 (MQ=255)
        cAGCGGACCAACCAGTTTGTAGAAGAAGGTTTCTTTCGCGCCGAGCATATTGCCGTGACGCACTGcgcg                                                     >  1:1676413/1‑69 (MQ=255)
        cAGCGGACCAACCAGTTTGTAGAAGAAGGTTTCTTTCGCGCCGAGCATATTACCGTGACg                                                              >  1:3006356/1‑60 (MQ=255)
        cAGCGGACCAACCAGTTTGTAGAAGAAGGTTTCTTTCGCGCCGAGCATATTACCGTGACGCACTgcgc                                                      >  1:2596216/1‑68 (MQ=255)
        cAGCGGACCAACCAGTTTGTAGAAGAAGGTTTCTTTCGCGCCGAGCATATTACCGTGACGCACTGcgcg                                                     >  1:476330/1‑69 (MQ=255)
                    caGTTTGTAGAAGAAGGTTTCTTTCGCGCCGAGCATATTACCGTGACGCACTGCGCGACGAATGATacg                                         >  1:2943213/1‑69 (MQ=255)
                    caGTTTGTAGAAGAAGGTTTCTTTCGCGCCGAGCATATTACCGTGACGCACTGCGCGACGAATGATacg                                         >  1:1033495/1‑69 (MQ=255)
                    caGTTTGTAGAAGAAGGTTTCTTTCGCGCCGAGCATATTACCGTGACGCACTGCGCGACGAATGATacg                                         >  1:2289938/1‑69 (MQ=255)
                    caGTTTGTAGAAGAAGGTTTCTTTCGCGCCGAGCATATTACCGTGACGCACTGCGCGACGAATGATacg                                         >  1:1682049/1‑69 (MQ=255)
                    caGTTTGTAGAAGAAGGTTTCTTTCGCGCCGAGCATATTACCGTGACGCACTGCGCGACGAATGATacg                                         >  1:1644390/1‑69 (MQ=255)
                    caGTTTGTAGAAGAAGGTTTCTTTCGCGCCGAGCATATTACCGTGACGCACTGCGCGACGAATGATacg                                         >  1:2280646/1‑69 (MQ=255)
                    caGTTTGTAGAAGAAGGTTTCTTTCGCGCCGAGCATATTACCGTGACGCACTGCGCGACGAATGATacg                                         >  1:136528/1‑69 (MQ=255)
                    caGTTTGTAGAAGAAGGTTTCTTTCGCGCCGAGCATATTACCGTGACGCACTGCGCGACGAATGATac                                          >  1:1607283/1‑68 (MQ=255)
                    caGTTTGTAGAAGAAGGTTTCTTTAGCGCCGAGCATATTACCGTGACGCACTGCGCGACGAATGATacg                                         >  1:2549722/1‑69 (MQ=255)
                    caGTTTGTAGAAGAAGGTTACTTTCGCGCCGAGCATATAACCGTCACGCACTGCGCGACGAATGATac                                          >  1:483539/1‑68 (MQ=255)
                            agaagaAGGTTTCTTTCGCGCCGAGCATATTACCGTGACGCACTGCGCGACGAATGATacgacg                                      <  1:1930359/64‑1 (MQ=255)
                               agaagGTTTCTTTCGCGCCGAGCATATTGCCGTGACGCACTGCGCGACGAATGATACGACGCAGTACAt                              <  1:510417/69‑1 (MQ=255)
                               agaagGTTTCTTTCGCGCCGAGCATATTGCCGTGACGCACTGCGCGACGAATGATACGACGCAGTACAt                              <  1:2070779/69‑1 (MQ=255)
                               agaagGTTTCTTTCGCGCCGAGCATATTGCCGTGACGCACTGCGCGACGAATGATACGACGCAGTACAt                              <  1:2725365/69‑1 (MQ=255)
                               agaagGTTTCTTTCGCGCCGAGCATATTGCCGTGACGCACTGCGCGACGAATGATACGACGCAGTACAt                              <  1:58925/69‑1 (MQ=255)
                               agaagGTTTCTTTCGCGCCGAGCATATTGCCGTGACGCACTGCGCGACGAATGATACGACGCAGTACAt                              <  1:822075/69‑1 (MQ=255)
                               agaagGTTTCTTTAGCGCCGAGCATATTGCCGTGACGCACTGCGCGACGAATGATACGACGCAGTACAt                              <  1:2151897/69‑1 (MQ=255)
                                   ggTTTCTTTCGCGCCGAGCATATTACCGTGACGCACTGCGCGACGAATGATACGACGCAGTACATAAcc                          >  1:529566/1‑69 (MQ=255)
                                   ggTTTCTTTCGCGCCGAGCATATTACCGTGACGCACTGCGCGACGAATGATACGACGCAGTACATAAcc                          >  1:232690/1‑69 (MQ=255)
                                   ggTTTCTTTCGCGCCGAGCATATTACCGTGACGCACTGCGCGACGAATGATACGACGCAGTACATAAcc                          >  1:3079133/1‑69 (MQ=255)
                                   ggTTTCTTTCGCGCCGAGCATATTACCGTGACGCACTGCGCGACGAATGATACGACGCAGTACATAAcc                          >  1:866660/1‑69 (MQ=255)
                                   ggTTTCTTTCGCGCCGAGCATATTACCGTGACGCACTGCGCGACGAATGATACGACGCAGTACATAAcc                          >  1:2592274/1‑69 (MQ=255)
                                   ggTTTCTTTCGCGCCGAGCATATTACCGTGACGCACTGCGCGACGAATGATACGACGCAGTACATAAcc                          >  1:2464392/1‑69 (MQ=255)
                                   ggTTTCTTTCGCGCCGAGCATATTACCGTGACGCACTGCGCGACGAATGATACGACGCAGTACATAAcc                          >  1:2479738/1‑69 (MQ=255)
                                   ggTTTCTTTCGCGCCGAGCATATTACCGTGACGCACTGCGCGACGAATGATACGACGCAGTACATAAcc                          >  1:2029541/1‑69 (MQ=255)
                                   ggTTTCTTTCGCGCCGAGCATATTACCGTGACGCACTGCGCGACGAATGATACGACGCAGTACATAAcc                          >  1:1946483/1‑69 (MQ=255)
                                   ggTTTCTTTCGCGCCGAGCATATTACCGTGACGCACTGCGCGACGAATGATACGACGCAGTACATAAc                           >  1:60123/1‑68 (MQ=255)
                                    gTTTCTTTCGCGCCGAGCATATTGCCGTGACGCACTGCGCGACGAATGATACGACGCAGTACATAACCa                         <  1:573731/69‑1 (MQ=255)
                                    gTTTCTTTCGCGCCGAGCATATTGCCGTGACGCACTGCGCGACGAATGATACGACGCAGTACATAACCa                         <  1:1888994/69‑1 (MQ=255)
                                    gTTTCTTTCGCGCCGAGCATATTGCCGTGACGCACTGCGCGACGAATGATACGACGCAGTACATAACCa                         <  1:840951/69‑1 (MQ=255)
                                    gTTTCTTTCGCGCCGAGCATATTGCCGTGACGCACTGCGCGACGAATGATACGACGCAGTACATAACCa                         <  1:1224134/69‑1 (MQ=255)
                                    gTTTCTTTCGCGCCGAGCATATTGCCGTGACGCACTGCGCGACGAATGATACGACGCAGTACATAACCa                         <  1:3033446/69‑1 (MQ=255)
                                    gTTTCTTTCGCGCCGAGCATATTGCCGTGACGCACTGCGCGACGAATGATACGACGCAGTACATAACCa                         <  1:2557889/69‑1 (MQ=255)
                                    gTTTCTTTCGCGCCGAGCATATTGCCGTGACGCACTGCGCGACGAATGATACGACGCAGTACATAACCa                         <  1:2818045/69‑1 (MQ=255)
                                    gTTTCTTTCGCGCCGAGCATATTGCCGTGACGCACTGCGCGACGAATGATACGACGCAGTACATAACCa                         <  1:2565679/69‑1 (MQ=255)
                                    gTTTCTTTCGCGCCGAGCATATTGCCGTGACGCACTGCGCGACGAATGATACGACGCAGTACATAACCa                         <  1:52752/69‑1 (MQ=255)
                                    gTTTCTTTCGCGCCGAGCATATTGCCGTGACGCACTGCGCGACGAATGATACGACGCAGTACATAACCa                         <  1:210604/69‑1 (MQ=255)
                                    gTTTCTTTCGCGCCGAGCATATTGCCGTGACGCACTGCGCGACGAATGATACGACGCAGTACATAACCa                         <  1:1521599/69‑1 (MQ=255)
                                    gTTTCTTTCGCGCCGAGCATATTGCCGTGACGCACTGCGCGACGAATGATACGACGCAGTACATAACCa                         <  1:352206/69‑1 (MQ=255)
                                    gTTTCTTTCGCGCCGAGCATATTGCCGTGACGCACTGCGCGACGAATGATACGACGCAGTACATAACCa                         <  1:1746405/69‑1 (MQ=255)
                                          ttcGCGCCGAGCATATTGCCGTGACGCACTGCGCGACGAATGATACGACGCAGTACATAACCa                         <  1:1240153/63‑1 (MQ=255)
                                                        atTACCGTGACGCACTGCGCGACGAATGATACGACGCAGTACATAACCACGGTTTTCATTGGACGGCAt     <  1:672190/69‑1 (MQ=255)
                                                        atTACCGTGACGCACTGCGCGACGAATGATACGACGCAGTACATAACCACGGTTTTCATTGGACGGCAt     <  1:739948/69‑1 (MQ=255)
                                                         ttACCGTGACGCACTGCGCGACGAATGATACGACGCAGTACATAACCACGGTTTTCATTg             >  1:2565052/1‑60 (MQ=255)
                                                         ttACCGTGACGCACTGCGCGACGAATGATACGACGCAGTACATAACCACGGTTTTCATTGGACGGCAtt    >  1:1921664/1‑69 (MQ=255)
                                                         ttACCGTGACGCACTGCGCGACGAATGATACGACGCAGTACATAACCACGGTTTTCATTGGACGGCAtt    >  1:1512698/1‑69 (MQ=255)
                                                         ttACCGTGACGCACTGCGCGACGAATGATACGACGCAGTACATAACCACGGTTTTCATTGGACGGCAtt    >  1:3225683/1‑69 (MQ=255)
                                                         ttACCGTGACGCACTGCGCGACGAATGATACGACGCAGTACATAACCACGGTTTTCATTGGACGGCAtt    >  1:166741/1‑69 (MQ=255)
                                                         ttACCGTGACGCACTGCGCGACGAATGATACGACGCAGTACATAACCACGGTTTTCATTGGACGGCAt     >  1:291032/1‑68 (MQ=255)
                                                           aCCGTGACGCACTGCGCGACGAATGATACGACGCAGTACATAACCACGGTTTTCATTGGACGGCATTa   >  1:820175/1‑68 (MQ=255)
                                                           aCCGTGACGCACTGCGCGACGAATGATACGACGCAGTACATAACCACGGTTTTCATTGGACGGCATTAc  >  1:1538386/1‑69 (MQ=255)
                                                           aCCGTGACGCACTGCGCGACGAATGATACGACGCAGTACATAACCACGGTTTTCATTGGACGGCATTAc  >  1:724396/1‑69 (MQ=255)
                                                           |                                                                    
ACGTCGATCAGCGGACCAACCAGTTTGTAGAAGAAGGTTTCTTTCGCGCCGAGCATATTACCGTGACGCACTGCGCGACGAATGATACGACGCAGTACATAACCACGGTTTTCATTGGACGGCATTAC  >  minE/1814155‑1814282

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑