Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	71,134	G→A	G250S (GGT→AGT)	araC →	DNA‑binding transcriptional dual regulator AraC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	71,134	0	G	A	100.0%	21.5 / NA	12	G250S (GGT→AGT)	araC	DNA‑binding transcriptional dual regulator AraC
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (8/4);  total (8/4)

CATCAGCGATCACCTGGCAGACAGCAATTTTGATATCGCCAGCGTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCGTCCATTTGACAGGCACAT  >  NC_000913/70941‑71349                                                                                                                                                                                                  |                                                                                                                                                                                                                        cATCAGCGATCACCTGGCAGACAGCAATTTTGATATCGCCAGCGTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCAGTCGCAATGTTGGTTTTGan                                                                                                                                                                                                     <  1:76015/214‑2 (MQ=255)    cAGCGATCACCTTGCAGACAGCAATTTTGATATCGCCAGCGTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCAGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTcg                                                                                                                                                                                   <  1:189439/229‑1 (MQ=255)            aCCTGGCAGACAGCAATTTTGATATCGCCAGCGTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCAGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGc                                                                                                                                                     >  1:213672/1‑251 (MQ=255)                   aGACAGCAATTTTGATATCGCCAGCGTCGCACAGCATGTTTACTTGTCGCCGTCGCGTCTGTCACCTCTTTTACGCCCGAAGTTAGGGATTAGAGTCTTAAGCTGGCGCGAGGACACACGCATTAGTCCGGCGAAGCTGCTTTTGAGCTCTACCCGGATGCATATCGCAACCGTCAGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCa                                                                                                                                                                >  1:243647/1‑233 (MQ=255)                                                           tGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCCGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGCCCAACGCATTAGTCAGGCGAAGCTGCTTTTGTGCACTACCCGGCTGCCTATCGCCACCGNCAGNNNNNAtnnnnnnnnnn                                                                                                                                                                                                       >  2:156235/1‑144 (MQ=255)                                                                                           ccGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGTCCAACGCATTAGTCAGGAGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCAGTCGCAATGTTGGTTTTGACGATAAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATc                                                                         >  1:98687/1‑248 (MQ=255)                                                                                                                          tGGCGCGAGGACCAACGCATTAGTCAGGCGACGCTGCTTTTGAGCACTACACGGATGCCTATCGCCACCGTCAGTCGCAATGTTGGTTTTGACGCTCACCTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGGGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACg                                                             >  1:217644/1‑229 (MQ=255)                                                                                                                                                                          gcccGGATGCCTCTCGCCACCGTCAGTCGCCATGTTGGTTTTGACGATCAACTCTATTTCTcg                                                                                                                                                                                   <  2:9951/62‑1 (MQ=255)                                                                                                                                                                          aCCCGGATGCCTATCGCCACCGTCAGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTcg                                                                                                                                                                                   >  1:9951/1‑63 (MQ=255)                                                                                                                                                                                  gCCTATCGCCACCGTCAGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCGTcaagccgtcaattgtct  >  1:196688/1‑216 (MQ=255)                                                                                                                                                                                   cttATCGCCCCCGTCATTCGCAATGTTGTTTTTGACTATCATCTCTATTTCTCGCGATTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTgaag                                                                                                                         <  1:183054/110‑1 (MQ=255)                                                                                                                                                                                         gCCACCGTCAGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTNGTCNNNNNTgnnnnnn                                                                             >  2:150279/1‑144 (MQ=255)                                                                                                                                                                                                  |                                                                                                                                                                                                                        CATCAGCGATCACCTGGCAGACAGCAATTTTGATATCGCCAGCGTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCGTCCATTTGACAGGCACAT  >  NC_000913/70941‑71349

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