Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_bme_pgi	3,486,574	C→A	A152E (GCA→GAA)	crp →	cAMP‑activated global transcription factor, mediator of catabolite repression

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_bme_pgi	3,486,574	0	C	A	100.0%	69.7 / NA	26	A152E (GCA→GAA)	crp	cAMP‑activated global transcription factor, mediator of catabolite repression
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (12/14);  total (12/14)

GCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGA  >  NC_000913_3_bme_pgi/3486430‑3486697                                                                                                                                                 |                                                                                                                            gCCAATTGATTCAGGTACACCCAGACACTCTGACGCGTTAGTCTACACAGAAGCCGCGTCGTCTACAAGTCACCTCAGAGAAAGTGGGCAACCTGCCGTTCCTCGACGTGACCGGCCGCATTGCACAGACTCTGCTGAATCTGGAaaaa                                                                                                                         >  2:290131/1‑149 (MQ=255)             aGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCGCAGACTCTGCTGAATCTGGAAAAACAACCAGACGCt                                                                                                             >  2:148631/1‑149 (MQ=255)             aGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCt                                                                                                             >  2:483143/1‑149 (MQ=255)               gTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTAt                                                                                                           <  2:385565/149‑1 (MQ=255)                        gACGTTCTGATGCGGTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCAccc                                                                                                  <  2:115604/149‑1 (MQ=255)                               tGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGt                                                                                           <  1:173763/149‑1 (MQ=255)                                gATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCATAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTa                                                                                          <  1:501784/149‑1 (MQ=255)                                      tttGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACCGTATGCaaa                                                                                    >  2:461214/1‑149 (MQ=255)                                            tGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGTCTCGCCCGGACGATATGCAAATCaaaa                                                                              >  1:194791/1‑149 (MQ=255)                                                caGATGGGGCCTCGTCTGCAAGTCACTTCCGAGAAGTTGGCCAACCTGGCGTTCGTTGACGTGACGGTCCGCATTGCACAGACTCTGCTTATTCTTTAAAAACAACCAGTCGCTATGACTCTCCCGGACTTTATGCAAATCAAAATTAc                                                                          <  2:121983/149‑1 (MQ=255)                                                           tcgtcTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGCATCTGGAAAAACAACCAGACGCTATGAATCACCCGGACGGTGTGCAAATCAAAATTACCAGTCAGGaaa                                                               >  1:451955/1‑149 (MQ=255)                                                                        cctgcAGAGGAAGGGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACGCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGt                                                  <  2:424661/145‑1 (MQ=255)                                                                          ctcctagaAAGTCGGCAACCGGTCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTAGGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCg                                                <  1:592030/144‑1 (MQ=255)                                                                             agagAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCt                                             >  2:402337/1‑149 (MQ=255)                                                                                   aGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCAGTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTtctc                                       <  2:490990/149‑1 (MQ=255)                                                                                            ccTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCt                                             <  2:263252/134‑1 (MQ=255)                                                                                            ccTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCt                                             >  1:263252/1‑134 (MQ=255)                                                                                             cTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGt                             <  1:334003/149‑1 (MQ=255)                                                                                                   ttCCTCGACGTGAAGGGCCGCAGTGCACAGACTCTGCTGCATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCCGTTCTCGTGAAACCGTGGGagg                       >  2:554535/1‑147 (MQ=255)                                                                                                           cGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGCCGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGa               >  1:143391/1‑149 (MQ=255)                                                                                                           cGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGa               >  1:76599/1‑149 (MQ=255)                                                                                                            gTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGGCTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGGAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCGTTCTGaa              >  1:336766/1‑149 (MQ=255)                                                                                                            gTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGAGGCATTCTGaa              <  1:620916/149‑1 (MQ=255)                                                                                                                 gggCAGCATTGCACAGACTTTTTTGAATCGTTAAAAACAACCAGACGCTATGGCTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGTGACGCATTCTGAAGATGc         <  2:451955/149‑1 (MQ=255)                                                                                                                    ccGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTgg      <  1:483143/149‑1 (MQ=255)                                                                                                                        aTTCCACAGACTCTGCTGAATCTGGAAAAACCCCCCAACAAATCCATGGAGCCTAACAGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTTGGACGCATTCTGAAGATGCTGGAAGa  <  1:211432/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                            GCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGA  >  NC_000913_3_bme_pgi/3486430‑3486697

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