Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,938,295	+T	coding (69/993 nt)	rbsR →	DNA‑binding transcriptional dual regulator RbsR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,938,295	1	.	T	100.0%	54.1 / NA	19	N23? (AAT→AAN)	rbsR	DNA‑binding transcriptional dual regulator RbsR
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base T (12/7);  total (12/7)

GGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTCCGTACCGTGGCGTGAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAAT‑AAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCA  >  NC_000913/3938091‑3938512                                                                                                                                                                                                              |                                                                                                                                                                                                                          ggAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTCCGTACCGTGGCGTGAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATTAAAGATc                                                                                                                                                                                                                    >  1:116795/1‑213 (MQ=255)          aCCATTGCCAGAGGCGATTCGTTTTACCCATGCTGCCGCTGCGATTGCCGTAACAAGTAAAACCGCACCACCTTCCGTACCGTGGCGTGAAGAGATCACCGCATTTTTAGACAGGCAGCGGAGCCGCTTGGCTACAATGAAATATGTTGACCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATTAAAGATCGCTTCGTCAGGGAAGCGCTTACCGCCAAAGTTGAAGCGGCGCTTaa                                                                                                                                                                      >  1:370669/1‑250 (MQ=255)                      ggNGATNCNNTNNGCCNATGCTGCCGCGGCGATTGCCGNNACACGTAAANGCGCACAACCTTCCGTACCGTGGCGTGAAGAGATCGACGCATTTTTAGACAGGNAGAGGTGACNNTTNGCTACNANNNANNNTGNTGCCCGCCTGGCGGGCGTTNCTNNCTCAACNGTTTCTCACGTTATCAATTAAAGATc                                                                                                                                                                                                                    <  2:116795/192‑1 (MQ=255)                          aTTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTCCGTACCGTGGCGTGAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGGTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCGCACGTTATCAATTAAAGATCGCTTCGGCAGTGTAGCGATTACCGCCAAAGGTGAAGCGGGTATTAAAGAATTCAATTACgcgc                                                                                                                                                     <  1:181924/251‑1 (MQ=255)                                    cccATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTCCGTACCGTGGCGTGAAGAGATCGACGCCTTTTTAGACAGGGAGGGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATTAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCaacgt                                                                                                                                                                                     >  1:243846/1‑206 (MQ=255)                                          ctgcCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTCCGTACCGTGGCGTGAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATTAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGcgcg                                                                                                                                     >  1:421085/1‑251 (MQ=255)                                                cTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTCCGTACCGTGGCGTGAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATTAAAGATCGCTTCGTCAGTGAAGc                                                                                                                                                                                                    <  1:35335/182‑1 (MQ=255)                                                                acGTAAAGGCGCACAACCTTCCGTACCGTGGCGTGAAGAGATCGACGCCTCTTTAGACAGGCCGGGGTGCCGCTTGGATACACTCAATGAAGTAGCCCCCCAGGCCGGCGTTTCTACCTCACCAGTTTCTCCCGTTCTCAATTAACGATCGCTTCGACAGTGAAGCGATTACCg                                                                                                                                                                                            >  1:248656/1‑174 (MQ=255)                                                                             naaCCTTCCGTACCGTGGCGTGAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGGTTGGCTACAATGAAAGATGTTGACCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATTAAAAATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCCAACTCAATCAAACACATACCATTGGc                                                                                                   >  1:99601/2‑250 (MQ=255)                                                                                               cGTGAAGAGCTCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATTAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCt                                                                                                                               >  1:115692/1‑204 (MQ=255)                                                                                               cGTGAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATTAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCt                                                                                                                               >  1:115700/1‑204 (MQ=255)                                                                                                         tCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATTAAAGATCGCTTCGTCAGtn                                                                                                                                                                                                        <  1:300480/121‑2 (MQ=255)                                                                                                         tCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGGTTCTACCTCAACAGTTTCTCACGTTATCAATTAAAGATCGCTTCGTCAGTg                                                                                                                                                                                                        >  2:300480/1‑121 (MQ=255)                                                                                                           gcngCANTNNTNNACANGCAGAGTTGACGCTTGGCTACNNTGAAAGATGNTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATTAAANATCGCTTCGNNAGNGAAGCNANNNCNNNCANAGTTGAAGCGGCGATTAAANAANNCAATTANGCGCCATCAGCTCTGGCGCGTAGCCt                                                                                                                               <  2:115700/189‑1 (MQ=255)                                                                                                           gANGCANTNNTNNACANGCAGAGGTGACGCTTGGCTACNNTGAAAGATGNTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATTAAANATCGCTTCGNNAGNGAAGCNANNNCNNNCANAGTTGAAGCGGCGATTAAANAANNCAATTANGCGCCATCAGCTCTGGCGCGTAGCCt                                                                                                                               <  2:115692/192‑1 (MQ=255)                                                                                                           gACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATTAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCCTTTCt                                                                    >  1:32039/1‑251 (MQ=255)                                                                                                                                                               cccGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATTAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAg                <  1:34313/251‑1 (MQ=255)                                                                                                                                                                                    ncTCAACAGTTTCTCACGTTATCAATTAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCa  >  1:96261/2‑244 (MQ=255)                                                                                                                                                                                                     gTTATCAATTAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGCTTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCATCACACTCAATCAAACACATACCATTGGCATGTTGCTCACTGCCAGTACCACTCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAgc                                    >  1:221947/1‑193 (MQ=255)                                                                                                                                                                                                         caatt‑AAAGATCGGTTGGTCAGTGACGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCCAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCCTTTCTATTACGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGt                     >  1:167102/5‑203 (MQ=255)                                                                                                                                                                                                              |                                                                                                                                                                                                                          GGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTCCGTACCGTGGCGTGAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAAT‑AAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCA  >  NC_000913/3938091‑3938512

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