Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,636,673	G→A	G223R (GGA→AGA)	creB →	response regulator in two‑component regulatory system with CreC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,636,673	0	G	A	100.0%	159.9 / NA	47	G223R (GGA→AGA)	creB	response regulator in two‑component regulatory system with CreC
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (24/23);  total (24/23)

CAACTGATGGATAGCGTATGGGAAGATGCGCAGGACACCTACGATCGCACCGTCGATACCCACATTAAAACGCTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAG  >  NC_000913/4636538‑4636802                                                                                                                                        |                                                                                                                                  cAACTGATGGATAGCGTATGGGAAGATGCGCAGGACACCTACGATCGCACCGTCGATACCCACATTAAAACGCTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGAGat                                                                                                                                >  2:496300/1‑139 (MQ=255)      gATGGATAGCGTATGGGAAGATGCGCAGGACACCTACGATCGCACCGTCGATACCCACATTAAAACGCTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGAGATATAGc                                                                                                                           <  1:676865/139‑1 (MQ=255)          gATAGCGTATGGGAAGATGCGCAGGACACCTACGATCGCACCGTCGATACCCACATTAAAACGCTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGAGATATAGCCTGa                                                                                                                       <  1:496300/139‑1 (MQ=255)               cGTATGGGAAGATGCGCAGGACACCTACGATCGCACCGTCGATACCCACATTAAAACGCTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGAGATATAGCCTGAGGGGc                                                                                                                  >  2:1183512/1‑139 (MQ=255)                      gAAGATGCGCAGGACACCTACGATCGCACCGTCGATACCCACATTAAAACGCTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGAGATATAGCCTGAGGGGCCTGTAAt                                                                                                           <  1:158329/139‑1 (MQ=255)                         gATGCGCAGGACACCTACGATCGCACCGTCGATACCCACATTAAAACGCTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGAGATATAGCCTGAGGGGCCTGTAATGCg                                                                                                        >  1:811371/1‑139 (MQ=255)                            gcgcAGGACACCTACGATCGCACCGTCGATACCCACATTAAAACGCTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGAGATATAGCCTGAGGGGCCTGTAATGCGTAt                                                                                                     <  2:338401/139‑1 (MQ=255)                                            aTCGCACCGTCGATACCCACATTAAAACGCTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGAGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGttgttg                                                                                     <  2:982319/139‑1 (MQ=255)                                              cGCACCGTCGATACCCACATTAAAACGCTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGAGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGTTGCt                                                                                   >  1:286633/1‑139 (MQ=255)                                                cACCGTCGATACCCACATTAAAACGCTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGAGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGTTGCTgg                                                                                 >  1:844972/1‑139 (MQ=255)                                                cACCGTCGATACCCACATTAAAACGCTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGAGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGTTGCTgg                                                                                 >  1:740193/1‑139 (MQ=255)                                                 aCCGTCGATACCCACATTAAAACGCTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGAGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGTTGCTggg                                                                                >  1:1122594/1‑139 (MQ=255)                                                 aCCGTCGATACCCACATTAAAACGCTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGAGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGTTGCTggg                                                                                <  1:246409/139‑1 (MQ=255)                                                   cGTCGATACCCACATTAAAACGCTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGAGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGTTGCTGGGCt                                                                              <  2:227206/139‑1 (MQ=255)                                                     tCGATACCCACATTAAAACGCTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGAGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTAt                                                                            >  2:1020994/1‑139 (MQ=255)                                                         tACCCACATTAAAACGCTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCCTGAGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGTTGCTGGGATAttttt                                                                        <  2:811371/139‑1 (MQ=255)                                                               caTTAAAACGCTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGAGATATAGCCTGAGGGGc                                                                                                                  >  2:274378/1‑91 (MQ=255)                                                               caTTAAAACGCTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGAGATATAGCCTGAGGGGc                                                                                                                  <  1:274378/91‑1 (MQ=255)                                                                          tGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGAGATATAGCCTGAGGGGc                                                                                                                  <  1:1179803/80‑1 (MQ=255)                                                                          tGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGAGATATAGCCTGAGGGGc                                                                                                                  >  2:1179803/1‑80 (MQ=255)                                                                            cGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGAGAtata                                                                                                                             >  2:446670/1‑67 (MQ=255)                                                                            cGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGAGAtata                                                                                                                             <  1:446670/67‑1 (MQ=255)                                                                            cGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGAGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGc                                                     <  1:1020994/139‑1 (MQ=255)                                                                            cGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGAGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGc                                                     <  2:1185799/139‑1 (MQ=255)                                                                            cGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGAGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGc                                                     >  1:177627/1‑139 (MQ=255)                                                                              tGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGAGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCt                                                   >  1:190676/1‑139 (MQ=255)                                                                                 cAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGAGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGt                                                <  2:1128635/139‑1 (MQ=255)                                                                                 cAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGAGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGt                                                <  1:382225/139‑1 (MQ=255)                                                                                      tGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGAGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTa                                           <  2:190676/139‑1 (MQ=255)                                                                                                       gATCTTTCACCGATTAATACTCATCGCGGCATGAGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGt                                                >  1:869623/1‑117 (MQ=255)                                                                                                       gATCTTTCACCGATTAATACTCATCGCGGCATGAGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGt                                                <  2:869623/117‑1 (MQ=255)                                                                                                              cACCGATTAATACTCATCGCGGCATGAGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGttgttg                                                                                     <  2:722814/73‑1 (MQ=255)                                                                                                              cACCGATTAATACTCATCGCGGCATGAGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGttgttg                                                                                     >  1:722814/1‑73 (MQ=255)                                                                                                              cACCGATTAATACTCATCGCGGCATGAGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGtt                   >  2:772049/1‑139 (MQ=255)                                                                                                                   aTTAATACTCATCGCGGCATGAGATATAGCCTGAggg                                                                                                                    <  1:742803/37‑1 (MQ=255)                                                                                                                   aTTAATACTCATCGCGGCATGAGATATAGCCTGAggg                                                                                                                    >  2:742803/1‑37 (MQ=255)                                                                                                                   aTTAATACTCATCGCGGCATGAGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGc                                                     <  2:248266/100‑1 (MQ=255)                                                                                                                   aTTAATACTCATCGCGGCATGAGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGc                                                     <  2:30759/100‑1 (MQ=255)                                                                                                                   aTTAATACTCATCGCGGCATGAGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGc                                                     >  1:248266/1‑100 (MQ=255)                                                                                                                   aTTAATACTCATCGCGGCATGAGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGc                                                     >  1:30759/1‑100 (MQ=255)                                                                                                                   aTTAATACTCATCGCGGCATGAGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAAcc              >  2:362797/1‑139 (MQ=255)                                                                                                                   aTTAATACTCATCGCGGCATGAGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAAcc              >  2:938885/1‑139 (MQ=255)                                                                                                                            cATCGCGGCATGAGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCg     <  1:362797/139‑1 (MQ=255)                                                                                                                             aTCGCGGCATGAGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGa    >  2:654920/1‑139 (MQ=255)                                                                                                                               cgcgGCATGAGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGTTGCTggg                                                                                <  2:1066109/61‑1 (MQ=255)                                                                                                                               cgcgGCATGAGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGTTGCTggg                                                                                >  1:1066109/1‑61 (MQ=255)                                                                                                                               cgcgGCATGAGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAag  >  2:26859/1‑139 (MQ=255)                                                                                                                                        |                                                                                                                                  CAACTGATGGATAGCGTATGGGAAGATGCGCAGGACACCTACGATCGCACCGTCGATACCCACATTAAAACGCTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAG  >  NC_000913/4636538‑4636802

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