Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	164,090	A→G	E662E (GAA→GAG)	hrpB →	putative ATP‑dependent RNA helicase HrpB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	164,090	0	A	G	100.0%	125.9 / NA	36	E662E (GAA→GAG)	hrpB	putative ATP‑dependent RNA helicase HrpB
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (15/21);  total (15/21)

AGGTTTAAGCGTGCTCAACTGGACGGCGGAAGCGGAACAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCATTACA  >  NC_000913/163961‑164226                                                                                                                                  |                                                                                                                                         aGGTTTAAGCGTGCTCAACTGGACGGCGGAAGCGGAACAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAGACGTGgctgc                                                                                                                                <  1:327336/140‑1 (MQ=255)   gTTTAAGCGTGCTCAACTGGACGGCGGAAGCGGAACAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAGACGTGGctgctg                                                                                                                              >  1:240836/1‑140 (MQ=255)   gTTTAAGCGTGCTCAACTGGACGGCGGAAGCGGAACAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAGACGTGGctgctg                                                                                                                              >  1:88034/1‑140 (MQ=255)    tttAAGCGTGCTCAACTGGACGGCGGAAGCGGAACAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAGACGTGGCtgctgc                                                                                                                             <  2:572341/140‑1 (MQ=255)    tttAAGCGTGCTCAACTGGACGGCGGAAGCGGAACAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAGACGTGGCtgctgc                                                                                                                             <  1:217470/140‑1 (MQ=255)      tAAGCGTGCTCAACTGGACGGCGGAATCGGAACAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAGACGTGGCTGCTGCca                                                                                                                           >  1:412628/1‑140 (MQ=255)             gCTCAACTGGACGGCGGAAGCGGAACAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAGACGTGGCTGCTGCCACATATGa                                                                                                                    <  1:511206/140‑1 (MQ=255)                      gACGGCGGAAGCGGAACAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAGACGTGGCTGCTGCCACATATGACTGGCGTAc                                                                                                           >  2:910204/1‑140 (MQ=255)                      gACGGCGGAAGCGGAACAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAGACGTGGCTGCTGCCACATATGACTGGCGTAc                                                                                                           >  2:459663/1‑140 (MQ=255)                                gcggaaCAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAGACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTAcg                                                                                                 <  2:88034/140‑1 (MQ=255)                                    aaCAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAGACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGc                                                                                             >  1:412822/1‑140 (MQ=255)                                     aCAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAGACGTGGCTGCTGCCACATATGACTGGCGTACAtt                                                                                                        <  2:657676/128‑1 (MQ=255)                                     aCAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAGACGTGGCTGCTGCCACATATGACTGGCGTACAtt                                                                                                        >  1:657676/1‑128 (MQ=255)                                     aCAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAGACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGcc                                                                                            >  2:59274/1‑140 (MQ=255)                                         cTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAGACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGaa                                                                                        >  1:449932/1‑140 (MQ=255)                                               ttGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGGCTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAGACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACt                                                                                  <  2:240836/140‑1 (MQ=255)                                                  cGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAGACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGa                                                                               >  2:595301/1‑140 (MQ=255)                                                      tGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAGACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACAtt                                                                           >  2:233938/1‑140 (MQ=255)                                                             cgccgcAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAGACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAgg                                                                    <  2:871086/140‑1 (MQ=255)                                                             cgccgcAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAGACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAgg                                                                    <  2:248119/140‑1 (MQ=255)                                                                   aaaGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAGACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTAc                                                              <  2:412628/140‑1 (MQ=255)                                                                    aaGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAGACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTAcg                                                             <  1:118625/140‑1 (MQ=255)                                                                    aaGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAGACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTAcg                                                             <  1:326187/140‑1 (MQ=255)                                                                             ccGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAGACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACt                                                    <  2:733082/140‑1 (MQ=255)                                                                               ggAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAGACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTg                                                  >  1:787633/1‑140 (MQ=255)                                                                                   tatGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAGACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTg                                              <  2:412822/140‑1 (MQ=255)                                                                                     tGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAGACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTggg                                            <  2:766167/140‑1 (MQ=255)                                                                                           gCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAGACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTggg                                            >  2:180382/1‑134 (MQ=255)                                                                                           gCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAGACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTggg                                            <  1:180382/134‑1 (MQ=255)                                                                                           gCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAGACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATgc                                      >  1:427681/1‑140 (MQ=255)                                                                                               gCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAGACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATgcagca                                  >  1:452785/1‑140 (MQ=255)                                                                                                   ttGATGATGAAAGTTTATTGGCAGCGCTGGAGACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGt                              <  1:910204/140‑1 (MQ=255)                                                                                                   ttGATGATGAAAGTTTATTGGCAGCGCTGGAGACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGt                              <  1:459663/140‑1 (MQ=255)                                                                                                               gTTTATTGGCAGCGCTGGAGACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGaa                  <  2:427681/140‑1 (MQ=255)                                                                                                                        cAGCGCTGGAGACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTgcg         <  1:591204/140‑1 (MQ=255)                                                                                                                        cAGCGCTGGAGACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTgcg         <  1:179307/140‑1 (MQ=255)                                                                                                                               ggAGACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCATTaca  <  2:452785/140‑1 (MQ=255)                                                                                                                                  |                                                                                                                                         AGGTTTAAGCGTGCTCAACTGGACGGCGGAAGCGGAACAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCATTACA  >  NC_000913/163961‑164226

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