Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,299,999	G→A	A12T (GCC→ACC)	yhbQ →	GIY‑YIG nuclease superfamily protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,299,999	0	G	A	100.0%	69.1 / NA	24	A12T (GCC→ACC)	yhbQ	GIY‑YIG nuclease superfamily protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (13/11);  total (13/11)

TTGCTGTTCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCGCCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGCGAAAGCACTGCGCGGGAAAGGAGAACTAACGCTGGTGTTTTCCGCGCCAGTGGGCGA  >  NC_000913/3299915‑3300132                                                                                     |                                                                                                                                      ttGCTGTTCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATcagc                                                                                 >  2:519518/1‑139 (MQ=255)   gCTGTTCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACa                                                                                                                                 >  1:156470/1‑89 (MQ=255)   gCTGTTCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACa                                                                                                                                 <  2:156470/89‑1 (MQ=255)                 taCTGTGGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCaa                                                                 >  2:569397/1‑139 (MQ=255)                      tgGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCggg                                                                                                                <  2:570987/87‑1 (MQ=255)                      tgGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCggg                                                                                                                >  1:570987/1‑87 (MQ=255)                      tgGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAAcgc                                                                                           >  2:362507/1‑108 (MQ=255)                      tgGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAAcgc                                                                                           <  1:362507/108‑1 (MQ=255)                        ggCCTGAGCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCAGCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCt                                                                                       >  2:548794/1‑110 (MQ=255)                        ggCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACTATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCt                                                                                       <  1:548794/110‑1 (MQ=255)                        ggCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGCg                                                          >  2:729465/1‑139 (MQ=255)                           cTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGCGaaa                                                       >  1:221767/1‑139 (MQ=255)                                 ttAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGCGAAAGCACTg                                                 <  2:221767/139‑1 (MQ=255)                                   aaCATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGCGAAAGCACTgcg                                               >  1:10264/1‑139 (MQ=255)                                         cGACGCAGACGATGACTCCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATcacca                                                                                                         <  2:451398/75‑1 (MQ=255)                                         cGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATcacca                                                                                                         >  1:451398/1‑75 (MQ=255)                                         cGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAAcg                                                                                            <  1:462504/88‑1 (MQ=255)                                         cGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAAcg                                                                                            >  2:462504/1‑88 (MQ=255)                                                      gACCCCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGCGAAAGCACTGCGCGGGAAAGGAGAACTAACg                            <  2:795554/139‑1 (MQ=255)                                                      gACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGCGAAAGCACTGCGCGGGAAAGGAGAACTAACg                            <  2:10264/139‑1 (MQ=255)                                                      gACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGCGAAAGCACTGCGCGGGAAAGGAGAACTAACg                            >  1:756308/1‑139 (MQ=255)                                                             ttgTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGTGATACACATGTAAGTCGAACGCCGATATCAGCAGCACCAAAGCGGCAAAGCGGCGAAAGCACTGCGCGGGAAAGGAGAACTAACGCTGgtgt                     >  2:746802/3‑139 (MQ=255)                                                                   cTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGCGAAAGCACTGCGCGGGAAAGGAGAACTAACGCTGGTGTTTTCcg               <  2:424397/139‑1 (MQ=255)                                                                                gTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGCGAAAGCACTGCGCGGGAAAGGAGAACTAACGCTGGTGTTTTCCGCGCCAGTGGGCGa  <  2:756308/139‑1 (MQ=255)                                                                                     |                                                                                                                                      TTGCTGTTCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCGCCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGCGAAAGCACTGCGCGGGAAAGGAGAACTAACGCTGGTGTTTTCCGCGCCAGTGGGCGA  >  NC_000913/3299915‑3300132

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