Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,203,428	A→G	S37G (AGC→GGC)	zraR →	fused DNA‑binding response regulator in two‑component regulatory system with ZraS: response regulator/sigma54 interaction protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,203,428	0	A	G	100.0%	164.3 / NA	49	S37G (AGC→GGC)	zraR	fused DNA‑binding response regulator in two‑component regulatory system with ZraS: response regulator/sigma54 interaction protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (23/26);  total (23/26)

GGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGACGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACAGCGGGCGACAGGCGTTGGAGCAGGTGCGGGAACAGGTTTTTGATCTTGTGCTTTGCGATGTGCGAATGGCGGAGATGGACGGCATCGCCACGCTGAAAGAGATCAAAGCGTTAAACCCGGCAATTCCGGTGCTGATTATGACT  >  NC_000913/4203283‑4203571                                                                                                                                                  |                                                                                                                                                ggCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGACGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAAcggcgg                                                                                                                                             <  2:593880/150‑1 (MQ=255) ggCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGACGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACGGCg                                                                                                                                              <  2:234232/149‑1 (MQ=255)  gCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGACGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACGGCggg                                                                                                                                            >  1:481327/1‑150 (MQ=255)  gCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGACGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACGGCGGGc                                                                                                                                           <  2:435958/151‑1 (MQ=255)             atatTACGCGTAAGGACCCACAAGGATGACGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACGGCGGGCGACAGGCGTTg                                                                                                                                <  2:403399/151‑1 (MQ=255)                  aCGCGTAAGGACCCACAAGGATGACGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACGGCGGGCGACAGGCGTTGGAGCa                                                                                                                           <  2:288733/151‑1 (MQ=255)                    gcgTAAGGACCCACAAGGATGACGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACGGCGGGCGACAGGCGTTGGAGCAgg                                                                                                                         <  2:109158/151‑1 (MQ=255)                         aGGACCCACAAGGATGACGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACGGCGGGCGACAGGCGTTGGAGCAGGTGCgg                                                                                                                    <  1:619867/151‑1 (MQ=255)                             cccACAAGGATGACGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACGGCGGGCGACAGGCGTTGGAGCAGGTGCGGGAAc                                                                                                                <  1:217588/151‑1 (MQ=255)                              ccACAAGGATGACGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACGGCGGGCGACAAGCGTTGGAGCAGGTGCGGGaa                                                                                                                 >  1:569833/1‑149 (MQ=255)                                  aaGGATGACGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACGGCGGGCGACAGGCGTTGGAGCAGGTGCGGGAACAGGtt                                                                                                           >  2:619824/1‑151 (MQ=255)                                   aGGATGACGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACGGCGGGCGACAGGCGTTGGAGCAGGTGCGGGAACAGGttt                                                                                                          <  2:580764/151‑1 (MQ=255)                                    ggATGACGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACGGCGGGCGACAGGCGTTGGAGCAGGTGCGGGAACAGGtttt                                                                                                         >  1:491707/1‑151 (MQ=255)                                    ggATGACGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACGGCGGGCCACAGGCGTTTGAGCAGGTGCGGGAACAGGtttt                                                                                                         >  1:722369/1‑151 (MQ=255)                                       tGACGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACGGCGGGCGACAGGCGTTGGAGCAGGTGCGGGAACAGGTTTTTg                                                                                                       >  1:663216/1‑150 (MQ=255)                                           gCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACGGCGGGCGACAGGCGTTGGAGCAGGTGCGGGAACAGGTTTTTGATCtt                                                                                                  <  2:352081/151‑1 (MQ=255)                                             aCGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACGGCGGGCGACAGGCGTTGGAGCAGGTGCGGGAACAGGTTTTTGATCTtgt                                                                                                <  1:110094/151‑1 (MQ=255)                                               gATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACGGCGGGCGACAGGCGTTGGAGCAGGTGCGGGAACAGGTTTTTGATCTTGTGc                                                                                              <  2:481363/151‑1 (MQ=255)                                                   atatCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACGGCGGGCGACAGGCGTTGGAGCAGGTGCGGGAACAGGTTTTTGATCTTGTGCTTTg                                                                                          >  2:794188/1‑151 (MQ=255)                                                       cGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACGGCGGTCGACAGGCGTTGGAGCAGGTGCGGGAACAGGTTTTTGATCTTGTGCTTTGCGAt                                                                                      >  2:88638/1‑151 (MQ=255)                                                         atatTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACGGCGGGCGACAGGCGTTGGAGCAGGTGCGGGAACAGGTTTTTGATCTTGTGCTTTGCGAtgt                                                                                    <  2:176498/151‑1 (MQ=255)                                                                ggtggtggATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACGGCGGGCGACAGGCGTTGGAGCAGGTGCGGGAACAGGTTTTTGATCTTGTGCTTTGCGATGTGCGAATg                                                                             <  1:441893/151‑1 (MQ=255)                                                                ggtggtggATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACGGCGGGCGACAGGCGTTGGAGCAGGTGCGGGAACAGGTTTTTGATCTTGTGCTTTGCGATGTGCGAATg                                                                             <  2:663258/151‑1 (MQ=255)                                                                    gtggATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACGGCGGGCGACAGGCGTTGGAGCAGGTGCGGGAACAGGTTTTTGATCTTGTGCTTTGCGATGTGCGAATGGCgg                                                                         <  1:843421/151‑1 (MQ=255)                                                                      ggATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACGGCGGGCGACAGGCGTTGGAGCAGGTGCGGGAACAGGTTTTTGATCTTGTGCTTTGCGATGTGCGAATGGCGga                                                                        <  2:708340/150‑1 (MQ=255)                                                                           atgaCATTAGACACTGCACTATTTTCCAGGCTTTACTGAGCGGCAGGGGCTATAACGTCGCGCTGGCGAACGGCGGGGGCCACGCGTTGGAGCAGGCGCGGGAACAGGTTTTTGGTCTTGTGCTTTGCGATGTTCGAATGGCGGAGATgg                                                                   >  1:698398/1‑150 (MQ=255)                                                                            tgaCATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACGGCGGGCGACAGGCGTTGGGGCAGGTGCGGGAACAGGTTTTTGATCTTGTGCTTTGCGATGTGCGAATGGCGGAGATgg                                                                   >  1:552225/1‑149 (MQ=255)                                                                             gaCATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACGGCGGGCGACAGGCGTTGGAGCAGGTGCGGGAACAGGTTTTTGATCTTGTGCTTTGCGATGTGCGAATGGCGGAGATGGACg                                                                >  2:72956/1‑151 (MQ=255)                                                                              aCATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACGGCGGGCGACAGGCGTTGGAGCAGGTGCGGGAACAGGTTTTTGATCTTGTGCTTTGCGATGTGCGAATGGCGGAGATGGACgg                                                               <  2:211328/151‑1 (MQ=255)                                                                                   aGCCCCTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACGGCGGGTGACAGGGGTTGGAGCAGGCGCGGGAACAGGTTTTTGATCTTGTGCTTTGCGGTGTGCGAACGGGGGAGATgg                                                                   >  1:724221/1‑142 (MQ=255)                                                                                   aGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACGGCGGGCGACAGGCGTTGGAGCAGGTGCGGGAACAGGTTTTTGATCTTGTGCTTTGCGATGTGCGAATGGCGGAGATGGACGGCATc                                                           <  2:266551/150‑1 (MQ=255)                                                                                     ccACTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACGGCGGGAGGCTGGCGTTGGAGGAGGTTGGGGATAAGGGTTTGGATCTTGAGCTTCGCGAGCTGCGAACGGCGGAGGGGGGCGGCGTc                                                           >  1:588018/1‑148 (MQ=255)                                                                                               aTTTTGCAGGCTTTACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACGGCGGGCGACAGGCGTTGGAGCAGGTGCGGGAACAGGTTTTTGATCTTGTGCTTTGCGATGTGCGAATGGCGGAGATGGACGGCATCGCCACGCTGAAag                                              <  1:849135/151‑1 (MQ=255)                                                                                               aTTTTGCAGGCTTTACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACGGCGGGCGACAGGCGTTGGAGCAGGTGCGGGAACAGGTTTTTGATCTTGTGCTTTGCGATGTGCGAATGGCGGAGATGGACGGCATCGCCACGCTGAAag                                              <  2:349079/151‑1 (MQ=255)                                                                                                ttttGCAGGCTTTACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACGGCGGGCGACAGGCGTTGGAGCAGGTGCGGGAACAGGTTTTTGCTCTTGTGCTTTGCGTTGTGCGAATGGCGGAGATGGACCGCATCACCACCCTTaaa                                               >  1:398320/1‑149 (MQ=255)                                                                                                 tttGCAGGCTTTACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACGGCGGGCGACAGGCGTTGGAGCAGGTGCGGGAACAGGTTTTTGATCTTGTGCTTTGCGATGTGCGAATGGCGGAGATGGACGGCATCGCCACGCTGAAag                                              >  1:136612/1‑149 (MQ=255)                                                                                                  ttGCAGGCTTTACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACGGCGGGCGACAGGCGTTGGAGCAGGTGCGGGAACAGGTTTTTGATCTTGTGCTTTGCGATGTGCGAATGGCGGAGATGGACGGCATCGCCACGCTGAAag                                              >  2:160479/1‑148 (MQ=255)                                                                                                   tGCAGGCTTTACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACGGCGGGCGACAGGCGTTGGAGCAGGTGCGGGAACAGGTTTTTGATCTTGTGCTTTGCGATGTGCGAATGGCGGAGATGGACGGCATCGCCACGCTGAAAgaga                                           <  1:192857/150‑1 (MQ=255)                                                                                                     cAGGCTTTACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACGGCGGGCGACAGGCGTTGGAGCAGGTGCGGGAACAGGTTTTTGATCTTGTGCTTTGCGATGTGCGAATGGCGGAGATGGACGGCATCGCCACGCTGAAAGAGATc                                         <  2:657461/150‑1 (MQ=255)                                                                                                            tACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACGGCGGGCGACAGGCGTTGGAGCAGGTGCGGGAACAGGTTTTTGATCTTGTGCTTTGCGATGTGCGAATGGCGGAGATGGACGGCATCGCCACGCTGAAAGAGATCAAAGCGtt                                 >  1:140929/1‑151 (MQ=255)                                                                                                            tACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACGGCGGGCGACAGGCGTTGGAGCAGGTGCGGGAACAGGTTTTTGATCTTGTGCTTTGCGATGTGCGAATGGCGGAGATGGACGGCATCGCCACGCTGAAAGAGATCAAAGCGtt                                 >  1:95248/1‑151 (MQ=255)                                                                                                                gcgcGGCTGGGGCTATAACGTCGCGCTGGCGAACGGCGGTCGACAGGCGTTGGAGCAGGTGCGGGAACAGGTTTTTGATCTTGTGCTTTGCGATGTGCGAATGGCGGAGATGGACGGCATCGCCACGCTGAAAGAGATCAAAGCGTTaa                               <  1:88615/149‑1 (MQ=255)                                                                                                                gcgcGGCTGGGGCTATAACGTCGCGCTGGCGAACGGCGGGCGACAGGCGTTGGAGCAGGTGCGGGAACAGGTTTTTGATCTTGTGCTTTGCGATGTGCGAATGGCGGAGATGGACGGCATCGCCACGCTGAAAGAGATCAAAGCGTTAAAc                             <  2:592975/151‑1 (MQ=255)                                                                                                                              ataACGTCGCGCTGGCGAACGGCGGGCGACAGGCGTTGGAGCAGGTGCGGGAACAGGTTTTTGATCTTGTGCTTTGCGATGTGCGAATGGCGGAGATGGACGGCATCGCCACGCTGAAAGAGATCAAAGCGTTa                                <  2:532164/134‑1 (MQ=255)                                                                                                                              ataACGTCGCGCTGGCGAACGGCGGGCGACAGGCGTTGGAGCAGGTGCGGGAACAGGTTTTTGATCTTGTGCTTTGCGATGTGCGAATGGCGGAGATGGACGGCATCGCCACGCTGAAAGAGATCAAAGCGTTAAACCCGGCAATTCCGGt               >  2:77444/1‑151 (MQ=255)                                                                                                                              ataACGTCGCGCTGGCGAACGGCGGGCGACAGGCGTTGGAGCAGGTGCGGGAACAGGGTTTTGATCTTGTGCTTTGCGATGTGCGAATGGCGGAGATGGACGGCATCGCCACGCTGAAAGAGATCAAAGCGTTa                                >  1:532127/1‑134 (MQ=255)                                                                                                                                 aCGTCGCGCTGGCGAACGGCGGGCGACAGGCGTTGGAGCAGGTGCGGGAACAGGTTTTTGATCTTGTGCTTTGCGATGTGCGAATGGCGGAGATGGACGGCATCGCCACGCTGGAAGAGATCAAAGCGTTAAACCCGGCAATTcc                  >  2:590534/1‑145 (MQ=255)                                                                                                                                 aCGTCGCGCTGGCGAACGGCGGGCGACAGGCGTTGGAGCAGGTGCGGGAACAGGTTTTTGATCTTGTGCTTTGCGATGTGCGAATGGCGGAGATGGACGGCATCGCCACGCTGAAAGAGATCAAAGCGTTAAACCCGGCAATTCCGGTg              >  1:698124/1‑149 (MQ=255)                                                                                                                                         cTGGCGAACGGCGGGCGACAGGCGTTGGAGCAGGTGCGGGAACAGGTTTTTGATCTTGTGCTTTGCGATGTGCGAATGGCGGAGATGGACGGCATCGCCACGCTGAAAGAGATCAAAGCGTTAAACCCGGCAATTCCGGTGCTGATTATg     >  2:314976/1‑150 (MQ=255)                                                                                                                                           ggCGAACGGCGGGCGACAGGCGTTGGAGCAGGTGCGGGAACAGGTTTTTGATCTTGTGCTTTGCGATGTGCGAATGGCGGAGATGGACGGCATCGCCACGCTGAAAGAGATCAAAGCGTTAAACCCGGCAATTCCGGTGCTGATTATGACt  <  2:455504/151‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                GGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGACGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACAGCGGGCGACAGGCGTTGGAGCAGGTGCGGGAACAGGTTTTTGATCTTGTGCTTTGCGATGTGCGAATGGCGGAGATGGACGGCATCGCCACGCTGAAAGAGATCAAAGCGTTAAACCCGGCAATTCCGGTGCTGATTATGACT  >  NC_000913/4203283‑4203571

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