Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	162,249	G→A	G49R (GGG→AGG)	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	162,249	0	G	A	100.0%	143.0 / NA	40	G49R (GGG→AGG)	hrpB	putative ATP‑dependent RNA helicase HrpB
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (19/21);  total (19/21)

TATTATTAAGTGCGCCGACCGGGGCCGGGAAATCAACCTGGCTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACGGGAAAATTATCCTGCTGGAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAAC  >  NC_000913/162172‑162323                                                                              |                                                                           tattatTAAGTGCGCCGACCGGGGCCGGGAAATCAACCTGGCTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACagg                                                                          >  1:35167/1‑80 (MQ=255)     attaAGTGCGCCGACCGGGGCCGGGAAATCAACCTGGCTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACAGGaaaa                                                                      <  2:35110/80‑1 (MQ=255)             cgcCGACCGGGGCCGGGAAATCAACCTGGCTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACAGGAAAATTATCctgc                                                             <  1:35111/81‑1 (MQ=255)                cGACCGGGGCCGGGAAATCAACCTGGCTGCCACTGCAACTGCTGGCGCATCCCGGCATTAACAGGAAAATTATCCTGCTgg                                                          <  2:35112/81‑1 (MQ=255)                 gACCGGGGCCGGGAAATCAACCTGGCTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACAGGAAAATTATCctgctg                                                           >  2:35164/1‑79 (MQ=255)                 gACCGGGGCCGGGAAATCAACCTGGCTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACAGGAAAATTATCctgctg                                                           >  1:35169/1‑79 (MQ=255)                   ccGGGGCCGGGAAATCAACCTGGCTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACAGGAAAATTATCCTGCTGGAGc                                                       <  2:35113/81‑1 (MQ=255)                    cGGGGCCGGGAAATCAACCTGGCTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACAGGAAAATTATCCTGCTGGAGcc                                                      <  1:35115/81‑1 (MQ=255)                    cGGGGCCGGGAAATCAACCTGGCTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACAGGAAAATTATCCTGCTGGAGcc                                                      <  2:35114/81‑1 (MQ=255)                     ggggCCGGGAAATCAACCTGGCCGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACAGGAAAATTATCCTGCTGGAGCcg                                                     >  1:35170/1‑81 (MQ=255)                      gggCCGGGAAATCAACCTGGCTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACAGGAAAATTATCCTGCTGGAGCcgc                                                    >  1:35173/1‑81 (MQ=255)                      gggCCGGGAAATCAACCTGGCTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACAGGAAAATTATCCTGCTGGAGCcgc                                                    >  1:35189/1‑81 (MQ=255)                      gggCCGGGAAATCAACCTGGCTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACAGGAAAATTATCCTGCTGGAGCcgc                                                    <  1:35116/81‑1 (MQ=255)                       ggCCGGGAAATCAACCTGGCTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACAGGAAAATTATCCTGCTGGAGCcgcg                                                   >  2:35157/1‑81 (MQ=255)                          cGGGAAATCAACCTGGCTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACAGGAAAATTATCCTGCTGGAGCCGcgtcg                                                <  1:35117/81‑1 (MQ=255)                           gggAAATCAACCTGGCTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACAGGAAAATTATCCTGCTGGAGCCGcgtcgt                                               >  2:35183/1‑81 (MQ=255)                            ggAAATCAACCTGGCTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACAGGAAAATTATCCTGCTGGAGCCGCgtcgtc                                              >  2:35151/1‑81 (MQ=255)                             gAAATCAACCTGGCTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACAGGAAAATTATCCTGCTGGAGCCGCGTCGTCt                                             >  1:35176/1‑81 (MQ=255)                              aaaTCAACCTGGCTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACAGGAAAATTATCCTGCTGGAGCCGCGTCGTCTg                                            <  2:35118/81‑1 (MQ=255)                               aaTCAACCTGGCTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACAGGAAAATTATCCTGCTGGAGCCGCGTCGTCTgg                                           <  1:35119/81‑1 (MQ=255)                                      cTGGCTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACAGGAAAATTATCCTGCTGGAGCCGCGTCGTCTGGCGGcgcg                                    >  2:35180/1‑81 (MQ=255)                                         gCTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACAGGAAAATTATCCTGCTGGAGCCGCGTCGTCTGGCGGCGCGt                                   <  1:35120/79‑1 (MQ=255)                                         gCTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACAGGAAAATTATCCTGCTGGAGCCGCGTCGTCTGGCGGCGCGTa                                  <  2:35121/80‑1 (MQ=255)                                          cTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACAGGAAAATTATCCTGCTGGAGCCGCGTCGTCTGGCGGCGCGTAAc                                >  2:35175/1‑81 (MQ=255)                                             ccGCTGCAACTGCTGGCGCATCCCGGCATTAACAGGAAAATTATCCTGCTGGAGCCGCGTCGTCTGGCGGCGCGTAACGTc                             >  2:35177/1‑81 (MQ=255)                                               gCTGCAACTGCTGGCGCATCCCGGCATTAACAGGAAAATTATCCTGCTGGAGCCGCGTCGTCTGGCGGCGCGTAACGTcgc                           >  2:35185/1‑81 (MQ=255)                                                   cAACTGCTGGCGCATCCCGGCATTAACAGGAAAATTATCCTGCTGGAGCCGCGTCGTCTGGCGGCGCGTAACGTCgcgc                         <  1:35122/79‑1 (MQ=255)                                                   cAACTGCTGGCGCATCCCGGCATTAACAGGAAAATTATCCTGCTGGAGCCGCGTCGTCTGGCGGCGCGTAACGTCgcgc                         <  2:35123/79‑1 (MQ=255)                                                   cAACTGCTGGCGCATCCCGGCATTAACAGGAAAATTATCCTGCTGGAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCaa                       >  2:35165/1‑81 (MQ=255)                                                     aCTGCTGGCGCATCCCGGCATTAACAGGAAAATTATCCTGCTGGAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACg                     <  2:35124/81‑1 (MQ=255)                                                     aCTGCTGGCGCATCCCGGCATTAACAGGAAAATTATCCTGCTGGAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACg                     >  2:35172/1‑81 (MQ=255)                                                        gctgGCGCATCCCGGCATTAACAGGAAAATTATCCTGCTGGAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCt                  >  2:35182/1‑81 (MQ=255)                                                         ctgGCGCATCCCGGCATTAACAGGAAAATTATCCTGCTGGAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCt                  <  2:35125/80‑1 (MQ=255)                                                         ctgGCGCATCCCGGCATTAACAGGAAAATTATCCTGCTGGAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCt                  <  1:35126/80‑1 (MQ=255)                                                            gcgcATCCCGGCATTAACAGGAAAATTATCCTGCTGGAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCg              >  1:35181/1‑81 (MQ=255)                                                              gcATCCCGGCATTGACAGGAAAATTATCCTGCTGGAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCg              >  2:35160/1‑79 (MQ=255)                                                                 tCCCGGCATTAACAGGAAAATTATCCTGCTGGAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAgct         <  1:35128/81‑1 (MQ=255)                                                                 tCCCGGCATTAACAGGAAAATTATCCTGCTGGAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAgct         <  1:35127/81‑1 (MQ=255)                                                                    cGGCATTAACAGGAAAATTATCCTGCTGGAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAgctgct      >  1:35208/1‑81 (MQ=255)                                                                     ggCATTAACAGGAAAATTATCCTGCTGGAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCtt     <  2:35129/81‑1 (MQ=255)                                                                        aTTAACAGGAAAATTATCCTGCTGGAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAAc  <  2:35130/81‑1 (MQ=255)                                                                              |                                                                           TATTATTAAGTGCGCCGACCGGGGCCGGGAAATCAACCTGGCTGCCGCTGCAACTGCTGGCGCATCCCGGCATTAACGGGAAAATTATCCTGCTGGAGCCGCGTCGTCTGGCGGCGCGTAACGTCGCGCAACGGCTGGCGGAGCTGCTTAAC  >  NC_000913/162172‑162323

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