Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,715,479	G→A	T19T (ACG→ACA)	eamB →	cysteine and O‑acetylserine exporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,715,479	0	G	A	100.0%	109.1 / NA	30	T19T (ACG→ACA)	eamB	cysteine and O‑acetylserine exporter
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (24/6);  total (24/6)

AAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTTAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGACGCCAGGACCGAACAATATTCTCGCCCTTAGCTCTGCTACGTCGCATGGATTTCGTCAAAGTACCCGCGTGCTGGCAGGGATGAGTCTGGGATTTTTGATTGTGATGTTACTGTGTGCGGGCATTTCATTTTCACTGGCAGTG  >  NC_000913/2715360‑2715620                                                                                                                        |                                                                                                                                              gaaCATTCATTTTTTTAAATGGTCCGTGGCGGGGACTGTCTACCAAAACCGAGGAGATAACAAGTGACACCGACCCTTTTAAGAGCTTTTTGGACTTACACCCTGATTACCGCTATGACACCAGGAACGAACAATATTCTCGACCTTAGct                                                                                                                >  2:590347/2‑151 (MQ=255)   aCATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTTAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGACACCAGGACCGAACAATATTCTCGCCCTTAGctct                                                                                                              >  1:259091/1‑151 (MQ=255)    cattcattTTTTTAAATGTTCCGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTTAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGACACCAGGACCGAACAATATTCTCGCCCTTAGCTCTg                                                                                                             >  1:275817/1‑151 (MQ=255)     attcattTTTTTAAATGTTCCGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTTAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGACACCAGGACCGAACAATATTCTCGCCCTTAGCTCTGc                                                                                                            >  1:165527/1‑151 (MQ=255)         attTTTTTAAATGTTCCGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTTAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGACACCAGGACCGAACAATATTCTCGCCCTTAGCTCTGCTACg                                                                                                        >  2:581956/1‑151 (MQ=255)            tttttAAATGTTCCGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTTAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGACACCAGGACCGAACAATATTCTCGCCCTTAGCTCTGCTACGTCg                                                                                                     >  2:583130/1‑151 (MQ=255)             ttttAAATGTTCCGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTTAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGACACCAGGACCGAACAATATTCTCGCCCTTAGCTCTGCTACGTCGc                                                                                                    >  1:250654/1‑151 (MQ=255)                tAAATGTTCCGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTTAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGACACCAGGACCGAACAATATTCTCGCCCTTAGCTCTGCTACGTCGCATg                                                                                                 >  2:393355/1‑151 (MQ=255)                   aTGTTCCGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTTAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGACACCAGGACCGAACAATATTCTCGCCCTTAGCTCTGCTACGTCGCATGGAt                                                                                              >  1:687549/1‑151 (MQ=255)                    tGTTCCGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTTAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGACACCAGGACCGAACAATATTCTCGCCCTTAGCTCTGCTACGTCGCATGGAtt                                                                                             >  1:189420/1‑151 (MQ=255)                      ttCCGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTTAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGACACCAGGACCGAACAATATTCTCGCCCTTAGCTCTGCTACGTCGCATGGATTTc                                                                                           >  1:636321/1‑151 (MQ=255)                      ttCCGTGTCGGGGACTGTCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTTAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGACACCAGGACCGAACAATATTCTCGCCCTTAGCTCTGCTACGTCGCATGGATTTc                                                                                           >  2:350497/1‑151 (MQ=255)                       tCCGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTTAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGACACCAGGACCGAACAATATTCTCGCCCTTAGCTCTGCTACGTCGCATGGATTTc                                                                                           <  1:278256/150‑1 (MQ=255)                              cGGGTACTGTCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTTAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGACACCAGGACCGAACAATATTCTCGCCCTTAGCTCTGCTACGTCGCATGGATTTCGTCAAAg                                                                                    <  2:248489/150‑1 (MQ=255)                                ggTACTGTCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTTAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGACACCAGGACCGAACAATATTCTCGCCCTTAGCTCTGCTACGTCGCATGGATTTCGTCAAAGTAc                                                                                 >  1:204763/1‑151 (MQ=255)                                                    aggagATAACAAGTGACACCGACCCTTTTAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGACACCAGGACCGAACAATATTCTCGCCCTTAGCTCTGCTACGTCGCATGGATTTCGTCAAAGTACCCGCGTGCTGGCAGGGATg                                                              >  1:582346/1‑150 (MQ=255)                                                           aaCAAGTGACACCGACCCTTTTAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGACACCAGGACCGAACAATATTCTCGCCCTTAGCTCTGCTACGTCGCATGGATTTCGTCAAAGTACCCGCGTGCTGGCAGGGATGAGTCTgg                                                       <  1:501130/150‑1 (MQ=255)                                                             cAAGTGACACCGACCCTTTTAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGACACCAGGACCGAACAATATTCTCGCCCTTAGCTCTGCTACGTCGCATGGATTTCGTCAAAGTACCCGCGTGCTGGCAGGGATGAGTCTGGGAt                                                    >  1:584876/1‑151 (MQ=255)                                                                    cacCGACCCTTTTAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGACACCAGGACCGAACAATATTCTCGCCCTTAGCTCTGCTACGTCGCATGGATTTCGTCAAAGTACCCGCGTGCTGGCAGGGATGAGTCTGGGATTTTTGAt                                             >  2:758410/1‑151 (MQ=255)                                                                     acCGACCCTTTTAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGACACCAGGACCGAACAATATTCTCGCCCTTAGCTCTGCTACGTCGCATGGATTTCGTCAAAGTACCCGCGTGCTGGCAGGGATGAGTCTGGGATTTTTGAtt                                            <  2:582310/151‑1 (MQ=255)                                                                       cGACCCTTTTAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGACACCAGGACCGAACAATATTCTCGCCCTTAGCTCTGCTACGTCGCATGGATTTCGTCAAAGTACCCGCGTGCTGGCAGGGATGAGTCTGGGATTTTTGATtgt                                          >  1:650728/1‑151 (MQ=255)                                                                             ttttATGTGCTGTTTGGACTTACACCCTGATTTCCGCTATGACACCAGGACCGGACAATATTCTCGCCCTTAGCTCTGCTACGTCGCATGGATTTCGTCACAGTACCCGCGTGCTGGCAGTGATGAGTCTTGGAATTTTGATTGAGCTGtt                                    >  2:124681/1‑151 (MQ=255)                                                                             ttttAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGACACCAGGACCGAACAATATTCTCGCCCTTAGCTCTGCTACGTCGCATGGATTTCGTCAAAGTGCCCGCGTGCTGGCAGGGATGAGTCTGGGATTTTTGATTGTGATGtt                                    >  2:630829/1‑151 (MQ=255)                                                                              tttAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGACACCAGGACCGAACAATATTCTCGCCCTTAGCTCTGCTACGTCGCATGGATTTCGTCAAAGTACCCGCGTGCTGGCAGGGATGAGTCTGGGATTTTTGATTGTGATGtt                                    >  2:382733/1‑150 (MQ=255)                                                                                 aaGTGCTTTTTGGACTTACACCCTGATTACCGCTATGACACCAGGACCGAACAATATTCTCGCCCTTAGCTCTGCTACGTCGCATGGATTTCGTCAAAGTACCCGCGTGCTGGCAGGGATGAGTCTGGGATTTTTGATTGTGATGTTACt                                 <  2:259080/150‑1 (MQ=255)                                                                                           tGGACTTACACCCTGATTACCGCTATGACACCAGGACCGAACAATATTCTCGCCCTTAGCTCTGCTACGTCGCATGGATTTCGTCAAAGTACCCGCGTGCTGGCAGGGATGAGTCTGGGATTTTTGATTGTGATGTTACTGTGTGCGGGc                       >  2:685943/1‑150 (MQ=255)                                                                                             gACTTACACCCTGATTACCGCTATGACACCAGGACCGAACAATATTCTCGCCCTTAGCTCTGCTACGTCGCATGGATTTCGTCAAAGTACCCGCGTGCTGGCAGGGATGAGTCTGGGATTTTTGATTGTGATGTTACTGTGTGCGGGcatt                    >  1:272957/1‑151 (MQ=255)                                                                                                              ccGCTATGACACCAGGACCGAACAATATTCTCGCCCTTAGCTCTGCTACGTCGCATGGATTTCGTCAAAGTACCCGCGTGCTGGCAGGGATGAGTCTGGGATTTTTGATTGTGATGTTACTGTGTGCGGGCATTTCATTTTCACTGGCAGt   >  1:596421/1‑151 (MQ=255)                                                                                                               cGCTATGACACCAGGACCGAACAATATTCTCGCCCTTAGCTCTGCTACGTCGCATGGATTTCGTCAAAGTACCCGCGTGCTGGCAGGGATGAGTCTGGGATTTTTGATTGTGATGTTACTGTGTGCGGGCATTTCATTTTCACTGGCAGTg  <  1:685985/151‑1 (MQ=255)                                                                                                                gCTATGACACCAGGACCGAACAATATTCTCGCCCTTAGCTCTGCTACGTCGCATGGATTTCGTCAAAGTACCCGCGTGCTGGCAGGGATGAGTCTGGGATTTTTGATTGTGATGTTACTGTGTGCGGGCATTTCATTTTCACTGGCAGTg  >  1:409384/1‑150 (MQ=255)                                                                                                                        |                                                                                                                                              AAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTTAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGACGCCAGGACCGAACAATATTCTCGCCCTTAGCTCTGCTACGTCGCATGGATTTCGTCAAAGTACCCGCGTGCTGGCAGGGATGAGTCTGGGATTTTTGATTGTGATGTTACTGTGTGCGGGCATTTCATTTTCACTGGCAGTG  >  NC_000913/2715360‑2715620

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