Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,970,940 C→G 8.0% E656Q (GAA→CAA)  speA ← biosynthetic arginine decarboxylase, PLP‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,970,9400CG 8.0% 330.5 / 12.6 111E656Q (GAA→CAA) speAbiosynthetic arginine decarboxylase, PLP‑binding
Reads supporting (aligned to +/- strand):  ref base C (59/43);  new base G (3/6);  total (62/49)
Fisher's exact test for biased strand distribution p-value = 1.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 120 total aligned reads displayed.

ACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAG  >  minE/1970871‑1971009
                                                                     |                                                                     
acaATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTc                                                                       >  1:341898/1‑70 (MQ=255)
  aaTAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTc                                                                       >  1:3803090/1‑68 (MQ=255)
    tAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGa                                                                   >  1:17758/1‑70 (MQ=255)
    tAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGAt                                                                  <  1:3122103/71‑1 (MQ=255)
    tAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGAt                                                                  <  1:200200/71‑1 (MQ=255)
        ttGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGt                                                              <  1:1362328/71‑1 (MQ=255)
        ttGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGt                                                              <  1:534074/71‑1 (MQ=255)
          gtgGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATaa                                                                >  1:1588477/1‑67 (MQ=255)
          gtgGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGta                                                             >  1:3676112/1‑70 (MQ=255)
          gtgGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGta                                                             >  1:44823/1‑70 (MQ=255)
          gtgGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGta                                                             >  1:627246/1‑70 (MQ=255)
           tgGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGtat                                                            <  1:1750548/70‑1 (MQ=255)
             gCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATaa                                                          <  1:2773499/70‑1 (MQ=255)
              cGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCaa                                                                     <  1:738448/58‑1 (MQ=255)
               ggATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAAcc                                                        <  1:529796/70‑1 (MQ=255)
                taTTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCg                                                       >  1:852155/2‑70 (MQ=255)
                gATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGAt                                                                  <  1:3535888/59‑1 (MQ=255)
                gATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAAcc                                                        >  1:52161/1‑69 (MQ=255)
                gATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAAc                                                         >  1:3070273/1‑68 (MQ=255)
                gATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCg                                                       >  1:2934206/1‑70 (MQ=255)
                 aTTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGAt                                                                  <  1:2480331/58‑1 (MQ=255)
                  ttATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTa                                                     >  1:2197282/1‑70 (MQ=255)
                  ttATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTa                                                     >  1:2137229/1‑70 (MQ=255)
                  ttATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTa                                                     >  1:2143892/1‑70 (MQ=255)
                  ttATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTa                                                     >  1:1169877/1‑70 (MQ=255)
                  ttATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTa                                                     >  1:862892/1‑70 (MQ=255)
                  ttATCACCGCTAAATTAAGCGGATTCAAGCAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTa                                                     >  1:2856310/1‑70 (MQ=255)
                     tCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTGAAGATAAGTATAACCGTTCaa                                                  <  1:2714076/70‑1 (MQ=255)
                     tCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGtat                                                            <  1:1468841/60‑1 (MQ=255)
                     tCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACaa                                                  <  1:969366/70‑1 (MQ=255)
                     tCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACaa                                                  <  1:1400826/70‑1 (MQ=255)
                     tCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACaa                                                  <  1:3733195/70‑1 (MQ=255)
                      cACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAAc                                                >  1:283978/1‑71 (MQ=255)
                       aCCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAAcc                                               <  1:3911257/71‑1 (MQ=255)
                       aCCGCTAAATTAAGCGGATTCAAGTAACAC‑‑GACTTACTCATCTTCAAGATaa                                                                >  1:2213445/1‑52 (MQ=255)
                        ccGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCt                                              <  1:871155/71‑1 (MQ=255)
                        ccGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCt                                              <  1:1579716/71‑1 (MQ=255)
                         cGCTAAATTAAGCGGATTCAAGTACCACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCt                                              >  1:3735059/1‑70 (MQ=255)
                          gCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTTATCTTGAAGAtga                                                                <  1:1070362/51‑3 (MQ=38)
                          gCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTTATCTTGAAGAtga                                                                <  1:1762797/51‑3 (MQ=38)
                          gCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACa                                                   <  1:3246306/64‑1 (MQ=255)
                          gCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGc                                            >  1:1408354/1‑71 (MQ=255)
                           cTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGcc                                           <  1:205824/71‑1 (MQ=255)
                             aaaTTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTc                                         >  1:1951648/1‑71 (MQ=255)
                              aaTTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAAc                                                         >  1:2150572/1‑54 (MQ=255)
                              aaTTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCg                                        <  1:2694898/71‑1 (MQ=255)
                              aaTTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCg                                        <  1:1519248/71‑1 (MQ=255)
                               aTTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTc                                         <  1:2864922/69‑1 (MQ=255)
                               aTTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCg                                        <  1:2711706/70‑1 (MQ=255)
                               aTTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCg                                        <  1:620801/70‑1 (MQ=255)
                                ttAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACaaa                                                 >  1:3868773/1‑60 (MQ=255)
                                ttAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACaaa                                                 >  1:1892324/1‑60 (MQ=255)
                                  aaGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAAct                                    >  1:456728/1‑71 (MQ=255)
                                    gCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAActct                                  <  1:2822201/71‑1 (MQ=255)
                                     cGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCtt                                 <  1:2958140/71‑1 (MQ=255)
                                     cGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCtt                                 <  1:1013174/71‑1 (MQ=255)
                                       gATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTc                                <  1:3810191/70‑1 (MQ=255)
                                       gATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTc                                <  1:3081021/70‑1 (MQ=255)
                                       gATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTc                                >  1:814379/1‑70 (MQ=255)
                                       gATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTCCCTCGAACTCTTc                                >  1:1184730/1‑70 (MQ=255)
                                        aTTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAActc                                   >  1:2655417/1‑66 (MQ=255)
                                        aTTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTc                                >  1:575903/1‑69 (MQ=255)
                                         ttCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAg                             <  1:1106049/71‑1 (MQ=255)
                                         ttCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAg                             <  1:3620488/71‑1 (MQ=255)
                                          tCAAGTAACACAGGACTTACTTATCTTGAAGATGAg                                                               >  1:4089161/1‑36 (MQ=18)
                                          tCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAg                             >  1:521913/1‑70 (MQ=255)
                                             aGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGcc                                           >  1:38910/1‑53 (MQ=255)
                                               tAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACtg                       >  1:1962623/1‑71 (MQ=255)
                                                aaCACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAAcc                                               >  1:3282876/1‑46 (MQ=255)
                                                 acacAGGACTTACTCATCTTCAAGATAAGTATAAc                                                         >  1:352935/1‑35 (MQ=255)
                                                 acacAGGACTTACTCATCTTCAAGATAAGTATAAc                                                         >  1:953380/1‑35 (MQ=255)
                                                 acacAGAACTTACTCATCTTCAAGATAAGTATAAc                                                         >  1:9717/1‑35 (MQ=255)
                                                     aGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGttgttg                 <  1:700902/71‑1 (MQ=255)
                                                     aGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGttgttg                 <  1:2521966/71‑1 (MQ=255)
                                                      ggACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGc                <  1:1967598/71‑1 (MQ=255)
                                                       gACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGttgttg                 >  1:3588070/1‑69 (MQ=255)
                                                       gACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGttgttg                 >  1:520884/1‑69 (MQ=255)
                                                        aCTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAg              <  1:592192/71‑1 (MQ=255)
                                                            aCTTATCTTGAAGATGAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGa                           <  1:31171/54‑1 (MQ=37)
                                                            aCTTATCTTGAAGATGAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGa                           <  1:402953/54‑1 (MQ=37)
                                                            aCTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACtgttg                    >  1:2883665/1‑61 (MQ=255)
                                                            aCTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACtgttg                    >  1:2583895/1‑61 (MQ=255)
                                                            aCTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGtt            >  1:280679/1‑69 (MQ=255)
                                                            aCTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGtt            >  1:4072493/1‑69 (MQ=255)
                                                            aCTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGtt            <  1:2499053/69‑1 (MQ=255)
                                                            aCTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGtt            >  1:1574777/1‑69 (MQ=255)
                                                              ttaTCTTGAAGATGAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAg         >  1:891039/3‑70 (MQ=255)
                                                                aTCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGaa                                      >  1:2336435/1‑39 (MQ=255)
                                                                aTCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGaa                                      >  1:3996891/1‑39 (MQ=255)
                                                                 tCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCAt      >  1:390647/1‑70 (MQ=255)
                                                                   ttCAAGATAAGTATAACCGTACAAACCTGCCTCGaa                                      >  1:651379/1‑36 (MQ=255)
                                                                   ttCAAGATAAGTATAACCGTACAAACCTGCCTCGaa                                      >  1:1385854/1‑36 (MQ=255)
                                                                   ttCAAGATAAGTATAACCGTACAAACCTGCCTCGaa                                      >  1:2854373/1‑36 (MQ=255)
                                                                   ttCAAGATAAGTATAACCGTACAAACCTGCCTCGaa                                      >  1:2138852/1‑36 (MQ=255)
                                                                   ttCAAGATAAGTATAACCGTACAAACCTGCCTCGaa                                      >  1:2747716/1‑36 (MQ=255)
                                                                     cAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATc     >  1:1805667/1‑67 (MQ=255)
                                                                     cAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATc     >  1:1396130/1‑67 (MQ=255)
                                                                     cAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATc     >  1:3247943/1‑67 (MQ=255)
                                                                     cAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAg  <  1:3490540/70‑1 (MQ=255)
                                                                     cAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAg  <  1:3394697/70‑1 (MQ=255)
                                                                     |                                                                     
ACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAG  >  minE/1970871‑1971009

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: