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

TAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGACCGGTAGCATCGATACCGTAAGGAACGATGTTGTACAGTTTACCGG  >  NC_000913/2684838‑2685101                                                                                                                                      |                                                                                                                                   taACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGa                                                                                                                               <  1:1253308/139‑1 (MQ=255)   acaacGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAat                                                                                                                             <  1:514209/139‑1 (MQ=255)    caacGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAata                                                                                                                            <  1:1195581/139‑1 (MQ=255)     aacGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAatat                                                                                                                           >  1:750817/1‑139 (MQ=255)        gTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCa                                                                                                                        >  1:766552/1‑139 (MQ=255)         tGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCag                                                                                                                       >  2:330885/1‑139 (MQ=255)             cATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAata                                                                                                                            <  1:366469/130‑1 (MQ=255)             cATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAata                                                                                                                            >  2:366469/1‑130 (MQ=255)               tGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAac                                                                                                                 <  2:1309814/139‑1 (MQ=255)                          gTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGTGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATc                                                                                                      >  2:1175236/1‑139 (MQ=255)                           ttCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCa                                                                                                     >  2:27818/1‑139 (MQ=255)                                gTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTc                                                                                                <  2:766552/139‑1 (MQ=255)                                   gACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGc                                                                                             <  2:750817/139‑1 (MQ=255)                                         agcagcAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTtg                                                                                          <  2:850880/136‑1 (MQ=255)                                         agaagcaACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTtg                                                                                          >  1:850880/4‑136 (MQ=255)                                                aCCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCg                                                                                               <  2:709430/124‑1 (MQ=255)                                                aCCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCg                                                                                               >  1:709430/1‑124 (MQ=255)                                                              gTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGTGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGa                                                                  <  1:1175236/139‑1 (MQ=255)                                                               tGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGAt                                                                 >  2:472145/1‑139 (MQ=255)                                                                 cgcCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCg                                                               <  1:472145/139‑1 (MQ=255)                                                                    cATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCg                                                            <  2:75427/139‑1 (MQ=255)                                                                                cAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGa                                                >  1:306745/1‑139 (MQ=255)                                                                                   gTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGt                                                                           >  1:14807/1‑109 (MQ=255)                                                                                   gTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGt                                                                           <  2:14807/109‑1 (MQ=255)                                                                                      aGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCg                                                                                                            >  2:1213149/1‑73 (MQ=255)                                                                                      aGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCg                                                                                                            <  1:1213149/73‑1 (MQ=255)                                                                                        cACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGa                                                                                                                               <  2:1064149/52‑1 (MQ=255)                                                                                        cACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGa                                                                                                                               >  1:1064149/1‑52 (MQ=255)                                                                                         aCCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAacc                                                                                                                <  2:505179/66‑1 (MQ=255)                                                                                         aCCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAacc                                                                                                                >  1:505179/1‑66 (MQ=255)                                                                                           cGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAat                                                                                                                             <  1:1186550/51‑1 (MQ=255)                                                                                           cGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAat                                                                                                                             >  2:1186550/1‑51 (MQ=255)                                                                                           cGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCa                                                                                                     >  1:264992/1‑75 (MQ=255)                                                                                           cGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCa                                                                                                     <  2:264992/75‑1 (MQ=255)                                                                                           cGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGGGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGACCGGTAGCATc                                     >  2:366531/1‑139 (MQ=255)                                                                                             aTGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTa                                                          >  2:589798/1‑116 (MQ=255)                                                                                             aTGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTa                                                          <  1:589798/116‑1 (MQ=255)                                                                                             aTGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGACCGGTa                                          <  2:62367/132‑1 (MQ=255)                                                                                             aTGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGACCGGTa                                          >  1:62367/1‑132 (MQ=255)                                                                                             aTGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGACCGGTAGCATCGa                                   >  1:1051139/1‑139 (MQ=255)                                                                                                  gTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGa                                                                                                                               <  2:122424/42‑1 (MQ=255)                                                                                                  gTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGa                                                                                                                               >  1:122424/1‑42 (MQ=255)                                                                                                  gTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAAt                                                                                                       <  2:1215271/66‑1 (MQ=255)                                                                                                  gTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAAt                                                                                                       >  1:1215271/1‑66 (MQ=255)                                                                                                  gTCAGCGATTTCACGCATTTTCGCCCACTCCACCACACCGGAATATGCAGGGAAACCAACGATAATAGTTTACGGCTTGTGCTCTTTGGCCTGTTTTTCCAGATCGGCGTAGTCGATATGACAGGTAGCATCGATACCg                              >  1:1393163/1‑139 (MQ=255)                                                                                                    cAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGACCGGTAGCATCGATACCGTa                            <  1:27818/139‑1 (MQ=255)                                                                                                     aGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGa                                                                                                                               >  1:66125/1‑39 (MQ=255)                                                                                                     aGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGa                                                                                                                               <  2:66125/39‑1 (MQ=255)                                                                                                         aTTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGACCGGTAGCATCGATACCGTAAGGaa                       <  2:306745/139‑1 (MQ=255)                                                                                                                gCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGACCGGTAGCATCGATACCGTAAGGAACGAtg                  >  2:237321/1‑137 (MQ=255)                                                                                                                gCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGACCGGTAGCATCGATACCGTAAGGAACGAtg                  <  1:237321/137‑1 (MQ=255)                                                                                                                   ttttCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGACCGGTAGCATCGATACCGTAAGGAACGATGTTGTa             <  1:972446/139‑1 (MQ=255)                                                                                                                              tCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTa                                                          <  1:747741/83‑1 (MQ=255)                                                                                                                              tCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTa                                                          >  2:747741/1‑83 (MQ=255)                                                                                                                              tCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGACCGGTAGCATCGATACCGTAAGGAACGATGTTGTACAGTTTACCgg  >  2:1476614/1‑139 (MQ=255)                                                                                                                              tCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGACCGGTAGCATCGATACCGTAAGGAACGATGTTGTACAGTTTACCgg  >  2:1106976/1‑139 (MQ=255)                                                                                                                                      |                                                                                                                                   TAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGACCGGTAGCATCGATACCGTAAGGAACGATGTTGTACAGTTTACCGG  >  NC_000913/2684838‑2685101

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