Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,604,277	A→G	T40A (ACC→GCC)	bglJ →	DNA‑binding transcriptional regulator BglJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,604,277	0	A	G	100.0%	69.0 / NA	20	T40A (ACC→GCC)	bglJ	DNA‑binding transcriptional regulator BglJ
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (10/10);  total (10/10)

CGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTT  >  NC_000913/4604154‑4604410                                                                                                                            |                                                                                                                                      cGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATGCCTATACCAGTCAGGa                                                                                                                       <  1:16686/140‑1 (MQ=255)   cAGAATGGAACCCAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGAATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATGCCTATACCAGTCAGGAGt                                                                                                                     >  2:50157/1‑140 (MQ=255)   cAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATGCCTATACCAGTCAGGAGt                                                                                                                     >  2:219516/1‑140 (MQ=255)          gAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATGCCTATACCAGTCAGGAGTCATTTCa                                                                                                              <  1:132684/140‑1 (MQ=255)                aGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATGCCTATACCAGTCAGGAGTCAt                                                                                                                  <  1:250790/130‑1 (MQ=255)                aGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATGCCTATACCAGTCAGGAGTCAt                                                                                                                  >  2:250790/1‑130 (MQ=255)                            aagaGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGATAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATGCCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGAt                                                                                            >  1:393223/1‑140 (MQ=255)                            aagaGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGATAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATGCCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGAt                                                                                            <  2:393223/140‑1 (MQ=255)                            aagaGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATGCCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGAt                                                                                            >  2:183421/1‑140 (MQ=255)                                 aaaTGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATGCCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGt                                                                                       >  2:285743/1‑140 (MQ=255)                                           cTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATGCCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGt                                                                             >  2:405936/1‑140 (MQ=255)                                              aTAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATGCCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTATGCGGCGGTctg                                                                          >  1:408907/1‑138 (MQ=255)                                                           tCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATGCCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTttttctttttct                                                             <  1:405936/140‑1 (MQ=255)                                                                               gagaGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATGCCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATGAGAAGTGAgcg                                            >  2:164878/1‑137 (MQ=255)                                                                               gagaGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATGCCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAgcg                                            <  1:164878/137‑1 (MQ=255)                                                                                 gagTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATGCCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCagag                                       <  1:142811/140‑1 (MQ=255)                                                                                 gagTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATGCCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCagag                                       <  1:369579/140‑1 (MQ=255)                                                                                       aTTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATGCCTATACCAGTCAGGAGTCAATTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCGTTTTTTTTTTTTTTTCCATGAGAAAGGAGCGCAGAGAggggg                                 >  2:388073/1‑138 (MQ=255)                                                                                            gAAAGTTTGCGGGTAACCCTTATAAGCTCCATGCCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGACGGTCATTTTTTTTTTTTTTTCCCTGAGAAGTGAGCGCAGAGAGGGGTTATCt                            >  2:260355/1‑140 (MQ=255)                                                                                                             ccTTATAAGCTCCATGCCTATACCAGTCAGGAGTCATTTCAGGATGCAATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGc           <  1:50157/140‑1 (MQ=255)                                                                                                                      cTCCATGCCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGtt  <  1:395423/140‑1 (MQ=255)                                                                                                                            |                                                                                                                                      CGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTT  >  NC_000913/4604154‑4604410

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