Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,693,818	C→G	V576L (GTG→CTG)	purL ←	phosphoribosylformyl‑glycineamide synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,693,818	0	C	G	100.0%	216.8 / NA	56	V576L (GTG→CTG)	purL	phosphoribosylformyl‑glycineamide synthetase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (31/25);  total (31/25)

GTCGCCTTTCGCTTTCAGCGTTTGTACATCGCGGGTCATCTTCGGCGTTTTACCAAGCAGGACGTCCAGCGGCAGATCGATCGGCTGATTATCAAAATGACGATCGTGCAGAGAAAGATGCAGTTCTTCGGTCGCTTCACCAATCACCGCGTAGGGTGCGCGCTCACGCTTACACAGTTCGTCAAACAGCGGTAATTGATCGGCAGCAACCGCCAGCACGTAGCGTTCCTGGGATTCGTTACACCAGATTTCCAGCGGGCTCATGCCCGGTTCGTCGCTTAGAATCTCGCGC  >  NC_000913/2693672‑2693963                                                                                                                                                   |                                                                                                                                                  gTCGCCTTTCGCTTTCAGCGTTTGTACATCGCGGGTCATCTTCGGCGTTTTACCAAGCAGGACGTCCAGCGGCAGATCGATCGGCTGATTATCAAAATGACGATCGTGCAGAGAAAGATGCAGTTCTTCGGTCGCTTCACCAATCAGcgcg                                                                                                                                               <  1:546231/151‑1 (MQ=255) gTCGCCTTTCGCTTTCAGCGTTTGTACATCGCGGGTCATCTTCGGCGTTTTACCAAGCAGGACGTCCAGCGGCAGATCGATCGGCTGATTATCAAAATGACGATCGTGCAGAGAAAGATGCAGTTCTTCGGTCGCTTCACCAATCAGcgcg                                                                                                                                               <  1:249655/151‑1 (MQ=255)        ttCGCTTTCAGCGTTTGTACATCGCGGGTCATCTTCGGCGTTTTACCAAGCAGGACGTCCAGCGGCAGATCGATCGGCTGATTATCAAAATGACGATCGTGCAGAGAAAGATGCAGTTCTTCGGTCGCTTCACCAATCAGCGCGTAGGGTg                                                                                                                                        >  2:779599/1‑151 (MQ=255)          cGCTTTCAGCGTTTGTACATCGCGGGTCATCTTCGGCGTTTTACCAAGCAGGACGTCCAGCGGCAGATCGATCGGCTGATTATCAAAATGACGATCGTGCAGAGAAAGATGCAGTTCTTCGGTCGCTTCACCAATCAGCGCGTAGGGTgcg                                                                                                                                      >  2:216983/1‑151 (MQ=255)           gCTTTCAGCGTTTGTACATCGCGGGTCATCTTCGGCGTTTTACCAAGCAGGACGTCCAGCGGCAGATCGATCGGCTGATTATCAAAATGACGATCGTGCAGAGAAAGATGCAGTTCTTCGGTCGCTTCACCAATCAGCGCGTAGGGTgcgc                                                                                                                                     <  1:77232/151‑1 (MQ=255)            cTTTCAGCGTTTGTACATCGCGGGTCATCTTCGGCGTTTTACCAAGCAGGACGTCCAGCGGCAGATCGATCGGCTGATTATCAAAATGACGATCGTGCAGAGAAAGATGCAGTTCTTCGGTCGCTTCACCAATCAGCGCGTAGGGTgcgcg                                                                                                                                    >  2:264034/1‑151 (MQ=255)            cTTTCAGCGTTTGTACATCGCGGGTCATCTTCGGCGTTTTACCAAGCAGGACGTCCAGCGGCAGATCGATCGGCTGATTATCAAAATGACGATCGTGCAGAGAAAGATGCAGTTCTTCGGTCGCTTCACCAATCAGCGCGTAGGGTgcgc                                                                                                                                     >  2:117709/1‑150 (MQ=255)              ttCAGCGTTTGTACATCGCGGGTCATCTTCGGCGTTTTACCAAGCAGGACGTCCAGCGGCAGATCGATCGGCTGATTATCAAAATGACGATCGTGCAGAGAAAGATGCAGTTCTTCGGTCGCTTCACCAATCAGCGCGTAGGGTgcgcg                                                                                                                                    <  2:663262/149‑1 (MQ=255)                 aGCGTTTGTACATCGCGGGTCATCTTCGGCGTTTTACCAAGCAGGACGTCCAGCGGCAGATCGATCGGCTGATTATCAAAATGACGATCGTGCAGAGAAAGATGCAGTTCTTCGGTCGCTTCACCAATCAGCGCGTAGGGTGCGCGCTCAc                                                                                                                               >  1:20527/1‑151 (MQ=255)                 aGCGTTTGTACATCGCGGGTCATCTTCGGCGTTTTACCAAGCAGGACGTCCAGCGGCAGATCGATCGGCTGATTATCAAAATGACGATCGTGCAGAGAAAGATGCAGTTCTTCGGTCGCTTCACCAATCAGCGCGTAGGGTGCGCGCTCAc                                                                                                                               >  1:89993/1‑151 (MQ=255)                   cGTTTGTACATCGCGGGTCATCTTCGGCGTTTTACCAAGCAGGACGTCCAGCGGCAGATCGATCGGCTGATTATCAAAATGACGATCGTGCAGAGAAAGATGCAGTTCTTCGGTCGCTTCACCAATCAGCGCGTAGGGTGCGCGCTCACGc                                                                                                                             <  2:227106/151‑1 (MQ=255)                   cGTTTGTACATCGCGGGTCATCTTCGGCGTTTTACCAAGCAGGACGTCCAGCGGCAGATCGATCGGCTGATTATCAAAATGACGATCGTGCAGAGAAAGATGCAGTTCTTCGGTCGCTTCACCAATCAGCGCGTAGGGTGCGCGCTCACGc                                                                                                                             >  2:152033/1‑151 (MQ=255)                   cGTTTGTACATCGCGGGTCATCTTCGGCGTTTTACCAAGCAGGACGTCCAGCGGCAGATCGATCGGCTGATTATCAAAATGACGATCGTGCAGAGAAAGATGCAGTTCTTCGGTCGCTTCACCAATCAGCGCGTAGGGTGCGCGCTCACGc                                                                                                                             >  1:760633/1‑151 (MQ=255)                        gTACATCGCGGGTCATCTTCGGCGTTTTACCAAGCAGGACGTCCAGCGGCAGATCGATCGGCTGATTATCAAAATGACGATCGTGCAGAGAAAGATGCAGTTCTTCGGTCGCTTCACCAATCAGCGCGTAGGGTGCGCGCTCACGCTTac                                                                                                                         <  2:626043/150‑1 (MQ=255)                            aTCGCGGGTCATCTTCGGCGTTTTACCAAGCAGGACGTCCAGCGGCAGATCGATCGGCTGATTATCAAAATGACGATCGTGCAGAGAAAGATGCAGTTCTTCGGTCGCTTCACCAATCAGCGCGTAGGGTGCGCGCTCACGCTTACACAGt                                                                                                                    <  2:355972/151‑1 (MQ=255)                              cgcgGGTCATCTTCGGCGTTTTACCAAGCAGGACGTCCAGCGGCAGATCGATCGGCTGATTATCAAAATGACGATCGTGCAGAGAAAGATGCAGTTCTTCGGTCGCTTCACCAATCAGCGCGTAGGGTGCGCGCTCACGCTTACACAGtt                                                                                                                   <  1:25226/150‑1 (MQ=255)                              cgcgGGTCATCTTCGGCGTTTTACCAAGCAGGACGTCCAGCGGCAGATCGATCGGCTGATTATCAAAATGACGATCGTGCAGAGAAAGATGCAGTTCTTCGGTCGCTTCACCAATCAGCGCGTAGGGTGCGCGCTCACGCTTACACAGTTc                                                                                                                  >  2:308223/1‑151 (MQ=255)                              cgcgGGTCATCTTCGGCGTTTTACCAAGCAGGACGTCCAGCGGCAGATCGATCGGCTGATTATCAAAATGACGATCGTGCAGAGAAAGATGCAGTTCTTCGGTCGCTTCACCAATCAGCGCGTAGGGTGCGCGCTCACGCTTACACAGTTc                                                                                                                  >  1:429242/1‑151 (MQ=255)                                          tCGGCGTTTTACCAAGCAGGACGTCCAGCGGCAGATCGATCGGCTGATTATCAAAATGACGATCGTGCAGAGAAAGATGCAGTTCTTCGGTCGCTTCACCAATCAGCGCGTAGGGTGCGCGCTCACGCTTACACAGTTCGTCAAACAg                                                                                                         <  2:172558/148‑1 (MQ=255)                                          tCGGCGTTTTACCAAGCAGGACGTCCAGCGGCAGATCGATCGGCTGATTATCAAAATGACGATCGTGCAGAGAAAGATGCAGTTCTTCGGTCGCTTCACCAATCAGCGCGTAGGGTGCGCGCTCACGCTTACACAGTTCGTCAAACAGc                                                                                                        <  2:576875/149‑1 (MQ=255)                                           cGGCGTTTTACCAAGCAGGACGTCCAGCGGCAGATCGATCGGCTGATTATCAAAATGACGATCGTGCAGAGAAAGATGCAGTTCTTCGGTCGCTTCACCAATCAGCGCGTAGGGTGCGCGCTCACGCTTACACAGTTCGTCAAACAGCGGt                                                                                                     <  1:143910/151‑1 (MQ=255)                                             gCGTTTTACCAAGCAGGACGTCCAGCGGCAGATCGATCGGCTGATTATCAAAATGACGATCGTGCAGAGAAAGATGCAGTTCTTCGGTCGCTTCACCAATCAGCGCGTAGGGTGCGCGCTCACGCTTACACAGTTCGTCAAACAGCGGTa                                                                                                    >  1:384731/1‑150 (MQ=255)                                                  ttACCAAGCAGGACGTCCAGCGGCAGATCGATCGGCTGATTATCAAAATGACGATCGTGCAGAGAAAGATGCAGTTCTTCGGTCGCTTCACCAATCAGCGCGTAGGGTGCGCGCTCACGCTTACACAGTTCGTCAAACAGCGGTAATTGAt                                                                                              >  2:658811/1‑151 (MQ=255)                                                      cAAGCAGGACGTCCAGCGGCAGATCGATCGGCTGATTATCAAAATGACGATCGTGCAGAGAAAGATGCAGTTCTTCGGTCGCTTCACCAATCAGCGCGTAGGGTGCGCGCTCACGCTTACACAGTTCGTCAAACAGCGGTAATTGATCGgc                                                                                          <  2:359779/151‑1 (MQ=255)                                                          cAGGACGTCCAGCGGCAGATCGATCGGCTGATTATCAAAATGACGATCGTGCAGAGAAAGATGCAGTTCTTCGGTCGCTTCACCAATCAGCGCGTAGGGTGCGCGCTCACGCTTACACAGTTCGTCAAACAGCGGTAATTGATCGgcagc                                                                                       >  1:489332/1‑150 (MQ=255)                                                               cGTCCAGCGGCAGATCGATCGGCTGATTATCAAAATGACGATCGTGCAGAGAAAGATGCAGTTCTTCGGTCGCTTCACCAATCAGCGCGTAGGGTGCGCGCTCACGCTTACACAGTTCGTCAAACAGCGGTAATTGATCGGCAGCAACCGc                                                                                 <  1:733303/151‑1 (MQ=255)                                                                   cAGCGGCAGATCGATCGGCTGATTATCAAAATGACGATCGTGCAGAGAAAGATGCAGTTCTTCGGTCGCTTCACCAATCAGCGCGTAGGGTGCGCGCTCACGCTTACACAGTTCGTCAAACAGCGGTAATTGATCGGCAGCAACCGCCAGc                                                                             >  2:537795/1‑151 (MQ=255)                                                                    aGCGGCAGATCGATCGGCTGATTATCAAAATGACGATCGTGCAGAGAAAGATGCAGTTCTTCGGTCGCTTCACCAATCAGCGCGTAGGGTGCGCGCTCACGCTTACACAGTTCGTCAAACAGCGGTAATTGATCGGCAGCAACCGCCAGc                                                                             >  1:369659/1‑150 (MQ=255)                                                                             tcgatcgGCTGATTATCAAAATGACGATCGTGCAGAGAAAGATGCAGTTCTTCGGTCGCTTCACCAATCAGCGCGTAGGGTGCGCGCTCACGCTTACACAGTTCGTCAAACAGCGGTAATTGATCGGCAGCAACCGCCAGCACGTAGCGtt                                                                   >  1:283260/1‑151 (MQ=255)                                                                             tcgatcgGCTGATTATCAAAATGACGATCGTGCAGAGAAAGATGCAGTTCTTCGGTCGCTTCACCAATCAGCGCGTAGGGTGCGCGCTCACGCTTACACAGTTCGTCAAACAGCGGTAATTGATCGGCAGCAACCGCCAGCACGTAGCGtt                                                                   >  1:504580/1‑151 (MQ=255)                                                                               gatcgGCTGATTATCAAAATGACGATCGTGCAGAGAAAGATGCAGTTCTTCGGTCGCTTCACCAATCAGCGCGTAGGGTGCGCGCTCACGCTTACACATTTCGTCAAACAGCGGTAATTGATCGGCAGCAACCGCCAGCACGTAGCGTTc                                                                  <  1:509749/150‑1 (MQ=255)                                                                                atcgGCTGATTATCAAAATGACGATCGTGCAGAGAAAGATGCAGTTCTTCGGTCGCTTCACCAATCAGCGCGTAGGGTGCGCGCTCACGCTTACACAGTTCGTCAAACAGCGGTAATTGATCGGCAGCAACCGCCAGCACGTAGCGTTCCt                                                                >  2:755910/1‑151 (MQ=255)                                                                                atcgGCTGATTATCAAAATGACGATCGTGCAGAGAAAGATGCAGTTCTTCGGTCGCTTCACCAATCAGCGCGTAGGGTGCGCGCTCACGCTTACACAGTTCGTCAAACAGCGGTAATTGATCGGCAGCAACCGCCAGCACGTAGCGTTCCt                                                                >  2:114430/1‑151 (MQ=255)                                                                                   ggCTGATTATCAAAATGACGATCGTGCAGAGAAAGATGCAGTTCTTCGGTCGCTTCACCAATCAGCGCGTAGGGTGCGCGCTCACGCTTACACAGTTCGTCAAACAGCGGTAATTGATCGGCAGCAACCGCCAGCACGTAGCGTTCCTggg                                                             <  2:122054/151‑1 (MQ=255)                                                                                    gCTGATTATCAAAATGACGATCGTGCAGAGAAAGATGCAGTTCTTCGGTCGCTTCACCAATCAGCGCGTAGGGTGCGCGCTCACGCTTACACAGTTCGTCAAACAGCGGTAATTGATCGGCAGCAACCGCCAGCACGTAGCGTTCCTggg                                                             >  2:258419/1‑150 (MQ=255)                                                                                      tGATTATCAAAATGACGATCGTGCAGAGAAAGATGCAGTTCTTCGGTCGCTTCACCAATCAGCGCGTAGGGTGCGCGCTCACGCTTACACAGTTCGTCAAACAGCGGTAATTGATCGGCAGCAACCGCCAGCACGTAGCGTTCCTGGGAtt                                                          >  1:789148/1‑151 (MQ=255)                                                                                      tGATTATCAAAATGACGATCGTGCAGAGAAAGATGCAGTTCTTCGGTCGCTTCACCAATCAGCGCGTAGGGTGCGCGCTCACGCTTACACAGTTCGTCAAACAGCGGTAATTGATCGGCAGCAACCGCCAGCACGTAGCGTTCCTGGGAtt                                                          >  2:353470/1‑151 (MQ=255)                                                                                             cAAAATGACGATCGTGCAGAGAAAGATGCAGTTCTTCGGTCGCTTCACCAATCAGCGCGTAGGGTGCGCGCTCACGCTTACACAGTTCGTCAAACAGCGGTAATTGATCGGCAGCAACCGCCAGCACGTAGCGTTCCTGGGATTCGTTa                                                     <  2:953761/149‑1 (MQ=255)                                                                                              aaaaTGACGATCGTGCAGAGAAAGATGCAGTTCTTCGGTCGCTTCACCAATCAGCGCGTAGGGTGCGCGCTCACGCTTACACAGTTCGTCAAACAGCGGTAATTGATCGGCAGCAACCGCCAGCACGTAGCGTTCCTGGGATTCGTTacac                                                  <  2:290852/151‑1 (MQ=255)                                                                                                   gACGATCGTGCAGAGAAAGATGCAGTTCTTCGGTCGCTTCACCAATCAGCGCGTAGGGTGCGCGCTCCCGCTTACACAGTTCGTCAAACAGCGGTAATTGATCGGCAGCAACCGCCAGCACGTAGCGTTCCTGGGATTCGTTACACCAGAt                                             <  2:865863/151‑1 (MQ=255)                                                                                                      gATCGTGCAGAGAAAGATGCAGTTCTTCGGTCGCTTCACCAATCAGCGCGTAGGGTGCGCGCTCACGCTTACACAGTTCGTCAAACAGCGGTAATTGATCGGCAGCAACCGCCAGCACGTAGCGTTCCTGGGATTCGTTACACCAGATTTc                                          >  2:709175/1‑151 (MQ=255)                                                                                                       atagTGCAGAGAAAAGTGAGGTGCTTCGTTCCATCCCCCAATCAGCGCGTGGGGTGCGCGCCCACGCTACCACCTTTCGTCAACCAGGGGTAATGGCCCGCAAGAAACCGCCAGCACGTAGCGTTCCTGGGATTCGTTACACCAGAt                                             <  1:788249/144‑1 (MQ=255)                                                                                                         cGTGCAGAGAAAGATGCAGTTCTTCGGTCGCTTCACCAATCAGCGCGTAGGGTGCGCGCTCACGCTTACACAGTTCGTCAAACAGCGTTAATTGATCGGCAGCAACCGCCAGCACGTAGCGTTCCTGGGATTCGTTACACCAGATTTCCa                                        <  1:18328/150‑1 (MQ=255)                                                                                                              agagAAAGATGCAGTTCTTCGGTCGCTTCACCAATCAGCGCGTAGGGTGCGCGCTCACGCTTACACAGTTCGTCAAACAGCGGTAATTGATCGGCAGCAACCGCCAGCACGTAGCGTTCCTGGGATTCGTTACACCAGATTTCCAGCGGGc                                  >  1:89756/1‑151 (MQ=255)                                                                                                              agagAAAGATGCAGTTCTTCGGTCGCTTCACCAATCAGCGCGTAGGGTGCGCGCTCACGCTTACACAGTCCGTCAAACAGCGGTAATTGATCGGCAGCAACCGCCAGCACGTAGCGTTCCTGGGATTCGTTACACCAGATTTCCAGCGGGc                                  >  2:903254/1‑151 (MQ=255)                                                                                                                  aaaGATGCAGTTCTTCGGTCGCTTCACCAATCAGCGCGTAGGGTGCGCGCTCACGCTTACACAGTTCGTCAAACAGCGGTAATTGATCGGCAGCAACCGCCAGCACGTAGCGTTCCTGGGATTCGTTACACCAGATTTCCAGCGGGCTCAt                              >  1:47361/1‑151 (MQ=255)                                                                                                                  aaaGATGCAGTTCTTCGGTCGCTTCACCAATCAGCGCGTAGGGTGCGCGCTCACGCTTACACAGTTCGTCAAACAGCGGTAATTGATCGGCAGCAACCGCCAGCACGTAGCGTTCCTGGGATTCGTTACACCAGATTTCCAGCGGGCTCAt                              >  1:57221/1‑151 (MQ=255)                                                                                                                  aaaGATGCAGTCCTTCGGTCGCTTCACCAATCAGCGCGTAGGGTGCGCGCCCACGCTTACACAGTTCGTCAAACAGCGGTAATTGATCGGCAGCAACCGCCAGCACGTAGCGTTCCTGGGATTCGTTACACCAGATTTCCAGCGGGCTCAt                              <  1:117707/151‑1 (MQ=255)                                                                                                                       tGCAGTTCTTCGGTCGCTTCACCAATCAGCGCGTAGGGTGCGCGCTCACGCTTACACAGTTCGTCAAACAGCGGTAATTGATCGGCAGCAACCGCCAGCACGTAGCGTTCCTGGGATTCGTTACACCAGATTTCCAGCGGGCTCATGccc                          <  2:802765/150‑1 (MQ=255)                                                                                                                        gCAGTTCTTCGGTCGCTTCACCAATCAGCGCGTAGGGTGCGCGCTCACGCTTACACAGTTCGTCAAACAGCGGTAATTGATCGGCAGCAACCGCCAGCACGTAGCGTTCCTGGGATTCGTTACACCAGATTTCCAGCGGGCTCATGCCCgg                        >  1:145583/1‑151 (MQ=255)                                                                                                                         cAGTTCTTCGGTCGCTTCACCAATCAGCGCGTAGGGTGCGCGCTCACGCTTACACAGTTCGTCAAACAGCGGTAATTGATCGGCAGCAACCGCCAGCACGTAGCGTTCCTGGGATTCGTTACACCAGATTTCCAGCGGGCTCATGCCCGGt                       <  2:917200/151‑1 (MQ=255)                                                                                                                                      gCTTCACCAATCAGCGCGTAGGGTGCGCGCTCACGCTTACACAGTTCGTCAAACAGCGGTAATTGATCGGCAGCAACCGCCAGCACGTAGCGTTCCTGGGATTCGTTACACCAGATTTCCAGCGGGCTCATGCCCGGTTCGTCGCTTAGa           >  2:671195/1‑150 (MQ=255)                                                                                                                                            ccAACCAGCGCGTAGGGTGCGCGCTCACGCTTACACATTTCGTCAAACAGCGGTAATTGATCGGCAGCAACCGCCAGCACGTAGCGTTCCTGGGATTCGTTACACCAGATTTCCAGCGGGCTCATGCCCGGTTCGTCGCTTAGAATCTcgc    <  1:907769/151‑1 (MQ=255)                                                                                                                                             cAATCAGCGCGTAGGGTGCGCGCTCACGCTTACACAGTTCGTCAAACAGCGGTAATTGATCGGCAGCAACCGCCAGCACGTAGCGTTCCTGGGATTCGTTACACCAGATTTCCAGCGGGCTCATGCCCGGTTCGTCGCTTAGAATCTcgcg   <  2:775010/151‑1 (MQ=255)                                                                                                                                             cAATCAGCGCGTAGGGTGCGCGCTCACGCTTACACAGTTCGTCAAACAGCGGTAATTGATCGGCAGCAACCGCCAGCACGTAGCGTTCCTGGGATTCGTTACACCAGATTTCCAGCGGGCTCATGCCCGGTTCGTCGCTTAGAATCTcgcg   >  1:931011/1‑151 (MQ=255)                                                                                                                                                tCAGCGCGTAGGGTGCGCGCTCACGCTTACACAGTTCGTCAAACAGCGGTAATTGATCGGCAGCAACCGCCAGCACGTAGCGTTCCTGGGATTCGTTACACCAGATTTCCAGCGGGCTCATGCCCGGTTCGTCGCTTAGAATCTCgcgc  >  1:673216/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                  GTCGCCTTTCGCTTTCAGCGTTTGTACATCGCGGGTCATCTTCGGCGTTTTACCAAGCAGGACGTCCAGCGGCAGATCGATCGGCTGATTATCAAAATGACGATCGTGCAGAGAAAGATGCAGTTCTTCGGTCGCTTCACCAATCACCGCGTAGGGTGCGCGCTCACGCTTACACAGTTCGTCAAACAGCGGTAATTGATCGGCAGCAACCGCCAGCACGTAGCGTTCCTGGGATTCGTTACACCAGATTTCCAGCGGGCTCATGCCCGGTTCGTCGCTTAGAATCTCGCGC  >  NC_000913/2693672‑2693963

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