Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	365,306	G→A	R334C (CGC→TGC)	lacZ ←	beta‑galactosidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	365,306	0	G	A	100.0%	27.5 / NA	12	R334C (CGC→TGC)	lacZ	beta‑galactosidase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (7/5);  total (7/5)

CAGGCCGTAGCGGTCGCACAGCGTGTACCACAGCGGATGGTTCGGATAATGCGAACAGCGCACGGCGTTAAAGTTGTTCTGCTTCATCAGCAGGATATCCTGCACCATCGTCTGCTCATCCATGACCTGACCATGCAGAGGATGATGCTCGTGACGGTTAACGCCTCGAATCAGCAACGGCTTGCCGTTCAGCAGCAGCAGACCATTTTCAATCCGCACCTCGCGGAAACCGACATCGCAGGCTTCTGCTTCAATCAGCGTGCCGTCGGCGGTGTGCAGTTCAACCACCGCACGATAGAGATTCGGGATTTCGGCGCTCCACAGTTTCGGGTTTTCGACGTTCAGACGTAGTGTGACGCGATCGGCATAACCACCACGCTCATCGATAATTTCACCGCCGAAAGGCGCGGTGCCGCTGGCGAC  >  NC_000913/365082‑365504                                                                                                                                                                                                                                 |                                                                                                                                                                                                       cAGGCCGTAGCGGTCGCACAGCGTGTACCACAGCGGATGGTTCGGATAATGCGAACAGCGCACGGCGTTAAAGTTGTTCTGCTTCATCAGCAGGATATCCTGCACCAGCGTCTGCTCATCCATGACCTGACCATGCAGAGGATGATGCTCGTGACGGTTAACGCCTCGAATCAGCAACGGCTTGCCGTTCAGCAGCAGCAGACCATTTTCAATCCGCTCCTCGCAGAAACCGACATCGCAGGCTTCTGCtt                                                                                                                                                                              <  1:231146/251‑1 (MQ=255)       gTAGCGGTCGCACAGCGTGTACCACAGCGGATGGTTCGGATAATGCGAACAGCGAAAGGCGTTAAAGTTGTTCTGCTTAATAAGCAGGATATCCTGCACCATCGTCTGCTCATCCATGACCTGACCATGCAGAGGATGATGATCGTGACGGTTAACGCCTCGAATCAGCAACGGCTTGCCGTTCAGCAGCAGCAGACCATTTTCAATCCGCACATCGCAGAAACCGACATCGCAGGCTTCTGCTTCAATc                                                                                                                                                                         >  1:328223/1‑250 (MQ=255)                         gtACCACAGCGGATGGTTCGGATATTGCAAACCGCGCAAGGCGTTACAGTTGTTCTGCTTAATCAGCAGGATCTACTGCACCATCGTCTGATCATCCATGCCCTGAACATGCAGAGGATGATGCTCGTGAAGGTTAACGCCTCGAATCAACAACGGCTTGCCGTTCAGCAGCAGCAGACCATTTTCAATCCGAACCTAGCAGAAACCGACATCGCAGGCTTCTGCTTCACTCAGCGTGCCGTAGGCGgtgt                                                                                                                                                      >  1:342495/1‑251 (MQ=255)                                                                                 ncTTCATCAGCAGGATATCCTGCACCATCGTCTGCTCATCCATGACCTGACCATGCAGAGGATGATGCTCGTGACGGTTAACGCCTCGAATCAGCAACGGCTTGCCGTTCAGCAGCAGCAGACCATTTTCAATCCGCACCTCGCAGAAACCGACATCGCAGGCTTCTGCTTCAATCAGCGTGCCGTCGGCGgt                                                                                                                                                        >  1:92427/2‑193 (MQ=255)                                                                                             ggATATCCTGCACCATCGTCTGCTCATCCATGACCTGACCATGCAGAGGATGATGCTCGTGACGGTTAACGCCTCGAATCAGCAACGGCTTGCCGTTCAGCAGCAGCAGACCATTTTCAATCCGCACCTCGCAGAAACCGACATCGCAGGCTTCTGCTTCAATCAGCGTGCCGTCGGCGGTGTGCAGTTCAACCACCGCACGATAGAGATTCGGGATTTCGGCGCTCCACAGTTTCGGGTTTTCGACGTTc                                                                                  <  1:96111/251‑1 (MQ=255)                                                                                                                                                gatgCTCGTGACGGTTAACGCCTCGAATCAGCAAAGCCTTGCCGTTCAGCAGCAGCAGACCATTTTCAATCCGCACCTCGCAGAAACCGACATCGCAGGCTTCTGCTTCAATCAGCGTGCCGTCGGCGGTGTGCAGTTCAACCACCGCACGATAGAGATTCGGGATTTCGGCGCTCCACCGTTTCGGGTTTTCGACGTTCAGACGTAGTGTGACGCGAACGGCATAACCACCACGCTCATCGATAATTTc                                >  1:295757/1‑250 (MQ=255)                                                                                                                                                         gACGGTTAACGCCTCGAATCAGCAACGGCTTGCCGTTCAGCAGCAGCAGACCATTTTCAATCCGCACCTCGCAGAAACCGACATCGCAGGCTTCTGCTTCAATCAGCGTGCCGTCGGCGGTGTGCAGTTCAACCACCGCACGATAGAGATTCGGGATTTCGGCGCTCCACAGTTTCGGGTTTTCGACGTTCAGACGTAGTGTGACGCGGTCGGCATAACCACCACGCTCATCGATAATTTCAccgcc                          >  1:122227/1‑247 (MQ=255)                                                                                                                                                                     cTCGAATCAGCAACGGCTTGCCGTTCAGCAGCAGCAGACCATTTTCAATCCGCACCTCGCAGAAACCGACATCGCAGGCTTCTGCTTCAATCAGCGTGCCGTCGGCGGtgtg                                                                                                                                                     <  1:115413/112‑1 (MQ=255)                                                                                                                                                                        gAATCAGCAACGGCTTGCCGTTCAGCAGCAGAAGACCATTTTCAATCCGCACCTCGCAGAAACCGACATCGCAGGCTTCTGCTTCAATCAGCGTGACGTCGGCGGTGTGCAGTTCAACCACCGCACGATAGAGATTCGGGATTTCGGCGCTCCACAGTTTCGGGTTTTCGACGTTCAGACGTAGTGTGACGCGATCGGCATAACCACCACGCTCATCGATAATTTCACCGCCGAAAGGCGCGGTGCCGCTg       >  1:321358/1‑251 (MQ=255)                                                                                                                                                                                    gCTTGCCGTTCAGCAGCAGCAGACCATTTTCAATCCGCACCTCGCAGAAACCGACATCGCAGGCTTCTGCTTCAATCAGCGTGCCGTCGGCGGTGTGCAGTTCAaccac                                                                                                                                         <  1:301550/109‑1 (MQ=255)                                                                                                                                                                                     cTTGCCGTTCAGCAGCAGCAGACCATTTTCAATCCGCACCGCGCAGAAACCGACATCGCAGGCTTCTGCTTCAATCAGCGTGCCGTGGGCGGTGGGCGGTGCAACCCCCGCAGGATAGAGATTCGGGATTGCGGCGCTGCACAGGTTCGGGTTGTCGAGGTTCAAGCGTAGTGTGACGCGTTCGGCATACCCCCCACGCTCATCGATAATTTCACCGCCGAAAGGCGCGGTGCCGCTGGCGAc  <  1:166165/243‑1 (MQ=255)                                                                                                                                                                                                                          gccTCGCAGAAACCGACATCGCAGGCNTCTGCTNNAANCAGCGTGCCGTCGGCGGTGNGCNNNTNNNCNCCCGCNCGNNAGCGATTCGNGATTTCGGCGCTCCACAGTTTCGGGTTTTCGACGTTCAGACGTAGTGTGACGCNACCGGCATANNCACAACGCTCATCGAAAATTTNCCCNNCNNANGgcntg                >  2:81137/2‑189 (MQ=255)                                                                                                                                                                                                                                 |                                                                                                                                                                                                       CAGGCCGTAGCGGTCGCACAGCGTGTACCACAGCGGATGGTTCGGATAATGCGAACAGCGCACGGCGTTAAAGTTGTTCTGCTTCATCAGCAGGATATCCTGCACCATCGTCTGCTCATCCATGACCTGACCATGCAGAGGATGATGCTCGTGACGGTTAACGCCTCGAATCAGCAACGGCTTGCCGTTCAGCAGCAGCAGACCATTTTCAATCCGCACCTCGCGGAAACCGACATCGCAGGCTTCTGCTTCAATCAGCGTGCCGTCGGCGGTGTGCAGTTCAACCACCGCACGATAGAGATTCGGGATTTCGGCGCTCCACAGTTTCGGGTTTTCGACGTTCAGACGTAGTGTGACGCGATCGGCATAACCACCACGCTCATCGATAATTTCACCGCCGAAAGGCGCGGTGCCGCTGGCGAC  >  NC_000913/365082‑365504

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