Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,619,325	G→A	G69E (GGA→GAA)	marR →	DNA‑binding transcriptional repressor MarR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,619,325	0	G	A	100.0%	114.3 / NA	31	G69E (GGA→GAA)	marR	DNA‑binding transcriptional repressor MarR
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (13/18);  total (13/18)

ATCCGCTGCGCGGCGTGTATTACTCCGGTTGAACTGAAAAAGGTATTGTCGGTCGACCTGGGAGCACTGACCCGTATGCTGGATCGCCTGGTCTGTAAAGGCTGGGTGGAAAGG  >  NC_000913/1619264‑1619377                                                              |                                                     aTCCGCTGCGCGGCGTGTATTACTCCGGTTGAACTGAAAAAGGTATTGTCGGTCGACCTGGAAGc                                                   <  1:1276521/65‑1 (MQ=255) aTCCGCTGCGCGGCGTGTATTACTCCGGTTGAACTGAAAAAGGTATTGTCGGTCGACCTGGAAGc                                                   >  1:70182/1‑65 (MQ=255) aTCCGCTGCGCGGCGTGTATTACTCCGGTTGAACTGAAAAAGGTATTGTCGGTCGACCTGGAAGc                                                   <  1:3799324/65‑1 (MQ=255) aTCCGCTGCGCGGCGTGTATTACTCCGGTTGAACTGAAAAAGGTATTGTCGGTCGACCTGGAAGc                                                   <  1:374835/65‑1 (MQ=255) aTCCGCTGCGCGGCGTGTATTACTCCGGTTGAACTGAAAAAGGTATTGTCGGTCGACCTGGAAGc                                                   <  1:1814085/65‑1 (MQ=255) aTCCGCTGCGCGGCGTGTATTACTCCGGTTGAACTGAAAAAGGTATTGTCGGTCGACCTGGAAGc                                                   <  1:2250195/65‑1 (MQ=255) aTCCGCTGCGCGGCGTGTATTACTCCGGTTGAACTGAAAAAGGTATTGTCGGTCGACCTGGAAGc                                                   <  1:2653959/65‑1 (MQ=255)  tCCGCTGCGCGGCGTGTATTACTCCGGTTGAACTGAAAAAGGTATTGTCGGTCGACCTGGAAGCa                                                  <  1:745232/65‑1 (MQ=255)   ccGCTGCGCGGCGTGTATTACTCCGGTTGAACTGAAAAAGGTATTGTCGGTCGACCTGGAAGCAc                                                 >  1:4009660/1‑65 (MQ=255)   ccGCTGCGCGGCGTGTATTACTCCGGTTGAACTGAAAAAGGTATTGTCGATCGACCTGGAAGCAc                                                 >  1:532806/1‑65 (MQ=255)      cTGCGCGGCGTGTATTACTCCGGTTGAACTGAAAAAGGTATTGTCGGTCGACCTGGAAGCACTGa                                              <  1:887784/65‑1 (MQ=255)      cTGCGCGGCGTGTATTACTCCGGTTGAACTGAAAAAGGTATTGTCGGTCGACCTGGAAGCACTGa                                              >  1:1206108/1‑65 (MQ=255)      cTGCGCGGCGTGTATTACTCCGGTTGAACTGAAAAAGGTATTGTCGGTCGACCTGGAAGCACTGa                                              >  1:3741762/1‑65 (MQ=255)      cTGCGCGGCGTGTATTACTCCGGTTGAACTGAAAAAGGTATTGTCGGTCGACCTGGAAGCACTGa                                              <  1:285199/65‑1 (MQ=255)       tGCGCGGCGTGTATTACTCCGGTTGAACTGAAAAAGGTATTGTCGGTCGACCTGGAAGCACTGAc                                             >  1:4312726/1‑65 (MQ=255)        gcgcGGCGTGTATTACTCCGGTTGAACTGAAAAAGGTATTGTCGGTCGACCTGGAAGCACTGAcc                                            <  1:239899/65‑1 (MQ=255)        gcgcGGCGTGTATTACTCCGGTTGAACTGAAAAAGGTATTGTCGGTCGACCTGGAAGCACTGAcc                                            <  1:1596621/65‑1 (MQ=255)        gcgcGGCGTGTATTACTCCGGTTGAACTGAAAAAGGTATTGTCGGTCGACCTGGAAGCACTGAcc                                            <  1:3632558/65‑1 (MQ=255)        gcgcGGCGTGTATTACTCCGGTTGAACTGAAAAAGGTATTGTCGGTCGACCTGGAAGCACTGAcc                                            >  1:1752549/1‑65 (MQ=255)             gcgTGTATTACTCCGGTTGAGCTGAAAAAGGTATTGTCGGTCGACCTGGAAGCACTGACCCGTAt                                       <  1:1290289/65‑1 (MQ=255)             gcgTGTATTACTCCGGTTGAACTGAAAAAGGTATTGTCGGTCGACCTGGAAGCACTGACCCGTAt                                       <  1:2371898/65‑1 (MQ=255)                  tATTACTCCGGTTGAACTGAAAAAGGTATTGTCGGTCGACCTGGAAGCACTGACTCGTATGccgt                                  >  1:1722008/1‑62 (MQ=255)                  tATTACTCCGGTTGAACTGAAAAAGGTATTGTCGGTCGACCTGGAAGCACTGACCCGTATGCTgg                                  <  1:4077347/65‑1 (MQ=255)                    ttACTCCGGTTGAACTGAAAAAGGTATTGTCGGTCGACCTGGAAGCACTGACCCGTATGCCggct                                >  1:1526179/1‑63 (MQ=255)                    ttACTCCGGTTGAACTGAAAAAGGTATTGTCGGGCGACCTGGAAGCACTGACCCGTATGCTGGAt                                >  1:2120826/1‑65 (MQ=255)                      aCTCCGGTTGAACTGAAAAAGGTATTGTCGGTCGACCTGGAAGCACTGACCCGTATCCTGGATCg                              <  1:3883610/65‑1 (MQ=255)                        tCCGGTTGAACTGAAAAAGGTATTGTCGGTCGACCTGGAAGCACTGACCCGTATGCTGGATCGcc                            >  1:1241302/1‑65 (MQ=255)                        tCCGGTTGAACTGAAAAAGGTATTGTCGGTCGACCTGGAAGCACTGACCCGTATGCTGGATCGcc                            >  1:612068/1‑65 (MQ=255)                        tCCGGTTGAACTGAAAAAGGTATTGTCGGTCGACCTGGAAGCACTGACCCGTATGCTGGATCGcc                            >  1:1659999/1‑65 (MQ=255)                                aaCTGAAAAAGGTATTGTCGGTCGACCTGGAAGCACTGACCCGTATGCTGGATCGCCTGGTCTGt                    <  1:4328234/65‑1 (MQ=255)                                         aGGTATTGTCGGTCGACCTGGAAGCACTGACCCGTATGCTGGATCGCCTGGTCTGTAAAGGCTgg           >  1:2445806/1‑65 (MQ=255)                                           gTATTGTCGGTCGACCTGGAAGCACTTGCCCGTATGCTGGGTCGCCTGGTCTGTAAAGGCTGGGt         <  1:766582/65‑1 (MQ=255)                                           gTATTGTCGGTCGACCTGGAAGCACTGACCCGTATGCTGGATCGCCTGGTCTGTAAAGGCTGGGt         <  1:2592460/65‑1 (MQ=255)                                           agaTTGTCGGTCGACCTGGAAGCACTGACCCGTATGCTGGATCGCCTGGTCTGTAAAGGCTGGGt         <  1:2607518/63‑1 (MQ=255)                                                gtcggtcgACCTGGAAGCACTGACCCGTATGCTGGATCGCCTGGTCTGTAAAGGCTGGGTGGaaa    >  1:329376/1‑65 (MQ=255)                                                  cggtcgACCTGGAAGCACTGACCCGTATGCTGGATCGCCTGGTCTGTAAAGGCTGGGTGGAAAgg  >  1:236728/1‑65 (MQ=255)                                                              |                                                     ATCCGCTGCGCGGCGTGTATTACTCCGGTTGAACTGAAAAAGGTATTGTCGGTCGACCTGGGAGCACTGACCCGTATGCTGGATCGCCTGGTCTGTAAAGGCTGGGTGGAAAGG  >  NC_000913/1619264‑1619377

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