Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	164,352	C→T	L750L (CTG→TTG)	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,352	0	C	T	100.0%	114.3 / NA	30	L750L (CTG→TTG)	hrpB	putative ATP‑dependent RNA helicase HrpB
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (14/16);  total (14/16)

GCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATCACACGAGATTTGAGCGACTTCTGGAAAGGAGC  >  NC_000913/164276‑164428                                                                             |                                                                             gcccgcGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGt                                                                          >  1:32780/1‑81 (MQ=255)   ccgcGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGc                                                                        >  1:32796/1‑81 (MQ=255)   ccgcGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGc                                                                        >  1:32791/1‑81 (MQ=255)    cgcgCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCt                                                                       >  1:32827/1‑81 (MQ=255)          ggCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGtt                                                                 >  2:32792/1‑81 (MQ=255)           gCGGTGAGAATGCGAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGtt                                                                 <  1:32753/80‑1 (MQ=255)               tGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCttt                                                            <  1:32754/81‑1 (MQ=255)                   aaTGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCAcc                                                        >  2:32789/1‑81 (MQ=255)                       cAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGcc                                                    <  2:491564/81‑1 (MQ=255)                       cAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGcc                                                    <  1:32755/81‑1 (MQ=255)                        aaGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGccc                                                   <  2:32756/81‑1 (MQ=255)                            gaTGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAg                                               <  2:32757/81‑1 (MQ=255)                                 ttGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCAt                                          >  1:32798/1‑81 (MQ=255)                                    gCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTAc                                       <  1:32758/81‑1 (MQ=255)                                     cGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTAc                                       >  1:32802/1‑80 (MQ=255)                                     cGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACa                                      <  2:32759/81‑1 (MQ=255)                                           caccaATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATcac                                <  1:32761/81‑1 (MQ=255)                                                  ccGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATCACACGAGAt                         >  1:32797/1‑81 (MQ=255)                                                    gacgaTCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATCACACGAGAttt                       >  2:32803/1‑81 (MQ=255)                                                     acgaTCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATCACACGAGATTTg                      <  1:32762/81‑1 (MQ=255)                                                      cgaTCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATCACACGAGATTTGa                     <  2:32763/81‑1 (MQ=255)                                                           gCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATCACACGAGATTTGAGCGAc                <  1:32764/81‑1 (MQ=255)                                                            cccAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATCACACGAGATTTGAGCGACt               >  2:32812/1‑81 (MQ=255)                                                              cAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATCACACGAGATTTGAGCGACTTc             <  1:32765/81‑1 (MQ=255)                                                              cAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATCACACGAGATTTGAGCGACTTc             >  2:32816/1‑81 (MQ=255)                                                               aGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATCACACGAGATTTGAGCGACTTCt            <  1:32766/81‑1 (MQ=255)                                                                  ggCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATCACACGAGATTTGAGCGACTTCTgg          >  1:32814/1‑80 (MQ=255)                                                                  ggCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATCACACGAGATTTGAGCGACTTCTGGa         <  1:32767/81‑1 (MQ=255)                                                                       gTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATCACACGAGATTTGAGCGACTTCTGGAAAg      >  2:32795/1‑79 (MQ=255)                                                                          ccGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATCACACGAGATTTGAGCGACTTCTGGAAAGGAGc  <  1:32768/80‑1 (MQ=255)                                                                             |                                                                             GCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATCACACGAGATTTGAGCGACTTCTGGAAAGGAGC  >  NC_000913/164276‑164428

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑