Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,684,971	G→A	G179G (GGC→GGT)	glyA ←	serine hydroxymethyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,684,971	0	G	A	100.0%	201.2 / NA	63	G179G (GGC→GGT)	glyA	serine hydroxymethyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (31/32);  total (31/32)

TAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGACCGGTAGCATCGATACCGTAAGGAACGATGTTGTACAGTTTACCGG  >  NC_000913/2684835‑2685101                                                                                                                                         |                                                                                                                                   tagtaACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACAcc                                                                                                                                  <  1:436278/139‑1 (MQ=255)    taACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGa                                                                                                                               <  2:1569019/139‑1 (MQ=255)          cGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGc                                                                                                                         >  2:1390217/1‑139 (MQ=255)          cGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGc                                                                                                                         >  1:82877/1‑139 (MQ=255)            tGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCag                                                                                                                       >  1:1251822/1‑139 (MQ=255)            tGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCag                                                                                                                       >  2:1433338/1‑139 (MQ=255)              aGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCagag                                                                                                                     >  1:747660/1‑139 (MQ=255)              aGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCagag                                                                                                                     >  2:48642/1‑139 (MQ=255)                  tGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAac                                                                                                                 >  1:778519/1‑139 (MQ=255)                    aGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAacca                                                                                                               <  1:1317012/139‑1 (MQ=255)                         ccGGGTTCGGGTAGACGCCAGCAGCAATCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGAt                                                                                                          >  2:691688/1‑139 (MQ=255)                           gggTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAacc                                                                                                                <  2:1351020/131‑1 (MQ=255)                           gggTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAacc                                                                                                                >  1:1351020/1‑131 (MQ=255)                           gggTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATaa                                                                                                        >  2:1078477/1‑139 (MQ=255)                                 gggTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCAttt                                                                                                  >  2:1454651/1‑139 (MQ=255)                                      gACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGc                                                                                             >  1:649840/1‑139 (MQ=255)                                      gACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGc                                                                                             >  2:406563/1‑139 (MQ=255)                                       aCGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAatat                                                                                                                           >  1:84223/1‑108 (MQ=255)                                       aCGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAatat                                                                                                                           <  2:84223/108‑1 (MQ=255)                                        cGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCtt                                                                                           <  2:649840/139‑1 (MQ=255)                                                      aGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCtt                                                                             <  1:691688/139‑1 (MQ=255)                                                             cAACGTGCGCCATATCAACGAACAGGTAAGCACCTATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAacca                                                                                                               <  2:88455/98‑1 (MQ=255)                                                             cAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAacca                                                                                                               >  1:88455/1‑98 (MQ=255)                                                             cAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAac                                                                                                                 <  1:605317/96‑1 (MQ=255)                                                             cAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAac                                                                                                                 >  2:605317/1‑96 (MQ=255)                                                             cAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTc                                                                      >  1:1159823/1‑139 (MQ=255)                                                                       cATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCg                                                            <  1:406563/139‑1 (MQ=255)                                                                       cATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCg                                                            <  1:1433338/139‑1 (MQ=255)                                                                            cAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTc                                                       >  2:83134/1‑139 (MQ=255)                                                                             aaCGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCg                                                      <  1:48642/139‑1 (MQ=255)                                                                               cGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGat                                                    <  1:1474477/139‑1 (MQ=255)                                                                                   cAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAac                                                                                                                 <  2:997716/74‑1 (MQ=255)                                                                                   cAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAac                                                                                                                 >  1:997716/1‑74 (MQ=255)                                                                                   cAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGa                                                >  1:1471812/1‑139 (MQ=255)                                                                                     ggTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCt                                                                              <  2:94945/107‑1 (MQ=255)                                                                                     ggTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCt                                                                              >  1:94945/1‑107 (MQ=255)                                                                                     ggTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGAcc                                              >  1:1583649/1‑139 (MQ=255)                                                                                      gTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTc                                                                                                <  2:343592/88‑1 (MQ=255)                                                                                      gTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTc                                                                                                >  1:343592/1‑88 (MQ=255)                                                                                        aaGCACCGATGCTGTCAGCGATTTCACGCATTTTCGTCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGACCGGt                                           <  2:778519/139‑1 (MQ=255)                                                                                        aaGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGACCGGt                                           <  2:747660/139‑1 (MQ=255)                                                                                             ccGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGACCGGTAGCAt                                      <  1:1078477/139‑1 (MQ=255)                                                                                             ccGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGACCGGTAGCAt                                      <  2:1251822/139‑1 (MQ=255)                                                                                               gATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGACCGGTAGCATCg                                    <  2:1243324/139‑1 (MQ=255)                                                                                                aTGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCgg                                                                                              >  2:21822/1‑80 (MQ=255)                                                                                                aTGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCgg                                                                                              <  1:21822/80‑1 (MQ=255)                                                                                                aTGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTgg                                                                                >  2:289306/1‑94 (MQ=255)                                                                                                aTGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTgg                                                                                <  1:289306/94‑1 (MQ=255)                                                                                                aTGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGACCGGTAGCATCGa                                   >  2:847099/1‑139 (MQ=255)                                                                                                aTGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGACCGGTAGCATCGa                                   <  2:1159823/139‑1 (MQ=255)                                                                                                  gCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGACCGGTAGCATCGATa                                 <  2:1360487/139‑1 (MQ=255)                                                                                                            aTTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGACCGGTAGCATCGATACCGTAAGGaa                       <  2:1471812/139‑1 (MQ=255)                                                                                                             tttCACGCATTTTCGCCCAGTCCACCACACCGGCATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGACCGGTAGCATCGATACCGTAAGGAAc                      >  1:207009/1‑139 (MQ=255)                                                                                                              ttCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGACCGGTAGCATCGATACTGTAAGGAACg                     <  2:82877/139‑1 (MQ=255)                                                                                                               tCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAAt                                                                                                       <  2:1134544/56‑1 (MQ=255)                                                                                                               tCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAAt                                                                                                       >  1:1134544/1‑56 (MQ=255)                                                                                                                      ttttCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGACCGGTAGCATCGATACCGTAAGGAACGATGTTGTa             <  1:83134/139‑1 (MQ=255)                                                                                                                      ttttCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGACCGGTAGCATCGATACCGTAAGGAACGATGTTGTa             <  2:1583649/139‑1 (MQ=255)                                                                                                                                 tCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCg                                                               >  1:601596/1‑78 (MQ=255)                                                                                                                                 tCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCg                                                               <  2:601596/78‑1 (MQ=255)                                                                                                                                 tCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGACCGGTAGCATCGATACCGTa                            >  1:1321870/1‑113 (MQ=255)                                                                                                                                 tCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGACCGGTAGCATCGATACCGTa                            <  2:1321870/113‑1 (MQ=255)                                                                                                                                 tCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGACCGGTAGCATCGATACCGTAAGGAACGATGTTGTACAGTTTACCgg  >  1:189949/1‑139 (MQ=255)                                                                                                                                         |                                                                                                                                   TAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGACCGGTAGCATCGATACCGTAAGGAACGATGTTGTACAGTTTACCGG  >  NC_000913/2684835‑2685101

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