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

AACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGACCGGTAGCATCGATACCGTAAGGAACGATGTTGTACAGTTTACCGG  >  NC_000913/2684839‑2685101                                                                                                                                     |                                                                                                                                   aacaacGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGaa                                                                                                                              <  1:359060/139‑1 (MQ=255)  acaacGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAat                                                                                                                             <  1:511186/139‑1 (MQ=255)    aacGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAatat                                                                                                                           <  1:1027557/139‑1 (MQ=255)    aacGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAatat                                                                                                                           >  1:604896/1‑139 (MQ=255)      cGTGAGCATGAGGAACCGGATTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGc                                                                                                                         >  1:862702/1‑139 (MQ=255)       gTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGa                                                                                                                               >  2:544830/1‑132 (MQ=255)       gTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGa                                                                                                                               <  1:544830/132‑1 (MQ=255)          aGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCagag                                                                                                                     >  1:1009745/1‑139 (MQ=255)                        ggTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAAt                                                                                                       >  2:537289/1‑139 (MQ=255)                          ttCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCa                                                                                                     >  1:175390/1‑139 (MQ=255)                             gggTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCAttt                                                                                                  >  1:925328/1‑139 (MQ=255)                               gTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTc                                                                                                <  2:925328/139‑1 (MQ=255)                                   aCGCCAGCAGCAACCAGGCCCGCAACGTGCGCAATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCt                                                                                            >  2:767665/1‑139 (MQ=255)                                    cGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCtt                                                                                           <  2:175390/139‑1 (MQ=255)                                     gCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTtg                                                                                          <  2:1009745/139‑1 (MQ=255)                                                             gTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGa                                                                  <  2:604896/139‑1 (MQ=255)                                                                      atCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGc                                                                                             >  1:354125/1‑103 (MQ=255)                                                                      atCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGc                                                                                             <  2:354125/103‑1 (MQ=255)                                                                      atCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCACTCCACCACACCGGAATATGCAGAGAAACCACCGATAATATTTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAg                                                         >  1:230065/1‑139 (MQ=255)                                                                               cAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGa                                                >  2:1051047/1‑139 (MQ=255)                                                                                 ggTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATa                                                                                                         <  2:852338/80‑1 (MQ=255)                                                                                 ggTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATa                                                                                                         >  1:852338/1‑80 (MQ=255)                                                                                    aaGGACCGATGCTGTCAGCGATTTCACGCATTTTCGACCAGTCCACCACACCGGGATATGCAGAGAACCCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGACCGGt                                           <  2:230065/139‑1 (MQ=255)                                                                                    aaGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGACCGGt                                           <  1:1051047/139‑1 (MQ=255)                                                                                            aTGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGa                                                                                                                               <  2:659144/47‑1 (MQ=255)                                                                                            aTGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGa                                                                                                                               >  1:659144/1‑47 (MQ=255)                                                                                            aTGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCTAGATCGGCGTAGTCGATATGACCGGTAGCATCGa                                   >  1:722068/1‑139 (MQ=255)                                                                                            aTGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGACCGGTAGCATCGa                                   >  1:1066427/1‑139 (MQ=255)                                                                                            aTGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGACCGGTAGCATCGa                                   >  2:346138/1‑139 (MQ=255)                                                                                            aTGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGACCGGTAGCATCGa                                   >  1:1084150/1‑139 (MQ=255)                                                                                                 gTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGACCGGTAGCATCGATACCg                              >  2:502313/1‑139 (MQ=255)                                                                                                  tCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGACCGGTAGCATCGATACCGt                             <  1:435887/139‑1 (MQ=255)                                                                                                   cAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGACCGGTAGCATCGATACCGTa                            <  1:537289/139‑1 (MQ=255)                                                                                                        aTTTCACGCATTTTCGCCCAGTCCACCACATCGGAATATGCAGAGAAACCACCGATAATCATTTTGGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGACCGGTAGCATCGATACCGTAAGGaa                       <  2:862702/139‑1 (MQ=255)                                                                                                           tCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGACCGGTAGCATCGATACCGTAAGGAACGa                    <  1:346138/139‑1 (MQ=255)                                                                                                                    ttCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGACCGGTAGCATCGATACCGTAAGGAACGATGTTGTACa           >  2:580089/1‑139 (MQ=255)                                                                                                                     tCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGACCGGTAGCATCGATACCGTAAGGAACGATGTTGTACAg          <  2:1084150/139‑1 (MQ=255)                                                                                                                         ccAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGACCGGTAGCATCGATACCGTAAGGAACGATGTTGTACAGTTTa      <  2:277320/139‑1 (MQ=255)                                                                                                                         ccAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGACCGGTAGCATCGATACCGTAAGGAACGATGTTGTACAGTTTa      <  1:616390/139‑1 (MQ=255)                                                                                                                             tCCACCACACCGGAATATTCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGACCGGTAGCATCGATACCGTAAGGAACGATGTTGTACAGTTTACCgg  <  1:767665/139‑1 (MQ=255)                                                                                                                             tCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTc                                                                      >  2:900915/1‑71 (MQ=255)                                                                                                                             tCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTc                                                                      <  1:900915/71‑1 (MQ=255)                                                                                                                             tCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGACCGGTAGCATCGATACCGTAAGGaa                       >  1:692479/1‑118 (MQ=255)                                                                                                                             tCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGACCGGTAGCATCGATACCGTAAGGaa                       <  2:692479/118‑1 (MQ=255)                                                                                                                             tCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGACCGGTAGCATCGATACCGTAAGGAACGATGTTGTACAGTTTACCgg  >  2:675428/1‑139 (MQ=255)                                                                                                                                     |                                                                                                                                   AACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGACCGGTAGCATCGATACCGTAAGGAACGATGTTGTACAGTTTACCGG  >  NC_000913/2684839‑2685101

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