Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,316,813 G→A A213T (GCA→ACA)  rcsB → DNA‑binding transcriptional activator RcsB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,316,8130GA88.0% 370.2 / 34.1 117A213T (GCA→ACA) rcsBDNA‑binding transcriptional activator RcsB
Reads supporting (aligned to +/- strand):  ref base G (5/9);  new base A (45/58);  total (50/67)
Fisher's exact test for biased strand distribution p-value = 7.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.84e-01

Only 100 of 121 total aligned reads displayed.

GAACGATATCGCCCTGCTGAATTATCTCTCTTCAGTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTT  >  NC_000913/2316764‑2316862
                                                 |                                                 
  aCGATATCGCCCTGCTGAATTATCTCTCTTCAGTGACCTTAAGTCCGACa                                                 <  1:4151890/50‑1 (MQ=255)
   cGATATCGCCCTGCTGAATTATCTCTCTTCAGTGACCTTAAGTCCGACAg                                                <  1:7542569/50‑1 (MQ=255)
   cGATATCGCCCTGCTGAATTATCTCTCTTCAGTGACCTTAAGTCCGACAg                                                <  1:3261603/50‑1 (MQ=255)
    gATATCGCCCTGCTGAATTATCTCTCTTCAGTGACCTTAAGTCCgg                                                   >  1:3872923/1‑46 (MQ=255)
     atatCGCCCTGCTGAATTATCTCTCTTCAGTGACCTTAAGTCCGGCAGa                                               >  1:6866716/1‑49 (MQ=255)
     atatCGCCCTGCTGAATTATCTCTCTTCAGTGACCTTAAGTCCGGCAGa                                               >  1:832651/1‑49 (MQ=255)
     atatCGCCCTGCTGAATTATCTCTCTTCAGTGACCTTAAGTCCGGCAGAt                                              <  1:1060423/50‑1 (MQ=255)
     atatCGCCCTGCTGAATTATCTCTCTTCAGTGACCTTAAGTCCGACAGa                                               >  1:2858521/1‑49 (MQ=255)
     atatCGCCCTGCTGAATTATCTCTCTTCAGTGACCTTAAGTCCGACAGAt                                              <  1:4210388/50‑1 (MQ=255)
      tatCGCCCTGCTGAATTATCTCTCTTCAGTGACCTTAAGTCCGGCAGATa                                             <  1:2209150/50‑1 (MQ=255)
      tatCGCCCTGCTGAATTATCTCTCTTCAGTGACCTTAAGTCCGACAGATa                                             <  1:798685/50‑1 (MQ=255)
        tCGCCCTGCTGAATTATCTCTCTTCAGTGACCTTAAGTCCGACAGATa                                             >  1:4711651/1‑48 (MQ=255)
          gCCCTGCTGAATTATCTCTCTTCAGTGACCTTAAGTCCGACAGATAAAg                                          <  1:1569839/49‑1 (MQ=255)
          gCCCTGCTGAATTATCTCTCTTCAGTGACCTTAAGTCCGACAGATAAAg                                          <  1:3646058/49‑1 (MQ=255)
          gCCCTGCTGAATTATCTCTCTTCAGTGACCTTAAGTCCGACAGATAAAGa                                         <  1:4955200/50‑1 (MQ=255)
          gCCCTGCTGAATTATCTCTCTTCAGTGACCTTAAGTCCGACAGATAAAGa                                         >  1:4467241/1‑50 (MQ=255)
          gCCCTGCTGAATTATCTCTCTTCAGTGACCTTAAGTCCGACAGATAAAGa                                         <  1:3147226/50‑1 (MQ=255)
          gCCCTGCTGAATTATCTCTCTTCAGTGACCTTAAGTCCGACAGATAAAGa                                         >  1:2952686/1‑50 (MQ=255)
           cccTGCTGAATTATCTCTCTTCAGTGACCTTAAGTCCGACAGATAAAGAc                                        >  1:53551/1‑50 (MQ=255)
            ccTGCTGAATTATCTCTCTTCAGTGACCTTAAGTCCGACAGATAAAGAc                                        >  1:4550154/1‑49 (MQ=255)
            ccTGCTGAATTATCTCTCTTCAGTGACCTTAAGTCCGACAGATAAAGAc                                        >  1:2416036/1‑49 (MQ=255)
            ccTGCTGAATTATCTCTCTTCAGTGACCTTAAGTCCGACAGATAAAGACt                                       <  1:13283/50‑1 (MQ=255)
            ccTGCTGAATTATCGCTCTTCAGTGACCTTAAGTCCGACAGATAAAGACt                                       <  1:3705475/50‑1 (MQ=255)
               gctgAATTATCTCTCTTCAGTGACCTTAAGTCCGACAGATAAAGACTAAt                                    >  1:7507394/1‑50 (MQ=255)
                ctgAATTATCTCTCTTCAGTGACCTTAAGTCCGGCAGATAAAGACTAAt                                    <  1:1413616/49‑1 (MQ=255)
                ctgAATTATCTCTCTTCAGTGACCTTAAGTCCGACAGATAAAGACTaa                                     <  1:7042407/48‑1 (MQ=255)
                ctgAATTATCTCTCTTCAGTGACCTTAAGTCCGACAGATAAAGACTaa                                     <  1:3996298/48‑1 (MQ=255)
                ctgAATTATCTCTCTTCAGTGACCTTAAGTCCGACAGATAAAGACTAAt                                    <  1:4296527/49‑1 (MQ=255)
                ctgAATTATCTCTCTTCAGTGACCTTAAGTCCGACAGATAAAGACTAAt                                    <  1:775191/49‑1 (MQ=255)
                ctgAATTATCTCTCTTCAGTGACCTTAAGTCCGACAGATAAAGACTAAt                                    >  1:5653699/1‑49 (MQ=255)
                ctgAATTATCTCTCTTCAGTGACCTTAAGTCCGACAGATAAAGACTAATc                                   <  1:1941971/50‑1 (MQ=255)
                ctgAATTATCTCTCTTCAGTGACCTTAAGTCCGACAGATAAAGACTAATc                                   <  1:7572630/50‑1 (MQ=255)
                  gAATTATCTCTCTTCAGTGACCTTAAGTCCGACAGATAAAGACTAATCAc                                 >  1:5486119/1‑50 (MQ=255)
                   aaTTATCTCTCTTCAGTGACCTTAAGTCCGGCAGATAAAGACTAATCa                                  <  1:5672655/48‑1 (MQ=255)
                   aaTTATCTCTCTTCAGTGACCTTAAGTCCGACAGATAAAGACTAATCAcc                                >  1:3212182/1‑50 (MQ=255)
                   aaTTATCTCTCTTCAGTGACCTTAAGTCCGACAGATAAAGACTAATCAc                                 <  1:5192779/49‑1 (MQ=255)
                   aaTTATCTCTCTTCAGTGACCTTAAGTCCGACAGATAAAGACTAATCAc                                 <  1:3025031/49‑1 (MQ=255)
                    aTTATCTCTCTTCAGTGACCTTAAGTCCGACAGATAAAGACTAATCACCt                               <  1:8495962/50‑1 (MQ=255)
                     ttATCTCTCTTCAGTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTg                              <  1:6739846/50‑1 (MQ=255)
                     ttATCTCTCTTCAGTGACCTTAAGTCCGACAGATAAAGACTAATCACCTg                              <  1:5491735/50‑1 (MQ=255)
                     ttATCTCTCTTCAGTGACCTTAAGTCCGACAGATAAAGACTAATCACCTg                              <  1:4871033/50‑1 (MQ=255)
                      tATCTCTCTTCAGTGACCTTAAGTCCGACAGATAAAGACTAATCACCTg                              >  1:1105211/1‑49 (MQ=255)
                       aTCTCTCTTCAGTGACCTTAAGTCCGACAGATAAAGACTAATCACCTg                              >  1:3434996/1‑48 (MQ=255)
                       aTCTCTCTTCAGTGACCTTAAGTCCGACAGATAAAGACTAATCACCTGt                             <  1:3348897/49‑1 (MQ=255)
                        tctctcTTCAGTGACCTTAAGTCCGACAGATAAAGACTAATCACCTGTa                            <  1:5329227/49‑1 (MQ=255)
                        tctctcTTCAGTGACCTTAAGTCCGACAGATAAAGACTAATCACCTGTa                            <  1:6567060/49‑1 (MQ=255)
                          tctctTCAGTGACCTTAAGTCCGACAGATAAAGACTAATCACCTGTAGGc                         >  1:1684353/1‑50 (MQ=255)
                          tctctTCAGTGACCTTAAGTCCGACAGATAAAGACTAATCACCTGTAGGc                         >  1:8320137/1‑50 (MQ=255)
                           ctctTCAGTGACCTTAAGTCCGACAGATAAAGACTAATCACCTGTAGGc                         >  1:455154/1‑49 (MQ=255)
                            tctTCAGTGACCTTAAGTCCTACAGATAAAGACTAATCACCTGTAg                           <  1:2097727/46‑1 (MQ=255)
                            tctTCAGTGACCTTAAGTCCGACAGATAAAGACTAATCACCTGTAgg                          <  1:6190568/47‑1 (MQ=255)
                            tctTCAGTGACCTTAAGTCCGACAGATAAAGACTAATCACCTGTAg                           <  1:890585/46‑1 (MQ=255)
                            tctTCAGTGACCTTAAGTCCGACAGATAAAGACTAATCACCTGTAg                           <  1:1071496/46‑1 (MQ=255)
                            tctTCAGTGACCTTAAGTCCGACAGATAAAGACTAATCACCTGTAGGc                         >  1:4320694/1‑48 (MQ=255)
                            tctTCAGTGACCTTAAGTCCGACAGATAAAGACTAATCACCTGTAGGc                         <  1:622577/48‑1 (MQ=255)
                            tctTCAGTGACCTTAAGTCCGACAGATAAAGACTAATCACCTGTAGGc                         >  1:5530264/1‑48 (MQ=255)
                            tctTCAGTGACCTTAAGTCCGACAGATAAAGACTAATCACCTGTAGGc                         <  1:3960226/48‑1 (MQ=255)
                            tctTCAGTGACCTTAAGACCGACAGATAAAGACTAATCACCTGTAGGc                         >  1:8254625/1‑48 (MQ=255)
                             ctTCAGTGACCTTAAGTCCGACAGATAAAGACTAATCACCTGTAGGc                         >  1:428801/1‑47 (MQ=255)
                             ctTCAGTGACCTTAAGTCCGACAGATAAAGACTAATCACCTGTAGGc                         >  1:5055571/1‑47 (MQ=255)
                              ttCAGTGACCTTAAGTCCGACAGATAAAGACTAATCACCTGTAGGcc                        <  1:7956035/47‑1 (MQ=255)
                                cAGTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGa                     <  1:7936691/48‑1 (MQ=255)
                                cAGTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGAt                    >  1:2731437/1‑49 (MQ=255)
                                cAGTGACCTTAAGTCCGACAGATAAAGACTAATCACCTGTAGGCCAg                      <  1:1359991/47‑1 (MQ=255)
                                cAGTGACCTTAAGTCCGACAGATAAAGACTAATCACCTGTAGGCCAGAt                    <  1:8532884/49‑1 (MQ=255)
                                cAGTGACCTTAAGTCCGACAGATAAAGACTAATCACCTGTAGGCCAGAt                    <  1:8338005/49‑1 (MQ=255)
                                cAGTGACCTTAAGTCCGACAGATAAAGACTAATCACCTGTAGGCCAGAt                    >  1:2066333/1‑49 (MQ=255)
                                cAGTGACCTTAAGTCCGACAGATAAAGACTAATCACCTGTAGGCCAGAt                    <  1:7003154/49‑1 (MQ=255)
                                 aGTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATaa                  <  1:5018105/50‑1 (MQ=255)
                                 aGTGACCTTAAGTCCGACAGATAAAGACTAATCACCTGTAGGCCAGAt                    >  1:7275893/1‑48 (MQ=255)
                                 aGTGACCTTAAGTCCGACAGATAAAGACTAATCACCTGTAGGCCAGATaa                  <  1:420611/50‑1 (MQ=255)
                                 aGTGACCTTAAGTCCGACAGATAAAGACTAATCACCTGTAGGCCAGATaa                  <  1:2370922/50‑1 (MQ=255)
                                 aGTGACCTTAAGTCCGACAGATAAAGACTAATCACCTGTAGGCCAGATaa                  <  1:6406345/50‑1 (MQ=255)
                                 aGTGACCTTAAGTCCGACAGATAAAGACTAATCACCTGTAGGCCAGATa                   <  1:8034194/49‑1 (MQ=255)
                                 aGTGACCTTAAGTCCGACAGATAAAGACTAATCACCTGTAGGCCAGATa                   <  1:3085728/49‑1 (MQ=255)
                                  gTGACCTTAAGTCCGACAGATAAAGACTAATCACCTGTAGGCCAGAt                    >  1:1172773/1‑47 (MQ=255)
                                   tGACCTTAAGTCCGACAGATAAAGACTAATCACCTGTAGGCCAGATAAGa                >  1:4785839/1‑50 (MQ=255)
                                   tGACCTTAAGTCCGACAGATAAAGACTAATCACCTGTAGGCCAGATAAGa                <  1:4012765/50‑1 (MQ=255)
                                   tGACCTTAAGTCCGACAGATAAAGACTAATCACCTGTAGGCCAGATAAGa                >  1:5400056/1‑50 (MQ=255)
                                    gACCTTAAGTCCGACAGATAAAGACTAATCACCTGTAGGCCAGATAAGAc               >  1:4120050/1‑50 (MQ=255)
                                     aCCTTAAGTCCGACAGATAAAGACTAATCACCTGTAGGCCAGATAAGAcg              >  1:5203466/1‑50 (MQ=255)
                                      ccTTAAGTCCGACAGATAAAGACTAATCACCTGTAGGCCAGATAAGAcgc             >  1:4900655/1‑50 (MQ=255)
                                      ccTTAAGTCCGACAGATAAAGACTAATCACCTGTAGGCCAGATAAGAcgc             <  1:2292847/50‑1 (MQ=255)
                                      ccTTAAGTCCGACAGATAAAGACTAATCACCTGTAGGCCAGATAAGAcgc             <  1:326946/50‑1 (MQ=255)
                                       cTTAAGTCCGACAGATAAAGACTAATCACCTGTAGGCCAGATAAGAcgcg            <  1:2469397/50‑1 (MQ=255)
                                         tAAGTCCGACAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGt           >  1:27103/1‑49 (MQ=255)
                                         tAAGTCCGACAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGt           >  1:3337008/1‑49 (MQ=255)
                                         tAAGTCCGACAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGt           >  1:5890511/1‑49 (MQ=255)
                                          aaGTCCGACAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGtt          <  1:79293/49‑1 (MQ=255)
                                          aaGTCCGACAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGtt          <  1:1218442/49‑1 (MQ=255)
                                          aaGTCCGACAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGt           >  1:1588313/1‑48 (MQ=255)
                                           aGTCCGACAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTa         <  1:3217636/49‑1 (MQ=255)
                                           aGTCCGACAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTa         <  1:2111802/49‑1 (MQ=255)
                                           aGTCCGACAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTa         <  1:3744137/49‑1 (MQ=255)
                                           aGTCCGACAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAg        >  1:5669247/1‑50 (MQ=255)
                                           aGTCCGACAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAg        >  1:8565516/1‑50 (MQ=255)
                                           aGTCCGACAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAg        <  1:2770105/50‑1 (MQ=255)
                                             tCCGACAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAgtg      <  1:3076667/50‑1 (MQ=255)
                                                ggCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCt   >  1:83826/1‑50 (MQ=255)
                                                 gCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCtt  <  1:586033/50‑1 (MQ=255)
                                                 |                                                 
GAACGATATCGCCCTGCTGAATTATCTCTCTTCAGTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTT  >  NC_000913/2316764‑2316862

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: