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%	165.3 / NA	52	G179G (GGC→GGT)	glyA	serine hydroxymethyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (27/25);  total (27/25)

CGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGACCGGTAGCATCGATACCGTAAGGAACGATGTTGTACAGTTTACCGG  >  NC_000913/2684844‑2685101                                                                                                                                |                                                                                                                                   cGTGAGCATGAGGAACCGAGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAacatct                                                                                                                         >  2:1101790/1‑137 (MQ=255)         tGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAac                                                                                                                 >  1:1618132/1‑139 (MQ=255)                       cGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGa                                                                                                                               <  1:358422/111‑1 (MQ=255)                       cGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGa                                                                                                                               >  2:358422/1‑111 (MQ=255)                       cGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCAtt                                                                                                   >  1:616347/1‑139 (MQ=255)                             gACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGc                                                                                             >  2:684754/1‑139 (MQ=255)                             gACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGc                                                                                             >  1:891533/1‑139 (MQ=255)                                          aCCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTgg                                                                                <  2:891533/139‑1 (MQ=255)                                               gCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTc                                                                                                <  2:59125/118‑1 (MQ=255)                                               gCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTc                                                                                                >  1:59125/1‑118 (MQ=255)                                               gCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGt                                                                           >  1:663023/1‑139 (MQ=255)                                                       cGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAg                                                                   <  2:1618132/139‑1 (MQ=255)                                                         tGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCAttt                                                                                                  >  2:483313/1‑106 (MQ=255)                                                         tGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCAttt                                                                                                  <  1:483313/106‑1 (MQ=255)                                                         tGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGAt                                                                 <  1:1477274/139‑1 (MQ=255)                                                         tGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGAt                                                                 >  2:1151378/1‑139 (MQ=255)                                                           cgcCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCg                                                               >  1:961521/1‑139 (MQ=255)                                                                 atCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAg                                                         >  2:1160150/1‑139 (MQ=255)                                                                   cAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTc                                                       >  2:337866/1‑139 (MQ=255)                                                                          cAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTc                                                                                                <  1:1338454/91‑1 (MQ=255)                                                                          cAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTc                                                                                                >  2:1338454/1‑91 (MQ=255)                                                                             gTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGc                                                                                             <  2:1028603/91‑1 (MQ=255)                                                                             gTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGc                                                                                             >  1:1028603/1‑91 (MQ=255)                                                                               aaGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGACCGGt                                           <  2:961521/139‑1 (MQ=255)                                                                               aaGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGACCGGt                                           <  2:663023/139‑1 (MQ=255)                                                                               aaGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGACCGGt                                           <  1:1266517/139‑1 (MQ=255)                                                                                  cACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATg                                                 >  1:1041008/1‑130 (MQ=255)                                                                                  cACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAATCCACCACACCGGAATATGCAGAGCAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATg                                                 <  2:1041008/130‑1 (MQ=255)                                                                                       aTGCTGTCAGCGATTTCATGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGACCGCTAGCATCGa                                   <  2:1121588/139‑1 (MQ=255)                                                                                       aTGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCag                                                                                                                       <  1:271495/55‑1 (MQ=255)                                                                                       aTGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCag                                                                                                                       >  2:271495/1‑55 (MQ=255)                                                                                       aTGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCg                                                            >  2:1107137/1‑114 (MQ=255)                                                                                       aTGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCg                                                            <  1:1107137/114‑1 (MQ=255)                                                                                       aTGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGACCGGTAGCATCGa                                   >  2:327773/1‑139 (MQ=255)                                                                                       aTGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGACCGGTAGCATCGa                                   >  2:725226/1‑139 (MQ=255)                                                                                       aTGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGACCGGTAGCATCGa                                   <  1:1160150/139‑1 (MQ=255)                                                                                            gTCAGCGATTTCACGCATTTTCGCCCAGTCCACCATACCGGAATATGCAgaga                                                                                                                    <  1:954690/53‑1 (MQ=255)                                                                                            gTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAgaga                                                                                                                    >  2:954690/1‑53 (MQ=255)                                                                                             tCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGACCGGTAg                                         <  1:92458/127‑1 (MQ=255)                                                                                             tCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGACCGGTAg                                         >  2:92458/1‑127 (MQ=255)                                                                                                 cGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCg                                                            <  2:280868/104‑1 (MQ=255)                                                                                                 cGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCg                                                            >  1:280868/1‑104 (MQ=255)                                                                                                     ttCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGACCGGTAGCATCGATACCGTAAGGAACg                     <  1:337866/139‑1 (MQ=255)                                                                                                       cACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGc                                                                                             <  2:288920/65‑1 (MQ=255)                                                                                                       cACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGc                                                                                             >  1:288920/1‑65 (MQ=255)                                                                                                       cACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGACCGGTAGCATCGATACCGTAAGGAACg                     >  1:714382/1‑137 (MQ=255)                                                                                                       cACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGACCGGTAGCATCGATACCGTAAGGAACg                     <  2:714382/137‑1 (MQ=255)                                                                                                        aCGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGACCGGTAGCATCGATACCGTAAGGAACGAtg                  <  1:327773/139‑1 (MQ=255)                                                                                                             ttttCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGACCGGTAGCATCGATACCGTAAGGAACGATGTTGTa             <  2:616347/139‑1 (MQ=255)                                                                                                                        tCCACCACACCGGAATATGCAGAGAAACCACCGATAATCAttt                                                                                                  <  2:756867/43‑1 (MQ=255)                                                                                                                        tCCACCACACCGGAATATGCAGAGAAACCACCGATAATCAttt                                                                                                  >  1:756867/1‑43 (MQ=255)                                                                                                                        tCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGACCGGTAGCATCGATACCGTAAGGAACGATGTTGTACAGTTTACCgg  >  1:468386/1‑139 (MQ=255)                                                                                                                                |                                                                                                                                   CGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGACCGGTAGCATCGATACCGTAAGGAACGATGTTGTACAGTTTACCGG  >  NC_000913/2684844‑2685101

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