Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,042,246 A→G S655S (AGT→AGC)  etk ← tyrosine‑protein kinase, role in O‑antigen capsule formation

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,042,2460AG100.0% 177.7 / NA 47S655S (AGT→AGC) etktyrosine‑protein kinase, role in O‑antigen capsule formation
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (26/21);  total (26/21)

GATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACGCCCCACGACCGCGGCATCACTCACCGCCAGCATCGGCGGCGTATCGACAATCACCAGATCGTAATGGTCGTTCGCCCATTCCAGTAATTGACGCATCCGATCGCGCATCAGCAGTTCAGACGGGTTAGGTGGCACCTGACCGCGAGTAAT  >  NC_000913/1042114‑1042377
                                                                                                                                    |                                                                                                                                   
gATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACGCCCCACGACCGCGGCATCGCTCAcc                                                                                                                               >  1:128209/1‑139 (MQ=255)
   aGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACGCCCCACGACCGCGGCATCGCTCACCGcc                                                                                                                            >  1:726128/1‑139 (MQ=255)
   aGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACGCCCCACGACCGCGGCATCGCTCACCGcc                                                                                                                            <  1:969094/139‑1 (MQ=255)
     cgcCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACGCCCCACGACCGCGGCATCGCTCACCGCCAg                                                                                                                          >  1:496634/1‑139 (MQ=255)
     cgcCTTTAATATTGACGCCTGCCTGTTCCAGACGCTACATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACGCCCCACGACCGCGGCATCGCTCACCGCCAg                                                                                                                          >  2:25271/1‑139 (MQ=255)
               atTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACGCCCCACGACCGCGGCATCGCTCACCGCCAGCATCggcggc                                                                                                                >  1:432551/1‑139 (MQ=255)
                  gACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACGCCCCACGACCGCGGCATCGCTCACCGCCAGCATCGGCGGCGTa                                                                                                             >  2:563174/1‑139 (MQ=255)
                  gACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACGCCCCACGACCGCGGCATCGCTCACCGCCAGCATCGGCGGCGTa                                                                                                             >  2:1064723/1‑139 (MQ=255)
                          cctgTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACGCCCCACGACCGCGGCATCGCTCACCGCCAGCATCGGCGGCGTATCGACAAt                                                                                                     >  2:306611/1‑139 (MQ=255)
                          cctgTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACGCCCCACGACCGCGGCATCGCTCACCGCCAGCATCGGCGGCGTATCGACAAt                                                                                                     <  1:306611/139‑1 (MQ=255)
                          cctgTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTACCAACAGAACGCCCCACGACCGCGGCATCGCTCACCGCCAGCATCGGCGGCGTATCGACAAt                                                                                                     >  2:208801/1‑139 (MQ=255)
                            tgTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACGCCCCACGACCGCGGCATCGCTCACCGCCAGCATCGGCGGCGTATCGACAATCa                                                                                                   <  2:726128/139‑1 (MQ=255)
                             gTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACGCCCCACGACCGCGGCATCGCTCACCGCCAGCATCGGCGGCGTATCGACAATCAc                                                                                                  <  1:1110220/139‑1 (MQ=255)
                              ttCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACGCCCCACGACCGCGGCATCGCTCACCGCCAGCATCGGCGGCGTATCGACAATCAcc                                                                                                 <  2:1091255/139‑1 (MQ=255)
                                        tGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACGCCCCACGACCGCGGCATCGCTCACCGCCAGCATCGGCGGCGTACCGACAATCACCAGATCGTAAt                                                                                       >  2:1076130/1‑139 (MQ=255)
                                         gCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACGCCCCACGACCGCGGCATCGCTCACCGCCAGCATCGGCGGCGTATCGACAATCACCAGATCGTAATg                                                                                      >  1:564411/1‑139 (MQ=255)
                                                   aCTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACGCCCCACGACCGCGGCATCGCTCACCGCCAGCATCGGCGGCGTATCGACAATCACCAGATCGTAATGGTCGTTCGcc                                                                            <  1:1064723/139‑1 (MQ=255)
                                                                 cGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACGCCCCACGACCGCGGCATCGCTCACCGCCAGCATCGGCGGCGTATCGACAATCACCAGATCGTAATGGTCGTTCGCCCATTCCAGTAATTg                                                              <  2:855358/139‑1 (MQ=255)
                                                                 cGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACGCCCCACGACCGCGGCATCGCTCACCGCCAGCATCGGCGGCGTATCGACAATCACCAGATCGTAATGGTCGTTCGCCCATTCCAGTAATTg                                                              >  1:519405/1‑139 (MQ=255)
                                                                      ttCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACGCCCCACGACCGCGGCATCGCTCACCGCCAGCATCGGCGGCGTATCGACAATCACCAGATCGTAATGGTCGTTCGCCCATTCCAGTAATTGACGCa                                                         <  2:898856/139‑1 (MQ=255)
                                                                      ttCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACGCCCCACGACCGCGACATCGCTCACCGCCAGCATCGGCGGCGTATCGACAATCACCAGATCGTAATGGTCGTTCGcc                                                                            <  2:541965/120‑1 (MQ=255)
                                                                      ttCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACGCCCCACGACCGCGACATCGCTCACCGCCAGCATCGGCGGCGTATCGACAATCACCAGATCGTAATGGTCGTTCGcc                                                                            >  1:541965/1‑120 (MQ=255)
                                                                        cAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACGCCCCACGACCGCGGCATCGCTCACCGCCAGCATCGGCGGCGTATCGACAATCACCAGATCGTAATGGTCGTTCGCCCATTCCAGTAATTGACGCATc                                                       >  2:970313/1‑139 (MQ=255)
                                                                                aaCGCGCAACCAGCAGGCTGGTGCCAACAGAACGCCCCACGACCGCGGCATCGCTCACCGCCAGCATCGGCGGCGTATCGACAATCACCAGATCGTAATGGTCGTTCGCCCATTCCAGTAATTGACGCATCCGATcgcg                                               >  2:488364/1‑139 (MQ=255)
                                                                                   gcgcAACCAGCAGGCTGGTGCCAACAGAACGCCCCACGACCGCGGCATCGCTCACCGCCAGCATCGGCGGCGTATCGACAATCACCAGATCGTAATGGTCGTTCGCCCATTCCAGTAATTGACGCGTCCGATCGCGCAt                                            <  2:432551/139‑1 (MQ=255)
                                                                                   gcgcAACCAGCAGGCTGGTGCCAACAGAACGCCCCACGACCGCGGCATCGCTCACCGCCAGCATCGGCGGCGTATCGACAATCACCAGATCGTAATGGTCGTTCGCCCATTCCAGTAATTGACGCATCCGATCGCGCAt                                            <  2:519405/139‑1 (MQ=255)
                                                                                                   ggTGCCAACAGAACGCCCCACGACCGCGGCATCGCTCACCGCCAGCATCGGCGGCGTATCGACAATCACCAGATCGTAATGGTCGTTCGCCCATTCCAGTAATTGACGCATCCGATCGCGCATCAGCAGTTCAGACggg                            <  1:970313/139‑1 (MQ=255)
                                                                                                          aCAGAACGCCCCACGACCGCGGCATCGCTCACCGCCAGCATCGGCGGCGTATCGACAATCACCAGATCGTAATGGTCGTTCGCCCATTCCAGTAATTGACGCATCCGATCGCGCATCAGCAGTTCAGACGGGTTAGGTg                     >  1:856150/1‑139 (MQ=255)
                                                                                                              aaCGCCCCACGACCGCGGCATCGCTCACCGCCAGCATCGGCGGCGTATCGACAATCACCAGATCGTAATGGTCGTTCGCCCATTCCAGTAATTGACGCATCCGATCGCGCATCAGCAGTTCAGACGGGTTAGGTGGCAc                 <  1:583321/139‑1 (MQ=255)
                                                                                                              aaCGCCCCACGACCGCGGCATCGCTCACCGCCAGCATCGGCGGCGTATCGACAATCACCAGATCGTAATGGTCGTTCGCCCATTCCAGTAATTGACGCATCCGATCGCGCATCAGCAGTTCAGACGGGTTAGGTGGCAc                 <  1:563174/139‑1 (MQ=255)
                                                                                                              aaCGCCCCACGACCGCGGCATCGCTCACCGCCAGCATCGGCGGCGTATCGACAATCACCAGATCGTAATGGTCGTTCGCCCATTCCAGTAATTGACGCATCCGATCGCGCATCAGCAGTTCAGACGGGTTAGGTGGCAc                 >  1:802704/1‑139 (MQ=255)
                                                                                                                   cccACGACCGCGGCATCGCTCACCGCCAGCATCGGCGGCGTATCGACAATCACCAGATCGTAATGGTCGTTCGCCCATTCCAGTAATTGACGCATCCGATCGCGCATCAGCAGTTCAGACGGGTTAGGTGGCACCTGAc            <  1:25271/139‑1 (MQ=255)
                                                                                                                    ccACGACCGCGGCATCGCTCACCGCCAGCATCGGCGGCGTATCGACAATCACCAGATCGTAATGGTCGTTCGCCCATTCCAGTAATTGACGCATCCGATCGCGCATCAGCAGTTCAGACGGGTTAGGTGGCACCTGAcc           <  2:564411/139‑1 (MQ=255)
                                                                                                                        gACCGCGGCATCGCTCACCGCCAGCATCGGCGGCGTATCGACAATCACCa                                                                                                >  2:758279/1‑50 (MQ=255)
                                                                                                                        gACCGCGGCATCGCTCACCGCCAGCATCGGCGGCGTATCGACAATCACCa                                                                                                <  1:758279/50‑1 (MQ=255)
                                                                                                                             cgGCATCGCTCACCGCCAGCATCGGCGGCGTATCGACAATCa                                                                                                   <  1:370705/42‑1 (MQ=255)
                                                                                                                             cgGCATCGCTCACCGCCAGCATCGGCGGCGTATCGACAATCa                                                                                                   >  2:370705/1‑42 (MQ=255)
                                                                                                                             cgGCATCGCTCACCGCCAGCATCGGCGGCGTATCGACAATCACCAGATCGTAATGGTCGTTCGcc                                                                            >  2:187537/1‑65 (MQ=255)
                                                                                                                             cgGCATCGCTCACCGCCAGCATCGGCGGCGTATCGACAATCACCAGATCGTAATGGTCGTTCGcc                                                                            <  1:187537/65‑1 (MQ=255)
                                                                                                                             cgGCATCGCTCACCGCCAGCATCGGCGGCGTATCGACAATCACCAGATCGTAATGGTCGTTCGCCCATTCCAGTAATTGACGCATCCGATCGCGCATCAGCAGTTCAGa                                >  2:1047429/1‑109 (MQ=255)
                                                                                                                             cgGCATCGCTCACCGCCAGCATCGGCGGCGTATCGACAATCACCAGATCGTAATGGTCGTTCGCCCATTCCAGTAATTGACGCATCCGATCGCGCATCAGCAGTTCAGa                                <  1:1047429/109‑1 (MQ=255)
                                                                                                                             cgGCATCGCTCACCGCCAGCATCGGCGGCGTATCGACAATCACCAGATCGTAATGGTCGTTCGCCCATTCCAGTAATTGACGCATCCGATCGCGCATCAGCAGTTCAGACGGGTTAGGTGGCACCTGTCCGCGAGTAAt  >  1:727253/1‑139 (MQ=255)
                                                                                                                             cgGCATCGCTCACCGCCAGCATCGGCGGCGTATCGACAATCACCAGATCGTAATGGTCGTTCGCCCATTCCAGTAATTGACGCATCCGATCGCGCATCAGCAGTTCAGACGGGTTAGGTGGCACCTGACCGCGAGTAAt  >  1:415749/1‑139 (MQ=255)
                                                                                                                             cgGCATCGCTCACCGCCAGCATCGGCGGCGTATCGACAATCACCAGATCGTAATGGTCGTTCGCCCATTCCAGTAATTGACGCATCCGATCGCGCATCAGCAGTTCAGACGGGTTAGGTGGCACCTGACCGCGAGTAAt  >  2:1079883/1‑139 (MQ=255)
                                                                                                                             cgGCATCGCTCACCGCCAGCATCGGCGGCGTATCGACAATCACCAGATCGTAATGGTCGTTCGCCCATTCCAGTAATTGACGCATCCGATCGCGCATCAGCAGTTCAGACGGGTTAGGTGGCACCTGACCGCGAGTAAt  >  1:917841/1‑139 (MQ=255)
                                                                                                                              ggCATCGCTCACCGCCAGCATCGGCGGCGTATCGa                                                                                                         <  1:228892/35‑1 (MQ=255)
                                                                                                                              ggCATCGCTCACCGCCAGCATCGGCGGCGTATCGa                                                                                                         >  2:228892/1‑35 (MQ=255)
                                                                                                                                    |                                                                                                                                   
GATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACGCCCCACGACCGCGGCATCACTCACCGCCAGCATCGGCGGCGTATCGACAATCACCAGATCGTAATGGTCGTTCGCCCATTCCAGTAATTGACGCATCCGATCGCGCATCAGCAGTTCAGACGGGTTAGGTGGCACCTGACCGCGAGTAAT  >  NC_000913/1042114‑1042377

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