Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,931,976	(G)5→6	coding (1096/1179 nt)	purT →	phosphoribosylglycinamide formyltransferase 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,931,971	1	.	G	100.0%	147.6 / NA	41	L364R (CTG→CGG)	purT	phosphoribosylglycinamide formyltransferase 2
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (20/21);  total (20/21)

CAGCTTCTGCCGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCT‑GGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCC  >  NC_000913/1931842‑1932088                                                                                                                                   |                                                                                                                      cAGCTTCTGCCGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCT‑ggggggtgg                                                                                                              >  2:156513/1‑135 (MQ=255)      tCTGCCGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCT‑gggggg                                                                                                                 >  1:1037275/1‑130 (MQ=255)      tCTGCCGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCT‑gggggg                                                                                                                 <  2:1037275/131‑2 (MQ=255)                          aCTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCa                                                                                                             <  1:1360760/116‑1 (MQ=255)                          aCTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCa                                                                                                             >  2:1360760/1‑116 (MQ=255)                             gACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGa                                                                                   >  2:358614/1‑139 (MQ=255)                                   tCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCAt                                                                             >  1:1131011/1‑139 (MQ=255)                                       aaTGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGaa                                                                         <  2:1131011/139‑1 (MQ=255)                                            cACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAAcgcgc                                                                    >  2:720226/1‑139 (MQ=255)                                             aCGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAAcgcgcg                                                                   >  2:1450984/1‑139 (MQ=255)                                                  tGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCa                                                              <  1:358614/139‑1 (MQ=255)                                                  tGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCa                                                              <  1:720226/139‑1 (MQ=255)                                                       aTGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCAcgccg                                                         <  1:1021062/139‑1 (MQ=255)                                                         gtgCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGa                                                                          <  1:1026800/120‑1 (MQ=255)                                                         gtgCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGa                                                                          >  2:1026800/1‑120 (MQ=255)                                                            cAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCa                                                              <  2:288359/129‑1 (MQ=255)                                                            cAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCa                                                              >  1:288359/1‑129 (MQ=255)                                                             aGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCa                                                              >  1:77889/1‑128 (MQ=255)                                                             aGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCa                                                              <  2:77889/128‑1 (MQ=255)                                                                  gCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTACt                                                                                                    >  1:1403935/1‑85 (MQ=255)                                                                  gCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTACt                                                                                                    <  2:1403935/85‑1 (MQ=255)                                                                  gCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGAc                                                                                  >  2:231389/1‑103 (MQ=255)                                                                  gCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGAc                                                                                  <  1:231389/103‑1 (MQ=255)                                                                   ccGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTaa                                             <  1:156513/139‑1 (MQ=255)                                                                    cGTAGGTGCAGATTTGCAGATTCGTTTATTTGGTAAGCC‑GAAATTGATGGCAGCCGTCGTCT‑g                                                                                                                      <  1:292561/63‑1 (MQ=255)                                                                    cGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTaaa                                            >  1:860334/1‑139 (MQ=255)                                                                    cGTAGGCGCAGATTTGCAGATTCGTTTATTTAGTAAGCC‑GAAATTGTTGGCAGCCGTCGTCT‑g                                                                                                                      >  2:292561/1‑63 (MQ=255)                                                                      tAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCgg                                                     >  2:491556/1‑128 (MQ=255)                                                                      tAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCgg                                                     <  1:491556/128‑1 (MQ=255)                                                                            cAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTAc                                       >  2:689141/1‑136 (MQ=255)                                                                            cAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTAc                                       <  1:689141/136‑1 (MQ=255)                                                                             aGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGCCTGGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAggg                                   >  2:856360/1‑139 (MQ=255)                                                                                                 ttGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGc                                                                                                        >  2:964585/1‑50 (MQ=255)                                                                                                 ttGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGc                                                                                                        <  1:964585/50‑1 (MQ=255)                                                                                                     tAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTACTGCa                                                                                                 >  2:435393/1‑53 (MQ=255)                                                                                                     tAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTACTGCa                                                                                                 <  1:435393/53‑1 (MQ=255)                                                                                                      aaGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAAcgcgcg                                                                   <  2:58522/82‑1 (MQ=255)                                                                                                      aaGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAAcgcgcg                                                                   >  1:58522/1‑82 (MQ=255)                                                                                                         ccGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTAc                                                                                                     >  2:346789/1‑45 (MQ=255)                                                                                                         ccGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTAc                                                                                                     <  1:346789/45‑1 (MQ=255)                                                                                                          cGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCaaaaaaa                 >  2:299823/1‑128 (MQ=255)                                                                                                          cGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCaaaaaaa                 <  1:299823/128‑1 (MQ=255)                                                                                                             aaaTTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCa                                                              >  1:908768/1‑80 (MQ=255)                                                                                                             aaaTTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCa                                                              <  2:908768/80‑1 (MQ=255)                                                                                                              aaTTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGcc  <  1:1450984/139‑1 (MQ=255)                                                                                                              aaTTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGcc  <  2:860334/139‑1 (MQ=255)                                                                                                                                   |                                                                                                                      CAGCTTCTGCCGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCT‑GGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCC  >  NC_000913/1931842‑1932088

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