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

CCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATCACACGAGATTTGAGCGACTTCTGGAAAGGAGCG  >  NC_000913/164274‑164429                                                                               |                                                                              ccgccgggGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGGCGATCGCCCAGGGGCGCGTGCCGTTg                                                                             <  1:39171/81‑1 (MQ=255)  cgcccgcGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTgg                                                                            >  2:39224/1‑81 (MQ=255)  cgcccgcGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTgg                                                                            >  1:39208/1‑81 (MQ=255)   gcccgcGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGt                                                                           <  2:39172/81‑1 (MQ=255)   gcccgcGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGt                                                                           >  1:39215/1‑81 (MQ=255)      cgcgCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCt                                                                        <  1:39173/81‑1 (MQ=255)       gcgcTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTg                                                                       >  2:39219/1‑81 (MQ=255)         gcTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTg                                                                       >  2:39210/1‑79 (MQ=255)          cTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAg                                                                    >  1:39212/1‑81 (MQ=255)               ggTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCt                                                               >  1:39214/1‑81 (MQ=255)                gTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCtt                                                              <  1:39174/81‑1 (MQ=255)                  gagaATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCttt                                                             >  1:39231/1‑80 (MQ=255)                  gagaATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCttt                                                             >  1:39228/1‑80 (MQ=255)                   agaATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCa                                                           >  1:39220/1‑81 (MQ=255)                      aTGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCt                                                        <  1:39175/81‑1 (MQ=255)                          aaGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGccc                                                    <  1:39176/81‑1 (MQ=255)                          aaGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGcc                                                     >  2:39233/1‑80 (MQ=255)                             agaTGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCaaa                                                 <  2:39177/81‑1 (MQ=255)                                tGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGc                                              >  1:39225/1‑81 (MQ=255)                                 gTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGc                                              >  2:39207/1‑80 (MQ=255)                                      gCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTAc                                        <  2:39178/81‑1 (MQ=255)                                       cGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACa                                       <  2:39179/81‑1 (MQ=255)                                        gAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACa                                       >  1:39218/1‑80 (MQ=255)                                             caccaATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATca                                  >  1:39216/1‑80 (MQ=255)                                              accaATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATcaca                                >  2:39222/1‑81 (MQ=255)                                               ccaATCCGACGATCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATCacac                               >  1:39221/1‑81 (MQ=255)                                                     cgacgaTCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATCACACGAGAtt                         >  1:39247/1‑81 (MQ=255)                                                      gacgaTCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATCACACGAGAttt                        >  2:39253/1‑81 (MQ=255)                                                      gacgaTCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATCACACGAGAttt                        >  1:39245/1‑81 (MQ=255)                                                      gacgaTCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATCACACGAGAttt                        >  1:39237/1‑81 (MQ=255)                                                         gaTCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATCACACGAGATTTGAg                     >  2:39243/1‑81 (MQ=255)                                                         gaTCGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATCACACGAGATTTGAg                     <  1:39180/81‑1 (MQ=255)                                                            cGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATCACACGAGATTTGAGCGa                  >  2:39255/1‑81 (MQ=255)                                                            cGCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATCACACGAGATTTGAGCGa                  >  2:39256/1‑81 (MQ=255)                                                             gCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATCACACGAGATTTGAGCGAc                 >  2:39227/1‑81 (MQ=255)                                                             gCCCAGGGGCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATCACACGAGATTTGAGCGAc                 >  1:39230/1‑81 (MQ=255)                                                                  ggggCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATCACACGAGATTTGAGCGACTTCTg            <  2:39181/81‑1 (MQ=255)                                                                   gggCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATCACACGAGATTTGAGCGACTTCTg            >  2:39238/1‑80 (MQ=255)                                                                    ggCGCGTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATCACACGAGATTTGAGCGACTTCTGGa          <  2:39182/81‑1 (MQ=255)                                                                       gcgTGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATCACACGAGATTTGAGCGACTTCTGGAAAg       <  2:39183/81‑1 (MQ=255)                                                                          tGCCGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATCACACGAGATTTGAGCGACTTCTGGAAAGGAg    >  2:39232/1‑81 (MQ=255)                                                                            ccGTTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATCACACGAGATTTGAGCGACTTCTGGAAAGGAGCg  >  1:39248/1‑81 (MQ=255)                                                                               |                                                                              CCGCCCGCGCTGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATCACACGAGATTTGAGCGACTTCTGGAAAGGAGCG  >  NC_000913/164274‑164429

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: ‑