Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,693,373	G→C	A724G (GCC→GGC)	purL ←	phosphoribosylformylglycinamide synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,693,373	0	G	C	100.0%	135.6 / NA	37	A724G (GCC→GGC)	purL	phosphoribosylformylglycinamide synthetase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (20/17);  total (20/17)

TCACCAGGGTGGCCTGCCGCCGCCATCCAGTTGGCGGAAAGTTTGATGCGTTTGATATCGCCAATTTGTGTTGCGGCGATGTTGGTTAACGCTTCACCGACCGCCAGACGGGCAGAGGCGGCGAAATCCAGCAGCGCAACCGGCGCACGCTC  >  NC_000913/2693298‑2693449                                                                            |                                                                             tCACCAGGGTGGCCTGCCGCCGCCATCCAGTTGGCGGAAAGTTTGATGCGTTTGATATCGCCAATTTGTGTTGCGCCGATg                                                                         >  1:447971/1‑81 (MQ=255)  cACCAGGGTGGCCTGCCGCCGCCATCCAGTTGGCGGAAAGTTTGATGCGTTTGATATCGCCAATTTGTGTTGCGCCGa                                                                           <  1:447933/78‑1 (MQ=255)  cACCAGGGTGGCCTGCCGCCGCCATCCAGTTGGCGGAAAGTTTGATGCGTTTGATATCGCCAATTTGTGTTGCGCCGAt                                                                          <  2:447932/79‑1 (MQ=255)  cACCAGGGTGGCCTGCCGCCGCCATCCAGTTGGCGGAAAGTTTGATGCGTTTGATATCGCCAATTTGTGTTGCGCCGAt                                                                          <  1:447931/79‑1 (MQ=255)   aCCAGGGTGGCCTGCCGCCGCCATCCAGTTGGCGGAAAGTTTGATGCGTTTGATATCGCCAATTTGTGTTGCGCCGATGtt                                                                       >  2:447977/1‑81 (MQ=255)   aCCAGGGTGGCCTGCCGCCGCCATCCAGTTGGCGGAAAGTTTGATGCGTTTGATATCGCCAATTTGTGTTGCGCCGATGt                                                                        <  2:447934/80‑1 (MQ=255)    ccAGGGTGGCCTGCCGCCGCCATCCAGTTGGCGGAAAGTTTGATGCGTTTGATATCGCCAATTTGTGTTGCGCCGATGTTg                                                                      <  1:447935/81‑1 (MQ=255)     cAGGGTGGCCTGCCGCCGCCATCCAGTTGGCGGAAAGTTTGATGCGTTTGATATCGCCAATTTGTGTTGCGCCGATGTTgg                                                                     <  1:447936/81‑1 (MQ=255)      aGGGTGGCCTGCCGCCGCCATCCAGTTGGCGGAAAGTTTGATGCGTTTGATATCGCCAATTTGTGTTGCGCCGATGTTGGt                                                                    >  2:447976/1‑81 (MQ=255)      aGGGTGGCCTGCCGCCGCCATCCAGTTGGCGGAAAGTTTGATGCGTTTGATATCGCCAATTTGTGTTGCGCCGATGTTGGt                                                                    >  1:447973/1‑81 (MQ=255)          tGGCCTGCCGCCGCCATCCAGTTGGCGGAAAGTTTGATGCGTTTGATATCGCCAATTTGTGTTGCGCCGATGTTGGTTAAc                                                                <  1:447937/81‑1 (MQ=255)                gccgccgccATCCAGTTGGCGGAAAGTTTGATGCGTTTGATATCGCCAATTTGTGTTGCGCCGATGTTGGTTAACGCTTCa                                                          >  1:447995/1‑81 (MQ=255)                  cgccgccATCCAGTTGGCGGAAAGTTTGATGCGTTTGATATCGCCAATTTGTGTTGCGCCGATGTTGGTTAACGCTTCacc                                                        >  1:447991/1‑81 (MQ=255)                     cgccATCCAGTTGGCGGAAAGTTTGATGCGTTTGATATCGCCAATTTGTGTTGCGCCGATGTTGGTTAACGCTTCaccgac                                                     >  1:447998/1‑81 (MQ=255)                      gccATCCAGTTGGCGGAAAGTTTGATGCGTTTGATATCGCCAATTTGTGTTGCGCCGATGTTGGTTAACGCTTCaccgacc                                                    >  1:447979/1‑81 (MQ=255)                      gccATCCAGTTGGCGGAAAGTTTGATGCGTTTGATATCGCCAATTTGTGTTGCGCCGATGTTGGTTAACGCTTCaccgacc                                                    >  1:447972/1‑81 (MQ=255)                           ccAGTTGGCGGAAAGTTTGATGCGTTTGATATCGCCAATTTGTGTTGCGCCGATGTTGGTTAACGCTTCACCGACCGc                                                  >  1:447982/1‑78 (MQ=255)                                   cGGAAAGTTTGATGCGTTTGATATCGCCAATTTGTGTTGCGCCGATGTTGGTTAACGCTTCACCGACCGCCAGACGGGCag                                       <  1:447938/81‑1 (MQ=255)                                            tGATGCGTTTGATATCGCCAATTTGTGTTGCGCCGATGTTGGTTAACGCTTCACCGACCGCCAGACGGGCAGAGgcggcg                               >  2:447999/1‑80 (MQ=255)                                            tGATGCGTTTGATATCGCCAATTTGTGTTGCGCCGATGTTGGTTAACGCTTCACCGACCGCCAGACGGGCAGAGgcggcg                               <  1:447939/80‑1 (MQ=255)                                            tGATGCGTTTGATATCGCCAATTTGTGTTGCGCCGATGTTGGTTAACGCTTCACCGACCGCCAGACGGGCAGAGgcggcg                               <  2:447941/80‑1 (MQ=255)                                            tGATGCGTTTGATATCGCCAATTTGTGTTGCGCCGATGTTGGTTAACGCTTCACCGACCGCCAGACGGGCAGAGgcggcg                               <  1:447940/80‑1 (MQ=255)                                             gATGCGTTTGATATCGCCAATTTGTGTTGCGCCGATGTTGGTTAACGCTTCACCGACCGCCAGACGGGCAGAGGCGGCGa                              <  1:447942/80‑1 (MQ=255)                                             gATGCGTTTGATATCGCCAATTTGTGTTGCGCCGATGTTGGTTAACGCTTCACCGACCGCCAGACGGGCAGAGGCGGCGa                              >  2:447983/1‑80 (MQ=255)                                              aTGCGTTTGATATCGCCAATTTGTGTTGCGCCGATGTTGGTTAACGCTTCACCGACCGCCAGACGGGCAGAGGCGGCGaaa                            >  1:447987/1‑81 (MQ=255)                                                gCGTTTGATATCGCCAATTTGTGTTGCGCCGATGTTGGTTAACGCTTCACCGACCGCCAGACGGGCAGAGGCGGCGAAATc                          >  1:447985/1‑81 (MQ=255)                                                       atatCGCCAATTTGTGTTGCGCCGATGTTGGTTAACGCTTCACCGACCGCCAGACGGGCAGAGGCGGCGAAATCCagcagc                   >  1:447996/1‑81 (MQ=255)                                                         atCGCCAATTTGTGTTGCGCCGATGTTGGTTAACGCTTCACCGACCGCCAGACGGGCAGAGGCGGCGAAATCCAGCAgcg                  >  1:447990/1‑80 (MQ=255)                                                          tCGCCAATTTGTGTTGCGCCGATGTTGGTTAACGCTTCACCGACCGCCAGACGGGCAGAGGCGGCGAAATCCAGCAgcgc                 <  2:447944/80‑1 (MQ=255)                                                          tCGCCAATTTGTGTTGCGCCGATGTTGGTTAACGCTTCACCGACCGCCAGACGGGCAGAGGCGGCGAAATCCAGCAGCGCa                <  2:447943/81‑1 (MQ=255)                                                            gCCAATTTGTGTTGCGCCGATGTTGGTTAACGCTTCACCGACCGCCAGACGGGCAGAGGCGGCGAAATCCAGCAGCGCAAc              >  1:448009/1‑81 (MQ=255)                                                            gCCAATTTGTGTTGCGCCGATGTTGGTTAACGCTTCACCGACCGCCAGACGGGCAGAGGCGGCGAAATCCAGCAGCGCAAc              <  2:447945/81‑1 (MQ=255)                                                              cAATTTGTGTTGCGCCGATGTTGGTTAACGCTTCACCGACCGCCAGACGGGCAGAGGCGGCGAAATCCAGCAGCGCAACCg            >  2:447993/1‑81 (MQ=255)                                                                aTTTGTGTTGCGCCGATGTTGGTTAACGCTTCACCGACCGCCAGACGGGCAGAGGCGGCGAAATCCAGCAGCGCAACCGgc          <  2:447947/81‑1 (MQ=255)                                                                     tgtTGCGCCGATGTTGGTTAACGCTTCACCGACCGCCAGACGGGCAGAGGCGGCGAAATCCAGCAGCGCAACCGGCGCACg     >  2:447992/1‑81 (MQ=255)                                                                       ttGCGCCGATGTTGGTTAACGCTTCACCGACCGCCAGACGGGCAGAGGCGGCGAAATCCAGCAGCGCAACCGGCGCACGCt   <  1:447948/81‑1 (MQ=255)                                                                        tGCGCCGATGTTGGTTAACGCTTCACCGACCGCCAGACGGGCAGAGGCGGCGAAATCCAGCAGCGCAACCGGCGCACGCTc  >  1:448005/1‑81 (MQ=255)                                                                            |                                                                             TCACCAGGGTGGCCTGCCGCCGCCATCCAGTTGGCGGAAAGTTTGATGCGTTTGATATCGCCAATTTGTGTTGCGGCGATGTTGGTTAACGCTTCACCGACCGCCAGACGGGCAGAGGCGGCGAAATCCAGCAGCGCAACCGGCGCACGCTC  >  NC_000913/2693298‑2693449

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