Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,429,936	G→T	H101Q (CAC→CAA)	purF ←	amidophosphoribosyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,429,936	0	G	T	100.0%	102.0 / NA	30	H101Q (CAC→CAA)	purF	amidophosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (14/16);  total (14/16)

CCAGCGGGTAGTGGCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTCTTCAAACAGTTTTTTACGCAACTCGTGAGCGTTGGTCAGATTGCCGTTGTGGGCAAGCGTAATGCCATACGGGGAGTTAACGTAAAACGGCTGCGCTTCAGAGGCGCTGGAGCTGCCAGCCGTGGGGTAACGCACATGACCAATGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAA  >  NC_000913/2429794‑2430073                                                                                                                                               |                                                                                                                                          ccAGCGGGTAGTGGCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTCTTCAAACAGTTTTTTACGCAACTCGTGAGCGTTGGTCAGATTGCCGTTTTGGGc                                                                                                                                      >  7:47195/1‑148 (MQ=255)   aGCGGGTAGTGGCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTCTTCAAACAGTTTTTTACGCAACTCGTGAGCGTTGGTCAGATTGCCGTTTTGGGCAAg                                                                                                                                   >  1:57059/1‑149 (MQ=255)              gCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGATTATCTTCAAAAAGTTTTTTACGCAACTCGTGAGCGTTGGTCAGATTGCCGTTTTGGGCAAGCGTAATGCCa                                                                                                                         >  4:52519/1‑148 (MQ=255)                ggAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTCTTCAAACAGTTTTTTACGCAACTCGTGAGCGTTGGTCAGATTGCCGTTTTGGGCAAGCGTAATGCCATAc                                                                                                                      >  8:62358/1‑149 (MQ=255)                    gTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTCTTCAAACAGTTTTTTACGCAACTCGTGAGCGTTGGTCAGATTGCCGTTTTGGGCAAGCGTAATGCCATACgggg                                                                                                                  >  7:84616/1‑149 (MQ=255)                       gTCCAGCGCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGGTGATGTGGCGGCGTTTTTCTTCAAACAGTTTTTTACGCAACTCGTGAGCGTTGGTCAGATTGCCGTTTTGGGCAAGCGTAATGCCATACGGGGAg                                                                                                                <  3:1119/148‑1 (MQ=255)                          cAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTCTTCAAACAGTTTTTTACGCAACTCGTGAGCGTTGGTCAGATTGCCGTTTTGGGCAAGCGTAATGCCATACGGGGAGTTa                                                                                                             <  2:80555/148‑1 (MQ=255)                                    gCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTCTTCAAACAGTTTTTTACGCAACTCGTGAGCGTTGGTCAGATTGCCGTTTTGGGCAAGCGTAATGCCATACGGGGAGTTAACGTAAAACgg                                                                                                  >  3:98488/1‑149 (MQ=255)                                      gAAGATATTAAGCAGAATTTCCGAGTCGGATGTGGGGTTGATGTGGCGGCGTTTTTCTTCAACCAGTTGTTTACGCCACTCGTGAGCGTTGGTCAGATTGCCGTTTTGGGCACGCGTAATGCCCGACGGGGAGGGAACGTAAAACGGCt                                                                                                >  7:57638/1‑149 (MQ=255)                                      gAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTCTTCAAACAGTTTTTTACGCAACTCGTGAGCGTTGGTCAGATTGCCGTTTTGGGCAAGCGTAATGCCATACGGGGAGTTAACGTAAAACGGCt                                                                                                >  6:70963/1‑149 (MQ=255)                                       aaGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTCTTCAAACAGTTTTTTACGCAACTCGTGAGCGTTGGTCAGATTGCCGTTTTGGGCAAGCGTAATGCCATACGGGGAGTTAACGTAAAACGGCTg                                                                                               >  2:60990/1‑149 (MQ=255)                                          atatTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTCTTCAAACAGTTTTTTACGCAACTCGTGAGCGTTGGTCAGATTGCCGTTTTGGGCAAGCGTAATGCCATACGGGGAGTTAACGTAAAACGGCTgcgc                                                                                            <  8:84616/149‑1 (MQ=255)                                                 gCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTCTTCAAACAGTTTTTTACGCAACTCGTGAGCGTTGGTCAGATTGCCGTTTTGGGCAAGCGTAATGCCATACGGGGAGTTAAAGTAAAACGGCTGCGCTTCagag                                                                                     >  5:98786/1‑149 (MQ=255)                                                           cGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTCTTCAAACAGTTTTTTACGCAACTCGTGAGCGTTGGTCAGATTGCCGTTTTGGGCAAGCGTAATGCCATACGGGGAGTTAACGTAAAACGGCTGCGCTTCAGAGGCGCTGGAg                                                                            <  5:116339/148‑1 (MQ=255)                                                           cGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTCTTCAAACAGTTTTTTACGCAACTCGTGAGCGTTGGTCAGATTGCCGTTTTGGGCAAGCGTAATGCCATACGGGGAGTTAACGTAAAACGGCTGCGCTTCAGAGGCGCTGGAGc                                                                           >  2:109358/1‑149 (MQ=255)                                                            gAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTCTTCAAACAGTTTTTTACGCAACTCGTGAGCGTTGGTCAGATTGCCGTTTTGGGCAAGCGTAATGCCATACGGGGAGTTAACGTAAAACGGCTGCGCTTCAGAGGCGCTGGAGCt                                                                          <  4:98488/149‑1 (MQ=255)                                                                  gAAGTGGTGTTGATGTGGCGGCGTTTTTCTTCAAACAGTTTTTTACGCAACTCGTGAGCGTTGGTCAGATTGCCGTTTTGGGCAAGCGTAATGCCATACGGGGAGTTAACGTAAAACGGCTGCGCTTCAGAGGCGCTGGAGCTGCCAGc                                                                    >  7:73401/1‑149 (MQ=255)                                                                      tggtgTTGATGTGGCGGCGTTTTTCTTCAAACAGTTTTTTACGCAACTCGTGAGCGTTGGTCAGATTGCCGTTTTGGGCAAGCGTAATGCCATACGGGGAGTTAACGTAAAACGGCTGCGCTTCAGAGGCGCTGGAGCTGCCAGCCGTg                                                                <  5:70963/149‑1 (MQ=255)                                                                        gtgtTGATGTGGCGGCGTTTTTCTTCAAACAGTTTTTTACGCAACTCGTGAGCGTTGGTCAGATTGCCGTTTTGGGCAAGCGTAATGCCATACGGGGAGTTAACGTAAAACGGCTGCGCTTCAGAGGCGCTGGAGCTGCCAGCCGTggg                                                              <  4:630/149‑1 (MQ=255)                                                                                     ggcgTTTTTCTTCAAACAGTTTTTTACGCAACTCGTGAGCGTTGGTCAGATTGCCGTTTTGGGCAAGCGTAATGCCATACGGGGAGTTAACGTAAAACGGCTGCGCTTCAGAGGCGCTGGAGCTGCCAGCCGTGGGGTAACGCACAtaa                                                 <  7:18281/149‑3 (MQ=255)                                                                                           tttCTTCAAACAGTTTTTTACGCAACTCGTGAGCGTTGGTCAGATTGCCGTTTTGGGCAAGCGTAATGCCATACGGGGAGTTAACGTAAAACGGCTGCGCTTCAGAGGCGCTGGAGCTGCCAGCCGTGGGGTAACGCACATGACCAAt                                            <  8:47195/148‑1 (MQ=255)                                                                                                   aaCAGTTTTTTACGCAACTCGTGAGCGTTGGTCAGATTGCCGTTTTGGGCAAGCGTAATGCCATACGGGGAGTTAACGTAAAACgg                                                                                                  >  5:67701/1‑86 (MQ=255)                                                                                                   aaCAGTTTTTTACGCAACTCGTGAGCGTTGGTCAGATTGCCGTTTTGGGCAAGCGTAATGCCATACGGGGAGTTAACGTAAAACgg                                                                                                  <  6:67701/86‑1 (MQ=255)                                                                                                      aGTTTTTTACGCAACTCTGGCGCGTTGGTCAGATTGCCGTTTTGGGCAAGCGTAATGCCATACGGGGAGTTAACGTAAAACGGCTGCGCTTCAGAGGCGCTGGAGCTGCCAGCCGTGGGGTAACGCACATGACCAATGCCCATATTGcc                                <  8:57638/149‑1 (MQ=255)                                                                                                      aGTTTTTTACGCAACTCGTGAGCGTTGGTCAGATTGCCGTTTTGGGCAAGCGTAATGCCATACGGGGAGTTAACGTAAAACGGCTGCGCTTCAGAGGCGCTGGAGCTGCCAGCCGTGGGGTAACGCACATGACCAATGCCCATATTGcc                                <  1:60990/149‑1 (MQ=255)                                                                                                           tttACGCAACTCGTGAGCGTTGGTCAGATTGCCGTTTTGGGCAAGCGTAATGCCATACGGGGAGTTAACGTAAAACGGCTGCGCTTCAGAGGCGCTGGAGCTGCCAGCCGTGGGGTAACGCACATGACCAATGCCCATATTGCCCTGc                            <  2:57059/148‑1 (MQ=255)                                                                                                                  aaCTCGTGAGCGTTGGTCAGATTGCCGTTTTGGGCAAGCGTAATGCCATACGGGGAGTTAACGTAAAACGGCTGCGCTTCAGAGGCGCTGGAGCTGCCAGCCGTGGGGTAACGCACATGACCAATGCCCATATTGCCCTGCAAACGCTg                    <  8:73401/149‑1 (MQ=255)                                                                                                                     tCGTGAGCGTTGGTCAGATTGCCGTTTTGGGCAAGCGTAATGCCATACGGGGAGTTAACGTAAAACGGCTGCGCTTCAGAGGCGCTGGAGCTGCCAGCCGTGGGGTAACGCACATGACCAATGCCCATATTGCCCTGCAAACGCTGCa                  >  4:104249/1‑148 (MQ=255)                                                                                                                                  tCAGATTGCCGTTTTGGGCAAGCGTAATGCCATACGGGGAGTTAACGTAAAACGGCTGCGCTTCAGAGGCGCTGGAGCTGCCAGCCGTGGGGTAACGCACATGACCAATGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCa    <  6:104368/149‑1 (MQ=255)                                                                                                                                    aGATTGCCGTTTTGGGCAAGCGTAATGCCATACGGGGAGTTAACGTAAAACGGCTGCGCTTCAGAGGCGCTGGAGCTGCCAGCCGTGGGGTAACGCACATGACCAATGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCaaa  <  7:62358/149‑1 (MQ=255)                                                                                                                                               |                                                                                                                                          CCAGCGGGTAGTGGCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTCTTCAAACAGTTTTTTACGCAACTCGTGAGCGTTGGTCAGATTGCCGTTGTGGGCAAGCGTAATGCCATACGGGGAGTTAACGTAAAACGGCTGCGCTTCAGAGGCGCTGGAGCTGCCAGCCGTGGGGTAACGCACATGACCAATGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAA  >  NC_000913/2429794‑2430073

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