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%	106.6 / NA	30	L364R (CTG→CGG)	purT	phosphoribosylglycinamide formyltransferase 2
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (13/17);  total (13/17)

CTGCAGCTTCTGCCGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCT‑GGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCC  >  NC_000913/1931839‑1932088                                                                                                                                      |                                                                                                                      cTGCAGCTTCTGCCGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTGTATTTGGTAAGCCGGAAATTGGTGGAAGCCGTCGTCT‑gggggg                                                                                                                 >  2:944117/1‑138 (MQ=255)      gCTTCTGCCGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCa                                                                                                             <  2:1133130/139‑1 (MQ=255)        ttCTGCCGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACt                                                                                                           <  2:32681/139‑1 (MQ=255)        ttCTGCCGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACt                                                                                                           <  1:865724/139‑1 (MQ=255)        ttCTGCCGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCT‑ggggggtggtact                                                                                                          >  1:238919/1‑131 (MQ=255)         tCTGCCGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTg                                                                                                          >  1:1287407/1‑139 (MQ=255)                                  ccAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACg                                                                                 >  1:292171/1‑139 (MQ=255)                                  ccAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCT‑gggggg                                                                                                                 <  1:270615/106‑2 (MQ=255)                                  ccAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCT‑gggggg                                                                                                                 >  2:270615/1‑105 (MQ=255)                                      tCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCT‑gggggg                                                                                                                 >  2:834856/1‑101 (MQ=255)                                      tCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCT‑gggggg                                                                                                                 <  1:834856/102‑2 (MQ=255)                                                aCGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAAcgcgcg                                                                   >  1:90346/1‑139 (MQ=255)                                                aCGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCT‑gggggg                                                                                                                 <  2:667043/92‑2 (MQ=255)                                                aCGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCT‑gggggg                                                                                                                 >  1:667043/1‑91 (MQ=255)                                                     tGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCa                                                              <  1:340430/139‑1 (MQ=255)                                                      gATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTACTGCAGAAAGTGTTGTTGACGCCATTGAACGCGCGAAGCa                                                              <  2:768186/138‑1 (MQ=255)                                                      gATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATCGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCTCGAAGCa                                                              >  1:768186/1‑138 (MQ=255)                                                            gtgCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCa                                                                                                             <  1:387782/85‑1 (MQ=255)                                                            gtgCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCa                                                                                                             >  2:387782/1‑85 (MQ=255)                                                            gtgCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGaa                                                                         >  2:307429/1‑121 (MQ=255)                                                            gtgCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGaa                                                                         <  1:307429/121‑1 (MQ=255)                                                                      ccGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCTGTCGTCTGGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTaa                                             <  2:184376/139‑1 (MQ=255)                                                                       cGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTaaa                                            >  2:881139/1‑139 (MQ=255)                                                                         tAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTAc                                                                                                     <  1:649346/80‑1 (MQ=255)                                                                         tAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTAc                                                                                                     >  2:649346/1‑80 (MQ=255)                                                                           ggCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGTACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTa                                        <  2:238919/139‑1 (MQ=255)                                                                                            tCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCTGTCGTCTGGGGGGTGGCa                                                                                                             <  2:553587/53‑1 (MQ=255)                                                                                            tCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCTGTCGTCTGGGGGGTGGCa                                                                                                             >  1:553587/1‑53 (MQ=255)                                                                                            tCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCt                                                                                                                        <  2:361321/42‑1 (MQ=255)                                                                                            tCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCt                                                                                                                        >  1:361321/1‑42 (MQ=255)                                                                                                ttATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCaaaaa                   <  2:90346/139‑1 (MQ=255)                                                                                                 tATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGTCaaaaaa                  >  1:80217/1‑139 (MQ=255)                                                                                                      ggTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCaaaaaaa                 >  1:1295732/1‑135 (MQ=255)                                                                                                      ggTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCaaaaaaa                 <  2:1295732/135‑1 (MQ=255)                                                                                                                 aaTTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGTCAAAAAAACGCTTCAAAAATGcc  <  2:80217/139‑1 (MQ=255)                                                                                                                 aaTTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGcc  <  1:881139/139‑1 (MQ=255)                                                                                                                          cAGCCGTCGTCTGGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGa                                                    <  1:37506/80‑1 (MQ=255)                                                                                                                          cAGCCGTCGTCTGGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGa                                                    >  2:37506/1‑80 (MQ=255)                                                                                                                              cgtcgtCTGGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACa                                                  >  1:1211200/1‑78 (MQ=255)                                                                                                                              cgtcgtCTGGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACa                                                  <  2:1211200/78‑1 (MQ=255)                                                                                                                                      |                                                                                                                      CTGCAGCTTCTGCCGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCT‑GGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCC  >  NC_000913/1931839‑1932088

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