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%	148.5 / NA	41	L750L (CTG→TTG)	hrpB	putative ATP‑dependent RNA helicase HrpB
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (16/25);  total (16/25)

CCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATCACACGAGATTTGAGCGACTTCTGGAAAGGAGC  >  NC_000913/164278‑164428                                                                           |                                                                             ccgcGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGc                                                                        <  1:37864/81‑1 (MQ=255)    cgcTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTg                                                                      >  1:37914/1‑80 (MQ=255)     gcTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTg                                                                      >  2:37924/1‑79 (MQ=255)      cTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAg                                                                   <  1:37865/81‑1 (MQ=255)         gCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGtt                                                                 <  2:37866/80‑1 (MQ=255)          cGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGc                                                               <  1:37867/81‑1 (MQ=255)           ggTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCt                                                              >  2:37930/1‑81 (MQ=255)           ggTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCt                                                              >  2:37918/1‑81 (MQ=255)           ggTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCt                                                              >  2:37916/1‑81 (MQ=255)            gTGAGAATGCAAGAGATGTTTGGCGAGGCAACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCtt                                                             <  1:37868/81‑1 (MQ=255)               agaATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTc                                                           <  1:37869/80‑1 (MQ=255)               agaATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTc                                                           <  1:37870/80‑1 (MQ=255)                  aTGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCt                                                       >  2:37933/1‑81 (MQ=255)                   tGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTg                                                      <  2:37871/81‑1 (MQ=255)                       agagATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCa                                                  <  1:37872/81‑1 (MQ=255)                          gaTGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAg                                               >  1:37925/1‑81 (MQ=255)                              tttGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGcc                                            <  2:37873/80‑1 (MQ=255)                              tttGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCa                                           <  2:37874/81‑1 (MQ=255)                                 ggCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCAtt                                         <  1:37875/80‑1 (MQ=255)                                   cGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACa                                      <  1:37876/81‑1 (MQ=255)                                    gAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACa                                      >  2:37929/1‑80 (MQ=255)                                     aGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACaaa                                    >  2:37923/1‑81 (MQ=255)                                      ggCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAAt                                   >  1:37927/1‑81 (MQ=255)                                      ggCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAAt                                   >  2:37920/1‑81 (MQ=255)                                         caccaATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATcac                                <  2:37877/81‑1 (MQ=255)                                            caATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATCacac                              >  1:453563/1‑80 (MQ=255)                                               tCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATCACACgaga                          <  1:37878/81‑1 (MQ=255)                                                  gacgaTCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATCACACGAGAttt                       >  1:37941/1‑81 (MQ=255)                                                  gacgaTCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATCACACGAGAttt                       <  2:37879/81‑1 (MQ=255)                                                   acgaTCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATCACACGAGATTTg                      <  1:37880/81‑1 (MQ=255)                                                    cgaTCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTACTTTCACCTGCCCAAAGGCCATTACAAATCACACGAGATTTGa                     <  2:37881/81‑1 (MQ=255)                                                     gaTCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATCACACGAGATTTGAg                    >  1:37935/1‑81 (MQ=255)                                                      aTCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCTAAGGCCATTACAAATCACACGAGATTTGAGc                   >  2:37926/1‑81 (MQ=255)                                                      aTCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATCACACGAGATTTGAGc                   >  2:37928/1‑81 (MQ=255)                                                         gCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATCACACGAGATTTGAGCGAc                <  2:37882/81‑1 (MQ=255)                                                            cAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATCACACGAGATTTGAGCGACTTc             <  2:37883/81‑1 (MQ=255)                                                                  cgcgTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATCACACGAGATTTGAGCGACTTCTGGaaa       <  1:37884/81‑1 (MQ=255)                                                                   gcgTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATCACACGAGATTTGAGCGACTTCTGGAAAg      <  1:37885/81‑1 (MQ=255)                                                                     gTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATCACACGAGATTTGAGCGACTTCTGGAAAGGa    <  2:37887/81‑1 (MQ=255)                                                                     gTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATCACACGAGATTTGAGCGACTTCTGGAAAGGa    <  1:37886/81‑1 (MQ=255)                                                                        ccGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATCACACGAGATTTGAGCGACTTCTGGAAAGGAGc  <  1:37888/80‑1 (MQ=255)                                                                           |                                                                             CCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATCACACGAGATTTGAGCGACTTCTGGAAAGGAGC  >  NC_000913/164278‑164428

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