Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_pae_tpiA	2,106,741	A→T	S162S (TCT→TCA)	wbbH ←	O‑antigen polymerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_pae_tpiA	2,106,741	0	A	T	100.0%	125.7 / NA	38	S162S (TCT→TCA)	wbbH	O‑antigen polymerase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (13/25);  total (13/25)

ACTCTATTAACACCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTAGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACTTAAC  >  NC_000913_3_pae_tpiA/2106595‑2106882                                                                                                                                                   |                                                                                                                                              aCTCTATTAACACCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGATTACCAATATCAAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTAATTAACTTTGTGTTAGTAAATTGTTTTGa                                                                                                                                             >  2:516060/1‑149 (MQ=255)     tATTAACACCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCaa                                                                                                                                         >  2:602375/1‑149 (MQ=255)     tATTAACACCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTAACCaa                                                                                                                                         >  2:454831/1‑149 (MQ=255)      ataaaCACCGACAATGAATGCATAAGAGATGATAACCTTAAAGACCATTTGCCTACCCGTATTCAGTATAATTGCAAAGAGGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCaaa                                                                                                                                        <  2:435433/146‑1 (MQ=255)          acacCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCaaataaa                                                                                                                                    >  1:622176/1‑149 (MQ=255)          acacCTACGATGAAGGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCaaataaa                                                                                                                                    <  2:219731/149‑1 (MQ=255)                          cATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGtt                                                                                                                    >  1:381718/1‑149 (MQ=255)                          cATAAGAGATGAAAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGtt                                                                                                                    <  1:621246/149‑1 (MQ=255)                           aTAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAg                                                                                                                               >  1:16669/1‑137 (MQ=255)                           aTAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAg                                                                                                                               <  2:16669/137‑1 (MQ=255)                           aTAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTa                                                                                                                   <  1:603511/149‑1 (MQ=255)                                 tatgatAACCATAAAGACAATTTGCTTACTGATATTCAGTATAAGTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGa                                                                                                             <  2:332756/148‑1 (MQ=255)                                  atgatAAGCATAAAGACAATTGTCGTACCAGTATGCAGTATAATTGAAAAGATTAATACAATAAAAACAAGTAAAGTCAATGTTTTACTTACCTTTGTATTAGTAAATTTGTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGAt                                                                                                            <  2:244004/149‑1 (MQ=255)                                   tgatAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGAtt                                                                                                           <  1:165315/149‑1 (MQ=255)                                     atAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTgg                                                                                                         <  1:331513/149‑1 (MQ=255)                                       aaCCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCt                                                                                                       <  2:10054/149‑1 (MQ=255)                                                   aaTTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAAtaaataaa                                                                                                                                <  1:400823/112‑1 (MQ=255)                                                   aaTTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAAtaaataaa                                                                                                                                >  2:400823/1‑112 (MQ=255)                                                                 tATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGAt                                                                             >  1:487047/1‑149 (MQ=255)                                                                     cAGTATAATTGCAAAGATGAATACAATAAAACCAAGTAAATTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTTATCTAATAAATAATGCAAAAGTAGTTAGAATGATTGTCTGCATGTATGGTGCGAAATTTCTTGAtgttt                                                                         >  2:301925/1‑147 (MQ=255)                                                                     cAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATgtgt                                                                         <  1:76387/149‑1 (MQ=255)                                                                                   aGATGAATACAATAAAAACAAATAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTTTCTTCAACATCAGCa                                                           >  2:492010/1‑149 (MQ=255)                                                                                       gAATACAATAAAAACAAGTAAAGTAAATTTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCAtctc                                                       <  2:495269/149‑1 (MQ=255)                                                                                          tACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTa                                                    <  1:107910/149‑1 (MQ=255)                                                                                           aCAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCCAATAAAGAAAGCAAAAGGAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTAt                                                   <  2:554426/149‑1 (MQ=255)                                                                                                aaaaaCAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGATGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAAt                                              >  1:176135/1‑149 (MQ=255)                                                                                                    aCAAGTAAAGTAAATGTTTTACTTACTTTTTGATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCat                                          <  2:103945/149‑1 (MQ=255)                                                                                                            agtaaaTGTTTTACTTACCTTTGGATTAGTAAATTTTTTTGACCAAATAAATCACGAAATAGTAGTTAGAATGATTTGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCtaa                                  <  1:311755/149‑1 (MQ=255)                                                                                                               aaaTGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCtaagta                               <  1:377190/149‑1 (MQ=255)                                                                                                                              cTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATATCTTATCAAATTCATCTAGCTAAGTAAGCTAGTCCCGAACt                >  1:221845/1‑149 (MQ=255)                                                                                                                                  gTATTAGTAAATTTTTTTGACCAAATAAACAAAGCAAAAGTAGTTAGATTGAGTGGCTGCATGTATGCTGAGAAATTTCTTGATGGGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTa            <  1:147800/149‑1 (MQ=255)                                                                                                                                   tATTAGTAAATTTTTTTGACCAAATAAATAAAGAACAAGTAGTTATAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAg           <  1:602375/149‑1 (MQ=255)                                                                                                                                    aTTAGTAAATTTTTTTGAACAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCTTATATCAAAGTAAGCTCGTCCCGAACTGGTaat          >  2:617340/1‑147 (MQ=255)                                                                                                                                      tAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTa        <  1:63391/149‑1 (MQ=255)                                                                                                                                      tAGTAAATTTGTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTa        <  2:223403/149‑1 (MQ=255)                                                                                                                                       aGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTAc       <  1:551807/149‑1 (MQ=255)                                                                                                                                            aTTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCCACCTAATCATCTCTTATCATATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACTTAAc  <  1:255321/149‑1 (MQ=255)                                                                                                                                            aTTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACTTAAc  <  1:454831/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                              ACTCTATTAACACCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTAGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACTTAAC  >  NC_000913_3_pae_tpiA/2106595‑2106882

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