Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,320,580 G→A 6.6% N52N (AAC→AAT rcsC ← hybrid sensory kinase in two‑component regulatory system with RcsB and YojN

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,320,5800GA 6.6% 276.4 / 6.1 91N52N (AAC→AATrcsChybrid sensory kinase in two‑component regulatory system with RcsB and YojN
Reads supporting (aligned to +/- strand):  ref base G (39/46);  new base A (4/2);  total (43/48)
Fisher's exact test for biased strand distribution p-value = 4.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

CGATTACCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTAAGTACATTAAGCGCATCGACGCCATCATTCG  >  W3110S.gb/2320510‑2320650
                                                                      |                                                                      
cGATTACCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATg                                                                        >  1:3006894/1‑71 (MQ=255)
cGATTACCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATg                                                                        >  1:3002256/1‑71 (MQ=255)
cGATTACCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATg                                                                        >  1:300787/1‑71 (MQ=255)
cGATTACCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATAATTGGCATg                                                                        >  1:3613723/1‑71 (MQ=255)
  aTTACCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGtt                                                                      <  1:2914393/71‑1 (MQ=255)
  aTTACCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGtt                                                                      <  1:563401/71‑1 (MQ=255)
  aTTACCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGtt                                                                      <  1:2433917/71‑1 (MQ=255)
      ccGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACg                                                                  <  1:3257805/71‑1 (MQ=255)
         gCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGtt                                                                      <  1:1181371/64‑1 (MQ=255)
         gCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTg                                                                     <  1:2219334/65‑1 (MQ=255)
          cAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGtt                                                                      <  1:1914328/63‑1 (MQ=255)
           aaCGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCt                                                             <  1:3207344/71‑1 (MQ=255)
            aCGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTa                                                            <  1:2047696/71‑1 (MQ=255)
            aCGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTa                                                            <  1:1712593/71‑1 (MQ=255)
            aCGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTa                                                            <  1:1402641/71‑1 (MQ=255)
              gTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAg                                                          >  1:3238409/1‑71 (MQ=255)
              gTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAg                                                          >  1:3598976/1‑71 (MQ=255)
               tCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGt                                                                 >  1:3011836/1‑63 (MQ=255)
               tCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGt                                                                 >  1:824352/1‑63 (MQ=255)
                 aGTCCCAACTGACGAATGCGTTGCGTCAAGCTGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCAc                                                       >  1:2078163/1‑71 (MQ=255)
                 aGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCAc                                                       >  1:2576933/1‑71 (MQ=255)
                 aGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCAc                                                       <  1:3806067/71‑1 (MQ=255)
                 aGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCAc                                                       <  1:3588095/71‑1 (MQ=255)
                   tCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCAc                                                       <  1:837939/69‑1 (MQ=255)
                   tCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGa                                                     <  1:600676/71‑1 (MQ=255)
                    cccAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGa                                                     >  1:3520598/1‑70 (MQ=255)
                    cccAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGa                                                     >  1:3131045/1‑70 (MQ=255)
                    cccAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGa                                                     <  1:3636444/70‑1 (MQ=255)
                     ccAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGata                                                   >  1:1942142/1‑71 (MQ=255)
                     caaaCTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGata                                                   >  1:2813521/3‑71 (MQ=255)
                     caaaCTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGata                                                   >  1:2268616/3‑71 (MQ=255)
                       aaCTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGatat                                                  >  1:646001/1‑70 (MQ=255)
                       aaCTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGatat                                                  >  1:1248832/1‑70 (MQ=255)
                       aaCTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGatat                                                  >  1:1871863/1‑70 (MQ=255)
                       aaCTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGatat                                                  >  1:2896424/1‑70 (MQ=255)
                       aaCTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATc                                                 <  1:1770606/71‑1 (MQ=255)
                       aaCTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATATTGACGTCTCTAAGCACGatat                                                  >  1:2248187/1‑70 (MQ=255)
                       aaCTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATATTGACGTCGCTAAGCACGatat                                                  >  1:3578924/1‑70 (MQ=255)
                       aaCTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATATTGACGTCGCTAAGCACGatat                                                  >  1:2632954/1‑70 (MQ=255)
                       aaCTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATATTGACGTCGCTAAGCACGatat                                                  >  1:257109/1‑70 (MQ=255)
                        aCTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGtt                                                                      <  1:3689811/49‑1 (MQ=255)
                         cTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCAc                                                       <  1:2637454/63‑1 (MQ=255)
                         cTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCAc                                                       <  1:1918874/63‑1 (MQ=255)
                         cTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCAc                                                       <  1:3729537/63‑1 (MQ=255)
                         cTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCaa                                               <  1:3075429/71‑1 (MQ=255)
                         cTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCaa                                               <  1:340821/71‑1 (MQ=255)
                         cTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCaa                                               >  1:2249295/1‑71 (MQ=255)
                         cTGACGAATGCGTTGAGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCaa                                               <  1:3833985/71‑1 (MQ=255)
                          tGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGa                                                                    >  1:239055/1‑49 (MQ=255)
                          tGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGa                                                                    >  1:3285497/1‑49 (MQ=255)
                           gACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGa                                                                    >  1:2816034/1‑48 (MQ=255)
                           gACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGACATCAata                                             >  1:819005/1‑71 (MQ=255)
                            aCGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAatat                                            >  1:1119624/1‑71 (MQ=255)
                               aaTGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGAt                                         <  1:64850/71‑1 (MQ=255)
                                aTGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCaa                                               >  1:2472436/1‑64 (MQ=255)
                                aTGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCaa                                               >  1:23575/1‑64 (MQ=255)
                                 tgcgttgcgtCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTc                                       >  1:1574412/1‑71 (MQ=255)
                                  gcgttgcgtCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCt                                      >  1:1733879/1‑71 (MQ=255)
                                   cgttgcgtCAAGCGGTAACCATCCAGATTTGGCATGTTGACGTCGCTaa                                                           >  1:2668457/1‑49 (MQ=255)
                                    gttgcgtCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTg                                    <  1:1174355/71‑1 (MQ=255)
                                     ttgcgtCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCAc                                                       <  1:17585/51‑1 (MQ=255)
                                       gcgtCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCt                                  <  1:240866/70‑1 (MQ=255)
                                       gcgtCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCt                                  <  1:353545/70‑1 (MQ=255)
                                        cgtCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCt                                  <  1:3282913/69‑1 (MQ=255)
                                         gtCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTAAg                               <  1:3729771/71‑1 (MQ=255)
                                           cAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTAAGt                              >  1:1993918/1‑70 (MQ=255)
                                            aaGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTAAGTa                             >  1:1023359/1‑70 (MQ=255)
                                             aGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTAAGTACa                           <  1:3042343/71‑1 (MQ=255)
                                               cGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTAAGTACAtt                         <  1:1219896/71‑1 (MQ=255)
                                               cGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTAAGTACAtt                         <  1:324870/71‑1 (MQ=255)
                                               cGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTAAGTACAt                          >  1:1997716/1‑70 (MQ=255)
                                               cGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTAAGTACAt                          >  1:3764664/1‑70 (MQ=255)
                                               cGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTAAGTACAt                          >  1:3644035/1‑70 (MQ=255)
                                               cGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTAAGTACAt                          >  1:959345/1‑70 (MQ=255)
                                                ggTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTAAGTACATTa                        <  1:1381637/71‑1 (MQ=255)
                                                 gTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGGTATCAATATGATTCTTGCTAAGTACATTag                       <  1:321060/71‑2 (MQ=255)
                                                 gTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTAAGTACATTaa                       <  1:2450783/71‑1 (MQ=255)
                                                     ccctccatATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTAAGTACATTAAgcgc                   <  1:3710976/68‑1 (MQ=255)
                                                     ccatccatATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTAAGTACATTAAgcgc                   <  1:731629/71‑1 (MQ=255)
                                                     ccatccatATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTAAGTACATTAAgcgc                   <  1:1840552/71‑1 (MQ=255)
                                                       atccatATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTAAGTACATTAAGCGCAt                 >  1:3024498/1‑71 (MQ=255)
                                                          catATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTAAGTACATTAAGCGCa                  >  1:459530/1‑67 (MQ=255)
                                                          catATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTAAGTACATTAAGCGCATCGa              >  1:1918223/1‑71 (MQ=255)
                                                           atatTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTAAGTACATTAAgcgc                   <  1:119622/65‑1 (MQ=255)
                                                              tttGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCtt                                     >  1:2455717/1‑44 (MQ=255)
                                                              tttGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTAAGTACATTAAGCGCATCGACGc           >  1:995447/1‑70 (MQ=255)
                                                              tttGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTAAGTACATTAAGCGCATCGACGc           <  1:1805602/70‑1 (MQ=255)
                                                               ttGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTAAGTACATTAAGCGCATCGACGc           <  1:3386222/69‑1 (MQ=255)
                                                                tGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTAAGTACATTAAGCGCATCGACGCcat        >  1:2353824/1‑71 (MQ=255)
                                                                 ggCATGTTGACGTCGCTAAGCACGACATCAATATGATTCTTGCTAAGTACATTAAGCGCATCGACGCcat        <  1:1182182/70‑1 (MQ=255)
                                                                   cATATTGACGTCGCTAAGCACGATATCAATATGATTCTTGTTAAGTACATTAAGCGCATCGACGCcatcat     <  1:3286474/71‑1 (MQ=255)
                                                                   cATATTGACGTCGCTAAGCACGATATCAATATGATTCTTGTTAAGTACATTAAGCGCATCGACGCcatcat     <  1:2598423/71‑1 (MQ=255)
                                                                     tGTTGACTTCGCTAAGCACGATATCAATATGATTCTTGCTAAGTACATTAAGCGCATCGACGCCATCATTc   >  1:615171/1‑71 (MQ=255)
                                                                     tGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTAAGTACATTAAGCGCATCGACGCCATCATTc   <  1:2639844/71‑1 (MQ=255)
                                                                     tGTTGACGTCGCTAGCACGATATCAATATGATTCTTGCt                                  <  1:1928456/39‑1 (MQ=37)
                                                                      gTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTAAGTAc                            >  1:3754670/1‑45 (MQ=255)
                                                                      gTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTAAGTACATTAAGCGCATCGACGCCATCATTcg  >  1:1212432/1‑71 (MQ=255)
                                                                      |                                                                      
CGATTACCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTAAGTACATTAAGCGCATCGACGCCATCATTCG  >  W3110S.gb/2320510‑2320650

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: