Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	957,879	G→A	G76D (GGT→GAT)	serC →	phosphoserine/phosphohydroxythreonine aminotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	957,879	0	G	A	100.0%	103.4 / NA	31	G76D (GGT→GAT)	serC	phosphoserine/phosphohydroxythreonine aminotransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (13/18);  total (13/18)

GGTGAGGGGAAATGGCTCAAATCTTCAATTTTAGTTCTGGTCCGGCAATGCTACCGGCAGAGGTGCTTAAACAGGCTCAACAGGAACTGCGCGACTGGAACGGTCTTGGTACGTCGGTGATGGAAGTGAGTCACCGTGGCAAAGAGTTCATTCAGGTTGCAGAGGAAGCCGAGAAGGATTTTCGCGATCTTCTTAATGTCCCCTCCAACTACAAGGTATTATTCTGCCATGGCGGTGGTC‑GCGGTCAGTTTGCTGCGGTACCGCTGAATATTCTCGGTGATAAAACCACCGCAGATTATGTTGATGCCGGTTACTGGGCGGCAAGTGCCATTAAAGAAGCGAAAAAATACTGCACGCCTAATGTCTTTGACGCCAAAGTGACTGTTGATGGTCTGCGCGCGGTTAAGCCAATGCGTGAATGGCAACTCTCTGATAATGCTGCTTATATGCATTATT  >  NC_000913/957642‑958097                                                                                                                                                                                                                                              |                                                                                                                                                                                                                            ggtgAGGGGAAATGGCTCAAATCTTCAATTTTAGCTCAGGTCCGGCAATGCTACCGGCAGAGGTGCTTAAACAGGCTCAACAGGAACTGCGCGACTGGAACGGTCTTGGTACGTCGGTGATGGAAGCGAGTCACCGTGGCAAAGAGTTCATTCAGGTTGCAGAGGAAGCCGAGAAGGATTTTCGCGATCTTCTTAATGTCCCCTCCAACTACAAGGTTTTATTCTGCCATGGCGGTGATC‑GCGGTCAGttt                                                                                                                                                                                                               >  1:177411/1‑251 (MQ=255)         gAAATGGCTCAAATCTTCAATTTTAGTTCTGGTCCGGCAATGCTACCGGCAGAGGTGCTTAAACAGGCTCAACAGGAACTGCGCGACTGGAACGGTCTTGGTACGTCGGTGATGGAAGTGAGTCACCGTGGCAAAGAGTTCATTCAGGTTGCAGAGGAAGCCGAGAAGGATTTTCGCGATCTTCTTAATGTCCCCTCCAACTACAAGGTATTATTCTGCCATGGCGGTGATC‑GCGGTCCGTTTGCTGCggg                                                                                                                                                                                                       >  1:115819/1‑250 (MQ=255)           aaTGGCTCAAATCTTCAATTTTAGTTCTGGTCCGGCAATGCTACCGGCAGAGGTGCTTAAACAGGCTCAACAGGAACTGCGCGACTGGAACGGTCTTGGTACGTCGGTGATGGAAGTGAGTCACCGTGGCAAAGAGTTCATTCAGGTTGCAGAGGAAGCCGAGAAGGATTTTCGCGATCTTCTTAATGTCCCCTCCAACTACAAGGTATTATTCTGCCATGGCGGTGATC‑GCGGTCCGTTTGCTGCGGTAc                                                                                                                                                                                                     >  1:87185/1‑251 (MQ=255)                  cAAATCTTCAATTTTAGTTCTGGTCCGGCAATGCTACCGGCAGAGGTGCTTAAACAGGCTCAACAGGAACTGCGCGACTGGAACGGTCTTGGTACGTCGGTGATGGAAGTGAGTCACCGTGGCAAAGAGTTCATTCAGGTTGCAGAGGAAGCCGAGAAGGATTTTCGCGATCTTCTTAATGTCCCCTCCAACTACAAGGTATTATTCTGCCATGGCGGTGATC‑GCGGTCAGTTTGCTGCGGTACCGCTg                                                                                                                                                                                                <  1:10372/249‑1 (MQ=255)                           aaTTTTAGTTCTGGTCCGGCAATGCTACCGGCAGAGGTGCTTAAACAGGCTCAACAGGAACTGCGCGACTGGAACGGTCTTGGTACGTCGGTGATGGAAGTGAGTCACCGTGGCAAAGAGTTCATTCAGGTTGCAGAGGAAGCCGAGAAGGATTTTCGCGATCTTCTTAATGTCCCCTCCAACTACAAGGTATTATTCTGCCATGGCGGTGATC‑GCGGTCAGTTTGCTGCGGTACCGCTGAATATTCTCgg                                                                                                                                                                                     >  2:319673/1‑251 (MQ=255)                              tttAGTTCTGGTCCGGCAATGCTACCGGCAGAGGTGCTTAAACAGGCTCAACAGGAACTGCGCGACTGGAACGGTCTTGGTACGTCGGTGATGGAAGTGAGTCACCGTGGCAAAGAGTTCATTCAGGTTGCAGAGGAAGCCGAGAAGGATTTTCGCGATCTTCTTAATGTCCCCTCCAACTGCAAGGTATTATTCTGCCATGGCGGTGATC‑GCGGTCAGTTTGCTGCGGTACCGCTGAATATTCTCGGTg                                                                                                                                                                                   <  2:123537/250‑1 (MQ=255)                                tAGTTCTGGTCCGGCAATGCTACCGGCAGAGGTGCTTAAACAGGCTCAACAGGAACTGCGCGACTGGAACGGTCTTGGTACGTCGGTGATGGAAGTGAGTCACCGTGGCAAAGAGTTCATTCAGGTTGCAGAGGAAGCCGAGAAGGATTTTCGCGATCTTCTTAATGTCCCCTCCAACTACAAGGTATTATTCTGCCATGGCGGTGATC‑GCGGTCAGTTTGCTGCGGTACCGCTGAATATTCTCGGTGAt                                                                                                                                                                                 <  1:164988/250‑1 (MQ=255)                                       ggTCCGGCAATGCTACCGGCAGAGGTGCTTAAACAGGCTCAACAGGAACTGCGCGACTGGAACGGTCTTGGTACGTCGGTGATGGAAGTGAGTCACCGTGGCAAAGAGTTCATTCAGGTTGCAGAGGAAGCCGAGAAGGATTTTCGCGATCTTCTTAATGTCCCCTCCAACTACAAGGTATTATTCTGCCAAGGCGGTGATC‑GCGGTCAGTTTGCTGCGGTACCGCTGAATATTCTCGGTGATAAAaccac                                                                                                                                                                         >  2:332334/1‑251 (MQ=255)                                        gaccGGCAATGCTACCGGCAGAGGTGCTTAAACAGGCTCAACAGGAACTGCGCGACTGGAACGGTCTTGGTACGTCGGTGATGGAAGTGAGTCACCGTGGCAAAGAGTTCATTCAGGTTGCAGAGGAAGCCGAGAAGGATTTTCGCGATCTTCTGAATGTCCCCTCCAACTACAAGGTATTATTCTGCCATGGCGGTGATC‑GCGGTCAGTGTGCTGCGGGACCGCTGAATATTCTCGGTGATAAAaccacc                                                                                                                                                                        <  2:185754/249‑1 (MQ=255)                                            ggCAATGCTACCGGCAGAGGTGCTTAAACAGGCTCAACAGGAACTGCGCGACTGGAACGGTCTTGGTACGTCGGTGATGGAAGTGAGTCACCGTGGCAAAGAGTTCATTCAGGTTGCAGAGGAAGCCGAGAAGGATTTTCGCGATCTTCTTAATGTCCCCTCCAACTACAAGGTATTATTCTGCCATGGCGGTGATC‑GCGGTCAGTTTGCTGCGGTACCGCTGAATATTCTCGGTGATAAAACCACCGCAg                                                                                                                                                                    <  1:258084/251‑1 (MQ=255)                                                  gcttccGGCAGAGGTGCTTAAACAGGCTCAACAGGAACTGCGCGACTGGAACGGTCTTGGTACGTCGGTGATGGAAGTGAGTCTCCGTGGCAAAGAGTTCATTCAGGTTGCAGAGGAAGCCGAGAAGGATTTTCGCGATCTTCTTAATGTCCCCGCCAACTACAAGGTATTATTCTGCCATGGCGGTGATC‑GCTGTCAGTTTGCTGCTGTACCGCTGAATATTCTCGGTGATAAAACCACCGCAGATTAt                                                                                                                                                               <  1:44184/246‑1 (MQ=255)                                                     aCCGGCAGAGGTGCGTAAACAGGCTCAACAGGAACTGCGCGACTGGAACGGTCTTGGTACGTCGGTGATGGAAGTGAGTCACCGTGGCAAAGAGTTCATTCAGGTTGCAGAGGAAGCCGAGAAGGATTTTCGCGATCTTCTTAATGTCCCCTCCAACTACAAGGTATTATTCTGCCATGGCGGTGATC‑GCGGTCAGTTTGCTGCGGTACCGCTGAATATTCTCGGTGATAAACCCACCGCAGATTATGtt                                                                                                                                                            <  2:201540/250‑1 (MQ=255)                                                          cAGAGGTGCTTAAACAGGCTCAACAGGAACTGCGCGACTGGAACGGTCTTGGTACGTCGGTGATGGAAGTGAGTCACCGTGGCAAAGAGATCATTCAGGTTGCAGAGGAAGCCGAGAAGGATTTTAGCGATCTTCTTAATGTCCCCTCCAACTACAAGGTATTATTCTGCCATGGCGGTGATC‑GCGGTCAGTTTGCTGCGGTACCGCTGAATATTCTCGGTGATAAAACCACCGCAGATTATGTTGATGcc                                                                                                                                                      >  1:252115/1‑251 (MQ=255)                                                           agagGTGCTTAAACAGGCTCAACAGGAACTGCGCGACTGGAACGGTCTTGGTACGTCGGTGATGGAAGTGAGTCACCGTGGCAAAGAGTTCATTCAGGTTGCAGAGGAAGCCGAGAAGGATTTTCGCGATCTTCTTAATGTCCCCTCCAACTACAAGGTATTATTCTGCCATGGCGGTGATC‑GCGGTCAGTTTGCTGCGGTACCGCTGAATATTCTCGGTGATAAAACCACCGCAGATTATGTTGATGCCg                                                                                                                                                     <  1:319673/251‑1 (MQ=255)                                                                              caacaaGAACTGCGCGACTGGAACGGGTTTGGTACGGCGGTGAAGGAAGTTAATCCCCGTGGCAAAGAGTTCATTCAGGTTGCAGAGGAAGCCGAGAAGGATTTTCGCGATCTTCTTAATGTCCCCTCCAACTACAAGGTATTATTCTGCCATGGCGGTGATC‑GCGGTCAGTTTGCTGCGGTACCGCTGAATATTCTCGGTGATAAAACCACCGCAGATTATGTTGATGCCGGTTACTGGGCGGCAAGTGc                                                                                                                                  <  1:24738/251‑1 (MQ=255)                                                                                     aaCTGCGCGACTGGAACGGTCTTGGTACGTCGGTGATGGAAGTGAGTCACCGTGGCAAAGAGTTCATTCAGGTTGCAGAGGAAGCCGAGAAGGATTTTCGCGATCTTCTTAATGTCCCCTCCAACTACAAGGTATTATTCTGCCATGGCGGTGATC‑GCGGTCAGTTTGCTGCGGTACCGCTGAATATTCTCGGTGATAAAACCACCGCAGATTATGTTGATGCCGGTTACTGGGCGGCAAGTGCCATTaa                                                                                                                            >  1:325064/1‑250 (MQ=255)                                                                                                    aCGGTCTTGGTACGTCAGTGATGGAAGTGAGTCACCGTGGCAAAGAGTTCATTCAGGGTGCAGAGGAAGCCGAGAAGGATTTTCGCGATCTTCTTAATGTCCCCTCCAACTACAAGGTATTATTCTGCCATGGCGGTGATC‑GCGGTCAGTTTGCTGCGGTACCGCTGAATATTCTCGGTGATAAAACCACCGCAGATTAGGTTGATGCCGGTTACTGGGCGGCAAGTGCCATTAAAGAAGCGAAAAAATAc                                                                                                            <  2:243687/251‑1 (MQ=255)                                                                                                     cggccTTGGTACGTCGGTGATGGAAGTGAGTCACCGTGGCAAAGAGTTCATTCAGGTTGCAGAGGAAGCCGAGAAGGATTTTCGCGATCTTCTTAATGTCCCCTCCAACTACAAGGTATTATTCTGCCATGGCGGTGATC‑GCGGGCAGTTTGCTGCGGTACCGCTGAATATTCTCGGTGATAAAACCACCGCCGATTATGTTGGTTCCGGTTACTGGGCGGCAAGTGCCATTAAAGAAGCGAAAAAATACt                                                                                                           >  2:123278/5‑251 (MQ=255)                                                                                                                   cGGTGATGGAAGTGAGTCACCGTGGCAAAGAGTTCATTCAGGTTGCAGAGGAAGCCGAGAAGGATTTTCGCGATCTTCTTAATGTCCCCTCCAACGACAAGGTATGATTCTGCCATGGCGGTGATC‑GCGGTCAGTTTGCTGCGGTACCGCGGAATATTCTCGGTGATAAAACCACCGCAGATTATGTTGATGCCGGTGACTGTGCGGCAAGAGCCATTAAAGAAGCGAAAAAATACTGCACGCCTAATGt                                                                                              <  2:230725/250‑1 (MQ=255)                                                                                                                      tGATGGAAGTGAGTCACCGTGGCAAAGAGTTCATTCAGGTTGCAGAGGAAGCCGAGAAGGATTTTCGCGATCTTCTTAATGTCCCCTCCAACTACAAGGTATTATTCTGCCATGGCGGTGATC‑GCGGTCAGGTTGCTGCGGTACCGCTGAATATTCTCGGTGATAAAACCACCGCAGATTATGTTGATGCCGGTTACTGGGCGGCAAGTGCCATTAAAGAAGCGAAAAAATACTGCACGCCTAATGGCttt                                                                                          >  2:290543/1‑251 (MQ=255)                                                                                                                       gATGGAAGTGAGTCACCGTGGCAAAGAGTTCATTCAGGTTGCAGAGGAAGCCGAGAAGGATTTTCGCGATCTTCTTAATGTCCCCTCCAACTACAAGGTATTATTCTGCCATGGCGGTGATC‑GCGGTCAGTTTGCTGCGGTACCGCTGAATATTCTCGGTGATAAAACCACCGCAGATTATGTTGATGCCGGTTACTGGGCGGCAAGTGCCATTAAAGAAGCGAAAAAATACTGCACGCCTAATGTCTTTg                                                                                         <  1:274954/251‑1 (MQ=255)                                                                                                                                             aaaGAGTTCATTCAGGTTGCAGAGGAAGCCGAGAAGGATTTTCGCGATCTTCTTAATGTCCCCTCCAACTACAAGGTATTATTCTGCCATGGCGGTGATC‑GCGGTCAGTTTGCTGCGGTACCGCTGAATATTCTCGGTGATAAAACCACCGCAGATTATGTTGATGCCGGTTACTGGGCGGCAAGTGCCATTAAAGAAGCg                                                                                                                     <  1:248395/201‑1 (MQ=255)                                                                                                                                             aaaGAGTTCATTCAGGTTGCAGAGGAAGCCGAGAAGGATTTTCGCGATCTTCTTAATGTCCCCTCCAACTACAAGGTATTATTCTGCCATGGCGGTGATC‑GCGGTCAGTTTGCTGCGGTACCGCTGAATATTCTCGGTGATAAAACCACCGCAGATTATGTTGATGCCGGTTACTGGGCGGCAAGTGCCATTAAAGAAGCg                                                                                                                     >  2:248395/1‑201 (MQ=255)                                                                                                                                                                         ccGAGAAGGATTTTCGCGATCTTCTTAATGTCCCCTCCAACTACAAGGTATTATTCTGCCATGGCGGTGATC‑GCGGTCAGTTTGCTGCGGTACCGCTGAATATTCTCGGTGATAAAACCACCGCAGATTATGTTGATGCCGGTTACTGGGCGGCAAGTGCCATTAAAGAAGCGAAAAAATACTGCACGCCTAATGTCTTTGACGCCAAAGTGACTGTTGATGGTCTGCGCGCGGTTAAGCCAATGCGTGaa                                       <  2:325064/251‑1 (MQ=255)                                                                                                                                                                         ccGAGAAGGATTTTCGCGATCTTCTTAATGTCCCCTCCAACTACAAGGTATTATTCTGCCATGGCGGTGATC‑GCGGTCAGTTTGCTGCGGTACCGCTGAATATTCTCGGTGATAAAACCACCGCAGATTATGTTGATGCCGGTTACTGGGCGGCAAGTGCCATTAAAGAAGCGAAAAAATACTGCACGCCTAATGTCTTTGACGCAAAAGTGACTGTTGATGGACTGCGCGCGGTTAAGCCAATGCGTGa                                        >  2:374807/1‑250 (MQ=255)                                                                                                                                                                                  aTTTTCGCGATCTTCTTAATGTCCCCTCCAACTACAAGGTATTATTCTGCCATGGCGGTGATC‑GCGGTCCGTTTGCTGCGGTACCGCTGAATATTCTCGGTGATTAAACCACCGCAGATTATGTTGATGCCGGTTACTGGGCGGCAAGTGCCATTAAAGAAGCGAAAAAATACTGCACGCCTAATGTCTTTGACGCCAAAGTGACTGGTGATGGTCTGCGCGCGGTTAAGCCAATGCGTGAATGGCAActc                              >  2:185142/1‑251 (MQ=255)                                                                                                                                                                                    tttCGCGATCTTCTTAATGTCCCCTCCAACTACAAGGTATTATTCTGCCATGGCGGTGATC‑GCGGTCAGTTTGCTGCGGTACCGCTGAATATTCTCGGTGATAAAACCACCGCAGATTATGTTGATGCCGGTTACTGGGCGGCAAGTGCCATTAAAGAAGCGAAAAAATACTGCACGCCTAATGTCTTTGACGCCAAAGTGACTGTTGATGGTCTGCGCGCGGTTAAGCCAATGCGTGAATGGCAActct                             <  2:328152/250‑1 (MQ=255)                                                                                                                                                                                               tcttAATGTCCCCTCCAACTACAAGGTATTATTCTGCCATGGCGGTGATC‑GCGGTCAGTTTGCTGCGGTACCGCTGAATATTCTCGGTGATAAAACCACCGCAGATTATGTTGATGCCGGTTACTGGGCGGCAAGTGCCATTAAAGAAGCGAAAAAATACTGCACGCCTAATGTCTTTGACGCCAAAGTGACTGTTGATGGTCTGCGCGCGGTTAAGCCAATGCGTGAATGGCAACTCTCTGATAAtgct                  <  1:337263/250‑1 (MQ=255)                                                                                                                                                                                               tcttAATGTCCCCTCCAACTACAAGGTATTATTCTGCCATGGCGGTGATC‑GCGGTCAGTTTGCTGCGGTACCGCTGAATATTCTCGGTGATAAAACCACCGCAGATTATGTTGATGCCGGTTACTGGGCGGCAAGTGCCATTAAAGAAGCGAAAAAATACTGCACGCCTAATGTCTTTGACGCCAAAGTGACTGTTGATGGTCTGCGCGCGGTTAAGCCAATGCGTGAATGGCAACTCTCTGATAAtgct                  <  1:201423/250‑1 (MQ=255)                                                                                                                                                                                                             ccAACTACAAGGTATTATTCTGCCATGGCGGTGATC‑GCGGTCAGTTTGCTGCGGTACCGCTGAATATTCTCGGTGATAAAACCACCGCAGATTATGTTGATGCCGGTTACTGGGCGGCAAGTGCCATTAAAGAAGCGAAAAAATACTGCACGCCTAATGTCTTTGACGCCAAAGTGACTGTTGATGGTCTGCGCGCGGTTAAGCCAATGCGTTAATGGCAACTCTCTGATAATGCTGCTTATATGCatta    <  1:364135/250‑1 (MQ=255)                                                                                                                                                                                                              cAACTACAAGGTATTATGCTGCCATGTCGGTGATCCTGGGTCGGTTTTC‑GCGGTACCGCTGCAAATTCACGGTGATCAAACCACCGCAAATTATGTTTATACCGGTTACTGGGCGGCAAGTGCCATTAAAGAAGGGAAAAAATAATGCACGCATAATGTGTTTGGCGCCAAAGAGACTGTTGATGGGCTGCGCGCGGTTAAGACAATGCGTGAATGGCAACTCCCTGTTAATGCTGCTTATATGCattatt  >  2:187123/1‑251 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                            GGTGAGGGGAAATGGCTCAAATCTTCAATTTTAGTTCTGGTCCGGCAATGCTACCGGCAGAGGTGCTTAAACAGGCTCAACAGGAACTGCGCGACTGGAACGGTCTTGGTACGTCGGTGATGGAAGTGAGTCACCGTGGCAAAGAGTTCATTCAGGTTGCAGAGGAAGCCGAGAAGGATTTTCGCGATCTTCTTAATGTCCCCTCCAACTACAAGGTATTATTCTGCCATGGCGGTGGTC‑GCGGTCAGTTTGCTGCGGTACCGCTGAATATTCTCGGTGATAAAACCACCGCAGATTATGTTGATGCCGGTTACTGGGCGGCAAGTGCCATTAAAGAAGCGAAAAAATACTGCACGCCTAATGTCTTTGACGCCAAAGTGACTGTTGATGGTCTGCGCGCGGTTAAGCCAATGCGTGAATGGCAACTCTCTGATAATGCTGCTTATATGCATTATT  >  NC_000913/957642‑958097

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