Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,970,934 C→T 8.0% E658K (GAG→AAG)  speA ← biosynthetic arginine decarboxylase, PLP‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,970,9340CT 8.0% 364.5 / 14.7 112E658K (GAG→AAG) speAbiosynthetic arginine decarboxylase, PLP‑binding
Reads supporting (aligned to +/- strand):  ref base C (47/56);  new base T (5/4);  total (52/60)
Fisher's exact test for biased strand distribution p-value = 7.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 114 total aligned reads displayed.

ATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCA  >  minE/1970864‑1971001
                                                                      |                                                                   
aTTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTc                                                                     <  1:2407047/71‑1 (MQ=255)
   tGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATc                                                                  >  1:2021219/1‑71 (MQ=255)
     tCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCt                                                                 <  1:342051/70‑1 (MQ=255)
      cacaATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCtt                                                                >  1:1506936/1‑70 (MQ=255)
       acaATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTc                                                               <  1:2988995/70‑1 (MQ=255)
       acaATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTc                                                               <  1:2340619/70‑1 (MQ=255)
       acaATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTc                                                               <  1:2662975/70‑1 (MQ=255)
         aaTAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTc                                                               >  1:790751/1‑68 (MQ=255)
          aTAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAg                                                            >  1:3044521/1‑70 (MQ=255)
          aTAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAg                                                            >  1:3696189/1‑70 (MQ=255)
          aTAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAg                                                            >  1:1469805/1‑70 (MQ=255)
           tAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTc                                                               <  1:3184821/66‑1 (MQ=255)
           tAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTc                                                               <  1:3551628/66‑1 (MQ=255)
           tAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGAt                                                          <  1:1281996/71‑1 (MQ=255)
             aaTTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGAt                                                          >  1:1675546/1‑69 (MQ=255)
              aTTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAg                                                       <  1:2242609/71‑1 (MQ=255)
                  tgGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGtat                                                    <  1:3644119/70‑1 (MQ=255)
                     cGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTc                                                               <  1:2420041/56‑1 (MQ=255)
                       gATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCtt                                                                <  1:3040808/53‑1 (MQ=255)
                       gATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCtt                                                                <  1:2117574/53‑1 (MQ=255)
                        aTTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTa                                             <  1:231975/71‑1 (MQ=255)
                        aTTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTa                                             <  1:1058271/71‑1 (MQ=255)
                          tATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCtt                                                                <  1:2000543/50‑1 (MQ=255)
                          tATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTAc                                            >  1:2046323/1‑70 (MQ=255)
                          tATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTAc                                            >  1:2579035/1‑70 (MQ=255)
                           aTCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACa                                           <  1:1896854/70‑1 (MQ=255)
                           aTCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACa                                           <  1:3018315/70‑1 (MQ=255)
                           aTCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACa                                           >  1:3382546/1‑70 (MQ=255)
                           aTCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACa                                           >  1:2296689/1‑70 (MQ=255)
                            tCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACaa                                          <  1:3530000/70‑1 (MQ=255)
                            tCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACaa                                          <  1:3255953/70‑1 (MQ=255)
                              aCCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTTATCTTGAAGATGAGtaa                                                    >  1:2332778/1‑57 (MQ=37)
                              aCCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTTATATTGAAGATGAGtaa                                                    >  1:3217821/1‑57 (MQ=37)
                              aCCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAAc                                        <  1:535740/70‑1 (MQ=255)
                               ccGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCt                                      <  1:3031858/71‑1 (MQ=255)
                                cGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCg                                               <  1:3192478/61‑1 (MQ=255)
                                 gCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGc                                    >  1:2295436/1‑71 (MQ=255)
                                    aaaTTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGa                                                           >  1:3379663/1‑45 (MQ=255)
                                     aaTTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATaa                                                  >  1:1348056/1‑53 (MQ=255)
                                     aaTTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATaa                                                  >  1:400234/1‑53 (MQ=255)
                                     aaTTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATaa                                                  >  1:341813/1‑53 (MQ=255)
                                     aaTTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTc                                 >  1:2152490/1‑70 (MQ=255)
                                     aaTTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTc                                 >  1:1006693/1‑70 (MQ=255)
                                     aaTTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGAAAAGTATAAACGTACAAACCTGCCt                                  >  1:1234867/1‑69 (MQ=255)
                                     aaTTAAGAGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTc                                 >  1:1445415/1‑70 (MQ=255)
                                        tAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAAc                             >  1:1497443/1‑71 (MQ=255)
                                        tAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAAc                             >  1:2719551/1‑71 (MQ=255)
                                         aaGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAAct                            >  1:3735861/1‑71 (MQ=255)
                                         aaGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAAct                            >  1:666345/1‑71 (MQ=255)
                                          aGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTACCTCGAActc                           <  1:623457/71‑1 (MQ=255)
                                           gCGGATTCAAGTAACACATGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAActct                          <  1:1861423/71‑1 (MQ=255)
                                           gCGGATTCAAGTAACACAGGACTTACTTATCTTGAAGAtga                                                        <  1:1312115/41‑3 (MQ=25)
                                           gCGGATTCAAGTAACACAGGACTTACTTATCTTGAAGAtga                                                        <  1:2310047/41‑3 (MQ=25)
                                           gCGGATTCAAGTAACACAGGACTTACTTATCTTGAAGAtga                                                        <  1:2306144/41‑3 (MQ=25)
                                           gCGGATTCAAGTAACACAGGACTTACTTATCTTGAAGAtga                                                        <  1:2283439/41‑3 (MQ=25)
                                           gCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCTAActct                          <  1:496470/71‑1 (MQ=255)
                                           gCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAActct                          <  1:1659772/71‑1 (MQ=255)
                                           gCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAActct                          <  1:3151293/71‑1 (MQ=255)
                                           gCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAActct                          >  1:3650840/1‑71 (MQ=255)
                                             ggATTCAAGTAACACAGGACTTACTTATCTTGAAGATGAGt                                                      >  1:2672116/1‑41 (MQ=21)
                                             ggATTCAAGTAACACAGGACTTACTTATCTTGAAGATGAGt                                                      >  1:1249756/1‑41 (MQ=21)
                                             ggATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTc                        <  1:461489/71‑1 (MQ=255)
                                             ggATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTc                        <  1:3001148/71‑1 (MQ=255)
                                             ggATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTc                        <  1:3807977/71‑1 (MQ=255)
                                             ggATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTc                        <  1:654640/71‑1 (MQ=255)
                                                ttCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAg                     >  1:3270314/1‑71 (MQ=255)
                                                ttCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAg                     >  1:3820288/1‑71 (MQ=255)
                                                ttCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAg                     >  1:3718727/1‑71 (MQ=255)
                                                ttCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAg                     >  1:2494551/1‑71 (MQ=255)
                                                ttCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAg                     >  1:1623354/1‑71 (MQ=255)
                                                   aaGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCaa                      >  1:2019348/1‑67 (MQ=255)
                                                   aaGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCaa                      >  1:824290/1‑67 (MQ=255)
                                                   aaGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCaa                      >  1:1432688/1‑67 (MQ=255)
                                                   aaGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCaa                      >  1:3567843/1‑67 (MQ=255)
                                                   aaGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGaa                  >  1:1202950/1‑71 (MQ=255)
                                                   aaGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGaa                  >  1:3811060/1‑71 (MQ=255)
                                                   aaGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGa                   <  1:237469/70‑1 (MQ=255)
                                                   aaGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGa                   <  1:3818586/70‑1 (MQ=255)
                                                    aGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAAc                 <  1:3643259/71‑1 (MQ=255)
                                                    aGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAAc                 <  1:2766523/71‑1 (MQ=255)
                                                    aGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAAc                 <  1:229165/71‑1 (MQ=255)
                                                     gTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCg                                <  1:3121767/55‑1 (MQ=255)
                                                      tAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACtg               >  1:201976/1‑71 (MQ=255)
                                                        acacAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACtgtt             >  1:2859732/1‑71 (MQ=255)
                                                        acacAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACtgt              >  1:1155473/1‑70 (MQ=255)
                                                          acaGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCaa                      <  1:2281834/60‑1 (MQ=255)
                                                          acaGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACtgttgt           <  1:2356353/71‑1 (MQ=255)
                                                             ggACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGTAACtgt              <  1:240071/65‑1 (MQ=255)
                                                             ggACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACtgt              <  1:2316500/65‑1 (MQ=255)
                                                             ggACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACtgt              <  1:1174824/65‑1 (MQ=255)
                                                             ggACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGttgttg         >  1:1523663/1‑70 (MQ=255)
                                                              gACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACtgtt             <  1:507394/65‑1 (MQ=255)
                                                                cTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGt     <  1:1584352/71‑1 (MQ=255)
                                                                cTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGt     <  1:3768153/71‑1 (MQ=255)
                                                                 ttACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGt     <  1:214826/70‑1 (MQ=255)
                                                                  tACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCg                                <  1:2809447/42‑1 (MQ=255)
                                                                  tACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCg                                <  1:1366926/42‑1 (MQ=255)
                                                                   aCTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGaa                  >  1:2675561/1‑55 (MQ=255)
                                                                    cTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCa  <  1:1960283/70‑1 (MQ=255)
                                                                      cATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCg                                >  1:1531341/1‑38 (MQ=255)
                                                                      |                                                                   
ATTTGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCA  >  minE/1970864‑1971001

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: