Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,695,419	T→C	D42G (GAC→GGC)	purL ←	phosphoribosylformylglycinamide synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,695,419	0	T	C	100.0%	175.7 / NA	50	D42G (GAC→GGC)	purL	phosphoribosylformylglycinamide synthetase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (29/21);  total (29/21)

GGGAGAGATGGTGCCAGGACGCGGGGTCACCAGCAGGAGTTTGCCTTGCGGGGCGTGGCTGGCGAGTGCCGGGCCATATTTCAGCAGGCGTTCAAGTTGTGCGTGCTCATCATCGTTTAACGGCGCATTGAGGTCAGCAAAATGGACATACTCGGCGTAAATATTGTGAACCGGGAGCCTGGCAGCCTGAAAACGTGCCAGCAGTTTGTTGATTCGGAATGCCGACAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTC  >  NC_000913/2695286‑2695552                                                                                                                                      |                                                                                                                                      gggAGAGATGGTGCCCGGACGCCGGATCACCAGCAGGCGTTTGCCTTGCGGGGCGTGGCTGGCGAGTGCCGGGCCATATTTCAGCAGGCGTTCAAGTTGTGCGTGCTCATCATCGTTTAACGGCGCATTGAGGCCAGCaa                                                                                                                                 <  1:283246/140‑1 (MQ=255) gggAGAGATGGTGCCAGGACGCGGGGTCACCAGCAGGAGTTTGCCTTGCGGGGCGTGGCTGGCGAGTGCCGGGCCATATTTCAGCAGGCGTTCAAGTTGTGCGTGCTCATCATCGTTTAACGGCGCATTGAGGCCAGCaa                                                                                                                                 >  1:46646/1‑140 (MQ=255) gggAGAGATGGTGCCAGGACGCGGGGTCACCAGCAGGAGTTTGCCTTGCGGGGCGTGGCTGGCGAGTGCCGGGCCATATTTCAGCAGGCGTTCAAGTTGTGCGTGCTCATCATCGTTTAACGGCGCATTGAGGCCAGCaa                                                                                                                                 <  1:1064/140‑1 (MQ=255)   gagagaTGGTGCCAGGACGCGGGGTCACCAGCAGGAGTTTGCCTTGCGGGGCGTGGCTGGCGAGTGCCGGGCCATATTTCAGCAGGCGTTCAAGTTGTGCGTGCTCATCATCGTTTAACGGCGCATTGAGGCCAGCaaaa                                                                                                                               <  2:714409/140‑1 (MQ=255)           gTGCCAGGACGCGGGGTCACCAGCAGGAGTTTGCCTTGCGGGGCGTGGCTGGCGAGTGCCGGGCCATATTTCAGCAGGCGTTCAAGTTGTGCGTGCTCATCATCGTTTAACGGCGCATTGAGGCCAGCAAAATGGACATa                                                                                                                       >  1:647141/1‑140 (MQ=255)            tGCCAGGACGCGGGGTCACCAGCAGGAGTTTGCCTTGCGGGGCGTGGCTGGCGAGTGCCGGGCCATATTTCAGCAGGCGTTCAAGTTGTGCGTGCTCATCATCGTTTAACGGCGCATTGAGGCCAGCAAAATGGACATAc                                                                                                                      >  2:849855/1‑140 (MQ=255)            tGCCAGGACGCGGGGTCACCAGCAGGAGTTTGCCTTGCGGGGCGTGGCTGGCGAGTGCCGGGCCATATTTCAGCAGGCGTTCAAGTTGTGCGTGCTCATCATCGTTTAACGGCGCATTGAGGCCAGCAAAATGGACATAc                                                                                                                      >  2:739709/1‑140 (MQ=255)            tGCCAGGACGCGGGGTCACCAGCAGGAGTTTGCCTTGCGGGGCGTGGCTGGCGAGTGCCGGGCCATATTTCAGCAGGCGTTCAAGTTGTGCGTGCTCATCATCGTTTAACGGCGCATTGAGGCCAGCAAAATGGACATAc                                                                                                                      >  2:700220/1‑140 (MQ=255)             gCCAGGACGCGGGGTCACCAGCTGGTGTTTGCCTTGCCGGGCGTGGCTGGCGAGTGCCGGGCCATATTTCAGCAGGCGTTCAAGTTGTGCGTGCTCATCATCGTTTAACGGCGCATTGAGGCCAGCAAAATGGACATACt                                                                                                                     <  1:395576/140‑1 (MQ=255)             gCCAGGACGCGGGGTCACCAGCAGGAGTTTGCCTTGCGGGGCGTGGCTGGCGAGTGCCGGGCCATATTTCAGCAGGCGTTCAAGTTGTGCGTGCTCATCATCGTTTAACGGCGCATTGAGGCCAGCAAAATGGACATACt                                                                                                                     <  2:947151/140‑1 (MQ=255)             gCCAGGACGCGGGGTCACCAGCAGGAGTTTGCCTTGCGGGGCGTGGCTGGCGAGTGCCGGGCCATATTTCAGCAGGCGTTCAAGTTGTGCGTGCTCATCATCGTTTAACGGCGCATTGAGGCCAGCAAAATGGACATACt                                                                                                                     <  1:739709/140‑1 (MQ=255)              ccAGGCCGCGGGGTCACCAGCAGGAGTTTGCCTTGCGGGGCGTGGCTGGCGAGTGCCGGGCCATATTTCAGCAGGCGTTCAAGTTGTGCGTGCTCATCATCGTTTAACGGCGCATTGAGGCCAGCAAAATGGACATACTc                                                                                                                    <  1:167375/140‑1 (MQ=255)              ccAGGACGCGGGGTCACCAGCAGGAGTTTGCCTTGCGGGGCGTGGCTGGCGAGTGCCGGGCCATATTTCAGCAGGCGTTCAAGTTGTGCGTGCTCATCATCGTTTAACGGCGCATTGAGGCCAGCAAAATGGACATACTc                                                                                                                    <  1:215540/140‑1 (MQ=255)              ccAGGACGCGGGGTCACCAGCAGGAGTTTGCCTTGCGGGGCGTGGCTGGCGAGTGCCGGGCCATATTTCAGCAGGCGTTCAAGTTGTGCGTGCTCATCATCGTTTAACGGCGCATTGAGGCCAGCAAAATGGACATACTc                                                                                                                    >  1:709978/1‑140 (MQ=255)              ccAGGACGCGGGGTCACCAGCAGGAGTTTGCACTGGGGGCCGTGGCTGGCGAGTGCCGGGCCATATTTCAGCAGGCGTTCAAGTTGTGCGTGCTCATCATCGTTTAACGGCGCATTGAGGCCAGCAAAATGGACATACTc                                                                                                                    <  2:388891/140‑1 (MQ=255)                      cgGGGTCACCAGCAGGAGTTTGCCTTGCGGGGCGTGGCTGGCGAGTGCCGGGCCATATTTCAGTAGGCGTTCAAGTTGTGCGTGCTCATCATCGTTTAACGGCGCATTGAGGCCAGCAAAATGGACATACTCGGCGTaaa                                                                                                            >  2:96831/1‑140 (MQ=255)                      cgGGGTCACCAGCAGGAGTTTGCCTTGCGGGGCGTGGCTGGCGAGTGCCGGGCCATATTTCAGCAGGCGTTCAAGTTGTGCGTGCTCATCATCGTTTAACGGCGCATTGAGGCCAGCAAAATGGACATACTCGGCGTaaa                                                                                                            >  1:867900/1‑140 (MQ=255)                          gTCACCAGCAGGAGTTTGCCTTGCGGGGCGTGGCTGGCGAGTGCCGGGCCATATTTCAGCAGGCGTTCAAGTTGTGCGTGCTCATCATCGTTTAACGGCGCATTGAGGCCAGCAAAATGGACATACTCGGCGTAAATAtt                                                                                                        <  2:647141/140‑1 (MQ=255)                                          tGCCTTGCGGGGCGTGGCTGGCGAGTGCCGGGCCATATTTCAGCAGGCGTTCAAGTTGTGCGTGCTCATCATCGTTTAACGGCGCATTGAGGCCAGCAAAATGGACATACTCGGGGTGAATATTATGAACCGGGAGCCTg                                                                                        >  1:284022/1‑140 (MQ=255)                                          tGCCTTGCGGGGCGTGGCTGGCGAGTGCCGGGCCATATTTCAGCAGGCGTTCAAGTTGTGCGTGCTCATCATCGTTTAACGGCGCATTGAGGCCAGCAAAATGGACATACTCGGCGTAAATATTGTGAACCGGGAGCCTg                                                                                        >  1:926606/1‑140 (MQ=255)                                           gCCTTGCGGGGCGTGGCTGGCGAGTGCCGGGCCATATTTCAGCAGGCGTTCAAGTTGTGCGTGCTCATCATCGTTTAACGGCGCATTGAGGCCAGCAAAATGGACATACTCGGCGTAAATATTGTGAACCGGGAGCCTgg                                                                                       <  2:867900/140‑1 (MQ=255)                                            ccTTGCGGGGCGTGGCTGGCGAGTGCCGGGCCATATTTCAGCAGGCGTTCAAGTTGTGCGTGCTCATCATCGTTTAACGGCGCATTGAGGCCAGCAAAATGGACATACTCGGCGTAAATATTGTGAACCGGGAGCCTGGc                                                                                      <  1:233335/140‑1 (MQ=255)                                               tGCGGGGCGTGGCTGGCGAGTGCCGGGCCATATTTCAGCAGGCGTTCAAGTTGTGCGTGCTCATCATCGTTTAACGGCGCATTGAGGCCAGCAAAATGGACATACTCGGCGTAAATATTGTGAACCGGGAGCCTGGCAGc                                                                                   >  1:766251/1‑140 (MQ=255)                                                gCGGGGCGTGGCTGGCGAGTGCCGGGCCATATTTCAGCAGGCGTTCAAGTTGTGCGTGCTCATCATCGTTTAACGGCGCATTGAGGCCAGCAAAATGGACATACTCGGCGTAAATATTGTGAACCGGGAGCCTGGCAGcc                                                                                  >  1:43719/1‑140 (MQ=255)                                                    ggCGTGGCTGGCGAGTGCCGGGCCATATTTCAGCAGGCGTTCAAGTTGTGCGTGCTCATCATCGTTTAACGGCGCATTGAGGCCAGCAAAATGGACATACTCGGCGTAAATATTGTGAACCGGGAGCCTGGCAGCCTGaa                                                                              <  1:894323/140‑1 (MQ=255)                                                    ggCGTGGCTGGCGAGAGCCGGGGCACATTGCAGCAGGCGTTCAAGTTGTGCGTGCTCATCATCGTTTAACGGCGCATTGAGGCCAGCAAAATGGACATACTCGGCGTAAATATTGTGAACCGGGAGCCTGGCAGCCTGaa                                                                              <  1:430901/140‑1 (MQ=255)                                                         ggctggcGAGTGCCGGGCCATATTTCAGCAGGCGTTCAAGTTGTGCGTGCTCATCATCGTTTAACGGCGCATTGAGGCCAGCAAAATGGACATACTCGGCGTTAATATTGTGAACCGGGAGCCTGGCAGCCTGAAAACGt                                                                         >  1:404801/1‑140 (MQ=255)                                                                             tatTTCAGCAGGCGTTCAAGTTGTGCGTGCTCATCATCGTTTAACGGCGCATTGAGGCCAGCAAAATGGACATACTCGGCGTAAATATTGTGAACCGGGAGCCTGGCAGCCTGAAAACGTGCCAGCAGTTTGTTGATTCg                                                     >  1:569354/1‑140 (MQ=255)                                                                                 tCAGCAGGCGTTCAAGTTGTGCGTGCTCATCATCGTTTAACGGCGCATTGAGGCCAGCAAAATGGACATACTCGGCGTAAATATTGTGAACCGGGAGCCTGGCAGCCTGAAAACGTGCCAGCAGTTTGTTGATTCGGAAt                                                 <  1:714499/140‑1 (MQ=255)                                                                                   agcagGCGTTCAAGTTGTGCGTGCTCATCATCGTTTAACGGCGCATTGAGGCCAGCAAAATGGACATACTCGGCGTAAATATTGTGAACCGGGAGCCTGGCAGCCTGAAAACGTGCCAGCAGTTTGTTGATTCGGAATGc                                               >  2:584844/1‑140 (MQ=255)                                                                                      agGCGTTCAAGTTGTACGTGCTGATCATCCTTTAACGGCGCATTGAGGCCAGCAAAATGGACATACTCGGCGTAAATATTGTGAACCGGGAGCCTGGCAGCCTGAAAACGTGCCAGCAGTTTGTTGATTCGGAATGCCGa                                            <  1:476510/140‑1 (MQ=255)                                                                                       ggCGTTCAAGTTGTGCGTGCTCATCATCGTTTAACGGCGCATTGAGGCCAGCAAAATGGACATACTCGGCGTAAATATTGTGAACCGGGAGCCTGGCAGCCTTAAAACGTACCATCAGTTTGTTCATTCGGAATGCCGAc                                           >  1:373166/1‑140 (MQ=255)                                                                                       ggCGTTCAAGTTGTGCGTGCTCATCATCGTTTAACGGCGCATTGAGGCCAGCAAAATGGACATACTCGGCGTAAATATTGTGAACCGGGAGCCTGGCAGCCTGAAAACGTGCCAGCAGTTTGTTGATTCGGAATGCCGAc                                           >  1:815226/1‑140 (MQ=255)                                                                                         cGTTCAAGTTGTGCGTGCTCATCATCGTTTAACGGCGCATTGAGGCCAGCAAAATGGACATACTCGGCGTAAATATTGTGAACCGGGAGCCTGGCAGCCTGACAACGTGCAAGCAGTTTGTTGAATCGGAATGCCGACAg                                         >  1:461112/1‑140 (MQ=255)                                                                                         cGTTCAAGTTGTGCGTGCTCATCATCGTTTAACGGCGCATTGAGGCCAGCAAAATGGACATACTCGGCGTAAATATTGTGAACCGGGAGCCTGGCAGCCTGAAAACGTGCCAGCAGTTTGTTGATTCGGAATGCCGACAg                                         >  2:625364/1‑140 (MQ=255)                                                                                          gTTCAAGTTGTGCGTGCTCATCATCGTTTAACGGCGCATTGAGGCCAGCAAAATGGACATACTCGGCGTAAATATTGTGAACCGGGAGCCTGGCAGCCTGAAAACGTGCCAGCAGTTTGTTGATTCGGAATGCCGACAGt                                        <  2:766251/140‑1 (MQ=255)                                                                                                  tgtgCGTGCTCATCATCGTTTAACGGCGCATTGAGGCCAGCAAAATGGACATACTCGGCGTAAATATTGTGAACCGGGAGCCTGGCAGCCTGAAAACGTGCCAGCAGTTTGTTGATTCGGAATGCCGACAGTGCAGGCGa                                >  1:968518/1‑140 (MQ=255)                                                                                                    tgcgtgcTCATCATCGTTTAACGGCGCATTGAGGCCAGCAAAATGGACATACTCGGCGTAAATATTGTGAACCGGGAGCCTGGCAGCCTGAAAACGTGCCAGCAGTTTGTTGATTCGGAATGCCGACAGTGCAGGCGAAc                              <  2:968518/140‑1 (MQ=255)                                                                                                    tgcgtgcTCATCATCGTTTAACGGCGCATTGAGGCCAGCAAAATGGACATACTCGGCGTAAATATTGTGAACCGGGAGCCTGGCAGCCTGAAAACGTGCCAGCAGTTTGTTGATTCGGAATGCCGACAGTGCAGGCGAAc                              >  2:1017815/1‑140 (MQ=255)                                                                                                      tgtgcTCATCATCGTTTAACGGCGCATTGAGGCCAGCAAAATGGACATACTCGGCGTAAATATTGTGAACCGGGAGCCTGGCAGCCTGAAAACGTGCCAGCAGTTTGTTGATTCGGAATGCCGACAGTGCAGGCGAACCa                            >  1:23341/2‑140 (MQ=255)                                                                                                      cgtgcTCATCATCGTTTAACGGCGCATTGAGGCCAGCAAAATGGACATACTCGGCGTAAATATTGTGAACCGGGAGCCTGGCAGCCTGAAAACGTGCCAGCAGTTTGTTGATTCGGAATGCCGACAGTGCAGGCGAACCa                            >  2:198300/1‑140 (MQ=255)                                                                                                      cgtgcTCATCATCGTTTAACGGCGCATTGAGGCCAGCAAAATGGACATACTCGGCGTAAATATTGTGAACCGGGAGCCTGGCAGCCTGAAAACGTGCCAGCAGTTTGTTGATTCGGAATGCCGACAGTGCAGGCGAACCa                            >  2:111146/1‑140 (MQ=255)                                                                                                          cTCATCATCGTTTAACGGCGCATTGAGGCCAGCAAAATGGACATACTCGGCGTAAATATTGTGAACCGGGAGCCTGGCAGCCTGAAAACGTGCCAGCAGTTTGTTGATTCGGAATGCCGACAGTGCAGGCGAACCACGCa                        >  1:427693/1‑140 (MQ=255)                                                                                                                  cGTTTAACGGCGCATTGGGGCCGGCAAAATGGACATACTAGGCGTAAATATTGTGAACCGGGAGCCTGGCAGCCTGAAAACGTGCCAGCAGTTTGTTGATTCGGAATGCCGACAGTGCAGGCGAACCACGCAGAATTTcc                <  2:404801/140‑1 (MQ=255)                                                                                                                         cGGCGCATTGAGGCCAGCAAAATGGACATACTCGGCGTAAATATTGTGAACCGGGAGCCTGGCAGCCTGAAAACGTGCCAGCAGTTTGTTGATTCGGAATGCCGACAGTGCAGGCGAACCACGCAGAATTTCCATCATaa         >  2:632077/1‑140 (MQ=255)                                                                                                                         cGGCGCATTGAGGCCAGCAAAATGGACATACTCGGCGTAAATATTGTGAACCGGGAGCCTGGCAGCCTGAAAACGTGCCAGCAGTTTGTTGATTCGGAATGCCGACAGTGCAGGCGAACCACGCAGAATTTCCATCATaa         >  1:796611/1‑140 (MQ=255)                                                                                                                            cgcATTGAGGCCAGCAAAATGGACATACTCGGCGTAAATATTGTGAACCGGGAGCCTGGCAGCCTGAAAACGTGCCAGCAGTTTGTTGATTCGGAATGCCGACAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGtc      >  2:714459/1‑140 (MQ=255)                                                                                                                            cgcATTGAGGCCAGCAAAATCGACATACTCGGCGTAAATATTGTGAACCGGGAGCCTGGCAGCCTGAAAACGTGCCAGCAGTTTGTTGATTCGGAATGCCGACACTACAGGCGAACCGCGCACAATTTCCATCATAAGtc      >  2:500167/1‑140 (MQ=255)                                                                                                                               aTTGAGGCCAGCAAAATGGACATACTCGGCGTAAATATTGTGAACCGGGAGCCTGGCAGCCTGAAAACGTGCCAGCAGTTTGTTGATTCGGAATGCCGACAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGtctct   <  1:215899/140‑1 (MQ=255)                                                                                                                                ttGAGGCCAGCAAAATGGACATACTCGGCGTAAATATTGTGAACCGGGAGCCTGGCAGCCTGAAAACGTGCCAGCAGTTTGTTGATTCGGAATGCCGACAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGtctctc  <  1:96831/140‑1 (MQ=255)                                                                                                                                      |                                                                                                                                      GGGAGAGATGGTGCCAGGACGCGGGGTCACCAGCAGGAGTTTGCCTTGCGGGGCGTGGCTGGCGAGTGCCGGGCCATATTTCAGCAGGCGTTCAAGTTGTGCGTGCTCATCATCGTTTAACGGCGCATTGAGGTCAGCAAAATGGACATACTCGGCGTAAATATTGTGAACCGGGAGCCTGGCAGCCTGAAAACGTGCCAGCAGTTTGTTGATTCGGAATGCCGACAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTC  >  NC_000913/2695286‑2695552

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