Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,533,152 G→A 36.5% L181L (CTG→TTG)  prmB ← N5‑glutamine methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,533,1520GA36.5% 127.4 / 112.5 183L181L (CTG→TTG) prmBN5‑glutamine methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (63/53);  new base A (46/21);  total (109/74)
Fisher's exact test for biased strand distribution p-value = 6.23e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 184 total aligned reads displayed.

CGTACTGCACTTTCGGCAAGTCGCGGAACAGATCGGAACGAATCGGAATGACGTTGTGGATCAGACCGTGTTCTTCGATGTTCTGTTCAGCAACCGCCAGCGCGTCTGGAGAGATATCCACCGCGTCGAC  >  minE/1533089‑1533218
                                                               |                                                                  
 gTACTGCACTTTCGGCAAGTCGCGGAACAGATCGGAACGAATCGGAATGACGTTGTGGATCAGACCgtg                                                              <  1:2593422/69‑1 (MQ=255)
 gTACTGCACTTTCGGCAAGTCGCGGAACAGATCGGAACGAATCGGAATGACGTTGTGGATCAGACCgtg                                                              <  1:568040/69‑1 (MQ=255)
 gTACTGCACTTTCGGCAAGTCGCGGAACAGATCGGAACGAATCGGAATGACGTTGTGGATCAGACCgtg                                                              <  1:1501175/69‑1 (MQ=255)
 gTACTGCACTTTCGGCAAGTCGCGGAACAGATCGGAACGAATCGGAATGACGTTGTGGATCAGACCgtg                                                              <  1:2907699/69‑1 (MQ=255)
  tACTGCACTTTCGGCAAGTCACGGAACAAGTCGGAACGAATCGGAATGACGTTGTGGATCAAACCgtgt                                                             >  1:2845556/1‑69 (MQ=255)
  tACTGCACTTTCGGCAAGTCACGGAACAAGTCGGAACGAATCGGAATGACGTTGTGGATCAAACCgtgt                                                             >  1:1951751/1‑69 (MQ=255)
  tACTGCACTTTCGGCAAGTCACGGAACAAGTCGGAACGAATCGGAATGACGTTGTGGATCAAACCgtgt                                                             >  1:1733787/1‑69 (MQ=255)
  tACTGCACTTTCGGCAAGTCACGGAACAAGTCGGAACGAATCGGAATGACGTTGTGGATCAAACCgtgt                                                             >  1:2660410/1‑69 (MQ=255)
  tACTGCACTTTCGGCAAGTCACGGAACAAGTCGGAACGAATCGGAATGACGTTGTGGATCAAACCgtgt                                                             >  1:1628321/1‑69 (MQ=255)
  tACTGCACTTTCGGCAAGTCACGGAACAAGTCGGAACGAATCGGAATGACGTTGTGGATCAAACCgtgt                                                             >  1:2757169/1‑69 (MQ=255)
  tACTGCACTTTCGGCAAGTCACGGAACAAGTCGGAACGAATCGGAATGACGTTGTGGATCAAACCgtgt                                                             >  1:1367318/1‑69 (MQ=255)
  tACTGCACTTTCGGCAAGTCACGGAACAAGTCGGAACGAATCGGAATGACGTTGTGGATCAAACCgtgt                                                             >  1:1330850/1‑69 (MQ=255)
  tACTGCACTTTCGGCAAGTCACGGAACAAGTCGGAACGAATCGGAATGACGTTGTGGATCAAACCgtgt                                                             >  1:40769/1‑69 (MQ=255)
  tACTGCACTTTCGGCAAGTCACGGAACAAGTCGGAACGAATCGGAATGACGTTGTGGATCAAACCgtgt                                                             >  1:647969/1‑69 (MQ=255)
  tACTGCACTTTCGGCAAGTCACGGAACAAGTCGGAACGAATCGGAATGACGTTGTGGATCAAACCgtgt                                                             >  1:1319675/1‑69 (MQ=255)
  tACTGCACTTTCGGCAAGTCACGGAACAAGTCGGAACGAATCGGAATGACGTTGTGGATCAAACCgtg                                                              >  1:724065/1‑68 (MQ=255)
    cTGCACTTTCGGCAAGTCGCGGAACAGATCGGAACGAATCGGAATGACGTTGTGGATCAGACCgtg                                                              >  1:1008711/1‑66 (MQ=255)
    cTGCACTTTCGGCAAGTCGCGGAACAGATCGGAACGAATCGGAATGACGTTGTGGATCAGACCgtg                                                              >  1:200870/1‑66 (MQ=255)
    cTGCACTTTCGGCAAGTCGCGGAACAGATCGGAACGAATCGGAATGACGTTGTGGATCAGACCgtg                                                              >  1:3059817/1‑66 (MQ=255)
    cTGCACTTTCGGCAAGTCGCGGAACAGATCGGAACGAATCGGAATGACGTTGTGGATCAGACCgtg                                                              >  1:2264980/1‑66 (MQ=255)
    cTGCACTTTCGGCAAGTCGCGGAACAGATCGGAACGAATCGGAATGACGTTGTGGATCAGACCgtg                                                              >  1:2314940/1‑66 (MQ=255)
    cTGCACTTTCGGCAAGTCGCGGAACAGATCGGAACGAATCGGAATGACGTTGTGGATCAGACCgtg                                                              >  1:1779703/1‑66 (MQ=255)
    cTGCACTTTCGGCAAGTCGCGGAACAGATCGGAACGAATCGGAATGACGTTGTGGATCAGACCgtg                                                              >  1:662171/1‑66 (MQ=255)
    cTGCACTTTCGGCAAGTCGCGGAACAGATCGGAACGAATCGGAATGACGTTGTGGATCAGACCgtg                                                              >  1:685225/1‑66 (MQ=255)
    cTGCACTTTCGGCAAGTCGCGGAACAGATCGGAACGAATCGGAATGACGTTGTGGATCAGACCGTGttc                                                           >  1:165129/1‑69 (MQ=255)
    cTGCACTTTCGGCAAGTCGCGGAACAGATCGGAACGAATCGGAATGACGTTGTGGATCAGACCGTGttc                                                           >  1:1922322/1‑69 (MQ=255)
    cTGCACTTTCGGCAAGTCGCGGAACAGATCGGAACGAATCGGAATGACGTTGTGGATCAGACCGTGttc                                                           >  1:2480535/1‑69 (MQ=255)
    cTGCACTTTCGGCAAGTCGCGGAACAGATCGGAACGAATCGGAATGACGTTGTGGATCAGACCGTGttc                                                           >  1:1716220/1‑69 (MQ=255)
    cTGCACTTTCGGCAAGTCGCGGAACAGATCGGAACGAATCGGAATGACGTTGTGGATCAGACCGTGttc                                                           >  1:1511857/1‑69 (MQ=255)
    cTGCACTTTCGGCAAGTCGCGGAACAGATCGGAACGAATCGGAATGACGTTGTGGATCAGACCGTGttc                                                           >  1:3016816/1‑69 (MQ=255)
    cTGCACTTTCGGCAAGTCGCGGAACAGATCGGAACGAATCGGAATGACGTTGTGGATCAGACCGTGtt                                                            >  1:1815499/1‑68 (MQ=255)
    cTGCACTTTCGGCAAGTCGCGGAACAGATCGGAACGAATCGGAATGACGTTGTGGATCAGACCGTGtt                                                            >  1:3286386/1‑68 (MQ=255)
    cTGCACTTTCGGAAGTCGCGGAACAGATCGGAACGAATCGGAATGACGTTGTGGATCAGACCgt                                                               >  1:553202/1‑64 (MQ=255)
       cACTTTCGGCAAGTCGCGGAACAGATCGGAACGAATCGGAATGACGTTGTGGATCAGACCGTGttcttc                                                        >  1:2805708/1‑69 (MQ=255)
       cACTTTCGGCAAGTCGCGGAACAGATCGGAACGAATCGGAATGACGTTGTGGATCAGACCGTGttcttc                                                        >  1:3245693/1‑69 (MQ=255)
         cTTTCGGCAAGTCGCGGAACAGATCGGAACGAATCGGAATGACGTTGTGGATCAGACCGTGttcttc                                                        <  1:954796/67‑1 (MQ=255)
         cTTTCGGCAAGTCGCGGAACAGATCGGAACGAATCGGAATGACGTTGTGGATCAGACCGTGttcttc                                                        <  1:2467392/67‑1 (MQ=255)
         cTTTCGGCAAGTCGCGGAACAGATCGGAACGAATCGGAATGACGTTGTGGATCAGACCGTGttcttc                                                        <  1:3316758/67‑1 (MQ=255)
         cTTTCGGCAAGTCGCGGAACAGATCGGAACGAATCGGAATGACGTTGTGGATCAGACCGTGttcttc                                                        <  1:1277913/67‑1 (MQ=255)
         cTTTCGGCAAGTCGCGGAACAGATCGGAACGAATCGGAATGACGTTGTGGATCAGACCGTGttcttc                                                        <  1:1139468/67‑1 (MQ=255)
         cTTTCCGCAAGTCGCGGAACAGATCGGAACGAATCGGAATGACGTTGTGGATCAGACCGTGttcttc                                                        <  1:2805414/67‑1 (MQ=255)
           ttCGGCAAGTCGCGGAACAGATCGGAACGAATCGGAATGACGTTGTGGATCAGACCGTGTTCTTCGAtg                                                    >  1:1909891/1‑69 (MQ=255)
           ttCGGCAAGTCGCGGAACAGATCGGAACGAATCGGAATGACGTTGTGGATCAGACCGTGTTCTTCGAtg                                                    >  1:1416907/1‑69 (MQ=255)
           ttCGGCAAGTCGCGGAACAGATCGGAACGAATCGGAATGACGTTGTGGATCAGACCGTGTTCTTCGAtg                                                    >  1:2689944/1‑69 (MQ=255)
           ttCGGCAAGTCGCGGAACAGATCGGAACGAATCGGAATGACGTTGTGGATCAGACCGTGTTCTTCGAtg                                                    >  1:297130/1‑69 (MQ=255)
           ttCGGCAAGTCGCGGAACAGATCGGAACGAATCGGAATGACGTTGTGGATCAGACCGTGTTCTTCGAtg                                                    >  1:2853453/1‑69 (MQ=255)
           ttCGGCAAGTCGCGGAACAGATCGGAACGAATCGGAATGACGTTGTGGATCAGACCGTGTTCTTCGAt                                                     >  1:2696370/1‑68 (MQ=255)
            tCGGCAAGTCACGGAACAAGTCGGAACGAATCGGAATGACGTTGTGGATCAAACCGTGTTCTTCGAtgt                                                   >  1:2285258/1‑69 (MQ=255)
            tCGGCAAGTCACGGAACAAGTCGGAACGAATCGGAATGACGTTGTGGATCAAACCGTGTTCTTCGAtgt                                                   >  1:1202990/1‑69 (MQ=255)
                     cgcgGAACAGATCGGAACGAATCGGAATGACGTTGTGGATCAGACCGTGTTCTTCGATGTTCTGTTCa                                           >  1:2172741/1‑68 (MQ=255)
                      acgGAACAAGTCGGAACGAATCGGAATGACGTTGTGGATCAAACCGTGTTCTTCGATGTTCTGTTCAg                                          >  1:1102662/2‑68 (MQ=255)
                      acgGAACAAGTCGGAACGAATCGGAATGACGTTGTGGATCAAACCGTGTTCTTCGATGTTCTGTTCAGc                                         >  1:2786847/2‑69 (MQ=255)
                      acgGAACAAGTCGGAACGAATCGGAATGACGTTGTGGATCAAACCGTGTTCTTCGATGTTCTGTTCAGc                                         >  1:767846/2‑69 (MQ=255)
                      acgGAACAAGTCGGAACGAATCGGAATGACGTTGTGGATCAAACCGTGTTCTTCGATGTTCTGTTCAGc                                         >  1:252536/2‑69 (MQ=255)
                      acgGAACAAGTCGGAACGAATCGGAATGACGTTGTGGATCAAACCGTGTTCTTCGATGTTCTGTTCAGc                                         >  1:2499757/2‑69 (MQ=255)
                      acgGAACAAGTCGGAACGAATCGGAATGACGTTGTGGATCAAACCGTGTTCTTCGATGTTCTGTTCAGc                                         >  1:598247/2‑69 (MQ=255)
                      acgGAACAAGTCGGAACGAATCGGAATGACGTTGTGGATCAAACCGTGTTCTTCGATGTTCTGTTCAGc                                         >  1:2339952/2‑69 (MQ=255)
                      acgGAACAAGTCGGAACGAATCGGAATGACGTTGTGGATCAAACCGTGTTCTTCGATGTTCTGTTCAGc                                         >  1:1208666/2‑69 (MQ=255)
                        ggAACAGATCGGAACGAATCGGAATGACGTTGTGGATCAGACCGTGTTCTTCGATGTTCTGTTCAGCa                                        <  1:3169433/68‑1 (MQ=255)
                        ggAACAGATCGGAACGAATCGGAATGACGTTGTGGATCAGACCGTGTTCTTCGATGTTCTGTTCAGCa                                        <  1:479611/68‑1 (MQ=255)
                        ggAACAGATCGGAACGAATCGGAATGACGTTGTGGATCAGACCGTGTTCTTCGATGTTCTGTTCAGCa                                        <  1:3333168/68‑1 (MQ=255)
                        ggAACAGATCGGAACGAATCGGAATGACGTTGTGGATCAGACCGTGTTCTTCGATGTTCTGTTCAGCa                                        <  1:339804/68‑1 (MQ=255)
                        ggAACAGATCGGAACGAATCGGAATGACGTTGTGGATCAGACCGTGTTCTTCGATGTTCTGTTCAGCa                                        <  1:1864424/68‑1 (MQ=255)
                        ggAACAGATCGGAACGAATCGGAATGACGTTGTGGATCAGACCGTGTTCTTCGATGTTCTGTTCAGCa                                        <  1:189935/68‑1 (MQ=255)
                           aCAGATCGGAACGAATCGGAATGACGTTGTGGATCAGACCGTGTTCTTCGATGTTCTGTTCAGCAGcc                                     >  1:53936/1‑68 (MQ=255)
                           aCAGATCGGAACGAATCGGAATGACGTTGTGGATCAGACCGTGTTCTTCGATGTTCTGTTCAGCAAcc                                     >  1:195910/1‑68 (MQ=255)
                               aTCGGAACGAATCGGAATGACGTTGTGGATCAGACCGTGTTCTTCGATGTTCTGTTCAGCAACCGCCAg                                <  1:724333/69‑1 (MQ=255)
                               aTCGGAACGAATCGGAATGACGTTGTGGATCAGACCGTGTTCTTCGATGTTCTGTTCAGCAACCGCCAg                                <  1:379388/69‑1 (MQ=255)
                               aTCGGAACGAATCGGAATGACGTTGTGGATCAGACCGTGTTCTTCGATGTTCTGTTCAGCAACCGCCAg                                <  1:2037400/69‑1 (MQ=255)
                                tCGGAACGAATCGGAATGACGTTGTGGATCAGACCGTGTTCTTCGATGTTCTGTTCAGCAACCGCCAg                                <  1:2350741/68‑1 (MQ=255)
                                   gAACGAATCGGAATGACGTTGTGGATCAGACCGTGTTCTTCGATGTTCTGTTCAGCAACCGCCAg                                <  1:2073527/65‑1 (MQ=255)
                                     aCGAATCGGAATGACGTTGTGGATCAGACCGTGTTCTTCGATGTTCTGTTCAGCAACCGCCAGCGCGTc                          >  1:3157003/1‑69 (MQ=255)
                                     aCGAATCGGAATGACGTTGTGGATCAGACCGTGTTCTTCGATGTTCTGTTCAGCAACCGCCAGCGCGTc                          >  1:1006955/1‑69 (MQ=255)
                                     aCGAATCGGAATGACGTTGTGGATCAGACCGTGTTCTTCGATGTTCTGTTCAGCAACCGCCAGCGCGTc                          >  1:799569/1‑69 (MQ=255)
                                     aCGAATCGGAATGACGTTGTGGATCAGACCGTGTTCTTCGATGTTCTGTTCAGCAACCGCCAGCGCGTc                          >  1:621499/1‑69 (MQ=255)
                                          tCGGAATGACGTTGTGGATCAAACCGTGTTCTTCGATGTTCTGTTCAGCAACCGCCAGCGCGTCTGgag                     <  1:866105/69‑1 (MQ=255)
                                          tCGGAATGACGTTGTGGATCAAACCGTGTTCTTCGATGTTCTGTTCAGCAACCGCCAGCGCGTCTGgag                     <  1:1113898/69‑1 (MQ=255)
                                          tCGGAATGACGTTGTGGATCAAACCGTGTTCTTCGATGTTCTGTTCAGCAACCGCCAGCGCGTCTGgag                     <  1:415687/69‑1 (MQ=255)
                                          tCGGAATGACGTTGTGGATCAAACCGTGTTCTTCGATGTTCTGTTCAGCAACCGCCAGCGCGTCTGgag                     <  1:1026492/69‑1 (MQ=255)
                                          tCGGAATGACGTTGTGGATCAAACCGTGTTCTTCGATGTTCTGTTCAGCAACCGCCAGCGCGTCTGgag                     <  1:95644/69‑1 (MQ=255)
                                          tCGGAATGACGTTGTGGATCAAACCGTGTTCTTCGATGTTCTGTTCAGCAACCGCCAGCGCGTCTGgag                     <  1:1268833/69‑1 (MQ=255)
                                          tCGGAATGACGTTGTGGATCAAACCGTGTTCTTCGATGTTCTGTTCAGCAACCGCCAGCGCGTCTGgag                     <  1:1385205/69‑1 (MQ=255)
                                          tCGGAATGACGTTGTGGATCAAACCGTGTTCTTCGATGTTCTGTTCAGCAACCGCCAGCGCGTCTGgag                     <  1:3014922/69‑1 (MQ=255)
                                          tCGGAATGACGTTGTGGATCAAACCGTGTTCTTCGATGTTCTGTTCAGCAACCGCCAGCGCGTCTGgag                     <  1:1435967/69‑1 (MQ=255)
                                          tCGGAATGACGTTGTGGATCAAACCGTGTTCTTCGATGTTCTGTTCAGCAACCGCCAGCGCGTCTGgag                     <  1:2741722/69‑1 (MQ=255)
                                          tCGGAATGACGTTGTGGATCAAACCGTGTTCTTCGATGTTCTGTTCAGCAACCGCCAGCGCGTCTGgag                     <  1:1569188/69‑1 (MQ=255)
                                              aaTGACGTTGTGGATCAGACCGTGTTCTTCGATGTTCTGTTCAGCAACCGCCAGCGCGTCTGGAGAGat                 <  1:3387266/69‑1 (MQ=255)
                                              aaTGACGTTGTGGATCAGACCGTGTTCTTCGATGTTCTGTTCAGCAACCGCCAGCGCGTCTGGAGAGat                 <  1:2179153/69‑1 (MQ=255)
                                              aaTGACGTTGTGGATCAGACCGTGTTCTTCGATGTTCTGTTCAGCAACCGCCAGCGCGTCTGGAGAGat                 <  1:3203512/69‑1 (MQ=255)
                                              aaTGACGTTGTGGATCAGACCGTGTTCTTCGATGTTCTGTTCAGCAACCGCCAGCGCGTCTGGAGAGat                 <  1:1679712/69‑1 (MQ=255)
                                              aaTGACGTTGTGGATCAGACCGTGTTCTTCGATGTTCTGTTCAGCAACCGCCAGCGCGTCTGGAGAGat                 <  1:1429492/69‑1 (MQ=255)
                                               aTGACGTTGTGGATCAGACCGTGTTCTTCGATGTTCTGTTCAGCAACCGCCAGCGCGTCTGGAGAGata                <  1:1930666/69‑1 (MQ=255)
                                               aTGACGTTGTGGATCAGACCGTGTTCTTCGATGTTCTGTTCAGCAACCGCCAGCGCGTCTGGAGAGata                <  1:1354699/69‑1 (MQ=255)
                                                 gACGTTGTGGATCAGACCGTGTTCTTCGATGTTCTGTTCAGCAACCGCCAGCGCGTCTGGAGAGATATc              <  1:1539753/69‑1 (MQ=255)
                                                  aCGTTGTGGATCAGACCGTGTTCTTCGATGTTCTGTTCAGCAACCGCCAGCGCGTCTGGAGAGATATcc             >  1:1039698/1‑69 (MQ=255)
                                                  aCGTTGTGGATCAGACCGTGTTCTTCGATGTTCTGTTCAGCAACCGCCAGCGCGTCTGGAGAGATATcc             >  1:2736332/1‑69 (MQ=255)
                                                  aCGTTGTGGATCAGACCGTGTTCTTCGATGTTCTGTTCAGCAACCGCCAGCGCGTCTGGAGAGATATcc             >  1:1695828/1‑69 (MQ=255)
                                                             cAGACCGTGTTCTTCGATGTTCTGTTCAGCAACCGCCAGCGCGTCTGGAGAGATATCCACCGCGTCGa   >  1:3238222/1‑68 (MQ=255)
                                                             cAGACCGTGTTCTTCGATGTTCTGTTCAGCAACCGCCAGCGCGTCTGGAGAGATATCCACCGCGTCGAc  >  1:1155707/1‑69 (MQ=255)
                                                               gACCGTGTTCTTCGATGTTCGGTTCAGCAACCGCCAGCGCGTCTGGAGAGATATCCAc           <  1:205505/58‑1 (MQ=255)
                                                               |                                                                  
CGTACTGCACTTTCGGCAAGTCGCGGAACAGATCGGAACGAATCGGAATGACGTTGTGGATCAGACCGTGTTCTTCGATGTTCTGTTCAGCAACCGCCAGCGCGTCTGGAGAGATATCCACCGCGTCGAC  >  minE/1533089‑1533218

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