Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,204,856	A→G	V359A (GTG→GCG)	purD ←	phosphoribosylamine‑‑glycine ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,204,856	0	A	G	100.0%	85.1 / NA	24	V359A (GTG→GCG)	purD	phosphoribosylamine‑‑glycine ligase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (13/11);  total (13/11)

TAGGCGCGTTTCTGCGCTTCTGCCACGGTATGACCCAGCGCGGTGACGCACAGTACGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGG  >  NC_000913/4204727‑4204991                                                                                                                                  |                                                                                                                                        tAGGCGCGTTTCTGCGCTTCTGCCACGGTGTGACCCAGCGCGGTGACGCACAGTACGCACCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCGCTTCTTCCAg                                                                                                                               >  2:29388/1‑140 (MQ=255)   ggCGCGTTTCTGCGCTTCTGCCACGGTATGACCCAGCGCGGTGACGCACAGTACGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCGCTTCTTCCAGCg                                                                                                                             >  2:10452/1‑140 (MQ=255)    gcgcGTTTCTGCGCTTCTGCCACGGTATGACCCAGCGCGGTGACGCACAGTACGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCGCTTCTTCCAGCgg                                                                                                                            >  2:99215/1‑140 (MQ=255)      gggTTTTTGCGCTTCTGCCACGGTATGCCCCGGCGGGGGGACGCGCAGTAGCCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCGCTTCTTCCAGCGGCa                                                                                                                          <  2:68135/138‑1 (MQ=255)                    cTGCCACGGTATGACCCAGCGCGGTGACGCACAGTACGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCGCTTCTTCCAGCGGCAGGCCGTGGATCACg                                                                                                            >  1:354582/1‑140 (MQ=255)                    cTGCCACGGTATGACCCAGCGCGGTGACGCACAGTACGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCGCTTCTTCCAGCGGCAGGCCGTGGATCACg                                                                                                            >  1:415378/1‑140 (MQ=255)                              aTGACCCAGCGCGGTGACGCACAGTACGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCGCTTCTTCCAGCGGCAGGCCGTGGATCACGTCACCGGTGc                                                                                                  >  1:202709/1‑140 (MQ=255)                              aTGACCCAGCGCGGTGACGCACAGTACGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCGCTTCTTCCAGCGGCAGGCCGTGGATCACGTCACCGGTGc                                                                                                  >  1:370537/1‑140 (MQ=255)                                      gcggGGTGCCGCACAGTCCGCCCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCGCTTCTTCCAGCGGCAGGCCGTGGATCACGTCACCGGTGCGGTAATCa                                                                                          <  2:202709/140‑1 (MQ=255)                                       cgcgGTGACGCACAGTACGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCATGCCGCTTCTTCCAGCGGCAGGCCGTGGATCACGTCACCGGTGCGGTAATCAc                                                                                         >  1:116827/1‑140 (MQ=255)                                                   caGTACGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCGCTTCTTCCAGCGGCAGGCCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCAc                                                                             <  2:354582/140‑1 (MQ=255)                                                     gTACGCTCCCGCCCTTGGTCCTTACCTCCCCGTCATCCGCCAGTTTTGCGCCCGCGTGGAACACTTTCCCGCCTGCCGCTTCTTCCAGCGGCAGGTCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCAccg                                                                           <  2:385110/140‑2 (MQ=255)                                                     gTACGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCGCTTCTTCCAGCGGCAGGCCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCAccc                                                                           <  2:415378/140‑1 (MQ=255)                                                      tACGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCATGCCGCTTCTTCCAGCGGCAGGCCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCg                                                                          <  2:116827/140‑1 (MQ=255)                                                        cgcgCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCGCTTCTTCCAGCGGCAGGCCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCg                                                                          >  2:105194/1‑138 (MQ=255)                                                        cgcgCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCGCTTCTTCCAGCGGCAGGCCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCg                                                                          <  1:105194/138‑1 (MQ=255)                                                          cgcccgccGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCGCTTCTTCCAGCGGCAGGCCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGc                                                                      <  2:370537/140‑1 (MQ=255)                                                                                                ttttGTGCCCGCGTGGAACACTTTGCCGCCTGCCGCTTCTTCCAGCGGCAGACCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCgg                                >  2:73894/1‑140 (MQ=255)                                                                                                         cgcgTGGAACACTTTGCCGCCTGCCGCTTCTTCCAGCGGCAGACCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTtct                       <  1:73894/140‑1 (MQ=255)                                                                                                                    cTTTGCCGCCTGCCGCTTCTTCCAGCGGCAGGCCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTg            >  1:349543/1‑140 (MQ=255)                                                                                                                             cTGCCGCTTCTTCCAGCGGCAGGCCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAg   >  2:320912/1‑140 (MQ=255)                                                                                                                              tGCCGCTTCTTCCAGCGGCAGGCCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTg            >  1:40347/1‑130 (MQ=255)                                                                                                                              tGCCGCTTCTTCCAGCGGCAGGCCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTg            <  2:40347/130‑1 (MQ=255)                                                                                                                              tGCCGCTTCTTCCAGCGGCAGGCCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAgg  <  2:118839/140‑1 (MQ=255)                                                                                                                              tGCCGCTTCTTCCAGCGGCAGGCCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAgg  <  1:320912/140‑1 (MQ=255)                                                                                                                                  |                                                                                                                                        TAGGCGCGTTTCTGCGCTTCTGCCACGGTATGACCCAGCGCGGTGACGCACAGTACGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGG  >  NC_000913/4204727‑4204991

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: ‑