Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	163,723	G→A	G540D (GGC→GAC)	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	163,723	0	G	A	100.0%	206.0 / NA	59	G540D (GGC→GAC)	hrpB	putative ATP‑dependent RNA helicase HrpB
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (23/36);  total (23/36)

GCGTTTGCCGATCGCATTGCTCGTCGCCGTGGGCAAGATGGACGCTATCAACTGGCAAACGGCATGGGAGCGATGCTCGATGCCAACGACGCGCTAAGCCGCCACGAATGGTTGATCGCACCGTTATTATTGCAGGGCAGCGCCTCGCCGGATGCGCGGATTTTACTGGCGCTGCTGGTCGATATTGATGAGTTAGTACAACGCTGCCCGCAGCTGGTACAGCAGTCTGACACTGTGGAGTGGGATGACGCGCAAGGTACGCTGAAAGCCT  >  NC_000913/163587‑163857                                                                                                                                         |                                                                                                                                       gCGTTTGCCGATCGCATTGCTCGTCGCCGTGGGCAAGATGGACGCTATCAACTGGCAAACGGCATGGGAGCGATGCTCGATGCCAACGACGCGCTAAGCCGCCACGAATGGTTGATCGCACCGTTATTATTGCAGGACAg                                                                                                                                     <  2:956201/140‑1 (MQ=255)  cGTTTGCCGATCGCATTGCTCGTCGCCGTGGGCAAGATGGACGCTATCAACTGGCAAACGGCATGGGAGCGATGCTCGATGCCAACGACGCGCTAAGCCGCCACGAATGGTTGATCGCACCGTTATTATTGCAGGACAgc                                                                                                                                    <  2:816931/140‑1 (MQ=255)      tGCCGATCGCATTGCTCGTCGCCGTGGGCAAGATGGACGCTATCAACTGGCAAACGGCATGGGAGCGATGCTCGATGCCAACGACGCGCTAAGCCGCCACGAATGGTTGATCGCACCGTTATTATTGCAGGACAGCGCCt                                                                                                                                >  2:83914/1‑140 (MQ=255)           aTCGCATTGCTCGTCGCCGTGGGCAAGATGGACGCTATCAACTGGCAAACGGCATGGGAGCGATGCTCGATGCCAACGACGCGCTAAGCCGCCACGAATGGTTGATCGCACCGTTATTATTGCAGGACAGCGCCTCGCCg                                                                                                                           <  1:83914/140‑1 (MQ=255)            tCGCATTGCTCGTCGCCGTGGGCAAGATGGACGCTATCAACTGGCAAACGGCATGGGAGCGATGCTCGATGCCAACGACGCGCTAAGCCGCCACGAATGGTTGATCGCACCGTTATTATTGCAGGACAGCGCCTCGCCgg                                                                                                                          <  2:1097065/140‑1 (MQ=255)                              tGGGCAAGATGGACGCTATCAACTGGCAAACGGCATGGGAGCGATGCTCGATGCCAACGACGCGCTAAGCCGCCACGAATGGTTGATCGCACCGTTATTATTGCAGGACAGCGCCTCGCCGGATGCGCGGATTTTACTgg                                                                                                        <  2:512704/140‑1 (MQ=255)                              tGGGCAAGATGGACGCTATCAACCGGCAAACGGCATGGGAGCGATGCTCGATGCCAACGACGCGCTAAGCCGCCACGAATGGTTGATCGCACCGTTATTATTGCAGGACAGCGCCTAGCCGGATGCGCGGATTTTACTgg                                                                                                        >  2:60204/1‑140 (MQ=255)                                ggCAAGATGGACGCTATCAACTGGCAAACGGCATGGGAGCGATGCTCGATGCCAACGACGCGCTAAGCCGCCACGAATGGTTGATCGCACCGTTATTATTGCAGGACAGCGCCTCGCCGGATGCGCGGATTTTACTGgcg                                                                                                      <  1:818543/140‑1 (MQ=255)                                ggCAAGATGGACGCTATCAACTGGCAAACGGCATGGGAGCGATGCTCGATGCCAACGACGCGCTAAGCCGCCACGAATGGTTGATCGCACCGTTATTATTGCAGGACAGCGCCTCGCCGGATGCGCGGATTTTACTGgcg                                                                                                      >  2:518739/1‑140 (MQ=255)                                   aaGATGGACGCTATCAACTGGCAAACGGCATGGGAGCGATGCTCGATGCCAACGACGCGCTAAGCCGCCACGAATGGTTGATCGCACCGTTATTATTGCAGGACAGCGCCTCGCCGGATGCGCGGATTTTACTGGCgctg                                                                                                   <  1:556010/140‑1 (MQ=255)                                   aaGATGGACGCTATCAACTGGCAAACGGCATGGGAGCGATGCTCGATGCCAACGACGCGCTAAGCCGCCACGAATGGTTGATCGCACCGTTATTATTGCAGGACAGCGCCTCGCCGGATGCGCGGATTTTACTGGCgctg                                                                                                   <  1:948394/140‑1 (MQ=255)                                     gATGGACGCTATCAACTGGCAAACGGCATGGGAGCGATGCTCGATGCCAACGACGCGCTAAGCCGCCACGAATGGTTGATCGCACCGTTATTATTGCAGGACAGCGCCTCGCCGGATGCGCGGATTTTACTGGCgctgct                                                                                                 >  2:357503/1‑140 (MQ=255)                                      aTGGACGCTATCAACTGGCAAACGGCATGGGAGCGATGCTCGATGCCAACGACGCGCTAAGCCGCCACGAATGGTTGATCGCACCGTTATTATTGCAGGACAGCGCCTCGCCGGATGCGCGGATTTTACTGGCGctgctg                                                                                                <  2:736780/140‑1 (MQ=255)                                        ggACGCTATCAACTGGCAAACGGCATGGGAGCGATGCTCGATGCCAACGACGCGCTAAGCCGCCACGAATGGTTGATCGCACCGTTATTATTGCAGGACAGCGCCTCGCCGGATGCGCGGATTTTACTGGCGCTGCTGGt                                                                                              <  2:224516/140‑1 (MQ=255)                                          aCGCTATCAACTGGCAAACGGCATGGGAGCGATGCTCGATGCCAACGACGCGCTAAGCCGCCACGAATGGTTGATCGCACCGTTATTATTGCAGGACAGCGCCTCGCCGGATGCGCGGATTTTACTGGCGCTGCTGGTCg                                                                                            <  2:28734/140‑1 (MQ=255)                                          aCGCTATCAACTGGCAAACGGCATGGGAGCGATGCTCGATGCCAACGACGCGCTAAGCCGCCACGAATGGTTGATCGCACCGTTATTATTGCAGGACAGCGCCTCGCCGGATGCGCGGATTTTACTGGCGCTGCTGGTCg                                                                                            >  2:551466/1‑140 (MQ=255)                                          aCGCTATCAACTGGCAAACGGCATGGGAGCGATGCTCGATGCCAACGACGCGCTAAGCCGCCACGAATGGTTGATCGCACCGTTATTATTGCAGGACAGCGCCTCGCCGGATGCGCGGATTTTACTGGCGCTGCTGGTCg                                                                                            >  1:771345/1‑140 (MQ=255)                                          aCGCTATCAACTGGCAAACGGCATGGGAGCGATGCTCGATGCCAACGACGCGCTAAGCCGCCACGAATGGTTGATCGCACCATTATTATTGCAGGACAGTGCCTCGCCGGATGCGCGGATTTTACTGGCGCTGCTGGTCg                                                                                            <  2:1018481/140‑1 (MQ=255)                                            gCTATCAACTGGCAAACGGCATGGGAGCGATGCTCGATGCCAACGACGCGCTAAGCCGCCACGAATGGTTGATCGCACCGTTATTATTGCAGGACAGCGCCTCGCCGGATGCGCGGATTTTACTGGCGCTGCTGGTCGat                                                                                          <  1:12813/140‑1 (MQ=255)                                            gCTATCAACTGGCAAACGGCATGGGAGCGATGCTCGATGCCAACGACGCGCTAAGCCGCCACGAATGGTTGATCGCACCGTTATTATTGCAGGACAGCGCCTCGCCGGATGCGCGGATTTTACTGGCGCTGCTGGTCGat                                                                                          <  1:949751/140‑1 (MQ=255)                                               aTCAACTGGCAAACGGCATGGGAGCGATGCTCGATGCCAACGACGCGCTAAGCCGCCACGAATGGTTGATCGCACCGTTATTATTGCAGGACAGCGCCTCGCCGGATGCGCGGATTTTACTGGCGCTGCTGGTCGATAtt                                                                                       >  2:1101070/1‑140 (MQ=255)                                                 cAACTGGCAAACGGCATGGGAGCGATGCTCGATGCCAACGACGCGCTAAGCCGCCACGAATGGTTGATCGCACCGTTATTATTGCAGGACAGCGCCTCGCCGGATGCGCGGATTTTACTGGCGCTGCTGGTCGATATtga                                                                                     <  2:135808/140‑1 (MQ=255)                                                 cAACTGGCAAACGGCATGGGAGCGATGCTCGATGCCAACGACGCGCTAAGCCGCCACGAATGGTTGATCGCACCGTTATTATTGCAGGACAGCGCCTCGCCGGATGCGCGGATTTTACTGGCGCTGCTGGTCGATATtga                                                                                     <  1:761290/140‑1 (MQ=255)                                                 cAACTGGCAAACGGCATGGGAGCGATGCTCGATGCCAACGACGCGCTAAGCCGCCACGAATGGTTGATCGCACCGTTATTATTGCAGGACAGCGCCTCGCCGGATGCGCGGATTTTACTGGCGCTGCTGGTCGATATtga                                                                                     <  1:357503/140‑1 (MQ=255)                                                 cAACTGGCAAACGGCATGGGAGCGATGCTCGATGCCAACGACGCGCTAAGCCGCCACGAATGGTTGATCGCACCGTTATTATTGCAGGACAGCGCCTCGCCGGATGCGCGGATTTTACTGGCGCTGCTGGTCGATATtga                                                                                     <  1:15608/140‑1 (MQ=255)                                                   aCTGGCAAACGGCATGGGAGCGATGCTCGATGCCAACGACGCGCTTAGCCGCCACGAATGGTTGATCGCACCGTTATTATTGCAGGACAGCGCCTCGCCGGATGTGCGGATTTTACTGGCGCTGCTGGTCGATATtgatg                                                                                   >  2:781748/1‑140 (MQ=255)                                                        cAAACGGCATGGGAGCGATGCTCGATGCCAACGACGCGCTAAGCCGCCACGAATGGTTGATCGCACCGTTATTATTGCAGGACAGCGCCTCGCCGGATGCGCGGATTTTACTGGCGCTGCTGGTCGATATTGATGAGTTa                                                                              >  1:233663/1‑140 (MQ=255)                                                            cGGCATGGGAGCGATGCTCGATGCCAACGACGCGCTAAGCCGCCACGAATGGTTGATCGCACCGTTATTATTGCAGGACAGCGCCTCGCCGGATGCGCGGATTTTACTGGCGCTGCTGGTCGATATtga                                                                                     >  1:1092169/1‑129 (MQ=255)                                                            cGGCATGGGAGCGATGCTCGATGCCAACGACGCGCTAAGCCGCCACGAATGGTTGATCGCACCGTTATTATTGCAGGACAGCGCCTCGCCGGATGCGCGGATTTTACTGGCGCTGCTGGTCGATATtga                                                                                     <  2:1092169/129‑1 (MQ=255)                                                            cGGCATGGGAGCGATGCTCGATGCCAACGACGCGCTAAGCCGCCACGAATGGTTGATCGCACCGTTATTATTGCAGGACAGCGCCTCGCCGGATGCGCGGATTTTACTGGCGCTGCTGGTCGATATTGATGAGt                                                                                >  2:138743/1‑134 (MQ=255)                                                            cGGCATGGGAGCGATGCTCGATGCCAACGACGCGCTAAGCCGCCACGAATGGTTGATCGCACCGTTATTATTGCAGGACAGCGCCTCGCCGGATGCGCGGATTTTACTGGCGCTGCTGGTCGATATTGATGAGt                                                                                <  1:138743/134‑1 (MQ=255)                                                                 tGGGAGCGATGCTCGATGCCAACGACGCGCTAAGCCGCCACGAATGGTTGATCGCACCGTTATTATTGCAGGACAGCGCCTCGCCGGATGCGCGGATTTTACTGGCGCTGCTGGTCGATATTGATGAGTTAGTACAACGc                                                                     >  1:438587/1‑140 (MQ=255)                                                                     aGCGATGCTCGATGCCAACGACGCGCTAAGCCGCCACGAATGGTTGATCGCACCGTTATTATTGCAGGACAGCGCCTAGCCGGATGCGCGGATTTTACTGGCGCTGCTGGTCGATATTGATGAGTTAGTACAACGCTGcc                                                                 <  1:60204/140‑1 (MQ=255)                                                                      gCGATGCTCGATGCCAACGACGCGCTAAGCCGCCACGAATGGTTGATCGCACCGTTATTATTGCAGGACAGCGCCTCGCCGGATGCGCGGATTTTACTGGCGCTGCTGGTCGATATTGATGAGTTAGTACAACGCTGccc                                                                <  2:771345/140‑1 (MQ=255)                                                                         aTGCTCGATGCCAACGACGCGCTAAGCCGCCACGAATGGTTGATCGCACCGTTATTATTGCAGGACAGCGCCTCGCCGGATGCGCGGATTTTACTGGCGCTGCTGGTCGATATTGATGAGTTAGTACAACGCTGccc                                                                <  1:88729/137‑1 (MQ=255)                                                                         aTGCTCGATGCCAACGACGCGCTAAGCCGCCACGAATGGTTGATCGCACCGTTATTATTGCAGGACAGCGCCTCGCCGGATGCGCGGATTTTACTGGCGCTGCTGGTCGATATTGATGAGTTAGTACAACGCTGccc                                                                >  2:88729/1‑137 (MQ=255)                                                                         aTGCTCGATGCCAACGACGCGCTAAGCCGCCACGAATGGTTGATCGCACCGTTATTATTGCAGGACAGCGCCTCGCCGGATGCGCGGATTTTACTGGCGCTGCTGGTCGATATTGATGAGTTAGTACAACGCTGCCCGCa                                                             <  1:518739/140‑1 (MQ=255)                                                                                   ccAACGGCGCGCTAAGCCGCCACGAATGGTTGATCGCACCGTTATTATTGCAGGACAGCGCCTCGCCGGATGCGCGGATTTTACTGGCGCTGCTGGTCGATATTGATGAGTTAGTACAACGCTGCCCGCAGCTGGTAcag                                                   <  2:942583/140‑1 (MQ=255)                                                                                   ccAACGACGCGCTAAGCCGCCACGAATGGTTGATCGCACCGTTATTATTGCAGGACAGCGCCTCGCCGGATGCGCGGATTTTACTGGCGCTGCTGGTCGATATTGATGAGTTAGTACAACGCTGCCCGCAGCTGGTAcag                                                   <  1:551466/140‑1 (MQ=255)                                                                                      acgacgCGCTAAGCCGCCACGAATGGTTGATCGCACCGTTATTATTGCAGGACAGCGCCTCGCCGGATGCGCGGATTTTACTGGCGCTGCTGGTCGATATTGATGAGTTAGTACAACGCTGCCCGCAGCTGGTAcagcag                                                <  1:586823/140‑1 (MQ=255)                                                                                        gacgCGCTAAGCCGCCACGAATGGTTGATCGCACCGTTATTATTGCAGGACAGCGCCTCGCCGGATGCGCGGATTTTACTGGCGCTGCTGGTCGATATTGATGAGTTAGTACAACGCTGCCCGCAGCTGGTACAGCAGTc                                              <  1:1101070/140‑1 (MQ=255)                                                                                          cgcgCTAAGCCGCCACGAATGGTTGATCGCACCGTTATTATTGCAGGACAGCGCCTCGCCGGATGCGCGGATTTTACTGGCGCTGCTGGTCGATATTGATGAGTTAGTACAACGCTGCCCGCAGCTGGTACAGCAGTCTg                                            <  2:1114983/140‑1 (MQ=255)                                                                                             gcTAAGCCGCCACGAATGGTTGATCGCACCGTTATTATTGCAGGACAGCGCCTCGCCGGATGCGCGGATTTTACTGGCGCTGCTGGTCGATATTGATGAGTTAGTACAACGCTGCCCGCAGCTGGTAc                                                     <  2:1006344/128‑1 (MQ=255)                                                                                             gcTAAGCCGCCACGAATGGTTGATCGCACCGTTATTATTGCAGGACAGCGCCTCGCCGGATGCGCGGATTTTACTGGCGCTGCTGGTCGATATTGATGAGTTAGTACAACGCTGCCCGCAGCTGGTAc                                                     >  1:1006344/1‑128 (MQ=255)                                                                                             gcTAAGCCGCCACGAATGGTTGATCGCACCGTTATTATTGCAGGACAGCGCCTCGCCGGATGCGCGGATTTTACTGGCGCTGCTGGTCGATATTGATGAGTTAGTACAACGCTGCCCGCAGCTGGTACAGCAGTCTGaca                                         <  2:233663/140‑1 (MQ=255)                                                                                                aaGCCGCCACGAATGGTTGATCGCACCGTTATTATTGCAGGACAGCGCCTCGCCGGATGCGCGGATTTTACTGGCGCTGCTGGTCGATATTGATGAGTTAGTACAACGCTGCCCGCAGCTGGTACAGCAGTCTGACACtg                                      <  1:452301/140‑1 (MQ=255)                                                                                                      ccACGAATGGTTGATCGCACCGTTATTATTGCAGGACAGCGCCTCGCCGGATGCGCGGATTTTACTGGCGCTGCTGGTCGATATTGATGAGTTAGTACAACGCTGCCCGCAGCTGGTACAGCAGTCTg                                            <  2:502191/128‑1 (MQ=255)                                                                                                      ccACGAATGGTTGATCGCACCGTTATTATTGCAGGACAGCGCCTCGCCGGATGCGCGGATTTTACTGGCGCTGCTGGTCGATATTGATGAGTTAGTACAACGCTGCCCGCAGCTGGTACAGCAGTCTg                                            >  1:502191/1‑128 (MQ=255)                                                                                                       cACGAATGGTTGATCGCACCGTTATTATTGCAGGACAGCGCCTCGCCGGATGCGCGGATTTTACTGGCGCTGCTGGTCGATATTGATGAGTTAGTACAACGCTGCCCGCAGCTGGTACAGCAGTCTGACACTGTGGAGTg                               >  1:421720/1‑140 (MQ=255)                                                                                                           aaTGGTTGATCGCACCGTTATTATTGCAGGACAGCGCCTCGCCGGATGCGCGGATTTTACTGGCGCTGCTGGTCGATATTGATGAGTTAGTACAACGCTGCCCGCAGCTGGTACAGCAGTCTGACACTGTGGAGTGGGAt                           >  1:725596/1‑140 (MQ=255)                                                                                                            aTGGTTGATCGCACCGTTATTATTGCAGGACAGCGCCTCGCCGGATGCGCGGATTTTACTGGCGCTGCTGGTCGATATTGATGAGTTAGTACAACGCTGCCCGCAGCTGGTACAGCAGTCTGACACTGTGGAGTGGGATg                          <  2:421720/140‑1 (MQ=255)                                                                                                                 tGATCGCACCGTTATTATTGCAGGACAGCGCCTCGCCGGATGCGCGGATTTTACTGGCGCTGCTGGTCGATATTGATGAGTTAGTACAACGCTGCCCGCAGCTGGTACAGCAGTCTGACACTGTGGAGTGGGATGAcgcg                     <  2:438587/140‑1 (MQ=255)                                                                                                                 tGATCGCACCGTTATTATTGCAGGACAGCGCCTCGCCGGATGCGCGGATTTTACTGGCGCTGCTGGTCGATATTGATGAGTTAGTACAACGCTGCCCGCAGCTGGTACAGCAGTCTGACACTGTGGAGTGGGATGAcgcg                     >  1:315062/1‑140 (MQ=255)                                                                                                                   aTCGCACCGTTATTATTGCAGGACAGCGCCTCGCCGGATGCGCGGATTTTACTGGCGCTGCTGGTCGATATTGATGAGTTAGTACAACGCTGCCCGCAGCTGGTACAGCAGTCTGACACTGTGGAGTGGGATGACGCGCa                   >  2:1136298/1‑140 (MQ=255)                                                                                                                   aTCGCACCGTTATTATTGCAGGACAGCGCCTCGCCGGATGCGCGGATTTTACTGGCGCTGCTGGTCGATATTGATGAGTTAGTACAACGCTGCCCGCAGCTGGTACAGCAGTCTGACACTGTGGAGTGGGATGACGCGCa                   >  1:537813/1‑140 (MQ=255)                                                                                                                        aCCGTTATTATTGCAGGACAGCGCCTCGCCGGATGCGCGGATTTTACTGGCGCTGCTGGTCGATATTGATGAGTTAGTACAACGCTGCCCGCAGCTGGTACAGCAGTCTGACACTGTGGAGTGGGATGACGCGCAAGGTa              >  1:246018/1‑140 (MQ=255)                                                                                                                        aCCGTTATTATTGCAGGACAGCGCCTCGCCGGATGCGCGGATTTTACTGGCGCTGCTGGTCGATATTGATGAGTTAGTACAACGCTGCCCGCAGCTGGTACAGCAGTCTGACACTGTGGAGTGGGATGACGCGCAAGGTa              >  1:546304/1‑140 (MQ=255)                                                                                                                             tattatTGCAGGACAGCGCCTCGCCGGATGCGCGGATTTTACTGGCGCTGCTGGTCGATATTGATGAGTTAGTACAACGCTGCCCGCAGCTGGTACAGCAGTCTGACACTGTGGAGTGGGATGACGCGCAAGGTACGCTg         <  2:725596/140‑1 (MQ=255)                                                                                                                                    gCAGGACAGCGCCTCGCCGGATGCGCGGATTTTACTGGCGCTGCTGGTCGATATTGATGAGTTAGTACAACGCTGCCCGCAGCTGGTACAGCAGTCTGACACTGTGGAGTGGGATGACGCGCAAGGTACGCTGAAAGCCt  >  2:1125696/1‑140 (MQ=255)                                                                                                                                    gCAGGACAGCGCCTCGCCGGATGCGCGGATTTTACTGGCGCTGCTGGTCGATATTGATGAGTTAGTACAACGCTGCCCGCAGCTGGTACAGCAGTCTGACACTGTGGAGTGGGATGACGCGCAAGGTACGCTGAAAGCCt  >  2:154539/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                                       GCGTTTGCCGATCGCATTGCTCGTCGCCGTGGGCAAGATGGACGCTATCAACTGGCAAACGGCATGGGAGCGATGCTCGATGCCAACGACGCGCTAAGCCGCCACGAATGGTTGATCGCACCGTTATTATTGCAGGGCAGCGCCTCGCCGGATGCGCGGATTTTACTGGCGCTGCTGGTCGATATTGATGAGTTAGTACAACGCTGCCCGCAGCTGGTACAGCAGTCTGACACTGTGGAGTGGGATGACGCGCAAGGTACGCTGAAAGCCT  >  NC_000913/163587‑163857

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑