Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	164,554	C→T	intergenic (+20/‑176)	hrpB → / → mrcB	putative ATP‑dependent RNA helicase HrpB/peptidoglycan glycosyltransferase/peptidoglycan DD‑transpeptidase MrcB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	164,554	0	C	T	100.0%	201.9 / NA	57	intergenic (+20/‑176)	hrpB/mrcB	putative ATP‑dependent RNA helicase HrpB/peptidoglycan glycosyltransferase/peptidoglycan DD‑transpeptidase MrcB
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (23/34);  total (23/34)

TGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTC  >  NC_000913/164418‑164668                                                                                                                                         |                                                                                                                   tGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCTCtg                                                                                                                 <  2:856793/140‑1 (MQ=255)      aGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCTCTGTGGga                                                                                                            <  2:411546/140‑1 (MQ=255)        gAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCTCTGTGGgaga                                                                                                          >  1:1110626/1‑140 (MQ=255)        gAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCTCTGTGGgaga                                                                                                          >  1:523741/1‑140 (MQ=255)        gAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCTCTGTGGgaga                                                                                                          >  1:270823/1‑140 (MQ=255)         aGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCTCTGTGGGagag                                                                                                         >  1:746136/1‑140 (MQ=255)           cGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCTCTGTGGGAGAggg                                                                                                       >  1:677447/1‑140 (MQ=255)           cGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCTCTGTGGGAGAggg                                                                                                       >  2:155217/1‑140 (MQ=255)            gTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCTCTGTGGGAGAGGGc                                                                                                      <  2:257601/140‑1 (MQ=255)             tACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCTCTGTGGGAGAGGGcc                                                                                                     >  2:719661/1‑140 (MQ=255)             tACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCTCTGTGGGAGAGGGcc                                                                                                     >  1:665673/1‑140 (MQ=255)                  tGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAGATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCTCTGTGGGAGAGGGCCGGGGt                                                                                                <  1:360251/140‑1 (MQ=255)                  tGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCTCTGTGGGAGAGGGCCGGGGt                                                                                                >  1:1005221/1‑140 (MQ=255)                    aGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCTCTGTGGGAGAGGGCCGGGGTGa                                                                                              >  1:523759/1‑140 (MQ=255)                         cAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCTCTGTGGGAGAGGGCCGGGGTGAGGGCa                                                                                         >  1:198130/1‑140 (MQ=255)                                agaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCTCTGTGGGAGAGGGCCGGGGTGAGGGCATCAgcgc                                                                                  >  2:114767/1‑140 (MQ=255)                                 gaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCTCTGTGGGAGAGGGCCGGGGTGAGGGCATCAgcgcg                                                                                 <  2:1005221/140‑1 (MQ=255)                                       aGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCTCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtt                                                                           <  2:270823/140‑1 (MQ=255)                                       aGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCTCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtt                                                                           <  1:306480/140‑1 (MQ=255)                                       aGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCTCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtt                                                                           <  1:719661/140‑1 (MQ=255)                                       aGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCTCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtt                                                                           <  2:890974/140‑1 (MQ=255)                                          gCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACG‑AAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCTCTGTGGGAGAGGGCCGGGGTGAGGGCATCAgcgcg                                                                                 >  2:686964/1‑130 (MQ=255)                                          gCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACG‑AAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCTCTGTGGGAGAGGGCCGGGGTGAGGGCATCAgcgcg                                                                                 <  1:686964/130‑1 (MQ=255)                                           cGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCTCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcacc                                                                       >  1:996685/1‑140 (MQ=255)                                             ttATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCTCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccct                                                                     <  2:523759/140‑1 (MQ=255)                                             ttATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCTCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccct                                                                     <  2:792223/140‑1 (MQ=255)                                             ttATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCTCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccct                                                                     >  2:966187/1‑140 (MQ=255)                                                   cAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCTCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTcaccctcaccct                                                               <  1:114767/140‑1 (MQ=255)                                                    aaaCATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCTCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTa                                                              >  2:940767/1‑140 (MQ=255)                                                       cATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCTCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAAcc                                                           <  2:996685/140‑1 (MQ=255)                                                        aTGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCTCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAAccc                                                          <  1:672936/140‑1 (MQ=255)                                                          gTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCTCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcgccct                                                                     <  2:1196120/127‑1 (MQ=255)                                                          gTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCTCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAAccc                                                          <  1:767383/138‑1 (MQ=255)                                                          gTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCTCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAAccc                                                          >  2:767383/1‑138 (MQ=255)                                                          gTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCTCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCctc                                                        <  1:885682/140‑1 (MQ=255)                                                          gTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCTCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCctc                                                        <  2:1209764/140‑1 (MQ=255)                                                          gTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGCATTCGTAAGGCTGGTTTTCTCCCTCTCTCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccct                                                                     >  1:1196120/1‑127 (MQ=255)                                                            cTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCTCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTcaccctcaccct                                                               >  1:16616/1‑131 (MQ=255)                                                            cTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCTCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTcaccctcaccct                                                               <  2:16616/131‑1 (MQ=255)                                                            cTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCTCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCtctc                                                      <  2:56044/140‑1 (MQ=255)                                                            cTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCTCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCtctc                                                      >  2:18782/1‑140 (MQ=255)                                                                 cGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCTCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCa                                                 <  1:940767/140‑1 (MQ=255)                                                                 cGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCTCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCa                                                 <  1:18782/140‑1 (MQ=255)                                                                  ggACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCTCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCCCCCTCTCCCTAACCCTCTCCCTCaa                                                <  1:352303/140‑1 (MQ=255)                                                                   gacgacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCTCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                               <  1:1175637/140‑1 (MQ=255)                                                                   gacgacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCTCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                               <  2:1110626/140‑1 (MQ=255)                                                                         ccGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCTCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgagaa                                         <  2:12565/140‑2 (MQ=255)                                                                                      cACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCTCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGAc                                    <  1:239532/132‑1 (MQ=255)                                                                                      cACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCTCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGAc                                    >  2:239532/1‑132 (MQ=255)                                                                                      cACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCTCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAg                            >  1:284631/1‑140 (MQ=255)                                                                                      cACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCTCTGTGGGAGAGGGCCGGGATGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAg                            >  1:534004/1‑140 (MQ=255)                                                                                                     cgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCTCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgagag             <  2:665673/140‑1 (MQ=255)                                                                                                     cgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCTCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgagag             <  2:284631/140‑1 (MQ=255)                                                                                                        aaaaGTATTCGTAAGGTTGGTTTTCTCCCTCTCTCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGata          <  2:198130/140‑1 (MQ=255)                                                                                                         aaaGTATTCGTAAGGTTGGTTTTCTCCCTCTCTCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGatat         <  2:523741/140‑1 (MQ=255)                                                                                                          aaGTATTCGTAAGGTTGGTTTTCTCCCTCTCTCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtc        <  2:1023813/140‑1 (MQ=255)                                                                                                               ttCGTAAGGTTGGTTTTCTCCCTCTCTCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttctt   <  2:1209588/140‑1 (MQ=255)                                                                                                                tCGTAAGGTTGGTTTTCTCCCTCTCTCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttcttc  <  1:1032132/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                   TGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTC  >  NC_000913/164418‑164668

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