Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,636,119	C→T	P38L (CCG→CTG)	creB →	DNA‑binding transcriptional regulator CreB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,636,119	0	C	T	100.0%	127.2 / NA	38	P38L (CCG→CTG)	creB	DNA‑binding transcriptional regulator CreB
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (20/18);  total (20/18)

GCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCAAGCAGGTACCCGACGTCATGATTCTCGATGTTGGTCTGCCGGATATTAGCGGCTTTGAATTGTGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGC  >  NC_000913/4636000‑4636254                                                                                                                        |                                                                                                                                        gCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCTGGTGCTGGATAAAGCTCGCa                                                                                                                     >  2:134930/1‑140 (MQ=255)    aTTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCTGGTGCTGGATAAAGCTCGCAAGc                                                                                                                  >  2:799811/1‑140 (MQ=255)            aaCGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCTGGTGCTGGATAAAGCTCGCAAGCAGGTAccc                                                                                                          <  2:178931/140‑1 (MQ=255)            aaCGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCTGGTGCTGGATAAAGCTCGCAAGCAGGTAccc                                                                                                          <  1:799811/140‑1 (MQ=255)                 gAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCTGGTGCTGGATAAAGCTCGCAAGCAGGTACCCGACGt                                                                                                     <  1:79026/140‑1 (MQ=255)                  aaaCGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCTGGTGCTGGATAAAGCTCGCAAGCAGGTACCCGACGTc                                                                                                    <  1:566514/140‑1 (MQ=255)                          tGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCTGGTGCTGGATAAAGCTCGCAAGCAGGTACCCGACGTCATGATtct                                                                                            >  2:70968/1‑140 (MQ=255)                                     aGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCTGGTGCTGGATAAAGCTCGCAAGCAGGTACCCGACGTCATGATTCTCGATGTTGGTc                                                                                 <  1:70968/140‑1 (MQ=255)                                      gATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCTGGTGCTGGATAAAGCTCGCAAGCAGGTACCCGACGTCATGATTCTCGATGTTGGt                                                                                  <  2:590831/138‑1 (MQ=255)                                      gATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCTGGTGCTGGATAAAGCTCGCAAGCAGGTACCCGACGTCATGATTCTCGATGTTGGt                                                                                  >  1:590831/1‑138 (MQ=255)                                      gATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCTGGTGCTGGATAAAGCTCGCAAGCAGGTACCCGACGTCATGATTCTCGATGTTGGTCt                                                                                >  2:814969/1‑140 (MQ=255)                                      gATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCTGGTGCTGGATAAAGCTCGCAAGCAGGTACCCGACGTCATGATTCTCGATGTTGGTCt                                                                                >  1:24696/1‑140 (MQ=255)                                      gATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCTGGCGCTGGATAAAGCTCGCAAGCAGGTACCCGACGTCATGATTCTCGATGTTGGTCt                                                                                >  1:599630/1‑140 (MQ=255)                                       aTGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCTGGTGCTGGATAAAGCTCGCAAGCAGGTACCCGACGTCATGATTCTCGATGTTGGTCTg                                                                               <  1:134930/140‑1 (MQ=255)                                           gCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCTGGTGCTGGATAAAGCTCGCAAGCAGGTACCCGACGTCATGATTCTCGATGTTGGTCTGCCgg                                                                           >  1:424792/1‑140 (MQ=255)                                              aGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCTGGTGCTGGATAAAGCTCGCAAGCAGGTACCCGACGCCATCATGCTCGATGTTGGTCTGCCGGata                                                                        >  1:328218/1‑140 (MQ=255)                                                 gATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCTGGTGCTGGATAAAGCTCGCAAGCAGGTACCCGACGTCATGATTCTCGATGTTGGTCTGCCGGATATTa                                                                     >  2:569682/1‑140 (MQ=255)                                                  aTAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCTGGTGCTGGATAAAGCTCGCAAGCAGGTACCCGACGTCATGATTCTCGATGTTGGTCTGCCGGATATTAg                                                                    <  1:814969/140‑1 (MQ=255)                                                       cGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCTGGTGCTGGATAAAGCTCGCAAGCAGGTACCCGACGTCATGATTCTCGATGTTGGTCTGCCGGATATTAGCGGCt                                                               <  1:83975/140‑1 (MQ=255)                                                                ggTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCTGGTGCTGGATAAAGCTCGCAAGCAGGTACCCGACGTCATGATTCTCGATGTTGGTCTGCCGGATATTAGCGGCTTTGAATtgt                                                      >  1:679296/1‑140 (MQ=255)                                                                       aTGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCTGGTGCTGGATAAAGCTCGCAAGCAGGTACCCGACGTCATGATTCTCGATGTTGGTCTGCCGGATATTAGCGGCTTTGAATTGTGCCGCCa                                               >  1:119373/1‑140 (MQ=255)                                                                           tGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCTGGTGCTGGATAAAGCTCGCAAGCAGGTACCCGACGTCATGATTCTCGATGTTGGTCTGCCGGATATTAGCGGCTTTGAATTGTGCCGCCAGTTa                                           >  1:750108/1‑140 (MQ=255)                                                                           tGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCTGGTGCTGGATAAAGCTCGCAAGCAGGTACCCGACGTCATGATTCTCGATGTTGGTCTGCCGGATATTAGCGGCTTTGAATTGTGCCGCCAGTTa                                           >  2:718288/1‑140 (MQ=255)                                                                           tGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCTGGTGCTGGATAAAGCTCGCAAGCAGGTACCCGACGTCATGATTCTCGATGTTGGTCTGCCGGATATTAGCGGCTTTGAATTGTGCCGCCAGTTa                                           >  1:527437/1‑140 (MQ=255)                                                                           tGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCTGGTGCTGGATAAAGCTCGCAAGCAGGTACCCGACGTCATGATTCTCGATGTTGGTCTGCCGGATATTAGCGGCTTTGAATTGTGCCGCCAGTTa                                           >  1:731433/1‑140 (MQ=255)                                                                           tGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCTGGTGCTGGATAAAGCTCGCAAGCAGGTACCCGACGTCATGATTCTCGATGTTGGTCAGCCGGATGTTAGCGACTTTGAGTTGTGCCGCCAGTTa                                           >  1:261886/1‑140 (MQ=255)                                                                                cagGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCTGGTGCTGGATAAAGCTCGCAAGCAGGTACCCGACGTCATGATTCTCGATGTTGGTCTGCCGGATATTAGCGGCTTTGAATTGTGCCGCCAGTTACTGgc                                      >  2:533225/1‑140 (MQ=255)                                                                                cagGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCTGGTGCTGGATAAAGCTCGCAAGCAGGTACCCGACGTCATGATTCTCGATGTTGGTCTGCCGGATATTAGCGGCTTTGAATTGTGCCGCCAGTTACTGgc                                      >  1:834575/1‑140 (MQ=255)                                                                                       gTTTTGCCGTCGAGGTCTTTGACCGAGGCTTGCTGCTGCGGGATAAAGCTCGCAAGCAGGTACCCGACGTCATGATTCTCGATGTTGGTCTGCCGGATATTAGCGGCTTTGAATTGTGCCGCCAGTTACTGGCGCTCCAt                               <  2:328218/140‑1 (MQ=255)                                                                                            gCCGTCGAGGTCTTTGAGCGAGGCTTGCTGGTGCTGGATAAAGCTCGCAAGCAGGTACCCGACGTCATGATTCTCGATGTTGGTCTGCCGGATATTAGCGGCTTTGAATTGTGCCGCCAGTTACTGGCGCTCCATCCGGc                          <  2:527437/140‑1 (MQ=255)                                                                                            gCCGTCGAGGTCTTTGAGCGAGGCTTGCTGGTGCTGGATAAAGCTCGCAAGCAGGTACCCGACGTCATGATTCTCGATGTTGGTCTGCCGGATATTAGCGGCTTTGAATTGTGCCGCCAGTTACTGGCGCTCCATCCGGc                          <  1:533225/140‑1 (MQ=255)                                                                                            gCCGTCGAGGTCTTTGAGCGAGGCTTGCTGGTGCTGGATAAAGCTCGCAAGCAGGTACCCGACGTCATGATTCTCGATGTTGGTCTGCCGGATATTAGCGGCTTTGAATTGTGCCGCCAGTTACTGGCGCTCCATCCGGc                          <  2:679296/140‑1 (MQ=255)                                                                                            gCCGTCGAGGTCTTTGAGCGAGGCTTGCTGGTGCTGGATAAAGCTCGCAAGCAGGTACCCGACGTCATGATTCTCGATGTTGGTCTGCCGGATATTAGCGGCTTTGAATTGTGCCGCCAGTTACTGGCGCTCCATCCGGc                          <  1:718288/140‑1 (MQ=255)                                                                                               gTCGAGGTCTTTGAGCGCGGCTCGCTGGTGCTGGATAAAACTCGCAAGAAGGTACCCGACGTCATGATTCTCGATGTTGGTCTGCCGGATATTAGCGGCTTTGAATTGTGCCGCCAGTTACTGGCGCTCCATCCGGCGtt                       <  2:261886/140‑1 (MQ=255)                                                                                               gTCGAGGTCTTTGAGCGAGGCTTGCTGGTGCTGGATAAAGCTCGCAAGCAGGTACCCGACGTCATGATTCTCGATGTTGGTCTGCCGGATATTAGCGGCTTTGAATTGTGCCGCCAGTTACTGGCGCTCCATCCGGCGtt                       <  2:750108/140‑1 (MQ=255)                                                                                                                aGGCTTGCTGGTGCTGGATAAAGCTCGCAAGCAGGTACCCGACGTCATGATTCTCGATGTTGGTCTGCCGGATATTAGCGGCTTTGAATTGTGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGa      >  2:725366/1‑140 (MQ=255)                                                                                                                    ttGCTGGTGCTGGATAAAGCTCGCAAGCAGGTACCCGACGTCATGATTCTCGATGTTGGTCTGCCGGATATTAGCGGCTTTGAATTGTGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGc  <  2:24696/140‑1 (MQ=255)                                                                                                                    ttGCTGGTGCTGGATAAAGCTCGCAAGCAGGTACCCGACGTCATGATTCTCGATGTTGGTCTGCCGGATATTAGCGGCTTTGAATTGTGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGc  <  2:119373/140‑1 (MQ=255)                                                                                                                        |                                                                                                                                        GCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCAAGCAGGTACCCGACGTCATGATTCTCGATGTTGGTCTGCCGGATATTAGCGGCTTTGAATTGTGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGC  >  NC_000913/4636000‑4636254

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