Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,204,874	C→T	G353D (GGC→GAC)	purD ←	phosphoribosylamine‑‑glycine ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,204,874	0	C	T	100.0%	154.7 / NA	48	G353D (GGC→GAC)	purD	phosphoribosylamine‑‑glycine ligase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (20/28);  total (20/28)

GCTTCTGCCACGGTATGACCCAGCGCGGTGACGCACAGTACGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAA  >  NC_000913/4204742‑4205009                                                                                                                                     |                                                                                                                                        ctttCTGCCCGGGTATGCCCCGGCGCGGTGCCGCACAGTCCGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGTCGTGGAt                                                                                                                                  <  2:164077/138‑1 (MQ=255)     cTGCCACGGTATGACCCACCGCGGTGACGCACAGTACGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGTCGTGGATCACg                                                                                                                              >  1:96674/1‑140 (MQ=255)       gCCACGGTATGACCCAGCGCGGTGACGCACAGTACGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGTCGTGGATCACGTc                                                                                                                            >  2:330002/1‑140 (MQ=255)         cagggTATGCCCCAGCGGGGTGCCGCACAGTACGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGTCGTGGATCACGTCAc                                                                                                                          <  1:486020/137‑1 (MQ=255)            ggTATGCCCCAGCGCGGTGCCGCACAGTCCGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCCCTTCTTCCAGCGGCAGGTCGTGGATCACGTCACCgg                                                                                                                       <  2:341174/140‑1 (MQ=255)                 gccccAGCGGGGTGACGCACAGTACGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGTCGTGGATCACGTCACCGGTGCgg                                                                                                                  <  1:84382/138‑1 (MQ=255)                      aGCGCGGTGACGCACAGTACGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGTCGTGGATCACGTCACCGGTGCGGTAATc                                                                                                             >  2:181150/1‑140 (MQ=255)                        cgcgGTGACGCACAGTACGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGTCGTGGATCACGTCACCGGTGCGGTAATCAc                                                                                                           <  1:161207/140‑1 (MQ=255)                        cgcgGTGACGCACAGTACGCGCCCGCCGTTGGTCACTACCTGCTCGTCACCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGTCGTGGATCACGTCACCGGTGCGGTAATCAc                                                                                                           >  1:376943/1‑140 (MQ=255)                          cgGTGACGCACAGTACGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGTCGTGGATCACGTCACCGGTGCGGTAATCAccc                                                                                                         >  2:65020/1‑140 (MQ=255)                                       tccgcgCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGTCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCg                                                                                            <  1:315313/138‑1 (MQ=255)                                       tccgcgCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGTCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCg                                                                                            <  1:436918/138‑1 (MQ=255)                                       tACGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGTCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCg                                                                                            >  2:298219/1‑140 (MQ=255)                                               cgccgTTGGTCACTTCCTGCTCGTAATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGTCGTGGATCACGTCGCCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATc                                                                                    <  2:96674/140‑1 (MQ=255)                                               cgccgTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGTCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATc                                                                                    >  2:1370/1‑140 (MQ=255)                                                 ccgTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGTCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCCCAGCCATcac                                                                                  >  2:521795/1‑140 (MQ=255)                                                 ccgTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCC‑CGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGTCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATcacc                                                                                 >  1:249920/1‑140 (MQ=255)                                                  cgTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGTCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATcacc                                                                                 <  1:298219/140‑1 (MQ=255)                                                  cgTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGTCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATcacc                                                                                 <  1:382392/140‑1 (MQ=255)                                                   gTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCC‑CGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGTCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATcaccac                                                                               <  2:249920/140‑1 (MQ=255)                                                            tACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGTCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCgagag                                                                       >  1:387475/1‑140 (MQ=255)                                                            tACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGTCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCgagag                                                                       >  1:405231/1‑140 (MQ=255)                                                                 gCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGTCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCga                                                                          <  1:506665/132‑1 (MQ=255)                                                                 gCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGTCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCga                                                                          >  2:506665/1‑132 (MQ=255)                                                                   tCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGTCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGcgcg                                                                <  2:533704/140‑1 (MQ=255)                                                                      tCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGTCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTc                                                             >  1:2952/1‑140 (MQ=255)                                                                      tCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGTCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTc                                                             >  2:6449/1‑140 (MQ=255)                                                                      tCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGTCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTc                                                             >  2:188077/1‑140 (MQ=255)                                                                      tCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGTCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTc                                                             >  1:158620/1‑140 (MQ=255)                                                                          ccGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGTCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATcc                                                         <  2:69382/140‑1 (MQ=255)                                                                                  tttGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGTCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGa                                                 <  1:181150/140‑1 (MQ=255)                                                                                      tgCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGTCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGtt                                             >  1:101158/1‑140 (MQ=255)                                                                                       gCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGTCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGttt                                            <  2:310248/140‑1 (MQ=255)                                                                                                 aaCACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGTCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAg                                  <  1:1370/140‑1 (MQ=255)                                                                                                           cgccTGCCACTTCTTCCAGCGGCAGGTCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTca                        >  1:405826/1‑140 (MQ=255)                                                                                                            gccTGCCACTTCTTCCAGCGGCAGGTCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTcac                       <  1:6449/140‑1 (MQ=255)                                                                                                            gccTGCCACTTCTTCCAGCGGCAGGTCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTcac                       <  1:188077/140‑1 (MQ=255)                                                                                                             ccTGCCACTTCTTCCAGCGGCAGGTCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTcaca                      <  2:405826/140‑1 (MQ=255)                                                                                                               tGCCACTTCTTCCAGCGGCAGGTCGTGGATCACGTCACCGGTGCGGTAATCCCCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAgg                    <  1:219515/140‑1 (MQ=255)                                                                                                               tGCCACTTCTTCCAGCGGCAGGTCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAgg                    <  2:158620/140‑1 (MQ=255)                                                                                                                gCCACTTCTTCCAGCGGCAGGTCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAg                                  >  2:274857/1‑125 (MQ=255)                                                                                                                gCCACTTCTTCCAGCGGCAGGTCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAg                                  <  1:274857/125‑1 (MQ=255)                                                                                                                       cttcCAGCGGCAGGTCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAgg            <  1:65020/140‑1 (MQ=255)                                                                                                                       cttcCAGCGGCAGGTCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAgg            <  1:268351/140‑1 (MQ=255)                                                                                                                        ttcCAGCGGCAGGTCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGc           <  2:352638/140‑1 (MQ=255)                                                                                                                          ccAGCGGCAGGTCGTGGATAACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCag         <  2:413937/140‑1 (MQ=255)                                                                                                                            aGCGGCAGGTCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCagag       <  2:376943/140‑1 (MQ=255)                                                                                                                              cGGCAGGTCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCagag       >  1:93195/1‑138 (MQ=255)                                                                                                                              cGGCAGGTCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCagag       <  2:93195/138‑1 (MQ=255)                                                                                                                                 cAGGTCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCaa  >  1:537316/1‑140 (MQ=255)                                                                                                                                 cAGGTCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCaa  <  2:538958/140‑1 (MQ=255)                                                                                                                                     |                                                                                                                                        GCTTCTGCCACGGTATGACCCAGCGCGGTGACGCACAGTACGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAA  >  NC_000913/4204742‑4205009

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