Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,658,425	T→G	D20E (GAT→GAG)	gadE →	gad regulon transcriptional activator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,658,425	0	T	G	100.0%	54.4 / NA	16	D20E (GAT→GAG)	gadE	gad regulon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (6/10);  total (6/10)

TTACTATATTGAACAACGATTCGGACAAGGATGTAAATAATGAAAAGGATGACATATTCGAAACGATAACGGCTAAGGAGCAAGTTATGATTTTTCTCATGACGAAAGATTCTTTTCTTTTACAGGGCTTTTGGCAGTTGAAAGATAATCACGAAATGATAAAAATCAATTCCCTGTCAGAGATCAAAAAAGTAGGCAATAAACCCTTCAAGGTTATCATTGATACCTATCACAATCATATCCTTGATGAAGAAGCGATTAAATTTCTGGAG  >  NC_000913/3658280‑3658551                                                                                                                                                  |                                                                                                                               ttACTATATTGAACAACGATTCGGACAAGGATGTAAATAATGAAAAGGATGACATATTCGAAACGATAACGGCTAAGGAGCAAGTTATGATTTTTCTCATGACGAAAGATTCTTTTCTTTTACAGGGCTTTTGGCAGTTGAAAGAGAATCa                                                                                                                           >  1:242502/1‑151 (MQ=255)      atatTGAACAACGATTCGGACAAGGATGTAAATAATGAAAAGGATGACATATTCGAAACGATAACGGCTAAGGAGCAAGTTATGATTTTTCTCATGACGAAAGATTCTTTTCTTTTACAGGGCTTTTGGCAGTTGAAAGAGAATCACGaaa                                                                                                                      >  1:28116/1‑151 (MQ=255)                  gATTCGGACAAGGATGTAAATAATGAAAAGGATGACATATTCGAAACGATAACGGCTAAGGAGCAAGTTATGATTTTTCTCATGACGAAAGATTCTTTTCTTTTACAGGGCTTTTGGCAGTTGAAAGAGAATCACGAAATGATAAAAATCa                                                                                                          >  2:99456/1‑151 (MQ=255)                                tGTAAATAATGAAAAGGATGACATATTCGAAACGATAACGGCTAAGGAGCAAGTTATGATTTTTCTCATGACGAAAGATTCTTTTCTTTTACAGGGCTTTTGGCAGTTGAAAGAGAATCACGAAATGATAAAAATCAATTCCCTGTCagag                                                                                            <  1:201178/151‑1 (MQ=255)                                          gAAAAGGATGACATATTCGAAACGATAACGGCTAAGGAGCAAGTTATGATTTTTCTCATGACGAAAGATTCTTTTCTTTTACAGGGCTTTTGGCAGTTGAAAGAGAATCACGAAATGATAAAAATCAATTCCCTGTCAGAGATCAAAAAAg                                                                                  >  1:233557/1‑151 (MQ=255)                                                               aCGATAACGGCTAAGGAGCAAGTTATGATTTTTCTCATGACGAAAGATTCTTTTCTTTTACAGGGCTTTTGGCAGTTGAAAGAGAATCACGAAATGATAAAAATCAATTCCCTGTCAGAGATCAAAAAAGTAGGCAATAAACCCTTCAAgg                                                             <  1:309594/151‑1 (MQ=255)                                                               aCGATAACGGCTAAGGAGCAAGTTATGATTTTTCTCATGACGAAAGATTCTTTTCTTTTACAGGGCTTTTGGCAGTTGAAAGAGAATCACGAAATGATAAAAATCAATTCCCTGTCAGAGATCAAAAAAGTAGGCAATAAACCCTTCAAgg                                                             <  2:76814/151‑1 (MQ=255)                                                                cGATAACGGCTAAGGAGCAAGTTATGATTTTTCTCATGACGAAAGATTCTTTTCTTTTACAGGGCTTTTGGCAGTTGAAAGAGAATCACGAAATGATAAAAATCAATTCCCTGTCAGAGATCAAAAAAGTAGGCAATAAACCCTTCAAGGt                                                            <  2:244245/151‑1 (MQ=255)                                                                        gCTAAGGAGCAAGTTATGATTTTTCTCATGACGAAAGATTCTTTTCTTTTACAGGGCTTTTGGCAGTTGAAAGAGAATCACGAAATGATAAAAATCAATTCCCTGTCAGAGATCAAAAAAGTAGGCAATAAACCCTTCAAGGTTATCATTg                                                    >  2:235067/1‑151 (MQ=255)                                                                              gAGCAAGTTATGATTTTTCTCATGACGAAAGATTCTTTTCTTTTACAGGGCTTTTGGCAGTTGAAAGAGAATCACGAAATGATAAAAATCAATTCCCTGTCAGAGATCAAAAAAGTAGGCAATAAACCCTTCAAGGTTATCATTGATAc                                                <  2:211498/149‑1 (MQ=255)                                                                              gAGCAAGTTATGATTTTTCTCATGACGAAAGATTCTTTTCTTTTACAGGGCTTTTGGCAGTTGAAAGAGAATCACGAAATGATAAAAATCAATTCCCTGTCAGAGATCAAAAAAGTAGGCAATAAACCCTTCAAGGTTATCATTGATAc                                                <  2:254702/149‑1 (MQ=255)                                                                                     ttATGATTTTTCTCATGACGAAAGATTCTTTTCTTTTACAGGGCTTTTGGCAGTTGAAAGAGAATCACGAAATGATAAAAATCAATTCCCTGTCAGAGATCAAAAAAGTAGGCAATAAACCCTTCAAGGTTATCATTGATACCTATcac                                         <  2:218973/149‑1 (MQ=255)                                                                                             tttCTCATGACGAAAGATTCTTTTCTTTTACAGGGCTTTTGGCAGTTGAAAGAGAATCACGAAATGATAAAAATCAATTCCCTGTCAGAGATCAAAAAAGTAGGCAATAAACCCTTCAAGGTTATCATTGATACCTATCACAATCATATcc                               <  1:45373/151‑1 (MQ=255)                                                                                                      aCGAAAGATTCTTTTCTTTTACAGGGCTTTTGGCAGTTGAAAGAGAATCACGAAATGATAAAAATCAATTCCCTGTCAGAGATCAAAAAAGTAGGCAATAAACCCTTCAAGGTTATCATTGATACCTATCACAATCATATCCTTGATgaag                      >  1:28845/1‑151 (MQ=255)                                                                                                         aaaGATTCTTTTCTTTTACAGGGCTTTTGGCAGTTGAAAGAGAATCACGAAATGATAAAAATCAATTCCCTGTCAGAGATCAAAAAAGTAGGCAATAAACCCTTCAAGGTTATCATTGATACCTATCACAATCATATCCTTGATGaagaag                   <  1:283874/151‑1 (MQ=255)                                                                                                                              ggCTTTTGGCAGTTGAAAGAGAATCACGAAATGATAAAAATCAATTCCCTGTCAGAGATCAAAAAAGTAGGCAATAAACCCTTCAAGGTTATCATTGATACCTATCACAATCATATCCTTGATGAAGAAGCGATTAAATTTCTGgag  <  2:111793/147‑1 (MQ=255)                                                                                                                                                  |                                                                                                                               TTACTATATTGAACAACGATTCGGACAAGGATGTAAATAATGAAAAGGATGACATATTCGAAACGATAACGGCTAAGGAGCAAGTTATGATTTTTCTCATGACGAAAGATTCTTTTCTTTTACAGGGCTTTTGGCAGTTGAAAGATAATCACGAAATGATAAAAATCAATTCCCTGTCAGAGATCAAAAAAGTAGGCAATAAACCCTTCAAGGTTATCATTGATACCTATCACAATCATATCCTTGATGAAGAAGCGATTAAATTTCTGGAG  >  NC_000913/3658280‑3658551

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