Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,931,785	T→C	L302P (CTG→CCG)	purT →	phosphoribosylglycinamide formyltransferase 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,931,785	0	T	C	100.0%	134.6 / NA	37	L302P (CTG→CCG)	purT	phosphoribosylglycinamide formyltransferase 2
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (18/19);  total (18/19)

GCACTGGGCGGTTATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCTCCCCTCGTCCACATGATACCGGGATGGTGACGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCTGCATGTACGTGCCTTCCTCGGACTTCCGGTTGGCGGGATCCGTCAGTATGGTCCTGCAGCTTCTGCCGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCG  >  NC_000913/1931649‑1931916                                                                                                                                         |                                                                                                                                    gCACTGGGCGGTTATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCTCCCCTCGTCCACATGATACCGGGATGGTGACGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCCGc                                                                                                                                   <  2:648098/139‑1 (MQ=255)                      tttGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCTCCCCTCGTCCACATGATACCGGGATGGTGACGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCCGCATGTACGTGCCTTCCTCGGAc                                                                                                              >  2:820203/1‑139 (MQ=255)                      tttGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCTCCCCTCGTCCACATGATACCGGGATGGTGACGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCCGCATGTACGTGCCTTCCTCGGAc                                                                                                              >  1:500341/1‑139 (MQ=255)                           tgtCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCTCCCCTCGTCCACATGATACCGGGATGGTGACGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCCGCATGTACGTGCCTTCCTCGGACTTCCg                                                                                                         >  1:109287/1‑139 (MQ=255)                            gtCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCTCCCCTCGTCCACATGATACCGGGATGGTGACGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCCGCATGTACGTGCCTTCCTCGGACTTCCgg                                                                                                        <  1:855366/139‑1 (MQ=255)                                aGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCTCCCCTCGTCCACATGATACCGGGATGGTGACGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCCGCATGTACGTGCCTTCCTCGGACTTCCGGTTgg                                                                                                    >  2:690288/1‑139 (MQ=255)                                     tttGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCTCCCCTCGTCCACATGATACCGGGATGGTGACGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCCGCATGTACGTGCCTTCCTCGGACTTCCGGTTGGCGGGa                                                                                               >  1:187955/1‑139 (MQ=255)                                     tttGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCTCCCCTCGTCCACATGATACCGGGATGGTGACGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCCGCATGTACGTGCCTTCCTCGGACTTCCGGTTGGCGGGa                                                                                               <  2:500341/139‑1 (MQ=255)                                     tttGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCTCCCCTCGTCCACATGATACCGGGATGGTGACGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCCGCATGTACGTGCCTTCCTCGGACTTCCGGTTGGCGGGa                                                                                               <  1:554929/139‑1 (MQ=255)                                          cTGTGGTGATGAGGTGATTTTCAGTGAGGTCTCCCCTCGTCCACATGATACCGGGATGGTGACGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCCGCATGTACGTGCCTTCCTCGGACTTCCGGTTGGCGGGATCCGt                                                                                          <  2:109287/139‑1 (MQ=255)                                               gtgATGAGGTGATTTTCAGTGAGGTCTCCCCTCGTCCACATGATACCGGGATGGTGACGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCCGCATGTACGTGCCTTCCTCGGACTTCCGGTTGGCGGGATCCGTCAGTa                                                                                     <  2:10861/139‑1 (MQ=255)                                                     aGGTGATTTTCAGTGAGGTCTCCCCTCGTCCACATGATACCGGGATGGTGACGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCCGCATGTACGTGCCTTCCTCGGACTTCCGGTTGGCGGGATCCGTCAGTATGGTcc                                                                               >  1:134757/1‑139 (MQ=255)                                                       gTGATTTTCAGTGAGGTCTCCCCTCGTCCACATGATACCGGGATGGTGACGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCCGCATGTACGTGCCTTCCTCGGACTTCCGGTTGGCGGGAt                                                                                              <  1:351635/122‑1 (MQ=255)                                                       gTGATTTTCAGTGAGGTCTCCCCTCGTCCACATGATACCGGGATGGTGACGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCCGCATGTACGTGCCTTCCTCGGACTTCCGGTTGGCGGGAt                                                                                              >  2:351635/1‑122 (MQ=255)                                                                  tGAGGTCTCCCCTCGTCCACATGATACCGGGATGGTGACGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCCGCATGTACGTGCCTTCCTCGGACtt                                                                                                            >  2:394712/1‑97 (MQ=255)                                                                  tGAGGTCTCCCCTCGTCCACATGATACCGGGATGGTGACGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCCGCATGTACGTGCCTTCCTCGGACtt                                                                                                            <  1:394712/97‑1 (MQ=255)                                                                   gAGGTCTCCCCTCGTCCACATGATACCGGGATGGTGACGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCCGCATGTACGTGCCTTCCTCGGACTTCCGGTTGGCGGGATCCGTCAGTATGGTCCTGCAGCTTCTGCCg                                                                 <  2:187955/139‑1 (MQ=255)                                                                                   cacaTGATACCGGGATGGTGACGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCCGCATGTACGTGCCTTCCTCGGACTTCCGGTTGGCGGGATCCGTCAGTATGGTCCTGCAGCTTCTGCCGTTATTCTGCCACAACt                                                 >  1:530070/1‑139 (MQ=255)                                                                                   cacaTGATACCGGGATGGTGACGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCCGCATGTACGTGCCTTCCTCGGACTTCCGGTTGGCGGGATCCGTCAGTATGGTCCTGCAGCTTCTGCCGTTATTCTGCCACAACt                                                 >  1:23961/1‑139 (MQ=255)                                                                                          tACCGGGATGGTGACGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCCGCATGTACGTGCCTTCCTCGGACTTCCGGTTGGCGGGATCCGTCAGTATGGTCCTGCAGCTTCTGCCGTTATTCTGCCACAACTGACCAGt                                          <  1:761731/139‑1 (MQ=255)                                                                                           aCCGGGATGGTGACGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCCGc                                                                                                                                   >  2:408683/1‑49 (MQ=255)                                                                                             cGGGATGGTGACGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCCGCATGTACGTGCCTTCCTCGGACTTCCGGTTGGCGGGATCCGTCAGTATGGTCCTGCAGCTTCTGCCGTTATTCTGCCACAACTGACCAGTCAg                                       >  2:483848/1‑139 (MQ=255)                                                                                                gATGGTGACGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCCGCATGTACGTGCCTTCCTCGGACTTCCGGTTGGCGGGATCCGt                                                                                          <  2:831043/85‑1 (MQ=255)                                                                                                gATGGTGACGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCCGCATGTACGTGCCTTCCTCGGACTTCCGGTTGGCGGGATCCGt                                                                                          >  1:831043/1‑85 (MQ=255)                                                                                                gATGGTGACGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCCGCATGTACGTGCCTTCCTCGGACTTCCGGTTGGCGGGATCCGTCAGTATGGTCCTGCAGCTTCTGCCGTTATTCTGCCACAACTGAc                                              <  2:244770/129‑1 (MQ=255)                                                                                                gATGGTGACGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCCGCATGTACGTGCCTTCCTCGGACTTCCGGTTGGCGGGATCCGTCAGTATGGTCCTGCAGCTTCTGCCGTTATTCTGCCACAACTGAc                                              >  1:244770/1‑129 (MQ=255)                                                                                                   ggTGACGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCCGCATGTACGTGCCTTCCTCGGACTTCCGGTTGGCGGGATc                                                                                             <  1:370203/79‑1 (MQ=255)                                                                                                   ggTGACGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCCGCATGTACGTGCCTTCCTCGGACTTCCGGTTGGCGGGATc                                                                                             >  2:370203/1‑79 (MQ=255)                                                                                                   ggTGACGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCCGCATGTACGTGCCTTCCTCGGACTTCCGGTTGGCGGGATCCGTCAGTATGGTCCTGCAGCTTCTGCCGTTATTCTGc                                                        >  2:531383/1‑116 (MQ=255)                                                                                                   ggTGACGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCCGCATGTACGTGCCTTCCTCGGACTTCCGGTTGGCGGGATCCGTCAGTATGGTCCTGCAGCTTCTGCCGTTATTCTGc                                                        <  1:531383/116‑1 (MQ=255)                                                                                                   ggTGACGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCCGCATGTACGTGCCTTCCTCGGACTTCCGGTTGGCGGGATCCGTCAGTATAGTCCTGCAGCTTCTGCCGTTATTCTGCCACAACTGACCAGTCAGAATGTc                                 <  1:820203/139‑1 (MQ=255)                                                                                                             aTTTCTCAAGATCTCTCAGAGTTTGCCCCGCATGTACGTGCCTTCCTCGGACTTCCGGTTGGCGGGATCCGTCAGTATGGTCCTGCAGCTTCTGCCGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATa                       <  2:530070/139‑1 (MQ=255)                                                                                                                  tcAAGATCTCTCAGAGTTTGCCCCGCATGTACGTGCCTTCCTCGGACTTCCGGTTGGCGGGATCCGTCAGTATGGTCCTGCAGCTTCTGCCGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAAtgtg                  <  1:483848/139‑1 (MQ=255)                                                                                                                   cAAGATCTCTCAGAGTTTGCCCCGCATGTACGTGCCTTCCTCGGACTTCCGGTTGGCGGGATCCGTCAGTATGGTCCTGCAGCTTCTGCCGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAAt                     <  1:294092/135‑1 (MQ=255)                                                                                                                   cAAGATCTCTCAGAGTTTGCCCCGCATGTACGTGCCTTCCTCGGACTTCCGGTTGGCGGGATCCGTCAGTATGGTCCTGCAGCTTCTGCCGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAAt                     >  2:294092/1‑135 (MQ=255)                                                                                                                           ctcAGAGTTTGCCCCGCATGTACGTGCCTTCCTCGGACTTCCGGTTGGCGGGATCCGTCAGTATGGTCCTGCAGCTTCTGCCGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGcc         >  2:870595/1‑139 (MQ=255)                                                                                                                                  tttGCCCCGCATGTACGTGCCTTCCTCGGACTTCCGGTTGGCGGGATCCGTCAGTATGGTCCTGCAGCTTCTGCCGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGgcg  <  1:690288/139‑1 (MQ=255)                                                                                                                                         |                                                                                                                                    GCACTGGGCGGTTATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCTCCCCTCGTCCACATGATACCGGGATGGTGACGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCTGCATGTACGTGCCTTCCTCGGACTTCCGGTTGGCGGGATCCGTCAGTATGGTCCTGCAGCTTCTGCCGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCG  >  NC_000913/1931649‑1931916

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