Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,486,205	C→A	T29K (ACG→AAG)	crp →	DNA‑binding transcriptional dual regulator CRP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,486,205	0	C	A	100.0%	114.7 / NA	34	T29K (ACG→AAG)	crp	DNA‑binding transcriptional dual regulator CRP
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (15/19);  total (15/19)

AGTAATACATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGC‑AAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGG  >  NC_000913/3485961‑3486444                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                aGTAATACATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGC‑AAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTAt                                                                                                                                                                                                                                           >  2:246168/1‑251 (MQ=255)        cATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGC‑AAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCAc                                                                                                                                                                                                                                       <  1:132160/248‑1 (MQ=255)        cATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGC‑AAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCAc                                                                                                                                                                                                                                       >  2:132160/1‑248 (MQ=255)             aTGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGC‑AAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTg                                                                                                                                                                                                                                >  2:183968/1‑250 (MQ=255)                          aagCAACGGACGTCACATTACAGTGCAGTACTGTTTATCGCCCCTTCCCAGGAAACGGGACGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGATCTGGAGCAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAAACGC‑AAACAGACTCGACTCTCGGATGGATCTTGTCTCATTGCCACATCAATAAGTCCCCACCCAAGAGCAAGCTTATTCTACAGGGTTAGAAAGCGGagac                                                                                                                                                                                                                   >  2:82688/3‑247 (MQ=255)                          aTGCAAAGGACGTCTCATTACCGTGCAGTACAGTTGTTAGCCCCTGCCCAGGTAGCGTGAAGCATATTTCGGCAATCCATAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTAAAACAGAGGATACCCGCGAATGTTGCTTGGCAAACCGC‑AAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCTATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAc                                                                                                                                                                                                                   <  1:249501/250‑1 (MQ=255)                              aaaGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGC‑AAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGt                                                                                                                                                                                                              >  2:272422/1‑251 (MQ=255)                                       cacaTTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGC‑AAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATc                                                                                                                                                                                                      <  1:51966/250‑1 (MQ=255)                                         caTTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGC‑AAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGtt                                                                                                                                                                                                   <  1:164044/251‑1 (MQ=255)                                             aCCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACAGCGCATGGTGCTTGGCAAACCGC‑AAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTaa                                                                                                                                                                                                 >  1:170198/1‑249 (MQ=255)                                                           gggATAGCCCCTTCCCAGGTATCGGGAAGCATATTTCGGCAATCGAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGC‑AAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTg                                                                                                                                                                                 <  2:286602/249‑1 (MQ=255)                                                             gATAGCCCCTTCCCAGGTAGTGGGAAGAATATTCCGGCAATCCTAAACCAGCGGCGTTATCTGGCTTTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTTCTTGGCAAACCGC‑AAACAGACCCGACTCTCGAATGGTTCTTTTCTCATTGCCACATTCATAAATACCCATCCAAGAGCAAGCTTATTCACCAGTATGAAAAAGCGGAAACTCTGTACTACATCTTTAAAGTCTCTGTGCCAGTGc                                                                                                                                                                                <  2:142775/250‑1 (MQ=255)                                                                 gCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGC‑AAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATc                                                                                                                                                                           <  1:178737/251‑1 (MQ=255)                                                                              tAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGC‑AAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAgg                                                                                                                                                              >  2:273386/1‑251 (MQ=255)                                                                                gCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGC‑AAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGt                                                                                                                                                            <  1:273386/251‑1 (MQ=255)                                                                                gCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGC‑AAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGt                                                                                                                                                            <  1:272422/251‑1 (MQ=255)                                                                                   ggAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGC‑AAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTaaa                                                                                                                                                         <  1:183968/251‑1 (MQ=255)                                                                                     aaGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGC‑AAACAGACCCGACTCTCGAATGGTTCTTGTCTCATGGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACAGCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGATGGTAAAGa                                                                                                                                                       <  2:305021/251‑1 (MQ=255)                                                                                        cATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGC‑AAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAAt                                                                                                                                                    >  2:358233/1‑251 (MQ=255)                                                                                             ttCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGC‑AAACAGACCCGACTCTCGAATGGTTCTTGTCTCATCGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCGCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATc                                                                                                                                                <  1:161058/250‑1 (MQ=255)                                                                                              tCGGCAATCCAGAGACAGCGGCGTGATCTGGCTCTGGAGAAAGCTTATAACATAGGATAACCGCGAAGGGTGCTTGGCAAACCGC‑AAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATcc                                                                                                                                               <  2:274011/250‑1 (MQ=255)                                                                                                                 ggcgTTATATGGCTCTGGAGCAAGCTTATAACAGCGGATAACAGCGCATGGTGCTTGGAAAAACGCTCAAC‑GACACGACTAACGCATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCCTCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTACAGGCTCTGTGGAGGTGATGATAAACGTCGAATAGGGACAAGAAATGATCATCTCCTATCTGactccgggg                                                                                                                           >  2:159295/1‑243 (MQ=255)                                                                                                                                                     gATAACCGCGCATGGTGCTGGGCAAACCGC‑AAACAGTCCCGACGCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGTGAAAGAAATGATCCTCTCCTATCTGAATCAGGGTTATTGTATTGGCGAACTGGGCCTGTTTGAAGAGGGc                                                                                       <  1:232705/251‑1 (MQ=255)                                                                                                                                                     gATAACCGAGCATGGTGCTTGGCAAACCGC‑AAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAATACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATAGTTAAAGGCTCTGTGGCAATGCTGATAAAAGACGAAGAGGGTAAAGAAATAATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAAAGGGc                                                                                       >  2:210287/1‑251 (MQ=255)                                                                                                                                                      ctAACCGCGCCTGGTGCTTGGCAAACCGC‑AAACAGACCCGACTCTCGAATGGCCATTGTCTCCTAGCCACATTCATAAGTACACATCCACGAGCAAGCTTATTCAACAGGGTGCAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTATGGGGCCGTGCTGATCAAAGACGACGAGGGTAAAGAAATGATCCTCTCCTATCTGAATAAGGGTGATTTTCTTGGCGAACTGGCCCTGTTTGAAGAGGGcc                                                                                      >  2:338604/2‑251 (MQ=255)                                                                                                                                                                         ggCAAACCGC‑AAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGt                                                                   <  2:215084/251‑1 (MQ=255)                                                                                                                                                                           cAAACCGC‑AAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTAc                                                                 <  1:291975/251‑1 (MQ=255)                                                                                                                                                                            aaaCCGC‑AAACAGCCCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACg                                                                <  2:226983/251‑1 (MQ=255)                                                                                                                                                                                       cAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATGCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATTGGTACGTGCGAAAAcc                                                      <  2:170198/251‑1 (MQ=255)                                                                                                                                                                                                               tGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGAACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGGGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCg                              >  1:13753/1‑251 (MQ=255)                                                                                                                                                                                                                          aacaTTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGGACAAAAAAtt                   >  2:194598/2‑251 (MQ=255)                                                                                                                                                                                                                                           aTCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGc                                                                                                      >  2:281329/1‑151 (MQ=255)                                                                                                                                                                                                                                           aTCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGc                                                                                                      <  1:281329/151‑1 (MQ=255)                                                                                                                                                                                                                                           aTCCAAGAGCAAGCTTATTCACCAGGGGGAAAAAGCGGAAAAGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAgg  >  2:274559/1‑251 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                AGTAATACATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGC‑AAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGG  >  NC_000913/3485961‑3486444

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