Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	882,731	C→A	E2* (GAA→TAA)	deoR ←	DNA‑binding transcriptional repressor DeoR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	882,731	0	C	A	100.0%	44.7 / NA	16	E2* (GAA→TAA)	deoR	DNA‑binding transcriptional repressor DeoR
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (7/9);  total (7/9)

TTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATGGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAAAACGACGGGCGCACTGTGGTTGTTCAGATCGCGACGAATCGTCATCTCCGAAACCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCGACGTGT‑TTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGACGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGCACTAAGGCAGCCAATCATACCTGCCAGCAATCCACCGAGCATATCCAGCGGC‑CAGTGGACGCCAAGATAAAC  >  NC_000913/882500‑882971                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                  tttCGTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATGGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAAAACGACGGGCGCACTGTGGTTGTTCAGATCGCGACGAATCGTCATCTCCGAAACCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCGACGTGT‑TTACATAATCCCTCTGAATAg                                                                                                                                                                                                                                 <  1:69120/250‑1 (MQ=255)                    tttttGATCGCTTAACAGGTAATGGCTGGCACTGCGCGGTTCCTGAATAATATAGCCGCCGAGCAAAAAGACGGGCGCACTGTGGTTGTTCAGATCGCGACGAATCGTCATCTCCGAAACCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCGACGTGT‑TTACATAATCCATCTGAATAGTTATTGAAGCGAGCCGCTc                                                                                                                                                                                                              >  2:205164/1‑250 (MQ=255)                                tAACAGGTAATGGCGGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAAAACGACGGGCGCACTGTGGTTGTTCAGATCGCGACGAATCGTCATCTCCGAAACCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTCTCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATTCGCTCTTCGCGACGTGT‑TTACATAATCCCTCTGAATAGTTATTGAAGCGTGCCGCTCAATACTACACtt                                                                                                                                                                                                  <  2:213235/250‑1 (MQ=255)                                    aGGTAATGGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAAAACGACGGGCGCACTGTGGTTGTTCAGATCGCGACGAATCGTCATCTCCGAAACCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCGACGTGT‑TTACATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACttttt                                                                                                                                                                                               >  2:11445/1‑249 (MQ=255)                                            gctggcACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAAAACGACGGGCGCACTGTGGTTGTTCAGATCGCGACGAATCGTCATCTCCGAAACCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCGACGTGT‑TTACATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATc                                                                                                                                                                                     <  1:90987/251‑1 (MQ=255)                                            gctggcACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAAAACGACGGGCGCACTGTGGTTGTGCAGAGCGCGACGAATCGTCATCTCCGAAACCCCAAGCAGGGCGGCGGTGTTTTTAAGATGGAACTTATCGCTGCGTTTTAATTCTGGCAGCAGCTGCCCGATACGCTCTTCGCGACGTGT‑TTACATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAgagagc                                                                                                                                                                                     <  2:259429/251‑3 (MQ=255)                                                                     atatAGCCGCCGAGCAAAACGACGGGCGCACTGTGGTTGTTCAGATCGCGACGAATCGTCATCTCCGAAACCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCGACGTGT‑TTACATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGCGGAt                                                                                                                                                            <  1:11445/251‑1 (MQ=255)                                                                         aGCCGCCGAGCAAAACGACGGGCGCACTGTGGTTGTTCAGATCGAGACGAATCGTCATCTCCGAAAGACCAAGAAGGGCGGCGGCGGCTTTAAGATGGAACTTATCGCTGCGTTTAAATTCTTGCAGAAGCTGAACGATACGCTCTTCGCGACGAGT‑TTACATAATCCAACTGAATAGTTATTGAAGAGAGACGCTCAAAACTACACTTTATAGCAGAGATCAGTCACGCACCCAGCCTTTGCGGATCGGt                                                                                                                                                        >  1:180943/1‑251 (MQ=255)                                                                           ccgccgAGCAAAACGACGGGCGCACTGTGGTTGTTCAGATCGCGCCGACTCGTCATCTCCGAAACCCCAAGCAGGGCGGCCGCGTCTTTAAGATGTCACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCGCCGTGT‑TTACATCATCCCTCTGAATAGTTATTGAAGAGAGCCGCTCAATACTACACTTTTTAGCAGAGATAAGTCACGCACCCAGCCTTTGCGGATCGGTa                                                                                                                                                       >  1:12971/1‑250 (MQ=255)                                                                                gcgCAAAACGACGGGCGCACTGTGGTTGTTCAGATCGCGACGAATCGTCATCTCCGAAACCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCGACGTGT‑TTACATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGCGGATCGGTAATGCa                                                                                                                                                  >  2:67791/3‑250 (MQ=255)                                                                                             ggCGCACTGTGGTTGTTCAGATCGCGACGAATCGTCATCTCCGAAACCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCGACGTGT‑TTACATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATa                                                                                                                                    <  1:224061/251‑1 (MQ=255)                                                                                                    tgtgGTTGTTCAGATCGCGACGAATCGTCATCTCCGAAACCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCGACGTGT‑TTACATAATCCCTCTGACTAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGAt                                                                                                                             >  2:176407/1‑251 (MQ=255)                                                                                                                                                    gggCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCGACGTGT‑TTACATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGACGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTgg                                                                             <  2:186657/251‑1 (MQ=255)                                                                                                                                                                                                                    cctctctTCCCCACGTGT‑TTACATAATCCCTCTGAATAGTTATTCAAGCGGGCCGCGCAATACTACACTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCTCGATTGCAGTCGTTGATAGAGTTTCTGCCCCATCGCTTGCCAGATACTCTGGTCTCTAAGGCGGCCAATCATACCTGCCAGCAATCCGCCGATCATTTCCAGCGGC‑CCTTGGACGc            <  2:180943/248‑1 (MQ=255)                                                                                                                                                                                                                     cgcgctTCGCGACG‑GTGTTACATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATTCCACGATTGCAGACGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGCACTAAGGCAGCCAATCATACCTTCCAGCAATCCACCGAGCATAGCCAGCGGC‑CAGTGGACGcc           <  2:141595/247‑1 (MQ=255)                                                                                                                                                                                                                               gACGGGT‑TTACATAATCACTCTGAATAGTTATTGAAGCAAGCCGCTCAATCCTACACTTTTTATCAGATATCAGTCACGCACCCAGCCTTTGCGGATCTGTAATGCAAAACAGACGCAATACCACGATTGCAGACGATGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGCACTAAGGCATCCACTCAAACCTGCCCTCAATCCACCGAGCATATCCCGCGACACA‑TGGACGCCAAGCTaacc  >  2:248787/1‑248 (MQ=255)                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                  TTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATGGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAAAACGACGGGCGCACTGTGGTTGTTCAGATCGCGACGAATCGTCATCTCCGAAACCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCGACGTGT‑TTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGACGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGCACTAAGGCAGCCAATCATACCTGCCAGCAATCCACCGAGCATATCCAGCGGC‑CAGTGGACGCCAAGATAAAC  >  NC_000913/882500‑882971

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