Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	164,162	A→G	A686A (GCA→GCG)	hrpB →	RNA‑dependent NTPase HrpB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	164,162	0	A	G	100.0%	97.4 / NA	29	A686A (GCA→GCG)	hrpB	RNA‑dependent NTPase HrpB
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (12/17);  total (12/17)

AGTTGCATCAGGCGATGCTTAATGGCATCCGTGATAAAGGTTTAAGCGTGCTCAACTGGACGGCGGAAGCGGAACAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATCGCCATTCGTTATCATGAAGATAACCCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCA  >  NC_000913/163924‑164387                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                  aGTTGCATCAGGCGATGCTTAATGGCATCCGTGATAAAGGTTTAAGCGTGCTCAACTGGACGGCGGAAGCGGAACAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTAGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCGCTACGCGGATTa                                                                                                                                                                                                                       <  2:255008/251‑1 (MQ=255)      cATCAGGCGATGCTTAATGGCATCCGTGATAAAGGTTTAAGCGTGCTCAACTGGACGGCGGAAGCGGAACAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCGCTACGCGGATTACTTg                                                                                                                                                                                                                   >  1:2669/1‑250 (MQ=255)      cATCAGGCGATGCTTAATGGCATCCGTGATAAAGGTTTAAGCGTGCTCAACTGGACGGCGGAAGCGGAACAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAACGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCGCTACGCGGATTACTTGa                                                                                                                                                                                                                  <  2:253940/251‑1 (MQ=255)                 gCTTAATGGCGTCCGTGATAAAGGTTTAAGCGTGCTCAACTGGACGGCGGAAGCGGAACAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCGCTACGCGGATTACTTGATTGGGGAATgc                                                                                                                                                                                                       <  1:102226/251‑1 (MQ=255)                             ccGTGATAAAGGTTTAAGCGTGCTCAACTGGACGGCGGAAGCGGCACCGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCGCTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTgg                                                                                                                                                                                           >  1:210389/1‑251 (MQ=255)                                     aaGGTTTAAGCGTGCTCAACTGGACGGCGGAAGCGGAACAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCGCTACGCGGATTACTGGATTGGGGAATGCAGCAACGTCTGGATAGTGaa                                                                                                                                                                                   <  1:174145/251‑1 (MQ=255)                                            aaGCGTGCTCAACTGGACGGCGGAAGCGGAACAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACGTATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCGCTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCt                                                                                                                                                                             <  1:26204/250‑1 (MQ=255)                                                  gCTCAACTGGACGGCGGAAGCGGAACAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCGCTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCAt                                                                                                                                                                       <  2:225219/250‑1 (MQ=255)                                                       aCTGGACGGCGGAAGCGGAACAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCGCTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCATTACACt                                                                                                                                                                 <  2:321109/251‑1 (MQ=255)                                                                 ggaagcggaaCAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCGCTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCATTACACTGTGCCGACgg                                                                                                                                                       >  1:339424/1‑251 (MQ=255)                                                                                    ttGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCGCTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATCGCCATTCg                                                                                                                                    >  1:349444/1‑251 (MQ=255)                                                                                    ttGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCGCTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATCGCCATTCg                                                                                                                                    >  1:349417/1‑251 (MQ=255)                                                                                     tGCGTATGTTATACGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGAGGTTTATTATGAAAGTTTAGTTGCGGCGCTGGAAACGGGGAGGCTGACACATATGACCGGCGTACATTCAATACGCGGCCTGAAATCACTAGACATTTATCAGTCGCTACGCGGATTACTTGATTTGGGAATGCAGCAACGTCTGGCTAGTGAATTGCCTGCGCATTACCCTGTTGCGACGGTAAGACGGATTGCCATTc                                                                                                                                     <  2:220837/249‑1 (MQ=255)                                                                                      gCGTTTGTTATGCGCCGAAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCGCTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAAGTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATCGCCATTCGtt                                                                                                                                  <  2:43343/251‑1 (MQ=255)                                                                                        gTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCGCTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATCGCCATTCGTTAt                                                                                                                                >  1:224857/1‑251 (MQ=255)                                                                                          ttGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCGCTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATCGCCATTCGTTATc                                                                                                                               <  1:156104/250‑1 (MQ=255)                                                                                                                  ccGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCGCTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCATTACACTGTGACGACGGGAAGCCGGATCGCCATTCGTTATCATGAAGATAACCCGCCCGCGCTGGc                                                                                                      <  1:219995/251‑1 (MQ=255)                                                                                                                                                       tATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCGCTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATCGCCATTCGTTATCATGAAGATAACCCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCaa                                                                  <  2:89730/250‑1 (MQ=255)                                                                                                                                                        aTTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCGCTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATCGCCATTCGTTATCATGAAGATAACCCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATc                                                                <  1:180852/251‑1 (MQ=255)                                                                                                                                                              aGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCGCTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATCGCCATTCGTTATCATGAAGATAGCCCGCCCGCGCt                                                                                                         <  1:279109/204‑1 (MQ=255)                                                                                                                                                              aGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCGCTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATCGCCATTCGTTATCATGAAGATAACCCGCCCGCGCt                                                                                                         >  2:279109/1‑204 (MQ=255)                                                                                                                                                                       aaCGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCGCTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATCGCCATTCGTTATCATGAAGATAACCCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAgg                                                 <  1:358946/251‑1 (MQ=255)                                                                                                                                                                             gctgctGCACCATATGACTGGCGTACATTCAATACGCGGCCTGAAATCACTCGACATTTATCAGGCGCTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATCGCCATTCGTTATCATGAAGATAACCCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGGTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGcgcg                                           >  2:209817/1‑251 (MQ=255)                                                                                                                                                                                    ccACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCGCTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATCGCCATTCGTTATCATGAAGATAACCCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCt                                    >  2:221084/1‑251 (MQ=255)                                                                                                                                                                                             aCTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCGCTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATCGCCATTCGTTATCATGAAGATAACCCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTgg                            >  1:255811/1‑250 (MQ=255)                                                                                                                                                                                                 gCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCGCTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATCGCCATTCGTTATCATGAAGATAACCCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTg                       <  1:370030/251‑1 (MQ=255)                                                                                                                                                                                                               cgcgGCCTGAAATCACTCGACATTTATCAGGCGCTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATCGCCATTCGTTATCATGAAGATAACCCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTGCTTTCACCTGccc          >  2:143404/1‑250 (MQ=255)                                                                                                                                                                                                                     tggAAATAACTCGACATTTATCAGGCGCTACGCGGATTCCTTGATTGGGGAATGCAGCAACTTCTGGATAGTGAATTGCCTGAGCATTACACTGTGCCGACGGGAAGCCGGATCGCCATTCGTTATCATGATGATAACCCGCCCGCGCTGGCGGTGAGAAGGCAATAGATGGTTGGCGAGGCCACCAATCCGACGATCGTCCCGGGGCGCGTGGCGCTGGTGCTGGGGGTGGTTTCACCTGCCCAAATGcc   >  2:233652/3‑251 (MQ=255)                                                                                                                                                                                                                       gAAATCACTCGACATTTATCAGGCGCTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATCGCCATTCGTTATCATGAAGATAACCCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCa  <  2:363857/250‑1 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                  AGTTGCATCAGGCGATGCTTAATGGCATCCGTGATAAAGGTTTAAGCGTGCTCAACTGGACGGCGGAAGCGGAACAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATCGCCATTCGTTATCATGAAGATAACCCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCA  >  NC_000913/163924‑164387

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