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%	105.1 / 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)

GCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATCACACGAGATTTGAGCGACTTCTGGAAAGGAG  >  NC_000913/164276‑164427                                                                             |                                                                            gcccgcGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGt                                                                         <  2:35644/81‑1 (MQ=255) gcccgcGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGt                                                                         >  1:35695/1‑81 (MQ=255)    cgcgCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCt                                                                      <  2:35645/81‑1 (MQ=255)      cgcTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTgg                                                                    <  2:35646/81‑1 (MQ=255)        cTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAg                                                                  <  1:35647/81‑1 (MQ=255)         tGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGt                                                                 >  2:35683/1‑81 (MQ=255)            cGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGc                                                              <  1:35648/81‑1 (MQ=255)             ggTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCt                                                             <  2:35650/81‑1 (MQ=255)             ggTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCt                                                             <  1:35649/81‑1 (MQ=255)                gagaATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCttt                                                           >  1:35692/1‑80 (MQ=255)                 agaATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCa                                                         >  2:35682/1‑81 (MQ=255)                        aaGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGcc                                                   >  2:35684/1‑80 (MQ=255)                            gaTGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAg                                              >  1:35697/1‑81 (MQ=255)                             aTGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAgg                                             >  1:35690/1‑81 (MQ=255)                             aTGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAgg                                             <  1:35651/81‑1 (MQ=255)                              tGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGc                                            <  2:35652/81‑1 (MQ=255)                                tttGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCa                                          <  2:35653/81‑1 (MQ=255)                                  tGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCAt                                         >  1:35698/1‑80 (MQ=255)                                   ggCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTa                                       >  2:35701/1‑81 (MQ=255)                                      gAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACaa                                    <  2:35654/81‑1 (MQ=255)                                       aGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACaaa                                   <  2:35655/81‑1 (MQ=255)                                         gCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAAt                                  >  2:35694/1‑80 (MQ=255)                                         gCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATc                                 <  2:35656/81‑1 (MQ=255)                                             ccaATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATCacac                             <  2:35657/81‑1 (MQ=255)                                                 tCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATCACACgaga                         <  1:35658/81‑1 (MQ=255)                                                           gCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATCACACGAGATTTGAGCGa                <  1:35659/80‑1 (MQ=255)                                                           gCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATCACACGAGATTTGAGCGAc               >  2:35702/1‑81 (MQ=255)                                                            cccAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATCACACGAGATTTGAGCGACt              >  2:35691/1‑81 (MQ=255)                                                                  ggCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATCACACGAGATTTGAGCGACTTCTgg         >  2:35706/1‑80 (MQ=255)                                                                        tGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATCACACGAGATTTGAGCGACTTCTGGAAAGGAg  >  1:35696/1‑81 (MQ=255)                                                                             |                                                                            GCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATCACACGAGATTTGAGCGACTTCTGGAAAGGAG  >  NC_000913/164276‑164427

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