Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,010,827 C→T 6.0% Y85Y (TAC→TAT 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,8270CT 6.0% 383.3 / 9.3 117Y85Y (TAC→TATrstBsensory histidine kinase in two‑component regulatory system with RstA
Reads supporting (aligned to +/- strand):  ref base C (55/55);  new base T (2/5);  total (57/60)
Fisher's exact test for biased strand distribution p-value = 4.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 120 total aligned reads displayed.

GGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTACACGT  >  minE/1010758‑1010897
                                                                     |                                                                      
gggTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATAc                                                                        >  1:2729867/1‑70 (MQ=255)
     aaaCTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCtt                                                                  >  1:2032726/1‑71 (MQ=255)
      aaCTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCtt                                                                  >  1:3580212/1‑70 (MQ=255)
        ctctGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATc                                                                    >  1:901395/1‑66 (MQ=255)
        ctctGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATc                                                                    >  1:2696896/1‑66 (MQ=255)
        ctctGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTtgat                                                               >  1:3282860/1‑71 (MQ=255)
        ctctGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTtgat                                                               >  1:2736398/1‑71 (MQ=255)
        ctctGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTtgat                                                               >  1:3890886/1‑71 (MQ=255)
          ctGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTtgatga                                                             >  1:4095527/1‑71 (MQ=255)
          ctGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTtgatga                                                             >  1:3804111/1‑71 (MQ=255)
           tGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCtt                                                                  <  1:1587459/65‑1 (MQ=255)
           tGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCtt                                                                  <  1:1960692/65‑1 (MQ=255)
           tGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCtt                                                                  <  1:3280574/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:2144940/1‑71 (MQ=255)
               agagATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATAttt                                                        >  1:2688380/1‑71 (MQ=255)
                gagaTGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTc                                                       <  1:2309035/71‑1 (MQ=255)
                 agaTGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATAttt                                                        >  1:2726941/1‑69 (MQ=255)
                   aTGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCtt                                                                  <  1:1522335/57‑1 (MQ=255)
                   aTGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCtt                                                                  <  1:399489/57‑1 (MQ=255)
                   aTGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCtt                                                                  <  1:3657591/57‑1 (MQ=255)
                   aTGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCtt                                                                  <  1:3528430/57‑1 (MQ=255)
                   aTGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCtt                                                                  <  1:1045858/57‑1 (MQ=255)
                   aTGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCAt                                                    <  1:4230779/71‑1 (MQ=255)
                      gATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCa                                                 >  1:4039593/1‑71 (MQ=255)
                      gATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCa                                                 >  1:2173630/1‑71 (MQ=255)
                      gATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCa                                                 >  1:3437959/1‑71 (MQ=255)
                      gATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCa                                                 <  1:1897110/71‑1 (MQ=255)
                      gATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCa                                                 <  1:426754/71‑1 (MQ=255)
                        tttAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCtt                                                                  <  1:1905379/52‑1 (MQ=255)
                        tttAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCa                                                 <  1:88789/69‑1 (MQ=255)
                        tttAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCa                                                 <  1:2891462/69‑1 (MQ=255)
                        tttAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCAcc                                               <  1:916102/71‑1 (MQ=255)
                        tttAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCAc                                                >  1:165319/1‑70 (MQ=255)
                         ttATATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATAttt                                                        <  1:1427359/61‑1 (MQ=255)
                         ttAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATAttt                                                        <  1:3916122/61‑1 (MQ=255)
                         ttAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATAttt                                                        <  1:2599944/61‑1 (MQ=255)
                         ttAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCAcc                                               >  1:1418454/1‑70 (MQ=255)
                          tAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATTCCATCTTGATGATATTTCCATGCACCg                                              <  1:3758567/70‑1 (MQ=255)
                               ctctctTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCAcc                                               >  1:1130767/1‑64 (MQ=255)
                               ctctctTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCAcc                                               >  1:3396374/1‑64 (MQ=255)
                                 ctctTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTg                                      >  1:4118/1‑71 (MQ=255)
                                 ctctTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTg                                      >  1:4133515/1‑71 (MQ=255)
                                  tctTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTgg                                     >  1:1815905/1‑71 (MQ=255)
                                   ctTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTggc                                    >  1:715824/1‑71 (MQ=255)
                                    tttCGATCTGCGTGTCGAGCCACTGAGTAAATATCATCAagatgg                                                             <  1:1277065/45‑2 (MQ=25)
                                    tttCGATCTGCGTGTCGAGCCACTGAGTAAATATCATCAagatgg                                                             <  1:1510157/45‑2 (MQ=25)
                                    tttCGATCTGCGTGTCGAGCCACTGAGTAAATATCATCAagatgg                                                             <  1:2101237/45‑2 (MQ=25)
                                    tttCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCtt                                                                  <  1:583491/40‑1 (MQ=255)
                                    tttCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCtt                                                                  <  1:3093674/40‑1 (MQ=255)
                                    tttCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCtt                                                                  <  1:4181868/40‑1 (MQ=255)
                                    tttCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCtt                                                                  <  1:3393719/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)
                                      tCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTtga                                                                >  1:201095/1‑40 (MQ=255)
                                       cGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTggcggc                                 >  1:2425583/1‑70 (MQ=255)
                                           cTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATAttt                                                        >  1:2211669/1‑43 (MQ=255)
                                           cTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCa                                                     >  1:637345/1‑46 (MQ=255)
                                           cTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCa                                                     >  1:3721180/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)
                                              cGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGa                               <  1:1852828/65‑1 (MQ=255)
                                              cGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAAtt                           <  1:1603708/69‑1 (MQ=255)
                                              cGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAAt                            <  1:1538791/68‑1 (MQ=255)
                                                    gAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCt                   >  1:2973218/1‑71 (MQ=255)
                                                     aGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTg                  <  1:3420125/71‑1 (MQ=255)
                                                     aGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTg                  <  1:2187551/71‑1 (MQ=255)
                                                     aGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTg                  <  1:3498856/71‑1 (MQ=255)
                                                     aGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTg                  <  1:2316142/71‑1 (MQ=255)
                                                     aGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTg                  <  1:4001468/71‑1 (MQ=255)
                                                     aGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTg                  <  1:2066546/71‑1 (MQ=255)
                                                     aGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTg                  <  1:3745728/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)
                                                          cTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGa             <  1:2315262/71‑1 (MQ=255)
                                                            gAGTAAATACCATCTTGATGATATTTCCATGCACCGa                                             >  1:1512230/1‑37 (MQ=255)
                                                            gAGTAAATACCATCTTGATGATATTTCCATGCACCGa                                             >  1:2262486/1‑37 (MQ=255)
                                                            gAGTAAATACCATCTTGATGATATTTCCATGCACCGa                                             >  1:2195790/1‑37 (MQ=255)
                                                            gAGTAAATACCATCTTGATGATATTTCCATGCACCGa                                             >  1:3683958/1‑37 (MQ=255)
                                                            gAGTAAATACCATCTTGATGATATTTCCATGCACCGa                                             >  1:3398054/1‑37 (MQ=255)
                                                             aGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGAt            >  1:1734534/1‑69 (MQ=255)
                                                             aGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGAt            >  1:1372783/1‑69 (MQ=255)
                                                             aGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATc           <  1:4143871/70‑1 (MQ=255)
                                                             aGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATc           <  1:160757/70‑1 (MQ=255)
                                                             aGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATc           <  1:695291/70‑1 (MQ=255)
                                                              gTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAg         >  1:2447617/1‑71 (MQ=255)
                                                                 aaTACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGc                      >  1:1320987/1‑55 (MQ=255)
                                                                 aaTACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTa       >  1:2825192/1‑70 (MQ=255)
                                                                 aaTACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTa       >  1:3092976/1‑70 (MQ=255)
                                                                 aaTACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTa       >  1:1161707/1‑70 (MQ=255)
                                                                 aaTACCATCTTGATGATATATCCATGCACCGACTGCGTg                                      >  1:2966078/1‑39 (MQ=255)
                                                                  aTACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTaca     <  1:186883/71‑1 (MQ=255)
                                                                     |                                                                      
GGGTAAAACTCTGAAAGAGATGGATTTAAATCTCTCTTTCGATCTGCGTGTCGAGCCACTGAGTAAATACCATCTTGATGATATTTCCATGCACCGACTGCGTGGCGGCGAAATTGTCGCCCTGGACGATCAGTACACGT  >  minE/1010758‑1010897

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: