Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,486,157	A→C	E13A (GAA→GCA)	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,157	0	A	C	100.0%	60.4 / NA	24	E13A (GAA→GCA)	crp	DNA‑binding transcriptional dual regulator CRP
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (14/10);  total (14/10)

CCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTA  >  NC_000913/3486026‑3486291                                                                                                                                    |                                                                                                                                       ccccTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGCATTGGTCTTGt                                                                                                                             >  1:93061/1‑143 (MQ=255) ccccTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGCATGGTTCTTGt                                                                                                                             <  1:576649/143‑1 (MQ=255)  cccTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCGTGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGCATGGTTCTTGtc                                                                                                                            >  1:326243/1‑143 (MQ=255)     ttCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTGATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGCATGGTTCTTGTCTCa                                                                                                                         >  2:21804/1‑143 (MQ=255)        ccgggTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGCATGGTTCTTGTCTCATTg                                                                                                                      <  2:326591/140‑1 (MQ=255)         cAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTGATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACACGACTCTCGCATGGTTCTTGTCTCATTGc                                                                                                                     >  2:107211/1‑143 (MQ=255)                        atatTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGCATGGTTCTTGTCTCATTGCCACATTCATAAGTAc                                                                                                      <  1:608932/143‑1 (MQ=255)                           tttCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGATTGGCAAACCGCAAACAGACCCGACTCTCGCATGGTTCTTGTCTCATTGCCACATTCATAAGTagcaa                                                                                                   >  1:117353/1‑139 (MQ=255)                            ttCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAGCCGCGCATGGTGCTTGGCTAACCGCAAACAGACCCGACTCTCGCATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCAt                                                                                                  <  2:250127/143‑1 (MQ=255)                                             aGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCCAACAGACCCGACTCTCGCATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATGCCAGAGCACGCGTAtt                                                                                 >  2:227111/1‑143 (MQ=255)                                                    ttATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGGAAACAGACCCGACTCTCGCATGGTTCTTGTATCATTGACACATTCATAAGTACCCAACCAAGAGCACGGTTATTCACCCgg                                                                          >  1:311481/1‑143 (MQ=255)                                                       tCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAGCCGCAAACAGACCCGACTCTCGCATGGTTCTTGTCTCAGTGGGACATTCGTAAGTACCCATCCAAGAGCACGCTTATTCAGCAggggg                                                                       >  1:327010/1‑141 (MQ=255)                                                            ctctCGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGCATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGaaaaa                                                                  <  1:22340/143‑1 (MQ=255)                                                                 gagaAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACACGACTTTCGCATGGTTCTTGTCTCATTGCCACATTCATAAGTCGCCATCCAAGAGCACGCTTATTCACCCGGGTGAGAAAGGGGa                                                             >  1:160234/1‑143 (MQ=255)                                                                     aaGCTTATAACATCTGATAACCTCGCACGGTGCTTGGCAATCCGCAAACTGACCCGTCGGTCGCGTTGTTCGTCTCTCATTGCCACATTCATAAGTACCCATCCAAGATGTCGTTTATTTACCAGGGTGAAAAAGCGGAAACg                                                         <  2:340745/143‑1 (MQ=255)                                                                      aGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGCATGGTCCTTGTCTCATTGCCACATTCATAAGCACACATCCAAGAGCACGCTTATCCACAAGGGAGAAAAAGCGGAAAGGc                                                        >  1:281240/1‑143 (MQ=255)                                                                                 gagGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGCATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCATCAGGGTGAAAAAGCGGAAACGCTGTACTACATc                                             <  2:344473/143‑1 (MQ=255)                                                                                        aCCGCGCATGGTGCTTGGCAAACCGGAAACAGACCCGACTCTCGCATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAg                                      <  1:210982/143‑1 (MQ=255)                                                                                                      ttGGCAAACCGCAAACAGACCCGACTCTCGCATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTCCATCGTTAAAAGCTCTGTGGCAGTg                        >  1:501825/1‑143 (MQ=255)                                                                                                       tGGCAAACCGCAAACAGACCCGACTCTCGCATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGc                       >  1:422576/1‑143 (MQ=255)                                                                                                           aaaCCGCAAACAGACCCGACTCTCGCATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATACAAGAGAACGCTTATTCAGCAGGGGGAAAAAGCGGAAAAGCTGTAGTACAGCGTTAAAGGATATGTGGAAGTGCTGAt                   >  1:144123/1‑143 (MQ=255)                                                                                                                 cAAACAGACCCGAGTCTCTCCTGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCCTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGa             <  2:93061/143‑1 (MQ=255)                                                                                                                         actcaCTCGCGCATTGTTCTTGTCTCATTGCCGCATTCATTATTACCCATCCAAGAGCTCGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAgg     <  2:351463/139‑1 (MQ=255)                                                                                                                            gACTCTCGCATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTAAAATATTGGTTTCTTGCCTCTGTGGCTGTGCGGATCAAAGACGAAGAGGGTa  >  1:620093/1‑143 (MQ=255)                                                                                                                                    |                                                                                                                                       CCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTA  >  NC_000913/3486026‑3486291

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