Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	500,360	A→G	E79G (GAA→GGA)	gsk →	inosine/guanosine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	500,360	0	A	G	100.0%	183.2 / NA	51	E79G (GAA→GGA)	gsk	inosine/guanosine kinase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (29/22);  total (29/22)

ATCGATCAAACGCTGGTCGATATTGAAGCGAAAGTGGATGATGAATTTATTGAGCGTTATGGATTAAGCGCCGGGCATTCACTGGTGATTGAGGATGATGTAGCCGAAGCGCTTTATCAGGAACTAAAACAGAAAAACCTGATTACCCATCAGTTTGCGGGTGGCACCATTGGTAACACCATGCACAACTACTCGGTGCTCGCGGACGACCGTTCGGTGCTGCTGGGCGTCATGTGCAGCAATATTGAAATTGG  >  NC_000913/500239‑500492                                                                                                                          |                                                                                                                                     aTCGATCAAACGCTGGTCGATATTGAAGCGAAAGTGGATGATGAATTTATTGAGCGTTATGGATTAAGCGCCGGGCATTCACTGGTGATTGAGGATGATGTAGCCGAAGCGCTTTATCAGGGACTAAAACAGAAAAAcc                                                                                                                     >  1:673934/1‑139 (MQ=255)          aCGCTGGTCGATATTGAAGCGAAAGTGGATGATGAATTTATTGAGCGTTATGGATTAAGCGCCGGGCATTCACTGGTGATTGAGGATGATGTAGCCGAAGCGCTTTATCAGGGACTAAAACAGAAAAACCTGATTAccc                                                                                                            >  2:6390/1‑139 (MQ=255)                gTCGATATTGAAGCGAAAGTGGATGATGAATTTATTGAGCGTTATGGATTAAGCGCCGGGCATTCACTGGTGATTGAGGATGATGTAGCCGAAGCGCTTTATCAGGGACTaaa                                                                                                                                >  1:170893/1‑113 (MQ=255)                gTCGATATTGAAGCGAAAGTGGATGATGAATTTATTGAGCGTTATGGATTAAGCGCCGGGCATTCACTGGTGATTGAGGATGATGTAGCCGAAGCGCTTTATCAGGGACTaaa                                                                                                                                <  2:170893/113‑1 (MQ=255)                gTCGATATTGAAGCGAAAGTGGATGATGAATTTATTGAGCGTTATGGATTAAGCGCCGGGCATTCACTGGTGATTGAGGATGATGTAGCCGAAGCGCTTTATCAGGGACTAAAACAGAAAAACCTGATTACCCATCAGt                                                                                                      >  1:466781/1‑139 (MQ=255)                gTCGATATTGAAGCGAAAGTGGATGATGAATTTATTGAGCGTTATGGATTAAGCGCCGGGCATTCACTGGTGATTGAGGATGATGTAGCCGAAGCGCTTTATCAGGGACTAAAACAGAAAAACCTGATTACCCATCAGt                                                                                                      >  1:834485/1‑139 (MQ=255)                gTCGATATTGAAGCGAAAGTGGATGATGAATTTATTGAGCGTTATGGATTAAGCGCCGGGCATTCACTGGTGATTGAGGATGATGTAGCCGAAGCGCTTTATCAGGGACTAAAACAGAAAAACCTGATTACCCATCAGt                                                                                                      >  2:429815/1‑139 (MQ=255)                            gCGAAAGTTGATGATGAATTTATTGAGCGTTATGGATTAAGCGCCGGGCATTCACTGGTGATTGAGGATGATGTAGCCGAAGCGCTTTATCAGGGACTAAAACAGAAAAACCTGATTACCCATCAGTTTGCGGGTGGCa                                                                                          <  2:680589/139‑1 (MQ=255)                              gAAAGTGGATGATGAATTTATTGAGCGTTATGGATTAAGCGCCGGGCATTCACTGGTGATTGAGGATGATGTAGCCGAAGCGCTTTATCAGGGACTAAAACAGAAAAACCTGATTACCCATCAGTTTGCGGGTGGCAcc                                                                                        <  1:429815/139‑1 (MQ=255)                              gAAAGTGGATGATGAATTTATTGAGCGTTATGGATTAAGCGCCGGGCATTCACTGGTGATTGAGGATGATGTAGCCGAAGCGCTTTATCAGGGACTAAAACAGAAAAACCTGATTACCCATCAGTTTGCGGGTGGCAcc                                                                                        <  2:466781/139‑1 (MQ=255)                              gAAAGTGGATGATGAATTTATTGAGCGTTATGGATTAAGCGCCGGGCATTCACTGGTGATTGAGGATGATGTAGCCGAAGCGCTTTATCAGGGACTAAAACAGAAAAACCTGATTACCCATCAGTTTGCGGGTGGCAcc                                                                                        <  1:504825/139‑1 (MQ=255)                               aaaGTGGATGATGAATTTATTGAGCGTTATGGATTAAGCGCCGGGCATTCACTGGTGATTGAGGATGATGTAGCCGAAGCGCTTTATCAGGGACTAAAACAGAAAAACCTGATTACCCATCAGTTTGCGGGTGGCACCa                                                                                       >  1:595430/1‑139 (MQ=255)                                                      gCGTTATGGATTAAGCGTCGGGCATTCACTGGTGATTGAGGATGATGTAGCCGAAGCGCTTTATCAGGGACTAAAACAGaa                                                                                                                          >  1:524654/1‑81 (MQ=255)                                                      gCGTTATGGATTAAGCGCCGGGCATTCACTGGTGATTGAGGATGATGTAGCCGAAGCGCTTTATCAGGGACTAAAACAGaa                                                                                                                          <  2:524654/81‑1 (MQ=255)                                                      gCGTTATGGATTAAGCGCCGGGCATTCACTGGTGATTGAGGATGATGTAGCCGAAGCGCTTTATCAGGGACTAAAACAGAAAAACCTGATTACCCa                                                                                                           >  1:6507/1‑96 (MQ=255)                                                      gCGTTATGGATTAAGCGCCGGGCATTCACTGGTGATTGAGGATGATGTAGCCGAAGCGCTTTATCAGGGACTAAAACAGAAAAACCTGATTACCCa                                                                                                           <  2:6507/96‑1 (MQ=255)                                                                    gcgcCGGGCATTCACTGGTGATTGAGGATGATGTAGCCGAAGCGCTTTATCAGGGACTAAAACAGAAAAACCTGATTACCCATCAGTTTGCGGGTGGCACCATTGGTAACACCATGCACAACTACTCGGTGCTCGCGga                                                  >  2:562636/1‑139 (MQ=255)                                                                    gcgcCGGGCATTCACTGGTGATTGAGGATGATGTAGCCGAAGCGCTTTATCAGGGACTAAAACAGAAAAACCTGATTACCCATCAGTTTGCGGGTGGCACCATTGGTAACACCATGCACAACTACTCGGTGCTCGCGga                                                  >  1:38212/1‑139 (MQ=255)                                                                    gcgcCGGGCATTCACTGGTGATTGAGGATGATGTAGCCGAAGCGCTTTATCAGGGACTAAAACAGAAAAACCTGATTACCCATCAGTTTGCGGGTGGCACCATTGGTAACACCATGCACAACTACTCGGTGCTCGCGga                                                  >  2:49147/1‑139 (MQ=255)                                                                          ggCATTCACTGGTGATTGAGGATGATGTAGCCGAAGCGCTTTATCAGGGACTAAAACAGAAAAACCTGATTACCCATCAGTTTGCGGATGGCACCATTGGTAACACCATGCACAACTACTCGGTGCTCGCGGACGACCg                                            >  2:587434/1‑139 (MQ=255)                                                                                    ggTGATTGAGGATGATGTAGCCGAAGCGCTTTATCAGGGACTAAAACAGAAAAACCTGATTACCCAt                                                                                                          >  1:581485/1‑67 (MQ=255)                                                                                    ggTGATTGAGGATGATGTAGCCGAAGCGCTTTATCAGGGACTAAAACAGAAAAACCTGATTACCCAt                                                                                                          <  2:581485/67‑1 (MQ=255)                                                                                     gTGATTGAGGATGATGTAGCCGAAGCGCTTTATCAGGGACTAAAACAGAAAAACCTGATTACCCATCAGTTTGCGGGTGGCACCATTGGTAaca                                                                              >  2:31213/1‑94 (MQ=255)                                                                                     gTGATTGAGGATGATGTAGCCGAAGCGCTTTATCAGGGACTAAAACAGAAAAACCTGATTACCCATCAGTTTGCGGGTGGCACCATTGGTAaca                                                                              <  1:31213/94‑1 (MQ=255)                                                                                        aTTGAGGATGATGTAGCCGAAGCGCTTTATCAGGGACTAAAACAGAAAAACCTGATTACCCATCAGTTTGCGGGTgg                                                                                            >  2:541368/1‑77 (MQ=255)                                                                                        aTTGAGGATGATGTAGCCGAAGCGCTTTATCAGGGACTAAAACAGAAAAACCTGATTACCCATCAGTTTGCGGGTgg                                                                                            <  1:541368/77‑1 (MQ=255)                                                                                               atgatgTAGCCGAAGCGCTTTATCAGGGACTAAAACAGAAAAACCTGATTACCCATCAGTTTGCGGGTg                                                                                             <  1:681636/69‑1 (MQ=255)                                                                                               atgatgTAGCCGAAGCGCTTTATCAGGGACTAAAACAGAAAAACCTGATTACCCATCAGTTTGCGGGTg                                                                                             >  2:681636/1‑69 (MQ=255)                                                                                               atgatgTAGCCGAAGCGCTTTATCAGGGACTAAAACAGAAAAACCTGATTACCCATCAGTTTGCGGGTGGCACCATTGGTAACACCATg                                                                         <  2:142187/89‑1 (MQ=255)                                                                                               atgatgTAGCCGAAGCGCTTTATCAGGGACTAAAACAGAAAAACCTGATTACCCATCAGTTTGCGGGTGGCACCATTGGTAACACCATg                                                                         >  1:142187/1‑89 (MQ=255)                                                                                               atgatgTAGCCGAAGCGCTTTATCAGGGACTAAAACAGAAAAACCTGATTACCCATCAGTTTGCGGGTGGCACCATTGGTAACACCATGcac                                                                      >  2:570743/1‑92 (MQ=255)                                                                                               atgatgTAGCCGAAGCGCTTTATCAGGGACTAAAACAGAAAAACCTGATTACCCATCAGTTTGCGGGTGGCACCATTGGTAACACCATGcac                                                                      <  1:570743/92‑1 (MQ=255)                                                                                                 gatgTAGCCGAAGCGCTTTATCAGGGACTAAAACAGAAAAACCTGATTACCCATCAGTTTGCGGGTGGCACCATTGGTAACACCATGCACAACTACTCGGTGCTCGCGGACGACCGTTCGGTGCTGCTGGGCGTCAtgt                     >  1:824454/1‑139 (MQ=255)                                                                                                  atgTAGCCGAAGCGCTTTATCAGGGACTAAAACAGAAAAACCTGATTACCCATCAGTTTGCGGGTGGCACCATTGGTAACACCATGCACAACTACTCGGTGCTCGCGGACGACCGTTCGGTGCTGCTGGGCGTCAtgtg                    <  2:834485/139‑1 (MQ=255)                                                                                                  atgTAGCCGAAGCGCTTTATCAGGGACTAAAACAGAAAAACCTGATTACCCATCAGTTTGCGGGTGGCACCATTGGTAACACCATGCACAACTACTCGGTGCTCGCGGACGACCGTTCGGTGCTGCTGGGCGTCAtgtg                    >  2:352416/1‑139 (MQ=255)                                                                                                   tgTAGCCGAAGCGCTTTATCAGGGACTAAAACAGAAAAACCTGATTACCCATCAGTTTGCGGGTGGCACCATTGGTAACACCATGCACAACTACTCGGTGCTCGCGgacga                                               <  2:185871/111‑1 (MQ=255)                                                                                                   tgTAGCCGAAGCGCTTTATCAGGGACTAAAACAGAAAAACCTGATTACCCATCAGTTTGCGGGTGGCACCATTGGTAACACCATGCACAACTACTCGGTGCTCGCGgacga                                               >  1:185871/1‑111 (MQ=255)                                                                                                     tAGCCGAAGCGCTTTATCAGGGACTAAAACAGAAAAACCTGATTACCCATCAGTTTGCGGGTGGCACCATTGGTAACACCATGCACAACTACTCGGTGCTCGCGGACGACCGTTCGGTGCTGCTGGGCGTCATGTgcag                 <  1:352416/139‑1 (MQ=255)                                                                                                     tAGCCGAAGCGCTTTATCAGGGACTAAAACAGAAAAACCTGATTACCCATCAGTTTGCGGGTGGCACCATTGGTAACACCATGCACAACTACTCGGTGCTCGCGGACGACCGTTCGGTGCTGCTGGGCGTCATGTgcag                 <  2:313766/139‑1 (MQ=255)                                                                                                         cGAAGCGCTTTATCAGGGACTAAAACAGAAAAACCTGATTACCCATCAGTTTGCGGGTGGCACCATTGGTAACACCATGCACAACTACTCGGTGCTCGCGGACGACCGTTCGGTGCTGCTGGGCGTCATGTGCAGCAat             >  1:318652/1‑139 (MQ=255)                                                                                                            aGCGCTTTATCAGGGACTAAAACAGAAAAACCTGATTACCCATCAGTTTGCGGGTGGCAcc                                                                                        >  2:233519/1‑61 (MQ=255)                                                                                                            aGCGCTTTATCAGGGACTAAAACAGAAAAACCTGATTACCCATCAGTTTGCGGGTGGCAcc                                                                                        <  1:233519/61‑1 (MQ=255)                                                                                                            aGCGCTTTATCAGGGACTAAAACAGAAAAACCTGATTACCCATCAGTTTGCGGGTGGCACCATTGGTAACAcc                                                                            <  1:32426/73‑1 (MQ=255)                                                                                                            aGCGCTTTATCAGGGACTAAAACAGAAAAACCTGATTACCCATCAGTTTGCGGGTGGCACCATTGGTAACAcc                                                                            >  2:32426/1‑73 (MQ=255)                                                                                                            aGCGCTTTATCAGGGACTAAAACAGAAAAACCTGATTACCCATCAGTTTGCGGGTGGCACCATTGGTAACACCATGCACAACTACTCGGTGCTCGCGGACGACCGTTCGGTGCTGCTGGGCGTCATGTGCAGCAATAtt          >  1:717743/1‑139 (MQ=255)                                                                                                             gcgcTTTATCAGGGACTAAAACAGAAAAACCTGATTACCCATCAGTTTGCGGGTGGCACCAt                                                                                      <  1:291225/62‑1 (MQ=255)                                                                                                             gcgcTTTATCAGGGACTAAAACAGAAAAACCTGATTACCCATCAGTTTGCGGGTGGCACCAt                                                                                      >  2:291225/1‑62 (MQ=255)                                                                                                             gcgcTTTATCAGGGACTAAAACAGAAAAACCTGATTACCCATCAGTTTGCGGGTGGCACCATTGGTAACACCATGCACAACTACTCGGTGCTCGCGGACGACCGTTCGGTGCTGCTGGGCGTCATGTGCAGCAATATTg         >  1:617856/1‑139 (MQ=255)                                                                                                             gcgcTTTATCAGGGACTAAAACAGAAAAACCTGATTACCCATCAGTTTGCGGGTGGCACCATTGGTAACACCATGCACAACTACTCGGTGCTCGCGGACGACCGTTCGGTGCTGCTGGGCGTCATGTGCAGCAATATTg         >  2:385658/1‑139 (MQ=255)                                                                                                                    aTCAGGGACTAAAACAGAAAAACCTGATTACCCATCAGTTTGCGGGTGGCACCATTGGTAACACCATGCACAACTACTCGGTGCTCGCGGACGACCGTTCGGTGCTGCTGGGCGTCATGTGCAGCAATATTGAAATTgg  <  2:717743/139‑1 (MQ=255)                                                                                                                    aTCAGGGACTAAAACAGAAAAACCTGATTACCCATCAGTTTGCGGGTGGCACCATTGGTAACACCATGCACAACTACTCGGTGCTCGCGGACGACCGTTCGGTGCTGCTGGGCGTCATGTGCAGCAATATTGAAATTgg  <  2:824454/139‑1 (MQ=255)                                                                                                                          |                                                                                                                                     ATCGATCAAACGCTGGTCGATATTGAAGCGAAAGTGGATGATGAATTTATTGAGCGTTATGGATTAAGCGCCGGGCATTCACTGGTGATTGAGGATGATGTAGCCGAAGCGCTTTATCAGGAACTAAAACAGAAAAACCTGATTACCCATCAGTTTGCGGGTGGCACCATTGGTAACACCATGCACAACTACTCGGTGCTCGCGGACGACCGTTCGGTGCTGCTGGGCGTCATGTGCAGCAATATTGAAATTGG  >  NC_000913/500239‑500492

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