Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_hsa_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_hsa_tpiA	2,106,741	0	A	T	100.0%	106.2 / NA	31	S162S (TCT→TCA)	wbbH	O‑antigen polymerase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (14/17);  total (14/17)

TCTATTAACACCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTAGACCA‑AATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACTTAACC  >  NC_000913_3_hsa_tpiA/2106597‑2106883                                                                                                                                                 |                                                                                                                                                tctATTAACACCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCa                                                                                                                                            >  2:224121/1‑150 (MQ=255)             tACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCA‑AATAAATAAAga                                                                                                                               >  1:136770/1‑149 (MQ=255)              aCGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTACATTTTTTTGACCA‑actaaatatatca                                                                                                                              >  2:72827/1‑145 (MQ=255)               cGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCA‑AATAAATAAAGCaa                                                                                                                             >  1:601464/1‑150 (MQ=255)                gATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCA‑AATAAATAAAGCaaa                                                                                                                            >  2:123737/1‑150 (MQ=255)                      tacATGAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTTTTCA‑ATTAATTAAAGCAAAtgtaat                                                                                                                      >  2:472823/3‑148 (MQ=255)                      tGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCA‑AATAAATAAAGCAAAagtagt                                                                                                                      <  1:454136/150‑1 (MQ=255)                        cATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCA‑AATAAATAAAGCAAAAGTAGTTa                                                                                                                    <  2:601464/150‑1 (MQ=255)                         aTAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCA‑AATAAATAAAGCAAAAGTAGTTAg                                                                                                                   >  1:205381/1‑150 (MQ=255)                            agagATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCA‑AATAAATAAAGCAAAAGTAGTTAGAAt                                                                                                                <  2:215129/150‑1 (MQ=255)                               gatgatAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCA‑AATAAATAAAGCAAAAGTAGTTAGAATGAt                                                                                                             >  2:318956/1‑150 (MQ=255)                                       ccATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCA‑AGTAAATAAAGTAAAAGTAGTTAGAATGATTGGCTGCa                                                                                                     >  2:339315/1‑150 (MQ=255)                                        cATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCA‑AATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCAt                                                                                                    >  2:698626/1‑150 (MQ=255)                                                 aaTTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCA‑AATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTga                                                                                           >  1:72587/1‑150 (MQ=255)                                                        ttACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAGGATGAATCTAG‑TATAGTGGTTCGACTAGTGGTCTGCATGTATGCTGAGAAAttt                                                                                    >  1:510499/1‑150 (MQ=255)                                                              gTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCA‑AATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGAt                                                                              <  1:281022/150‑1 (MQ=255)                                                                aTTCAGTATAATTGCAAAGATGAATAAAATAAAAAAAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCA‑AATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGAtgt                                                                            <  2:257667/150‑1 (MQ=255)                                                                              cAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCA‑AATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCa                                                               <  1:415719/149‑1 (MQ=255)                                                                              cAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCA‑AATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCa                                                               >  2:415719/1‑149 (MQ=255)                                                                              cAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCA‑AATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAg                                                              <  2:72587/150‑1 (MQ=255)                                                                                  aaTGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCA‑AATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCAtc                                                          <  1:353485/149‑1 (MQ=255)                                                                                        tACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCA‑AATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTAt                                                    <  2:136770/150‑1 (MQ=255)                                                                                           aaTAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCA‑AATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCaa                                                 >  1:219373/1‑150 (MQ=255)                                                                                            aTAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCA‑AATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCaaa                                                <  2:205381/150‑1 (MQ=255)                                                                                                   cAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCA‑AATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCatat                                         <  2:219373/150‑1 (MQ=255)                                                                                                     agtaaagtaaATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCA‑AATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCGTATCAAATTCATATAg                                       <  1:556626/150‑1 (MQ=255)                                                                                                           gtaaaTGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCA‑AATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCtaagt                                 <  1:319415/150‑1 (MQ=255)                                                                                                                 gTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCA‑AATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTa                           <  1:698626/150‑1 (MQ=255)                                                                                                                 gTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCA‑AATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTa                           <  1:363309/150‑1 (MQ=255)                                                                                                                           ccTTTGTATTAGTAAATTTTTTTGACCA‑AATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACACCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACt                 <  1:123737/150‑1 (MQ=255)                                                                                                                                          aTTTTTTTGACCA‑AATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACTTAAcc  <  2:398777/150‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                TCTATTAACACCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTAGACCA‑AATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACTTAACC  >  NC_000913_3_hsa_tpiA/2106597‑2106883

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