Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,622,640	T→C	V136A (GTG→GCG)	purN →	phosphoribosylglycinamide formyltransferase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,622,640	0	T	C	100.0%	212.3 / NA	58	V136A (GTG→GCG)	purN	phosphoribosylglycinamide formyltransferase 1
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (29/29);  total (29/29)

TCTCAGCCCGGCGTTTGTCTCCCACTATGCCGGGCGTTTGCTGAACATTCACCCTTCTCTGCTGCCGAAATATCCCGGATTACACACCCATCGTCAGGCGCTGGAAAATGGCGATGAAGAGCACGGTACATCGGTGCATTTCGTCACCGATGAACTGGACGGTGGCCCGGTTATTTTACAGGCGAAAGTCCCGGTATTTGCTGGTGATTCGGAAGATGACATCACCGCCCGCGTGCAAACCCAGGAACACGCCATTTATCCACTGGTGA  >  NC_000913/2622506‑2622774                                                                                                                                       |                                                                                                                                       tctcAGCCCGGCGTTTGTCTCCCACTATGCCGGGCGTTTGCTGAACATTCACCCTTCTCTGCTGCCGAAATATCCCGGATTACACACCCATCGTCAGGCGCTGGAAAATGGCGATGAAGAGCACGGTACATCGGCGCAtt                                                                                                                                   >  1:1032259/1‑140 (MQ=255)    cAGCCCGGCGTTTGTCTCCCACTATGCCGGGCGTTTGCTGAACATTCACCCTTCTCTGCTGCCGAAATATCCCGGATTACACACCCATCGTCAGGCGCTGGAAAATGGCGATGAAGAGCACGGTACATCGGCGCATTTCg                                                                                                                                >  2:435539/1‑140 (MQ=255)     aGCCCGGCGTTTGTCTCCCACTATGCCGGGCGTTTGCTGAACATTCACCCTTCTCTGCTGCCGAAATATCCCGGATTACACACCCATCGTCAGGCGCTGGAAAATGGCGATGAAGAGCACGGTACATCGGCGCATTTCGt                                                                                                                               <  2:729659/140‑1 (MQ=255)          ggCGTTTGTCTCCCACTATGCCGGGCGTTTGCTGAACATTCACCCTTCTCTGCTGCCGAAATATCCCGGATTACACACCCATCGTCAGGCGCTGGAAAATGGCGATGAAGAGCACGGTACATCGGCGCATTTCGTCACCg                                                                                                                          <  2:119372/140‑1 (MQ=255)           gCGTTTGTCTCCCACTATGCCGGGCGTTTGCTGAACATTCACCCTTCTCTGCTGCCGAAATATCCCGGATTACACACCCATCGTCAGGCGCTGGAAAATGGCGATGAAGAGCACGGTACATCGGCGCATTTCGTCACCGa                                                                                                                         >  1:985171/1‑140 (MQ=255)           gCGTTTGTCTCCCACTATGCCGGGCGTTTGCTGAACATTCACCCTTCTCTGCTGCCGAAATATCCCGGATTACACACCCATCGTCAGGCGCTGGAAAATGGCGATGAAGAGCACGGTACATCGGCGCATTTCGTCACCGa                                                                                                                         >  2:134637/1‑140 (MQ=255)                 gTCTCCCACTATGCCGGGCGTTTGCTGAACATTCACCCTTCTCTGCTGCCGAAATATCCCGGATTACACACCCATCGTCAGGCGCTGGAAAATGGCGATGAAGAGCACGGTACATCGGCGCATTTCGTCACCGATGAACt                                                                                                                   <  1:797784/140‑1 (MQ=255)                      ccACTATGCCGGGCGTTTGCTGAACATTCACCCTTCTCTGCTGCCGAAATATCCCGGATTACACACCCATCGTCAGGCGCTGGAAAATGGCGATGAAGAGCACGGTACATCGGCGCATTTCGTc                                                                                                                              >  1:309466/1‑124 (MQ=255)                      ccACTATGCCGGGCGTTTGCTGAACATTCACCCTTCTCTGCTGCCGAAATATCCCGGATTACACACCCATCGTCAGGCGCTGGAAAATGGCGATGAAGAGCACGGTACATCGGCGCATTTCGTc                                                                                                                              <  2:309466/124‑1 (MQ=255)                        aCTATGCCGGGCGTTTGCTGAACATTCACCCTTCTCTGCTGCCGAAATATCCCGGATTACACACCCATCGTCAGGCGCTGGAAAATGGCGATGAAGAGCACGGTACATCGGCGCATTTCGTCACCGATGAACTGGACGGt                                                                                                            <  1:134637/140‑1 (MQ=255)                         cTATGCCGGGCGTTTGCTGAACATTCACCCTTCTCTGCTGCCGAAATATCCCGGATTACACACCCATCGTCAGGCGCTGGAAAATGGCGATGAAGAGCACGGTACATCGGCGCATTTCGTCACCGATGAACTGGACGGTg                                                                                                           >  1:533561/1‑140 (MQ=255)                          tATGCCGGGCGTTTGCTGAACATTCACCCTTCTCTGCTGCCGAAATATCCCGGATTACACACCCATCGTCAGGCGCTGGAAAATGGCGATGAAGAGCACGGTACATCGGCGCATTTCGTCACCGATGAACTGGACGGTgg                                                                                                          <  2:533561/140‑1 (MQ=255)                           aTGCCGGGCGTTTGCTGAACATTCACCCTTCTCTGCTGCCGAAATATCCCGGATTACACACCCATCGTCAGGCGCTGGAAAATGGCGATGAAGAGCACGGTACATCGGCGCATTTCGTCACCGATGAACTGGACGGTGGc                                                                                                         <  1:327383/140‑1 (MQ=255)                            tGCCGGGCGTTTGCTGAACATTCACCCTTCTCTGCTGCCGAAATATCCCGGATTACACACCCATCGTCAGGCGCTGGAAAATGGCGATGAAGAGCACGGTACATCGGCGCATTTCGTCACCGATGAACTGGACGGTGGcc                                                                                                        >  2:1036413/1‑140 (MQ=255)                             gCCGGGCGTTTGCTGAACATTCACCCTTCTCTGCTGCCGAAATATCCCGGATTACACACCCATCGTCAGGCGCTGGAAAATGGCGATGAAGAGCACGGTACATCGGCGCATTTCGTCACCGATGAACTGGACGGTGGccc                                                                                                       <  1:1036413/140‑1 (MQ=255)                             gCCGGGCGTTTGCTGAACATTCACCCTTCTCTGCTGCCGAAATATCCCGGATTACACACCCATCGTCAGGCGCTGGAAAATGGCGATGAAGAGCACGGTACATCGGCGCATTTCGTCACCGATGAACTGGACGGTGGccc                                                                                                       <  1:505711/140‑1 (MQ=255)                                 ggCGTTTGCTGAACATTCACCCTTCTCTGCTGCCGAAATATCCCGGATTACACACCCATCGTCAGGCGCTGGAAAATGGCGATGAAGAGCACGGTACATCGGCGCATTTCGTCACCGATGAACTGGACGGTGGCCCGGtt                                                                                                   >  1:865692/1‑140 (MQ=255)                                 ggCGTTTGCTGAACATTCACCCTTCTCTGCTGCCGAAATATCCCGGATTACACACCCATCGTCAGGCGCTGGAAAATGGCGATGAAGAGCACGGTACATCGGCGCATTTCGTCACCGATGAACTGGACGGTGGCCCGGt                                                                                                    <  2:738431/139‑1 (MQ=255)                                 ggCGTTTGCTGAACATTCACCCTTCTCTGCTGCCGAAATATCCCGGATTACACACCCATCGTCAGGCGCTGGAAAATGGCGATGAAGAGCACGGTACATCGGCGCATTTCGTCACCGATGAACTGGACGGTGGCCCGGt                                                                                                    >  1:738431/1‑139 (MQ=255)                                    gTTTGCTGAACATTCACCCTTCTCTGCTGCCGAAATATCCCGGATTACACACCCATCGTCAGGCGCTGGAAAATGGCGATGAAGAGCACGGTACATCGGCGCATTTCGTCACCGATGAACTGGACGGTGGCCCGGTTAtt                                                                                                >  2:499722/1‑140 (MQ=255)                                    gTTTGCTGAACATTCACCCTTCTCTGCTGCCGAAATATCCCGGATTACACACCCATCGTCAGGCGCTGGAAAATGGCGATGAAGAGCACGGTACATCGGCGCATTTCGTCACCGATGAACTGGACGGTGGCCCGGTTAtt                                                                                                >  1:929117/1‑140 (MQ=255)                                       tGCTGAACATTCACCCTTCTCTGCTGCCGAAATATCCCGGATTACACACCCATCGTCAGGCGCTGGAAAATGGCGATGAAGAGCACGGTACATCGGCGCATTTCGTCACCGATGAACTGGACGGTGGCCCGGTTATTTTa                                                                                             >  1:1091760/1‑140 (MQ=255)                                             aCATTCACCCTTCTCTGCTGCCGAAATATCCCGGATTACACACCCATCGTCAGGCGCTGGAAAATGGCGATGAAGAGCACGGTACATCGGCGCATTTCGTCACCGATGAACTGGACGGTGGCCCGGTTATTTTACAGGCg                                                                                       <  2:1091760/140‑1 (MQ=255)                                             aCATTCACCCTTCTCTGCTGCCGAAATATCCCGGATTACACACCCATCGTCAGGCGCTGGAAAATGGCGATGAAGAGCACGGTACATCGGCGCATTTCGTCACCGATGAACTGGACGGTGGCCCGGTTATTTTACAGGCg                                                                                       <  2:985171/140‑1 (MQ=255)                                               aTTCACCCTTCTCTGCTGCCGAAATATCCCGGATTACACACCCATCGTCAGGCGCTGGAAAATGGCGATGAAGAGCACGGTACATCGGCGCATTTCGTCACCGATGAACTGGACGGTGGCCCGGTTATTTTACAGGCGaa                                                                                     <  2:1032259/140‑1 (MQ=255)                                                 tCACCCTTCTCTGCTGCCGAAATATCCCGGATTACACACCCATCGTCAGGCGCTGGAAAATGGCGATGAAGAGCACGGTACATCGGCGCATTTCGTCACCGATGAACTGGACGGTGGCCCGGTTATTTTACAGGCGAAAg                                                                                   >  1:941183/1‑140 (MQ=255)                                                   ccccTTCTCTGCTGCCGAAATATCCCGGATTACACACCCATCGTCAGGCGCTGGAAAATGGCGATGAAGAGCACGGTACATCGGCGCATTTCGTCACCGATGAACTGGACGGTGGCCCGGTTATTTTACAGGCGAAAGTc                                                                                 <  1:542042/139‑1 (MQ=255)                                                        tgaatgctggCGAAATATCCCGGATTACGCTCCCAGCGGCAGGCGCTGGAAAATGGCGATGAAGAGCACGGTACATCGGCGCATTTCGTCACCGATGAACTGGACGGTGGCCCGGTTATTTTACAGGCGAAAGTCCCGGt                                                                            <  2:1001333/136‑1 (MQ=255)                                                         ctctGCTGCCGAAATATCCCGGATTACACACCCATCGTCAGGCGCTGGAAAATGGCGATGAAGAGCACGGTACATCGGCGCATTTCGTCACCGATGAACTGGACGGTGGCCCGGTTATTTTACAGGCGAAAGTCCCGGTa                                                                           >  2:363015/1‑140 (MQ=255)                                                         ctctGCTGCCGAAATATCCCGGATTACACACCCATCGTCAGGCGCTGGAAAATGGCGATGAAGAGCACGGTACATCGGCGCATTTCGTCACCGATGAACTGGACGGTGGCCCGGTTATTTTACAGGCGAAAGTCCCGGTa                                                                           >  1:393547/1‑140 (MQ=255)                                                              ctgcCGAAATATCCCGGATTACACACCCATCGTCAGGCGCTGGAAAATGGCGATGAAGAGCACGGTACATCGGCGCATTTCGTCACCGATGAACTGGACGGTGGCCCGGTTATTTTACAGGCGAAAGTCCCGGTATTTGc                                                                      >  2:364101/1‑140 (MQ=255)                                                               tgcCGAAATATCCCGGATTACACACCCATCGTCAGGCGCTGGAAAATGGCGATGAAGAGCACGGTACATCGGCGCATTTCGTCACCGATGAACTGGACGGTGGCCCGGTTATTTTACAGGCGAAAGTCCCGGTATTTGCt                                                                     >  2:713736/1‑140 (MQ=255)                                                                 ccGAAATATCCCGGATTACACACCCATCGTCAGGCGCTGGAAAATGGCGATGAAGAGCACGGTACATCGGCGCATTTCGTCACCGATGAACTGGACGGTGGCCCGGTTATTTTACAGGCGAAAGTCCCGGTATTTGCTgg                                                                   >  2:870062/1‑140 (MQ=255)                                                                  cGAAATATCCCGGATTACACACCCATCGTCAGGCGCTGGAAAATGGCGATGAAGAGCACGGTACATCGGCGCATTTCGTCACCGATGAACTGGACGGTGGCCCGGTTATTTTACAGGCGAAAGTCCCGGTATTTGCTGGt                                                                  <  1:640183/140‑1 (MQ=255)                                                                   gAAATATCCCGGATTACACACCCATCGTCAGGCGCTGGAAAATGGCGATGAAGAGCACGGTACATCGGCGCATTTCGTCACCGATGAACTGGACGGTGGCCCGGTTATTTTACAGGCGAAAGTCCCGGTATTTGCTGGTg                                                                 >  2:902574/1‑140 (MQ=255)                                                                           ccGGATTACACACCCATCGTCAGGCGCTGGAAAATGGCGATGAAGAGCACGGTACATCGGCGCATTTCGTCACCGATGAACTGGACGGTGGCCCGGTTATTTTACAGGCGAAAGTCCCGGTATTTGCTGGTGATTCGGaa                                                         >  1:1127467/1‑140 (MQ=255)                                                                                        ccATCGTCAGGCGCTGGAAAATGGCGATGAAGAGCACGGTACATCGGCGCATTTCGTCACCGATGAACTGGACGGTGGCCCGGTTATTTTACAGGCGAAAGTCCCGGTATTTGCTGGTGATTCGGAAGATGACATCAccg                                            >  1:13613/1‑140 (MQ=255)                                                                                         cATCGTCAGGCGCTGGAAAATGGCGATGAAGAGCACGGTACATCGGCGCATTTCGTCACCGATGAACTGGACGGTGGCCCGGTTATTTTACAGGCGAAAGTCCCGGTATTTGCTGGTGATTCGGAAGATGACATCAccgc                                           >  1:233497/1‑140 (MQ=255)                                                                                         cATCGTCAGGCGCTGGAAAATGGCGATGAAGAGCACGGTACATCGGCGCATTTCGTCACCGATGAACTGGACGGTGGCCCGGTTATTTTACAGGCGAAAGTCCCGGTATTTGCTGGTGATTCGGAAGATGACATCAccgc                                           >  1:111476/1‑140 (MQ=255)                                                                                           tCGTCAGGCGCTGGAAAATGGCGATGAAGAGCACGGTACATCGGCGCATTTCGTCACCGATGAACTGGACGGTGGCCCGGTTATTTTACAGGCGAAAGTCCCGGTATTTGCTGGTGATTCGGAAGATGACATCAccgccc                                         >  1:204913/1‑140 (MQ=255)                                                                                           tCGTCAGGCGCTGGAAAATGGCGATGAAGAGCACGGTACATCGGCGCATTTCGTCACCGATGAACTGGACGGTGGCCCGGTTATTTTACAGGCGAAAGTCCCGGTATTTGCTGGTGATTCGGAAGATGACATCAccgccc                                         <  1:713736/140‑1 (MQ=255)                                                                                           tCGTCAGGCGCTGGAAAATGGCGATGAAGAGCACGGTACATCGGCGCATTTCGTCACCGATGAACTGGACGGTGGCCCGGTTATTTTACAGGCGAAAGTCCCGGTATTTGCTGGTGATTCGGAAGATGACATCAccgcca                                         <  2:393547/140‑2 (MQ=255)                                                                                               cAGGCGCTGGAAAATGGCGATGAAGAGCACGGTACATCGGCGCATTTCGTCACCGATGAACTGGACGGTGGCCCGGTTATTTTACAGGCGAAAGTCCCGGTATTTGCTGGTGATTCGGAAGATGACATCACCGCCCGCGt                                     <  2:111476/140‑1 (MQ=255)                                                                                                  gcgcTGGAAAATGGCGATGAAGAGCACGGTACATCGGCGCATTTCGTCACCGATGAACTGGACGGTGGCCCGGTTATTTTACAGGCGAAAGTCCCGGTATTTGCTGGTGATTCGGAAGATGACATCACCGCCCGCGTGCa                                  <  1:870062/140‑1 (MQ=255)                                                                                                  gcgcTGGAAAATGGCGATGAAGAGCACGGTACATCGGCGCATTTCGTCACCGATGAACTGGACGGTGGCCCGGTTATTTTACAGGCGAAAGTCCCGGTATTTGCTGGTGATTCGGAAGATGACATCACCGCCCGCGTGCa                                  <  2:13613/140‑1 (MQ=255)                                                                                                     cTGGAAAATGGCGATGAAGAGCACGGTACATCGGCGCATTTCGTCACCGATGAACTGGACGGTGGCCCGGTTATTTTACAGGCGAAAGTCCCGGTATTTGCTGGTGATTCGGAAGATGACATCACCGCCCGCGTGCAAAc                               >  2:312979/1‑140 (MQ=255)                                                                                                       ggAAAATGGCGATGAAGAGCACGGTACATCGGCGCATTTCGTCACCGATGAACTGGACGGTGGCCCGGTTATTTTACAGGCGAAAGTCCCGGTATTTGCTGGTGATTCGGAAGATGACATCACCGCCCGCGTGCAAAccc                             <  2:19574/140‑1 (MQ=255)                                                                                                             tGGCGATGAAGAGCACGGTACATCGGCGCATTTCGTCACCGATGAACTGGACGGTGGCCCGGTTATTTTACAGGCGAAAGTCCCGGTATTTGCTGGTGATTCGGAAGATGACATCACCGCCCGCGTGCAAACCCAGGAac                       >  1:379756/1‑140 (MQ=255)                                                                                                             tGGCGATGAAGAGCACGGTACATCGGCGCATTTCGTCACCGATGAACTGGACGGTGGCCCGGTTATTTTACAGGCGAAAGTCCCGGTATTTGCTGGTGATTCGGAAGATGACATCACCGCCCGCGTGCAAACCCAGGAac                       <  2:941183/140‑1 (MQ=255)                                                                                                              ggCGATGAAGAGCACGGTACATCGGCGCATTTCGTCACCGATGAACTGGACGGTGGCCCGGTTATTTTACAGGCGAAAGTCCCGGTATTTGCTGGTGATTCGGAAGATGACATCACCGCCCGCGTGCAAACCCAGGAaca                      <  1:435539/140‑1 (MQ=255)                                                                                                              ggCGATGAAGAGCACGGTACATCGGCGCATTTCGTCACCGATGAACTGGACGGTGGCCCGGTTATTTTACAGGCGAAAGTCCCGGTATTTGCTGGTGATTCGGAAGATGACATCACCGCCCGCGTGCAAACCCAGGAaca                      <  1:499722/140‑1 (MQ=255)                                                                                                               gCGATGAAGAGCACGGTACATCGGCGCATTTCGTCACCGATGAACTGGACGGTGGCCCGGTTATTTTACAGGCGAAAGTCCCGGTATTTGCTGGTGATTCGGAAGATGACATCACCGCCCGCGTGCAAACCCAGGAacac                     <  2:204913/140‑1 (MQ=255)                                                                                                                 gATGAAGAGCACGGTACATCGGCGCATTTCGTCACCGATGAACTGGACGGTGGCCCGGTTATTTTACAGGCGAAAGTCCCGGTATTTGCTGGTGATTCGGAAGATGACATCACCGCCCGCGTGCAAACCCAGGAACACGc                   <  2:929117/140‑1 (MQ=255)                                                                                                                  aTGAAGAGCACGGTACATCGGCGCATTTCGTCACCGATGAACTGGACGGTGGCCCGGTTATTTTACAGGCGAAAGTCCCGGTATTTGCTGGTGATTCGGAAGATGACATCACCGCCCGCGTGCAAACCCAGGAACACGcc                  >  2:631946/1‑140 (MQ=255)                                                                                                                         gCACGGTACATCGGCGCATTTCGTCACCGATGAACTGGACGGTGGCCCGGTTATTTTACAGGCGAAAGTCCCGGTATTTGCTGGTGATTCGGAAGATGACATCACCGCCCGCGTGCAAACCCAGGAACACGCCATTTATc           <  2:233497/140‑1 (MQ=255)                                                                                                                         gCACGGTACATCGGCGCATTTCGTCACCGATGAACTGGACGGTGGCCCGGTTATTTTACAGGCGAAAGTCCCGGTATTTGCTGGTGATTCGGAAGATGACATCACCGCCCGCGTGCAAACCCAGGAACACGCCATTTATc           <  2:1127467/140‑1 (MQ=255)                                                                                                                           aCGGTACATCGGCGCATTTCGTCACCGATGAACTGGACGGTGGCCCGGTTATTTTACAGGCGAAAGTCCCGGTATTTGCTGGTGATTCGGAAGATGACATCACCGCCCGCGTGCAAACCCAGGAACTCGCCATTTATCCa         <  1:312979/140‑1 (MQ=255)                                                                                                                               tACATCGGCGCATTTCGTCACCGATGAACTGGACGGTGGCCCGGTTATTTTACAGGCGAAAGTCCCGGTATTTGCTGGTGATTCGGAAGATGACATCACCGCCCGCGTGCAAACCAAGGAACACGCCATTTATCCACTgg     >  2:834850/1‑140 (MQ=255)                                                                                                                               tACATCGGCGCATTTCGTCACCGATGAACTGGACGGTGGCCCGGTTATTTTACAGGCGAAAGTCCCGGTATTTGCTGGTGATTCGGAAGATGACATCACCGCCCGCGTGCAAACCAAGGAACACGCCATTTATCCACTgg     <  1:834850/140‑1 (MQ=255)                                                                                                                                  aTCGGCGCATTTCGTCACCGATGAACTGGACGGTGGCCCGGTTATTTTACAGGCGAAAGTCCCGGTATTTGCTGGTGATTCGGAAGATGACATCACCGCCCGCGTGCAAACCCAGGAACACGCCATTTATCCACTGGTGa  >  2:561306/1‑140 (MQ=255)                                                                                                                                       |                                                                                                                                       TCTCAGCCCGGCGTTTGTCTCCCACTATGCCGGGCGTTTGCTGAACATTCACCCTTCTCTGCTGCCGAAATATCCCGGATTACACACCCATCGTCAGGCGCTGGAAAATGGCGATGAAGAGCACGGTACATCGGTGCATTTCGTCACCGATGAACTGGACGGTGGCCCGGTTATTTTACAGGCGAAAGTCCCGGTATTTGCTGGTGATTCGGAAGATGACATCACCGCCCGCGTGCAAACCCAGGAACACGCCATTTATCCACTGGTGA  >  NC_000913/2622506‑2622774

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