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%	114.8 / NA	37	A12T (GCC→ACC)	yhbQ	GIY‑YIG nuclease superfamily protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (22/15);  total (22/15)

AACATGTTGTTTTGCCAGCCAACGGCTGATGGCGATGAGTGTTTCCATTGCTGTTCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCGCCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGCGAAAGCACTGCGCGGGAAAGGAGAACTAACGCTGGTGTTTTCCGCGCCAGTGGGC  >  NC_000913/3299868‑3300130                                                                                                                                    |                                                                                                                                    aacaTGTTGTTTTGCCAGCCAACGGCTGATGGCGATGAGTGTTTCCATTGCTGTTCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACaa                                                                                                                              >  2:303342/1‑139 (MQ=255)     tgttgtTTTGCCAGCCAACGGCTGATGGCGATGAGTGTTTCCATTGCTGTTCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAg                                                                                                                          <  1:303342/139‑1 (MQ=255)             tGCCAGCCAACGGCTGATGGCGATGAGTGTTTCCATTGCTGTTCTCTTTTATACTGTGGGTCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTtat                                                                                                                    <  1:872022/137‑1 (MQ=255)             tGCCAGCCAACGGCTGATGGCGATGAGTGTTTCCATTGCTGTTCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTtat                                                                                                                    >  2:872022/1‑137 (MQ=255)              gccagccaACGGCTGATGGCGATGAGTGTTTCCATTGCTGTTCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACTGACAATAAGCTTTATAcc                                                                                                                 >  2:147192/1‑139 (MQ=255)                             aTGGCGATGAGTGTTTCCATTGCTGTTCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCggg                                                                                                              >  1:568306/1‑127 (MQ=255)                             aTGGCGATGAGTGTTTCCATTGCTGTTCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCggg                                                                                                              <  2:568306/127‑1 (MQ=255)                                     gAGTGTTTCCATTGCTGTTCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAAcg                                                                                          <  1:274101/139‑1 (MQ=255)                                                           ttttATACTGTGGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCgg                                                                                                               >  1:55926/1‑96 (MQ=255)                                                           ttttATACTGTGGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCgg                                                                                                               <  2:55926/96‑1 (MQ=255)                                                                taCTGTGGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGAGAATAAGCTTTATACAGGGATCACCACAGATGTCGAACGCCGCTATCAGCAGCCCCAAAGCGGCaa                                                               >  1:494190/1‑139 (MQ=255)                                                                     tgGGCCTGACCTTAACATGCGACGCAGACGATGACTCCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCaaa                                                              <  1:1049409/135‑1 (MQ=255)                                                                     tgGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAAcgccg                                                                                       <  1:547064/110‑1 (MQ=255)                                                                     tgGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAAcgccg                                                                                       >  2:547064/1‑110 (MQ=255)                                                                     tgGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCaaa                                                              >  2:1049409/1‑135 (MQ=255)                                                                     tgGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAgggg                                                          >  2:1202615/1‑139 (MQ=255)                                                                       ggCCTGATCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCAACGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGCg                                                        >  1:798371/1‑139 (MQ=255)                                                                       ggCCTGACCTTAACATGCGACGCAGACGCTGACACCCTGGTTTCTTTACTTGATCCGTACCACCTACAATAAGCTTTATAGCGGGATCACCACGGATGTCGAACTCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGCg                                                        >  1:1196312/1‑139 (MQ=255)                                                                       ggCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGCg                                                        >  1:95655/1‑139 (MQ=255)                                                                       ggCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGCg                                                        >  1:611965/1‑139 (MQ=255)                                                                          cTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAAt                                                                                                                             >  1:825507/1‑67 (MQ=255)                                                                          cTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAAt                                                                                                                             <  2:825507/67‑1 (MQ=255)                                                                             aCCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGCGAAAGCa                                                  >  2:686020/1‑139 (MQ=255)                                                                                  aaCATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGCGAAAGCACTgcg                                             <  1:686020/139‑1 (MQ=255)                                                                                          aCGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGCGAAAGCACTGCGCGGGAAAg                                     >  2:78518/1‑139 (MQ=255)                                                                                          aCGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGCGAAAGCACTGCGCGGGAAAg                                     <  2:95655/139‑1 (MQ=255)                                                                                             cAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCAcc                                                                         >  1:866139/1‑100 (MQ=255)                                                                                             cAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCAcc                                                                         <  2:866139/100‑1 (MQ=255)                                                                                                      acacCCTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGCGAAAGCACTGCGCGGGAAAGGAGAACTAACGc                         <  1:78518/139‑1 (MQ=255)                                                                                                          ccTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGa                                                                                             <  1:738142/67‑1 (MQ=255)                                                                                                          ccTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGa                                                                                             >  2:738142/1‑67 (MQ=255)                                                                                                          ccTGGTTTCTTTACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGCGAAAGCACTGCGCGGGAAAGGAGAACTAACGCTGgt                     >  2:256630/1‑139 (MQ=255)                                                                                                                     tACTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGCGAAAGCACTGCGCGGGAAAGGAGAACTAACGCTGGTGTTTTCCgcgc          <  1:621109/139‑1 (MQ=255)                                                                                                                       cTTGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGCGAAAGCACTGCGCGGGAAAGGAGAACTAACGCTGGTGTTTTCCGCGCCa        <  1:940963/139‑1 (MQ=255)                                                                                                                        ttGATCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGCGAAAGCACTGCGCGGGAAAGGAGAACTAACGCTGGTGTTTTCCGCGCCAg       >  1:566113/1‑139 (MQ=255)                                                                                                                           aTCCGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGCGAAAGCACTGCGCGGGAAAGGAGAACTAACGCTGGTGTTTTCCGCGCCAGTgg    >  1:189400/1‑139 (MQ=255)                                                                                                                             ccGTACCACCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGCGAAAGCACTGCGCGGGAAAGGAGAACTAACGCTGGTGTTTTCCGCGCCAGTGGGc  >  1:516610/1‑139 (MQ=255)                                                                                                                                    |                                                                                                                                    AACATGTTGTTTTGCCAGCCAACGGCTGATGGCGATGAGTGTTTCCATTGCTGTTCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCGCCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGCGAAAGCACTGCGCGGGAAAGGAGAACTAACGCTGGTGTTTTCCGCGCCAGTGGGC  >  NC_000913/3299868‑3300130

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