Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,043,741	A→G	V157A (GTT→GCT)	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_000913	1,043,741	0	A	G	100.0%	150.8 / NA	43	V157A (GTT→GCT)	etk	tyrosine‑protein kinase, role in O‑antigen capsule formation
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (23/20);  total (23/20)

GTACGCTGGCTCAGGACAAACTGTGTTCCTGGTTTGGCCTTAATGTCCGCGATAGTCAGCGCAACGCCATCTTTTTCCAGACGCTGGCCGACCATACCATTGACGGTGAACTCTTCACCTTCCAGTGTATAGTGGCCGTTTTCCCCAACCGTGAGTGTCAGTTGCTGATCCTGACCATTCAGTTGTGGAATATGCATCCAGCTGATCGCCAGCTCACCTGGTTTTTCTTTGGTTAATCTCGCCCAGCCGCGACCCACAATCGGAAAATACTTCTGCTCAACTATGTCGCGC  >  NC_000913/1043594‑1043884                                                                                                                                                    |                                                                                                                                                gTACGCTGGCTCAGGACAAACTGTGTTCCTGGTTTGGCCTTAATGTCCGCGATAGTCAGCGCAACGCCATCTTTTTCCAGACGCTGGCCGACCATACCATTGACGGTGAACTCTTCACCTTCCAGTGTATAGTGGCCGTTTTCCCCAGCCg                                                                                                                                              >  1:675164/1‑151 (MQ=255)   aCGCTGGCTCAGGACAAACTGTGTTCCTGGTTTGGCCTTAATGTCCGCGATAGTCAGCGCAACGCCATCTTTTTCCAGACGCTGGCCGACCATACCATTGACGGTGAACTCTTCACCTTCCAGTGTATAGTGGCCGTTTTCCCCAGCCGt                                                                                                                                             <  2:439068/150‑1 (MQ=255)          cTCAGGACAAACTGTGTTCCTGGTTTGGCCTTAATGTCCGCGATAGTCAGCGCAACGCCATCTTTTTCCAGACGCTGGCCGACCATACCATTGACGGTGAACTCTTCACCTTCCAGTGTATAGTGGCCGTTTTCCCCAGCCGTGAGTGTc                                                                                                                                      >  1:270155/1‑150 (MQ=255)            cAGGACAAACTGTGTTCCTGGTTTGGCCTTAATGTCCGCGATAGTCAGCGCAACGCCATCTTTTTCCAGACGCTGGCCGACCATACCATTGACGGTGAACTCTTCACCTTCCAGTGTATAGTGGCCGTTTTCCCCAGCCGTGAGTGTCAg                                                                                                                                    >  1:406706/1‑150 (MQ=255)            cAGGACAAACTGTGTTCCTGGTTTGGCCTTAATGTCCGCGATAGTCAGCGCAACGCCATCTTTTTCCAGACGCTGGCCGACCATACCATTGACGGTGAACTCTTCACCTTCCAGTGTATAGTGGCCGTTTTCCCCAGCCGTGAGTGTCAGt                                                                                                                                   >  2:276331/1‑151 (MQ=255)                  aaaCTGTGTTCCTGGTTTGGCCTTAATGTCCGCGATAGTCAGCGCAACGCCATCTTTTTCCAGACGCTGGCCGACCATACCATTGACGGTGAACTCTTCACCTTCCAGTGTATAGTGGCCGTTTTCCCCAGCCGTGAGTGTCAGTTGCTg                                                                                                                              >  2:197093/1‑150 (MQ=255)                  aaaCTGTGTTCCTGGTTTGGCCTTAATGTCCGCGATAGTCAGCGCAACGCCATCTTTTTCCAGACGCTGGCCGACCATACCATTGACGGTGAACTCTTCACCTTCCAGTGTATAGTGGCCGTTTTCCCCAGCCGTGAGTGTCAGTTGCTGa                                                                                                                             <  2:219135/151‑1 (MQ=255)                    aCTGTGTTCCTGGTTTGGCCTTAATGTCCGCGATAGTCAGCGCAACGCCATCTTTTTCCAGACGCTGGCCGACCATACCATTGACGGTGAACTCTTCACCTTCCAGTGTATAGTGGCCGTTTTCCCCAGCCGTGAGTGTCAGTTGCTGATc                                                                                                                           <  2:444280/151‑1 (MQ=255)                          ttCCTGGTTTGGCCTTAATGTCCGCGATAGTCAGCGCAACGCCATCTTTTTCCAGACGCTGGCCGACCATACCATTGACGGTGAACTCTTCACCTTCCAGTGTATAGTGGCCGTTTTCCCCAGCCGTGAGTGTCAGTTGCTGATCCTGAcc                                                                                                                     >  1:393882/1‑151 (MQ=255)                               ggTTTGGCCTTAATGTCCGCGATAGTCAGCGCAACGCCATCTTTTTCCAGACGCTGGCCGACCATACCATTGACGGTGAACTCTTCACCTTCCAGTGTATAGTGGCCGTTTTCCCCAGCCGTGAGTGTCAGTTGCTGATCCTGACCAtt                                                                                                                  <  1:641728/149‑1 (MQ=255)                                gTTTGGCCTTAATGTCCGCGATAGTCAGCGCAACGCCATCTTTTTCCAGACGCTGGCCGACCATACCATTGACGGTGAACTCTTCACCTTCCAGTGTATAGTGGCCGTTTTCCCCAGCCGTGAGTGTCAGTTGCTGATCCTGACCATTc                                                                                                                 >  2:611219/1‑149 (MQ=255)                                        ttAATGTCCGCGATAGTCAGCGCAACGCCATCTTTTTCCAGACGCTGGCCGACCATACCATTGACGGTGAACTCTTCACCTTCCAGTGTATAGTGGCCGTTTTCCCCAGCCGTGAGTGTCAGTTGCTGATCCTGACCATTCAGTTGTGGa                                                                                                        >  1:65880/1‑150 (MQ=255)                                                cgcgATAGTCAGCGCAACGCCATCTTTTTCCAGACGCTGGCCGACCATACCATTGACGGTGAACTCTTCACCTTCCAGTGTATAGTGGCCGTTTTCCCCAGCCGTGAGTGTCAGTTGCTGATCCTGACCATTCAGTTGTGGAATATGCATc                                                                                               <  2:376969/151‑1 (MQ=255)                                                      aGTCAGCGCAACGCCATCTTTTTCCAGACGCTGGCCGACCATACCATTGACGGTGAACTCTTCACCTTCCAGTGTATAGTGGCCGTTTTCCCCAGCCGTGAGTGTCAGTTGCTGATCCTGACCATTCAGTTGTGGAATATGCATCCAGCt                                                                                          <  1:551224/150‑1 (MQ=255)                                                               aaCGCCATCTTTTTCCAGACGCTGGCCGACCATACCATTGACGGTGAACTCTTCACCTTCCAGTGTATAGTGGCCGTTTTCCCCAGCCGTGAGTGTCAGTTGCTGATCCTGACCATTCAGTTGTGGAATATGCATCCAGCTGATCGCCAGc                                                                                >  1:36193/1‑151 (MQ=255)                                                               aaCGCCATCTTTTTCCAGACGCTGGCCGACCATACCATTGACGGTGAACTCTTCACCTTCCAGTGTATAGTGGCCGTTTTCCCCAGCCGTGAGTGTCAGTTGCTGATCCTGACCATTCAGTTGTGGAATATGCATCCAGCTGATCGCCAGc                                                                                >  2:145695/1‑151 (MQ=255)                                                                aCGCCATCTTTTTCCAGACGCTGGCCGACCATACCATTGACGGTGAACTCTTCACCTTCCAGTGTATAGTGGCCGTTTTCCCCAGCCGTGAGTGTCAGTTGCTGATCCTGACCATTCAGTTGTGGAATATGCATCCAGCTGATCGCCAGCt                                                                               <  2:8599/151‑1 (MQ=255)                                                                  gCCATCTTTTTCCAGACGCTGGCCGACCATACCATTGACGGTGAACTCTTCACCTTCCAGTGTATAGTGGCCGTTTTCCCCAGCCGTGAGTGTCAGTTGCTGATCCTGACCATTCAGTTGTGGAATATGCATCCAGCTGATCGCCAGCTc                                                                              >  2:350596/1‑150 (MQ=255)                                                                        tttttCCAGACGCTGGCCGACCATACCATTGACGGTGAACTCTTCACCTTCCAGTGTATAGTGGCCGTTTTCCCCAGCCGTGAGTGTCAGTTGCTGATCCTGACCATTCAGTTGTGGAATATGCATCCAGCTGATCGCCAGCTCACCTGGt                                                                       >  2:495849/1‑151 (MQ=255)                                                                        tttttCCAGACGCTGGCCGACCATACCATTGACGGTGAACTCTTCACCTTCCAGTGTATAGTGGCCGTTTTCCCCAGCCGTGAGTGTCAGTTGCTGATCCTGACCATTCAGTTGTGGAATATGCATCCAGCTGATCGCCAGCTCACCTGGt                                                                       >  1:113559/1‑151 (MQ=255)                                                                         ttttCCAGACGCTGGCCGACCATACCATTGACGGTGAACTCTTCACCTTCCAGTGTATAGTGGCCGTTTTCCCCAGCCGTGAGTGTCAGTTGCTGATCCTGACCATTCAGTTGTGGAATATGCATCCAGCTGATCGCCAGCTCACCt                                                                          <  2:634766/147‑1 (MQ=255)                                                                          tttCCAGACGCTGGCCGACCATACCATTGACGGTGAACTCTTCACCTTCCAGTGTATAGTGGCCGTTTTCCCCAGCCGTGAGTGTCAGTTGCTGATCCTGACCATTCAGTTGTGGAATATGCATCCAGCTGATCGCCAGCTCACCTGGtt                                                                      <  2:259085/150‑1 (MQ=255)                                                                           ttCCAGACGCTGGCCGACCATACCATTGACGGTGAACTCTTCACCTTCCAGTGTATAGTGGCCGTTTTCCCCAGCCGTGAGTGTCAGTTGCTGATCCTGACCATTCAGTTGTGGAATATGCATCCAGCTGATCGCCAGCTCACCTGGttt                                                                     <  2:358860/150‑1 (MQ=255)                                                                            tCCAGACGCTGGCCGACCATACCATTGACGGTGAACTCTTCACCTTCCAGTGTATAGTGGCTGTTTTCCCCAGCCGTGAGTGTCAGTTGCTGATCCTGACCATTCAGTTGTGGAATATGCATCCAGCTGATCGCCAGCTCACCTGGtttt                                                                    <  1:303462/150‑1 (MQ=255)                                                                            tCCAGACGCTGGCCGACCATACCATTGACGGTGAACTCTTCACCTTCCAGTGTATAGTGGCCGTTTTCCCCAGCCGGGAGTGTCAGTTGCTGATCCTGACCATTCAGTTGTGGAATATGCATCCAGCTGATCGCCAGCTCACCTGGttttt                                                                   >  1:488952/1‑151 (MQ=255)                                                                             ccAGACGCTGGCCGACCATACCATTGACGGTGAACTCTTCACCTTCCAGTGTATAGTGGCCGTTTTCCCCAGCCGTGAGTGTCAGTTGCTGATCCTGACCATTCAGTTGTGGAATATGCATCCAGCTGATCGCCAGCTCACCTGGttt                                                                     <  1:32075/148‑1 (MQ=255)                                                                                 aCGCTGGCCGACCATACCATTGACGGTGAACTCTTCACCTTCCAGTGTATAGTGGCCGTTTTCCCCAGCCGTGAGTGTCAGTTGCTGATCCTGACCATTCAGTTGTGGAATATGCATCCAGCTGATCGCCAGCTCACCTGGTTTTTCTTTg                                                              <  2:554800/151‑1 (MQ=255)                                                                                       gCCGACCATACCATTGACGGTGAACTCTTCACCTTCCAGTGTATAGTGGCCGTTTTCCCCAGCCGTGAGTGTCAGTTGCTGATCCTGACCATTCAGTTGTGGAATATGCATCCAGCTGATCGCCAGCTCACCTGGTTTTTCTTTGGTTAAt                                                        <  1:35291/151‑1 (MQ=255)                                                                                         cGACCATACCATTGACGGTGAACTCTTCACCTTCCAGTGTATAGTGGCCGTTTTCCCCAGCCGTGAGTGTCAGTTGCTGATCCTGACCATTCAGTTGTGGAATATGCATCCAGCTGATCGCCAGCTCACCTGGTTTTTCTTTGGTTAAtct                                                      <  2:642238/151‑1 (MQ=255)                                                                                          gACCATACCATTGACGGTGAACTCTTCACCTTCCAGTGTATAGTGGCCGTTTTCCCCAGCCGTGAGTGTCAGTTGCTGATCCTGACCATTCAGTTGTGGAATATGCATCCAGCTGATCGCCAGCTCACCTGGTTTTTCTTTGGTTAAtctc                                                     <  1:61089/151‑1 (MQ=255)                                                                                                   atTGACGGTGAACTCTTCACCTTCCAGTGTATAGTGGCCGTTTTCCCCAGCCGTGAGTGTCAGTTGCTGATCCTGACCATTCAGTTGTGGAATATGCATCCAGCTGATCGCCAGCTCACCTGGTTTTTCTTTGGTTAATCTCGCCCAGCcg                                            >  2:186890/1‑151 (MQ=255)                                                                                                   atTGACGGTGAACTCTTCACCTTCCAGTGTATAGTGGCCGTTTTCCCCAGCCGTGAGTGTCAGTTGCTGATCCTGACCATTCAGTTGTGGAATATGCATCCAGCTGATCGCCAGCTCACCTGGTTTTTCTTTGGTTAATCTCGCCCAGCcg                                            <  1:656780/151‑1 (MQ=255)                                                                                                    ttGACGGTGAACTCTTCACCTTCCAGTGTATAGTGGCCGTTTTCCCCAGCCGTGAGTGTCAGTTGCTGATCCTGACCATTCAGTTGTGGAATATGCATCCAGCTGATCGCCAGCTCACCTGGTTTTTCTTTGGTTAATCTCGCCCAGCcgc                                           >  2:212432/1‑151 (MQ=255)                                                                                                     tGACGGTGAACTCTTCACCTTCCAGTGTATAGTGGCCGTTTTCCCCAGCCGTGAGTGTCAGTTGCTGATCCTGACCATTCAGTTGTGGAATATGCATCCAGCTGATCGCCAGCTCACCTGGTTTTTCTTTGGTTAATCTCGCCCAGCcgcg                                          >  1:278315/1‑151 (MQ=255)                                                                                                      gACGGTGAACTCTTCACCTTCCAGTGTATAGTGGCCGTTTTCCCCAGCCGTGAGTGTCAGTTGCTGATCCTGACCATTCAGTTGTGGAATATGCATCCAGCTGATCGCCAGCTCACCTGGTTTTTCTTTGGTTAATCTCGCCCAGCCGCGa                                         <  2:396610/151‑1 (MQ=255)                                                                                                       aCGGTGAACTCTTCACCTTCCAGTGTATAGTGGCCGTTTTCCCCAGCCGTGAGTGTCAGTTGCTGATCCTGACCATTCAGTTGTGGAATATGCATCCAGCTGATCGCCAGCTCACCTGGTTTTTCTTTGGTTAATCTCGCCCAGCCGCGAc                                        >  1:308376/1‑151 (MQ=255)                                                                                                                    cACCTTCCAGTGTATAGTGGCCGTTTTCCCCAGCCGTGAGTGTCAGTTGCTGATCCTGACCATTCAGTTGTGGAATATGCATCCAGCTGATCGCCAGCTCACCTGGTTTTTCTTTGGTTAATCTCGCCCAGCCGCGACCCACAATCGGaa                            >  2:15246/1‑150 (MQ=255)                                                                                                                            aGTGTATAGTGGCCGTTTTCCCCAGCCGTGAGTGTCAGTTGCTGATCCTGACCATTCAGTTGTGGAATATGCATCCAGCTGATCGCCAGCTCACCTGGTTTTTCTTTGGTTAATCTCGCCCAGCCGCGACCCACAATCGGAAAATACTTCt                   >  1:110643/1‑151 (MQ=255)                                                                                                                                    gTGGCCGTTTTCCCCAGCCGTGAGTGTCAGTTGCTGATCCTGACCATTCAGTTGTGGAATATGCATCCAGCTGATCGCCAGCTCACCTGGTTTTTCTTTGGTTAATCTCGCCCAGCCGCGACCCACAATCGGAAAATACTTCTGCTCAACt           <  1:279307/151‑1 (MQ=255)                                                                                                                                          gTTTTCCCCAGCCGTGAGTGTCAGTTGCTGATCCTGACCATTCAGTTGTGGAATATGCATCCAGCTGATCGCCAGCTCACCTGGTTTTTCTTTGGTTAATCTCGCCCAGCCGCGACCCACAATCGGAAAATACTTCTGCTCAACTATGTc      >  1:198116/1‑150 (MQ=255)                                                                                                                                              tCCCCAGCCGTGAGTGTCAGTTGCTGATCCTGACCATTCAGTTGTGGAATATGCATCCAGCTGATCGCCAGCTCACCTGGTTTTTCTTTGGTTAATCTCg                                                    <  2:422738/100‑1 (MQ=255)                                                                                                                                              tCCCCAGCCGTGAGTGTCAGTTGCTGATCCTGACCATTCAGTTGTGGAATATGCATCCAGCTGATCGCCAGCTCACCTGGTTTTTCTTTGGTTAATCTCg                                                    >  1:422725/1‑100 (MQ=255)                                                                                                                                               ccccAGCCGTGAGTGTCAGTTGCTGATCCTGACCATTCAGTTGTGGAATATGCATCCAGCTGATCGCCAGCTCACCTGGTTTTTCTTTGGTTAATCTCGCCCAGCCGCGACCCACAATCGGAAAATACTTCTGCTCAACTATGTCgcgc  >  2:358375/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                GTACGCTGGCTCAGGACAAACTGTGTTCCTGGTTTGGCCTTAATGTCCGCGATAGTCAGCGCAACGCCATCTTTTTCCAGACGCTGGCCGACCATACCATTGACGGTGAACTCTTCACCTTCCAGTGTATAGTGGCCGTTTTCCCCAACCGTGAGTGTCAGTTGCTGATCCTGACCATTCAGTTGTGGAATATGCATCCAGCTGATCGCCAGCTCACCTGGTTTTTCTTTGGTTAATCTCGCCCAGCCGCGACCCACAATCGGAAAATACTTCTGCTCAACTATGTCGCGC  >  NC_000913/1043594‑1043884

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