Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,010,832 T→A 7.1% L87H (CTT→CAT)  rstB → sensory histidine kinase in two‑component regulatory system with RstA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,010,8320TA 7.1% 279.4 / 11.9 98L87H (CTT→CAT) rstBsensory histidine kinase in two‑component regulatory system with RstA
Reads supporting (aligned to +/- strand):  ref base T (52/39);  new base A (2/5);  total (54/44)
Fisher's exact test for biased strand distribution p-value = 2.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 118 total aligned reads displayed.

AAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTACACGTTTTTG  >  minE/1010763‑1010902
                                                                     |                                                                      
aaaCTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCtt                                                                       >  1:2032726/1‑71 (MQ=255)
 aaCTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCtt                                                                       >  1:3580212/1‑70 (MQ=255)
   ctctGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTtgat                                                                    >  1:1577806/1‑71 (MQ=255)
   ctctGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTtgat                                                                    >  1:3890886/1‑71 (MQ=255)
   ctctGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTtgat                                                                    >  1:2736398/1‑71 (MQ=255)
     ctGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTtgatga                                                                  >  1:4095527/1‑71 (MQ=255)
     ctGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTtgatga                                                                  >  1:3804111/1‑71 (MQ=255)
      tGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCtt                                                                       <  1:3280574/65‑1 (MQ=255)
      tGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCtt                                                                       <  1:1960692/65‑1 (MQ=255)
      tGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCtt                                                                       <  1:1587459/65‑1 (MQ=255)
      tGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTtgatga                                                                  >  1:2949320/1‑70 (MQ=255)
         aaGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTgatgat                                                                 <  1:2604526/68‑1 (MQ=255)
         aaGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATAtt                                                              <  1:1244997/71‑1 (MQ=255)
          agagATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATAttt                                                             >  1:2688380/1‑71 (MQ=255)
          agagATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATAttt                                                             >  1:4134216/1‑71 (MQ=255)
          agagATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATAttt                                                             >  1:2144940/1‑71 (MQ=255)
           gagaTGGTTTTAAATCTCTCTTTCGATCTGATTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTc                                                            <  1:1266613/71‑1 (MQ=255)
           gagaTGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTc                                                            <  1:1437397/71‑1 (MQ=255)
           gagaTGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTc                                                            <  1:2309035/71‑1 (MQ=255)
            agaTGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATAttt                                                             >  1:2726941/1‑69 (MQ=255)
              aTGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCtt                                                                       <  1:3657591/57‑1 (MQ=255)
              aTGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCtt                                                                       <  1:399489/57‑1 (MQ=255)
              aTGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCtt                                                                       <  1:3528430/57‑1 (MQ=255)
              aTGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCtt                                                                       <  1:1045858/57‑1 (MQ=255)
              aTGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCtt                                                                       <  1:1522335/57‑1 (MQ=255)
              aTGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCAt                                                         <  1:437486/71‑1 (MQ=255)
                 gATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCa                                                      <  1:1897110/71‑1 (MQ=255)
                 gATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCa                                                      >  1:4039593/1‑71 (MQ=255)
                 gATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCa                                                      >  1:2173630/1‑71 (MQ=255)
                 gATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCa                                                      <  1:426754/71‑1 (MQ=255)
                   tttAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCtt                                                                       <  1:1905379/52‑1 (MQ=255)
                   tttAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCa                                                      <  1:2891462/69‑1 (MQ=255)
                   tttAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCa                                                      <  1:88789/69‑1 (MQ=255)
                   tttAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCAcc                                                    <  1:916102/71‑1 (MQ=255)
                   tttAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCAcc                                                    <  1:3597638/71‑1 (MQ=255)
                   tttAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCAc                                                     >  1:165319/1‑70 (MQ=255)
                    ttATATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATAttt                                                             <  1:1427359/61‑1 (MQ=255)
                    ttAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCAcc                                                    >  1:1418454/1‑70 (MQ=255)
                     tAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATTCCATCTTGATGATATTTCCATGCACCg                                                   <  1:3758567/70‑1 (MQ=255)
                     tAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCg                                                   <  1:3399667/70‑1 (MQ=255)
                          ctctctTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCAcc                                                    >  1:1130767/1‑64 (MQ=255)
                          ctctctTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCAcc                                                    >  1:3396374/1‑64 (MQ=255)
                          ctctctTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCAcc                                                    >  1:794715/1‑64 (MQ=255)
                          ctctctTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCAcc                                                    >  1:1900703/1‑64 (MQ=255)
                            ctctTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTg                                           >  1:4133515/1‑71 (MQ=255)
                             tctTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTgg                                          >  1:1815905/1‑71 (MQ=255)
                               tttCGATCTGCGTGTCGAGCCACTGAGTAAATATCATCAagatgg                                                                  <  1:2375507/45‑2 (MQ=25)
                               tttCGATCTGCGTGTCGAGCCACTGAGTAAATATCATCAagatgg                                                                  <  1:2101237/45‑2 (MQ=25)
                               tttCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCtt                                                                       <  1:2210164/40‑1 (MQ=255)
                               tttCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCtt                                                                       <  1:3393719/40‑1 (MQ=255)
                               tttCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCtt                                                                       <  1:3093674/40‑1 (MQ=255)
                               tttCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCtt                                                                       <  1:4181868/40‑1 (MQ=255)
                               tttCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCtt                                                                       <  1:583491/40‑1 (MQ=255)
                               tttCGATCTGCGTGTCGAGACACTGAGTAAATACCATCtt                                                                       <  1:4193140/40‑1 (MQ=255)
                               tttCGATCTCCGTGTCGAGCCACTGAGTAAATACCATCtt                                                                       <  1:1630160/40‑1 (MQ=255)
                                ttCGATCTGCGTGTCGAGCCACTGAGTAAATATCATCAagatgg                                                                  <  1:988371/44‑2 (MQ=25)
                                 tCGGTCTGCGTGTCGAGCCACTGAGTAAATATCATCAAGATGgta                                                                >  1:2855237/1‑45 (MQ=18)
                                 tCGATCTGCGTGTCGAGCCACTGAGTAAATATCATCAAGATGgta                                                                >  1:2390632/1‑45 (MQ=21)
                                  cGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTggcggc                                      >  1:2425583/1‑70 (MQ=255)
                                      cTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATAttt                                                             >  1:2211669/1‑43 (MQ=255)
                                      cTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCa                                                          >  1:637345/1‑46 (MQ=255)
                                      cTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCa                                                          >  1:1504285/1‑46 (MQ=255)
                                      cTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCa                                                          >  1:68821/1‑46 (MQ=255)
                                      cTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAAt                                 >  1:1127982/1‑71 (MQ=255)
                                         cGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAAt                                 <  1:1538791/68‑1 (MQ=255)
                                               gAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCt                        >  1:403791/1‑71 (MQ=255)
                                               gAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCt                        >  1:2973218/1‑71 (MQ=255)
                                                aGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTg                       <  1:3498856/71‑1 (MQ=255)
                                                aGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTg                       <  1:3420125/71‑1 (MQ=255)
                                                aGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTg                       <  1:2066546/71‑1 (MQ=255)
                                                aGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTg                       <  1:3745728/71‑1 (MQ=255)
                                                aGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTg                       <  1:2316142/71‑1 (MQ=255)
                                                aGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTg                       <  1:4001468/71‑1 (MQ=255)
                                                 gCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTg                       >  1:840732/1‑70 (MQ=255)
                                                 gCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTg                       >  1:4104180/1‑70 (MQ=255)
                                                    aCTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGAc                    >  1:2828165/1‑70 (MQ=255)
                                                     cTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGa                  <  1:2899699/71‑1 (MQ=255)
                                                       gAGTAAATACCATCTTGATGATATTTCCATGCACCGa                                                  >  1:2086650/1‑37 (MQ=255)
                                                       gAGTAAATACCATCTTGATGATATTTCCATGCACCGa                                                  >  1:2195790/1‑37 (MQ=255)
                                                       gAGTAAATACCATCTTGATGATATTTCCATGCACCGa                                                  >  1:2262486/1‑37 (MQ=255)
                                                       gAGTAAATACCATCTTGATGATATTTCCATGCACCGa                                                  >  1:3398054/1‑37 (MQ=255)
                                                       gAGTAAATACCATCTTGATGATATTTCCATGCACCGa                                                  >  1:1512230/1‑37 (MQ=255)
                                                       gAGTAAATACCATCTTGATGATATTTCCATGCACCGa                                                  >  1:3683958/1‑37 (MQ=255)
                                                        aGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGAt                 <  1:3777532/69‑1 (MQ=255)
                                                        aGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGAt                 >  1:1372783/1‑69 (MQ=255)
                                                        aGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGAt                 >  1:1734534/1‑69 (MQ=255)
                                                        aGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATc                <  1:2681424/70‑1 (MQ=255)
                                                        aGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATc                <  1:4143871/70‑1 (MQ=255)
                                                        aGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATc                <  1:695291/70‑1 (MQ=255)
                                                        aGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATc                <  1:160757/70‑1 (MQ=255)
                                                         gTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAg              >  1:2447617/1‑71 (MQ=255)
                                                            aaTACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGc                           >  1:1320987/1‑55 (MQ=255)
                                                            aaTACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTa            >  1:1161707/1‑70 (MQ=255)
                                                            aaTACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTa            >  1:2825192/1‑70 (MQ=255)
                                                            aaTACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTa            >  1:3092976/1‑70 (MQ=255)
                                                            aaTACCATCTTGATGATATATCCATGCACCGACTGCGTg                                           >  1:2966078/1‑39 (MQ=255)
                                                             aTACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTaca          <  1:186883/71‑1 (MQ=255)
                                                                ccATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTACACGt       >  1:1687085/1‑71 (MQ=255)
                                                                     ttGATGATATTTCCATGCACCGACTGCGTGGCGGCGaa                                   >  1:2515489/1‑38 (MQ=255)
                                                                     ttGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTACACGTTTTTg  >  1:276199/1‑71 (MQ=255)
                                                                     |                                                                      
AAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTACACGTTTTTG  >  minE/1010763‑1010902

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: