Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,684,867	G→T	P214Q (CCG→CAG)	glyA ←	serine hydroxymethyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,684,867	0	G	T	100.0%	100.2 / NA	29	P214Q (CCG→CAG)	glyA	serine hydroxymethyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (16/13);  total (16/13)

ACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAATATGCAGAGA  >  NC_000913/2684731‑2684987                                                                                                                                         |                                                                                                                         accaGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCTGGt                                                                                                                       >  2:460812/1‑140 (MQ=255)                 gagTTCAGTTTTTTGTACAGCTCTTCGCTACCCCCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCTGGTAGACGCCAGCAGCAAc                                                                                                       <  2:162460/140‑1 (MQ=255)                         ttttttGTACAGCTCTTCGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCTGGTAGACGCCAGCAGCAACCAGGCCCg                                                                                               <  1:53569/140‑1 (MQ=255)                            tttGTACAGCTCTTCGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCTGGTAGACGCCAGCAGCAACCAGGCCCGCaa                                                                                            >  2:10004/1‑140 (MQ=255)                              tGTACAGCTCTTCGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCTGGTAGACGCCAGCAGCAACCAGGCCCGCAACg                                                                                          >  2:282355/1‑140 (MQ=255)                                     ctctTGGCTGCCTCTTTTCGCCGGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCTGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCa                                                                                   <  2:54372/140‑1 (MQ=255)                                     ctctTCGCTCCCCCCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCTGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCa                                                                                   <  1:198469/140‑1 (MQ=255)                                     ctctTCGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCTGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCa                                                                                   >  2:103512/1‑140 (MQ=255)                                        ttCGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCTGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCatat                                                                                <  1:103512/140‑1 (MQ=255)                                          cGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCTGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCa                                                                              >  2:143929/1‑140 (MQ=255)                                               ccaccTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCTGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGaa                                                                         <  1:282355/140‑1 (MQ=255)                                               ccaccTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCTGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGaa                                                                         <  1:460812/140‑1 (MQ=255)                                                  ccTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCTGGTAGACGCCAGTAGCAACCAGGCCCGCAACGTGCGCCATATCAACGaa                                                                         <  2:62489/137‑1 (MQ=255)                                                  ccTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCTGGTAGACGCCAGTAGCAACCAGGCCCGCAACGTGCGCCATATCAACGaa                                                                         >  1:62489/1‑137 (MQ=255)                                                                        cgcgcgGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCTGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCg                                                >  2:611092/1‑140 (MQ=255)                                                                        cgcgcgGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCTGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCg                                                >  2:237435/1‑140 (MQ=255)                                                                                                  gtggtgGTAGTAACAACGTGAGCATGAGGAACCGGGTTCTGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTccac                      <  1:611092/140‑1 (MQ=255)                                                                                                  gtggtgGTAGTAACAACGTGAGCATGAGGAACCGGGTTCTGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTccac                      <  2:534299/140‑1 (MQ=255)                                                                                                  gtggtgGTAGTAACAACGTGAGCATGAGGAACCGGGTTCTGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTccac                      >  1:534299/1‑140 (MQ=255)                                                                                                   tggtggTAGTAACAACGTGAGCATGAGGAACCGGGTTCTGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTccacc                     <  1:502417/140‑1 (MQ=255)                                                                                                    ggtggtAGTAACAACGTGAGCATGAGGAACCGGGTTCTGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTccacca                    <  1:256319/140‑1 (MQ=255)                                                                                                    ggtgggAGTAACAACGTGAGCATGAGGAACCGGCTTCTGGTAGACCCCAGCAGCAACCAGGCCCGCAACGCGCGCCATATCAACGAACCACGAAGCACCGAAGCTGTCACCGACTTCACGCATTTTCGCCGAGTcaaccc                    >  2:181888/1‑139 (MQ=255)                                                                                                      tggtAGTAACAACGTGAGCATGAGGAACCGGGTTCTGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACg                  >  2:427759/1‑140 (MQ=255)                                                                                                              acaacGTGAGCATGAGGAACCGGGTTCTGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACg                  >  2:508305/1‑132 (MQ=255)                                                                                                              acaacGTGAGCATGAGGAACCGGGTTCTGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACg                  <  1:508305/132‑1 (MQ=255)                                                                                                               caacGTGAGCATGAGGAACCGGGTTCTGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAatat         >  1:533706/1‑140 (MQ=255)                                                                                                                 acGTGAGCATGAGGAACCGGGTTCTGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAATATGc       >  1:234958/1‑140 (MQ=255)                                                                                                                    tGAGCATGAGGAACCGGGTTCTGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAATATGCaga    >  2:407958/1‑140 (MQ=255)                                                                                                                      aGCATGAGGAACCGGGTTCTGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCATGCATTTTCGCCCAGTCCACCACGCCGGAATATGCAgaga  >  2:242976/1‑140 (MQ=255)                                                                                                                      aGCATGAGGAACCGGGTTCTGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAATATGCAgaga  >  2:151070/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                         ACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAATATGCAGAGA  >  NC_000913/2684731‑2684987

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