Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,970,946 A→G 5.4% Y654H (TAT→CAT)  speA ← biosynthetic arginine decarboxylase, PLP‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,970,9460AG 5.4% 316.5 / 8.6 110Y654H (TAT→CAT) speAbiosynthetic arginine decarboxylase, PLP‑binding
Reads supporting (aligned to +/- strand):  ref base A (61/43);  new base G (3/3);  total (64/46)
Fisher's exact test for biased strand distribution p-value = 6.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 112 total aligned reads displayed.

TTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATC  >  minE/1970879‑1971012
                                                                   |                                                                  
ttGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGt                                                                 <  1:534074/71‑1 (MQ=255)
ttGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGt                                                                 <  1:1362328/71‑1 (MQ=255)
  gtgGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATaa                                                                   >  1:1588477/1‑67 (MQ=255)
  gtgGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGta                                                                >  1:3676112/1‑70 (MQ=255)
  gtgGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGta                                                                >  1:627246/1‑70 (MQ=255)
  gtgGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGta                                                                >  1:44823/1‑70 (MQ=255)
     gCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATaa                                                             <  1:2773499/70‑1 (MQ=255)
       ggATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAAcc                                                           <  1:529796/70‑1 (MQ=255)
        taTTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCg                                                          >  1:852155/2‑70 (MQ=255)
        gATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAAcc                                                           >  1:52161/1‑69 (MQ=255)
        gATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAAc                                                            >  1:3070273/1‑68 (MQ=255)
        gATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCg                                                          >  1:2934206/1‑70 (MQ=255)
          ttATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTa                                                        >  1:862892/1‑70 (MQ=255)
          ttATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTa                                                        >  1:1169877/1‑70 (MQ=255)
          ttATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTa                                                        >  1:2197282/1‑70 (MQ=255)
          ttATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTa                                                        >  1:2143892/1‑70 (MQ=255)
          ttATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTa                                                        >  1:2137229/1‑70 (MQ=255)
          ttATCACCGCTAAATTAAGCGGATTCAAGCAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTa                                                        >  1:2856310/1‑70 (MQ=255)
           tATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTAc                                                       >  1:516027/1‑70 (MQ=255)
             tCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTGAAGATAAGTATAACCGTTCaa                                                     <  1:2714076/70‑1 (MQ=255)
             tCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGtat                                                               <  1:1468841/60‑1 (MQ=255)
             tCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACaa                                                     <  1:3733195/70‑1 (MQ=255)
             tCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACaa                                                     <  1:969366/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)
                 cGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTg                                                >  1:401566/1‑71 (MQ=255)
                  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)
                      aaTTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTa                                                        >  1:1711650/1‑58 (MQ=255)
                      aaTTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTc                                            >  1:677719/1‑70 (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)
                       aTTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGa                                          <  1:366178/71‑1 (MQ=255)
                        ttAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACaaa                                                    >  1:3868773/1‑60 (MQ=255)
                        ttAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACaaa                                                    >  1:1892324/1‑60 (MQ=255)
                          aaGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAAct                                       >  1:456728/1‑71 (MQ=255)
                             cGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCtt                                    <  1:1013174/71‑1 (MQ=255)
                              ggATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCATACAAACCTGCCTCGAACTCTTc                                   <  1:2205852/71‑1 (MQ=255)
                               gATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTc                                   <  1:3081021/70‑1 (MQ=255)
                               gATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTc                                   <  1:3810191/70‑1 (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)
                                     aGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAAc                            >  1:868974/1‑71 (MQ=255)
                                     aGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAAc                            >  1:2948369/1‑71 (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)
                                             aGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGATCTCTTCAAGGAACTGttgttg                    <  1:3459809/71‑1 (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:2741255/54‑1 (MQ=37)
                                                    aCTTATCTTGAAGATGAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGa                              <  1:402953/54‑1 (MQ=37)
                                                    aCTTATCTTGAAGATGAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGa                              <  1:31171/54‑1 (MQ=37)
                                                    aCTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACtgttg                       >  1:2883665/1‑61 (MQ=255)
                                                    aCTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGtt               <  1:2499053/69‑1 (MQ=255)
                                                    aCTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGtt               >  1:2852962/1‑69 (MQ=255)
                                                    aCTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGtt               >  1:280679/1‑69 (MQ=255)
                                                    aCTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGtt               >  1:1574777/1‑69 (MQ=255)
                                                      ttaTCTTGAAGATGAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAg            >  1:891039/3‑70 (MQ=255)
                                                      ttaTCTTGAAGATGAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAg            >  1:2555796/3‑70 (MQ=255)
                                                        aTCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGaa                                         >  1:3996891/1‑39 (MQ=255)
                                                        aTCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGaa                                         >  1:2336435/1‑39 (MQ=255)
                                                         tCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCAt         >  1:390647/1‑70 (MQ=255)
                                                           ttCAAGATAAGTATAACCGTACAAACCTGCCTCGaa                                         >  1:2854373/1‑36 (MQ=255)
                                                           ttCAAGATAAGTATAACCGTACAAACCTGCCTCGaa                                         >  1:1385854/1‑36 (MQ=255)
                                                           ttCAAGATAAGTATAACCGTACAAACCTGCCTCGaa                                         >  1:651379/1‑36 (MQ=255)
                                                           ttCAAGATAAGTATAACCGTACAAACCTGCCTCGaa                                         >  1:2747716/1‑36 (MQ=255)
                                                           ttCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATc        >  1:1949459/1‑69 (MQ=255)
                                                             cAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGTAACTGTTGTTGCTGTTCAGCATc        >  1:1225576/1‑67 (MQ=255)
                                                             cAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATc        >  1:3239865/1‑67 (MQ=255)
                                                             cAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATc        >  1:1805667/1‑67 (MQ=255)
                                                             cAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATc        >  1:3247943/1‑67 (MQ=255)
                                                             cAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATc        >  1:1396130/1‑67 (MQ=255)
                                                             cAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAg     <  1:3394697/70‑1 (MQ=255)
                                                               aGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATc  <  1:674229/71‑1 (MQ=255)
                                                                gATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATc  <  1:1097951/70‑1 (MQ=255)
                                                                 aTAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGa                              <  1:4013195/41‑1 (MQ=255)
                                                                 aTAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGa                              <  1:1582215/41‑1 (MQ=255)
                                                                   |                                                                  
TTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATCAAGATC  >  minE/1970879‑1971012

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: