Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,154,832	G→A	P6S (CCG→TCG)	argE ←	acetylornithine deacetylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,154,832	0	G	A	100.0%	26.0 / NA	11	P6S (CCG→TCG)	argE	acetylornithine deacetylase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (6/5);  total (6/5)

CCGCCAGCAACAAGCCGCCAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGTTGAGCGTAATAAAACCCACCAGCCGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATAAGAAGGTGAATAGCCCCGATGTTGAATACGCTGATTGTGGGTGCCAGCGGCTACGCTGGCGCA  >  NC_000913/4154613‑4155045                                                                                                                                                                                                                            |                                                                                                                                                                                                                      ccGCCAGCAACAAGCCGCCAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCATAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGAGAAATGGCGATAATTTGTTTTTCATTGTTGACACACCTCTg                                                                                                                                                                                        <  1:28435/251‑1 (MQ=255)        caacaaGCCGCCAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGATAATTTGTTTTTCATTGTTGACACACCTCTGGTc                                                                                                                                                                                     >  1:9126/1‑247 (MQ=255)               ccgccAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCAGCAGCAGAGTGAGTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAATGTGTGGCGATCATAGCGCGGTAAATCTCGATAAATGGCGATAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATc                                                                                                                                                                          <  1:173656/251‑1 (MQ=255)                                                                           ggCTGCACTTCCACATTGAAGCCCAAATCTTTGAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGATAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATCTTCATGCAGTATTTATGAATAACAATACACTAACGTTGCGCGTAATAAAACccacc                                                                                                              >  1:89300/1‑251 (MQ=255)                                                                                         aTTGACGCCCAAATCTTTAAACCAGTNCGCCAGNNGANTGATTAAATCTGCATTGCTNTGNNNGNNNGNCTCTTNCGNNGCGCTTATTNAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGATAATTTGTTTNTCATTGTTGNNACACCTCTGGTCATGATAGTANCAANNTNNANGCANTan                                                                                                                                                          >  2:39721/1‑192 (MQ=255)                                                                                                                                gATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGATAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGTTGAGCGTAATAAAACCCACCAGCCGTAAGGTGAATGTTTTACGTTTAGCCTGGCAACCAGACATAAGAAGGTg                                                         <  1:39721/251‑1 (MQ=255)                                                                                                                                                    gATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGATAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCCATATTCATGCAGTATTTATGAATAAAAATACACTAACGTTGAGCGTAATCAAACCCACCAGCCGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATAAGAAGTTGAATCGCCCCGATGTTGact                                      >  1:120297/1‑248 (MQ=255)                                                                                                                                                            gNNTCTTCCNNGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGATAATTTGTTTTTCATTGTTGACACTCCTCTGGTc                                                                                                                                                                                     <  2:9126/99‑1 (MQ=255)                                                                                                                                                                                       gtgGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGATAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGTTGAGCGTAATAAAACCCACCAGCCGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATAAGAAGGTGAATAGCCCCGATGTTGAATACGCTGATTGTGGGTGCCAGCGGCTACGCTGgcgcc  >  1:36047/1‑250 (MQ=255)                                                                                                                                                                                                      cgGTAAATCTCGATAAATGGCGATAATTTGTTTTTCATTGTTGACACCCCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGaa                                                                                                                                                   <  1:2760/91‑1 (MQ=255)                                                                                                                                                                                                      cgGTAAATCTCGATAAATGGCGATAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATNna                                                                                                                                                   >  2:2760/1‑91 (MQ=255)                                                                                                                                                                                                                            |                                                                                                                                                                                                                      CCGCCAGCAACAAGCCGCCAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGTTGAGCGTAATAAAACCCACCAGCCGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATAAGAAGGTGAATAGCCCCGATGTTGAATACGCTGATTGTGGGTGCCAGCGGCTACGCTGGCGCA  >  NC_000913/4154613‑4155045

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