Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,341,410 T→C 22.6% V198V (GTA→GTG wcaA ← predicted glycosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,341,4100TC22.6% 158.1 / 38.2 120V198V (GTA→GTGwcaApredicted glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base T (76/16);  new base C (16/11);  total (93/27)
Fisher's exact test for biased strand distribution p-value = 1.75e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.95e-01

Only 100 of 121 total aligned reads displayed.

TTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCC  >  minE/1341346‑1341464
                                                                |                                                      
ttttCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACtt                                                    >  1:2264907/1‑69 (MQ=255)
ttttCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACtt                                                    >  1:43653/1‑69 (MQ=255)
            aGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGAGTCGCCTCTTCTACTTTCCACGGTTCGc                                        <  1:3201888/69‑1 (MQ=255)
              gTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCg                                      >  1:2491763/1‑69 (MQ=255)
              gTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCg                                      >  1:177184/1‑69 (MQ=255)
              gTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCg                                      >  1:1641636/1‑69 (MQ=255)
              gTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCg                                      >  1:2976091/1‑69 (MQ=255)
              gTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCg                                      >  1:1406661/1‑69 (MQ=255)
              gTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCg                                      >  1:443295/1‑69 (MQ=255)
              gTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCg                                      >  1:5475/1‑69 (MQ=255)
                gATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTGCACGGTTCGCCGTa                                    >  1:3423517/1‑69 (MQ=255)
                gATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGt                                     >  1:581454/1‑68 (MQ=255)
                gATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGt                                     >  1:786539/1‑68 (MQ=255)
                gATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTa                                    >  1:976208/1‑69 (MQ=255)
                gATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTa                                    >  1:250933/1‑69 (MQ=255)
                gATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTa                                    >  1:1124606/1‑69 (MQ=255)
                gATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTa                                    >  1:711856/1‑69 (MQ=255)
                gATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTa                                    >  1:635984/1‑69 (MQ=255)
                gATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTa                                    >  1:603062/1‑69 (MQ=255)
                gATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTAGCCTCTTCTACTTTCCACGGTTCGCCGTa                                    >  1:1975037/1‑69 (MQ=255)
                 aTCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTc                                                  >  1:3261401/1‑54 (MQ=255)
                 aTCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGc                                        >  1:1103799/1‑64 (MQ=255)
                 aTCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTa                                    >  1:2061288/1‑68 (MQ=255)
                 aTCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTAc                                   >  1:1207761/1‑69 (MQ=255)
                 aTCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTAc                                   >  1:494606/1‑69 (MQ=255)
                 aTCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTAc                                   >  1:2002160/1‑69 (MQ=255)
                 aTCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTAc                                   >  1:2777104/1‑69 (MQ=255)
                 aTCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTAc                                   >  1:2705143/1‑69 (MQ=255)
                 aTCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTAc                                   >  1:2610203/1‑69 (MQ=255)
                   cTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCtcttct                                                        >  1:1317789/1‑46 (MQ=255)
                   cTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCg                                      >  1:1826691/1‑64 (MQ=255)
                   cTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCg                                      >  1:2955262/1‑64 (MQ=255)
                   cTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCg                                      >  1:635823/1‑64 (MQ=255)
                   cTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCg                                      >  1:2066583/1‑64 (MQ=255)
                   cTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCg                                      >  1:3151705/1‑64 (MQ=255)
                   cTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCg                                      >  1:431753/1‑64 (MQ=255)
                   cTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACt                                  >  1:266610/1‑68 (MQ=255)
                   cTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTc                                 >  1:2194388/1‑69 (MQ=255)
                   cTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTc                                 >  1:293719/1‑69 (MQ=255)
                   cTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTc                                 >  1:2980338/1‑69 (MQ=255)
                   cTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTc                                 >  1:1845306/1‑69 (MQ=255)
                   cTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTc                                 >  1:1561117/1‑69 (MQ=255)
                   cTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTc                                 >  1:2672746/1‑69 (MQ=255)
                   cTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTc                                 >  1:3389231/1‑69 (MQ=255)
                   cTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTc                                 >  1:371430/1‑69 (MQ=255)
                   cTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTc                                 >  1:2513498/1‑69 (MQ=255)
                   cTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTc                                 >  1:518112/1‑69 (MQ=255)
                   cTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTc                                 >  1:640588/1‑69 (MQ=255)
                     gCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTc                                                  >  1:760392/1‑50 (MQ=255)
                     gCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTcc                                >  1:25371/1‑68 (MQ=255)
                     gCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTcc                                >  1:2538770/1‑68 (MQ=255)
                     gCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTcc                                >  1:2602484/1‑68 (MQ=255)
                      cATTTCGCCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCCAc                                                      >  1:1492146/1‑45 (MQ=38)
                      cATTTCGCCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCCACTTTCCACGGTTCGCCg                                      >  1:1429117/1‑61 (MQ=255)
                      cATTTCGCCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCCACTTTCCACGGTTCGCCg                                      >  1:692060/1‑61 (MQ=255)
                      cATTTCGCCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCCACTTTCCACGGTTCGCCg                                      >  1:2230988/1‑61 (MQ=255)
                      cATTTCGCCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCCACTTTCCACGGTTCGCCg                                      >  1:1624842/1‑61 (MQ=255)
                      cATTTCGCCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCCACTTTCCACGGTTCGCCg                                      >  1:1252019/1‑61 (MQ=255)
                      cATTTCGCCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCCACTTTCCACGGTTCGCCg                                      >  1:1236393/1‑61 (MQ=255)
                      cATTTCGCCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCCACTTTCCACGGTTCGCCg                                      >  1:973828/1‑61 (MQ=255)
                      cATTTCGCCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCCACTTTCCACGGTTCGCCg                                      >  1:2811214/1‑61 (MQ=255)
                      cATTTCGCCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCCACTTTCCACGGTTCGCCg                                      >  1:2338075/1‑61 (MQ=255)
                      cATTTCGCCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCCACTTTCCACGGTTCGCCg                                      >  1:1086678/1‑61 (MQ=255)
                      cATTTCGCCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCCACTTTCCACGGTTCGCCg                                      >  1:772460/1‑61 (MQ=255)
                      cATTTCGCCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCCACTTTCCAAGGTTCGCCg                                      >  1:1639953/1‑61 (MQ=38)
                        tttCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTccacc                             >  1:1247072/1‑68 (MQ=255)
                          tCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACgg                                             >  1:767352/1‑50 (MQ=255)
                          tCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTccaccacc                          >  1:2745658/1‑69 (MQ=255)
                          tCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTccaccacc                          >  1:734359/1‑69 (MQ=255)
                          tCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTccaccacc                          >  1:694611/1‑69 (MQ=255)
                          tCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTccaccacc                          >  1:1928415/1‑69 (MQ=255)
                                    tgatATGCAGGATCTGCGTCGCCTCTTCCACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGa                <  1:1847997/69‑1 (MQ=255)
                                    tgatATGCAGGATCTGCGTCGCCTCTTCCACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGa                <  1:1171814/69‑1 (MQ=255)
                                    tgatATGCAGGATCTGCGTCGCCTCTTCCACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGa                <  1:3188167/69‑1 (MQ=255)
                                    tgatATGCAGGATCTGCGTCGCCTCTTCCACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGa                <  1:2074222/69‑1 (MQ=255)
                                    tgatATGCAGGATCTGCGTCGCCTCTTCCACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGa                <  1:1695686/69‑1 (MQ=255)
                                    tgatATGCAGGATCTGCGTCGCCTCTTCCACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGa                <  1:911972/69‑1 (MQ=255)
                                    tgatATGCAGGATCTGCGTCGCCTCTTCCACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGa                <  1:2747667/69‑1 (MQ=255)
                                    tgatATGCAGGATCTGCGTCGCCTCTTCCACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGa                <  1:651213/69‑1 (MQ=255)
                                     gatATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCa                   >  1:885175/1‑65 (MQ=255)
                                     gatATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGaa               <  1:2378183/69‑1 (MQ=255)
                                     gatATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGaa               <  1:692339/69‑1 (MQ=255)
                                     gatATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGaa               <  1:1594794/69‑1 (MQ=255)
                                     gatATGCAGGATCTGCGTCGCCTCTTCCACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGa                <  1:3302098/68‑1 (MQ=255)
                                     gatATGCAGGATCTGCGTCGCCTCTTCCACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGa                <  1:1517668/68‑1 (MQ=255)
                                      atatGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGaa               <  1:110039/68‑1 (MQ=255)
                                                tCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCAt                        >  1:2724643/1‑49 (MQ=255)
                                                tCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTaa     >  1:780945/1‑68 (MQ=255)
                                                tCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAAt    >  1:1282013/1‑69 (MQ=255)
                                                tCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAAt    >  1:3353692/1‑69 (MQ=255)
                                                tCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAAt    >  1:3188524/1‑69 (MQ=255)
                                                tCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAAt    >  1:2517670/1‑69 (MQ=255)
                                                tCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAAt    >  1:2287828/1‑69 (MQ=255)
                                                tCTGCCTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAAt    >  1:3234961/1‑69 (MQ=255)
                                                      tCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGaa               <  1:1339541/52‑1 (MQ=255)
                                                      tCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATcc  <  1:1348034/65‑1 (MQ=255)
                                                      tCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATcc  <  1:1372739/65‑1 (MQ=255)
                                                      tCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATcc  <  1:336750/65‑1 (MQ=255)
                                                      tCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATcc  <  1:1122629/65‑1 (MQ=255)
                                                      tCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATcc  <  1:2097443/65‑1 (MQ=255)
                                                                |                                                      
TTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCC  >  minE/1341346‑1341464

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: