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%	149.0 / NA	43	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 (21/22);  total (21/22)

TGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGACGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACAGCGGGCGACAGGCGTTGGAGCAGGTGCGGGAACAGGTTTTTGATCTTGTGCTTTGCGATGTGCGAATGGCGGAGATGGACGGCATCGCCACGCTGAAAGAGATCAAAGCGTTAAACCCGGCAATTCCGGTGC  >  NC_000913/4203282‑4203560                                                                                                                                                   |                                                                                                                                     tGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGACGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAAcggcgg                                                                                                                                  >  2:554570/1‑151 (MQ=255)   gCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGACGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACGGCGGGc                                                                                                                                <  1:173954/151‑1 (MQ=255)       ccGGTCAATATTACGCGTAAGGACCCACAAGGATGACGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACGGCGGGCGACa                                                                                                                            <  1:66126/151‑1 (MQ=255)        cGGTCAATATTACGCGTAAGGACCCACAAGGATGACGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACGGCGGGCGACAg                                                                                                                           <  1:494113/151‑1 (MQ=255)                         aaGGACCCACAAGGATGACGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACGGCGGGCGACAGGCGTTGGAGCAGGTGCg                                                                                                          <  1:714278/151‑1 (MQ=255)                           ggACCCACAAGGATGACGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACGGCGGGCGACAGGCGTTGGAGCAGGTGCggg                                                                                                        >  2:235860/1‑151 (MQ=255)                           ggACCCACAAGGATGACGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACGGCGGGCGACAGGCGTTGGAGCAGGTGCggg                                                                                                        <  2:523724/151‑1 (MQ=255)                           ggACCCACAAGGATGACGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACGGCGGGCGACAGGCGTTGGAGCAGGTGCggg                                                                                                        <  1:744372/151‑1 (MQ=255)                             aCCCACAAGGATGACGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACGGCGGGCGACAGGCGTTGGAGCAGGTGCGGGaa                                                                                                      <  1:544941/151‑1 (MQ=255)                                      gATGACGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACGGCGGGCGACAGGCGTTGGAGCAGGTGCGGGAACAGGttttt                                                                                             >  2:528679/1‑151 (MQ=255)                                      gATGACGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACGGCGGGCGACAGGCGTTGGAGCAGGTGCGGGAACAGGttttt                                                                                             <  2:78111/151‑1 (MQ=255)                                           cGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACGGCGGGCGACAGGCGTTGGAGCAGGTGCGGGAACAGGTTTTTGATc                                                                                         >  1:484235/1‑150 (MQ=255)                                           cGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACGGCGGGCGACAGGCGTTGGAGCAGGTGCGGGAACAGGTTTTTGATCt                                                                                        <  2:356908/151‑1 (MQ=255)                                            gCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACGGCGGGCGACAGGCGTTGGAGCAGGTGCGGGAACAGGTTTTTGATCtt                                                                                       <  1:650007/151‑1 (MQ=255)                                            gCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACGGCGGGCGACAGGCGTTGGAGCAGGTGCGGGAACAGGTTTTTGATCtt                                                                                       <  1:845851/151‑1 (MQ=255)                                             cACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACGGCGGGCGACAGGCGTTGGAGCAGGTGCGGGAACAGGTTTTTGATCTtg                                                                                      <  1:340575/151‑1 (MQ=255)                                             cACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACGGCGGGCGACAGGCGTTGGAGCAGGTGCGGGAACAGGTTTTTGATCTtg                                                                                      >  1:307711/1‑151 (MQ=255)                                                gATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACGGCGGGCGACAGGCGTTGGAGCAGGTGCGGGAACAGGTTTTTGATCTTGTGc                                                                                   <  2:765256/151‑1 (MQ=255)                                                 ataataTCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACGGCGGGCGACAGGCGTTGGAGCAGGTGCGGGAACAGGTTTTTGATCTTGTGCt                                                                                  >  1:766080/1‑151 (MQ=255)                                                    atatCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACGGCGGGCGACAGGCGTTGGAGCAGGTGCGGGAACAGGTTTTTGATCTTGTGCTTTg                                                                               <  1:467837/151‑1 (MQ=255)                                                     tatCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACGGCGGGCGACAGGCGTTGGAGCAGGTGCGGGAACAGGTTTTTGATCTTGTGCTTTGc                                                                              >  2:131658/1‑151 (MQ=255)                                                      atCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACGGCGGGCGACAGGCGTTGGAGCAGGTGCGGGAACAGGTTTTTGATCTTGTGCTTTGc                                                                              <  2:871580/150‑1 (MQ=255)                                                        cGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACGGCGGGCGACAGGCGTTGGAGCAGGTGCGGGAACAGGTTTTTGATCTTGTGCTTTGCGAt                                                                           >  2:156744/1‑151 (MQ=255)                                                         gATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACGGCGGGCGACAGGCGTTGGAGCAGGTGCGGGAACAGGTTTTTGATCTTGTGCTTTGCGAtg                                                                          >  1:258960/1‑151 (MQ=255)                                                                        gatgatGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACGGCGGGCGACAGGCGTTGGAGCAGGTGCGGGAACAGGTTTTTGATCTTGTGCTTTGCGATGTGCGAATGGCgg                                                              >  1:753798/1‑148 (MQ=255)                                                                        gatgatGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACGGCGGGCGACAGGCGTTGGAGCAGGTGCGGGAACAGGTTTTTGATCTTGTGCTTTGCGATGTGCGAATGGCGgaga                                                           <  2:15935/151‑1 (MQ=255)                                                                           gatgaCATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACGGCGGGCGACAGGCGTTGGAGCAGGTGCGGGAACAGGTTTTTGATCTTGTGCTTTGCGATGTGCGAATgg                                                                 <  1:470069/142‑1 (MQ=255)                                                                           gatgaCATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACGGCGGGCGACAGGCGTTGGAGCAGGTGCGGGAACAGGTTTTTGATCTTGTGCTTTGCGATGTGCGAATgg                                                                 >  2:470022/1‑142 (MQ=255)                                                                              gaCATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACGGCGGGCGACAGGCGTTGGAGCAGGTGCGGGAACAGGTTTTTGATCTTGTGCTTTGCGATGTGCGAATGGCGGAGATGGACg                                                     <  2:183566/151‑1 (MQ=255)                                                                                  ttAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACGGCGGGCGACAGGCGTTGGAGCAGGTGCGGGAACAGGTTTTTGATCTTGTGCTTTGCGATGTGCGAATGGCGGAGATGGACGGCAt                                                 >  2:618462/1‑151 (MQ=255)                                                                                   tAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACGGCGGGCGACAGGCGTTGGAGCAGGTGCGGGAACAGGTTTTTGATCTTGTGCTTTGCGATGTGCGAATGGCGGAGATGGACGGCAt                                                 >  1:746056/1‑150 (MQ=255)                                                                                       cACTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACGGCGGGCGACAGGCGTTGGAGCAGGTGCGGGAACAGGTTTTTGATCTTGTGCTTTGCGATGTGCGAATGGCGGAGATGGACGGCATCGcc                                             >  2:41780/1‑150 (MQ=255)                                                                                        aCTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACGGCGGGCGACAGGCGTTGGAGCAGGTGCGGGAACAGGTTTTTGATCTTGTGCTTTGCGATGTGCGAATGGCGGAGATGGACGGCATCGCCAc                                           >  1:861609/1‑151 (MQ=255)                                                                                        aCTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACGGCGGGCGACAGGCGTTGGAGCAGGTGCGGGAACAGGTTTTTGATCTTGTGCTTTGCGATGTGCGAATGGCGGAGATGGACGGCATCGCCAc                                           >  1:410664/1‑151 (MQ=255)                                                                                           gCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACGGCGGGCGACAGGCGTTGGAGCAGGTGCGGGAACAGGTTTTTGATCTTGTGCTTTGCGATGTGCGAATGGCGGAGATGGACGGCATCGCCACGCt                                        <  2:245244/151‑1 (MQ=255)                                                                                            cACTATTTTCCAGGCTTTTCTGCGCGGCTGAGGCTATAACGTCGCGCTGGCGAACGGCGGGCGACAGGCGTTGGAGCAGGTGTGGGAACAGGGTTTTGATCTTGTGCTTTGCGATGTGCGAATGGCTGAGATGGACGGCATCGCCACGCTg                                       <  2:534912/151‑1 (MQ=255)                                                                                                         gCTTTACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACGGCGGGCGACAGGCGTTGGAGCAGGTGCGGGAACAGGTTTTTGATCTTGTGCTTTGCGATGTGCGAATGGCGGAGATGGACGGCATCGCCACGCTGAAAGAGATCaa                            >  1:125881/1‑149 (MQ=255)                                                                                                            ttACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACGGCGGGCGACAGGCGTTGGAGCAGGTGCGGGAACAGGTTTTTGATCTTGTGCTTTGCGATGTGCGAATGGCGGAGATGGACGGCATCGCCACGCTGAAAGAGATCAAAGCg                        >  2:436599/1‑150 (MQ=255)                                                                                                               cTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACGGCGGGCGACAGGCGTTGGAGCAGGTGCGGGAACAGGTTTTTGATCTTGTGCTTTGCGATGTGCGAATGGCGGAGATGGACGGCATCGCCACGCTGAAAGAGATCAAAGCGTTa                     >  2:525041/1‑150 (MQ=255)                                                                                                                           ggCTATAACGTCGCGCTGGCGAACGGCGGGCGACAGGCGTTGGAGCAGGTGCGGGAACAGGTTTTTGATCTTGTGCTTTGCGATGTGCGAATGGCGGAGATGGACGGCATCGCCACGCTGAAAGAGATCAAAGCGTTAAACCCGGCAATTc        >  1:808727/1‑151 (MQ=255)                                                                                                                            gCTATAACGTCGCGCTGGCGAACGGCGGGCGACAGGCGTTGGAGCAGGTGCGGGAACAGGTTTTTGATCTTGTGCTTTGCGATGTGCGAATGGCGGAGATGGACGGCATCGCCACGCTGAAAGAGATCAAAGCGTTAAACCCGGCAATTcc       <  2:620831/151‑1 (MQ=255)                                                                                                                            gCTATAACGTCGCGCTGGCGAACGGCGGGCGACAGGCGTTGGAGCAGGTGCGGGAACAGGTTTTTGATCTTGTGCTTTGCGATGTGCGAATGGCGGAGATGGACGGCATCGCCACGCTGAAAGAGATCAAAGCGTTAAACCCGGCAATTcc       <  2:785632/151‑1 (MQ=255)                                                                                                                              tataACGTCGCGCTGGCGAACGGCGGGCGACAGGCGTTGGAGCAGGTGCGGGAACAGGTTTTTGATCTTGTGCTTTGCGATGTGCGAATGGCGGAGATGGACGGCATCGCCACGCTGAAAGAGATCAAAGCGTTAAACCCGGCAATTCCgg     >  1:904724/1‑151 (MQ=255)                                                                                                                                 aaCGTCGCGCTGGCGAACGGCGGGCGACAGGCGTTGGAGCAGGTGCGGGAACAGGTTTTTGATCTTGTGCTTTGCGATGTGCGAATGGCGGAGATGGACGGCATCGCCACGCTGAAAGAGATCAAAGCGTTAAACCCGGCAATTCCGGTGc  >  1:390281/1‑151 (MQ=255)                                                                                                                                                   |                                                                                                                                     TGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGACGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCTATAACGTCGCGCTGGCGAACAGCGGGCGACAGGCGTTGGAGCAGGTGCGGGAACAGGTTTTTGATCTTGTGCTTTGCGATGTGCGAATGGCGGAGATGGACGGCATCGCCACGCTGAAAGAGATCAAAGCGTTAAACCCGGCAATTCCGGTGC  >  NC_000913/4203282‑4203560

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