Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,931,125	A→G	E82G (GAG→GGG)	purT →	phosphoribosylglycinamide formyltransferase 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,931,125	0	A	G	100.0%	146.7 / NA	41	E82G (GAG→GGG)	purT	phosphoribosylglycinamide formyltransferase 2
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (22/19);  total (22/19)

AGGTGATTGCCGTCGATCGCTATGCCGACGCACCAGCCATGCATGTCGCGCATCGCTCCCATGTCATTAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATAT‑CGTGCCGGAGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAACGAAATTAACGATGAATCGCGAGGGTATCCGTCGCCTGGCGGCAGAAGAGCTGCAG  >  NC_000913/1930990‑1931258                                                                                                                                         |                                                                                                                                      agGTGATTGCCGTCGATCGCTATGCCGACGCACCAGCCATGCATGTCGCGCATCGCTCCCATGTCATTAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATAT‑CGTGCCGGGGAt                                                                                                                                    >  2:406742/1‑139 (MQ=255)         gCCGTCGATCGCTATGCCGACGCACCAGCCATGCATGTCGCGCATCGCTCCCATGTCATTAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATAT‑CGTGCCGGGGATCGAAGCTa                                                                                                                            <  2:26071/139‑1 (MQ=255)           cGTCGATCGCTATGCCGACGCACCAGCCATGCATGTCGCGCATCGCTCCCATGTCATTAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATAT‑CGTGCCGGGGATCGAAGCTAtt                                                                                                                          >  2:790380/1‑139 (MQ=255)           cGTCGATCGCTATGCCGACGCACCAGCCATGCATGTCGCGCATCGCTCCCATGTCATTAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATAT‑CGTGCCGGGGATCGAAGCTAtt                                                                                                                          >  2:277540/1‑139 (MQ=255)             tCGATCGCTATGCCGACGCACCAGCCATGCATGTCGCGCATCGCTCCCATGTCATTAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATAT‑CGTGCCGGGGATCGAAGCTATTGc                                                                                                                        >  2:406192/1‑139 (MQ=255)                  cGCTATGCCGACGCACCAGCCATGCATGTCGCGCATCGCTCCCATGTCATTAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATAT‑CGTGCCGGGGATCGAAGCTATTGCCACCg                                                                                                                   <  1:406192/139‑1 (MQ=255)                            aCGCACCAGCCATGCATGTCGCGCATCGCTCCCATGTCATTAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATAT‑CGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGAt                                                                                                         <  1:632852/139‑1 (MQ=255)                               cACCAGCCATGCATGTCGCGCATCGCTCCCATGTCATTAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATAT‑CGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGAt                                                                                                         <  1:364617/136‑1 (MQ=255)                               cACCAGCCATGCATGTCGCGCATCGCTCCCATGTCATTAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATAT‑CGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGAt                                                                                                         >  2:364617/1‑136 (MQ=255)                                            tgTCGCGCATCGCTCCCATGTCATTAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGagga                                                                                          <  1:72194/139‑1 (MQ=255)                                                gcgcATCGCTCCCATGTCATTAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATAT‑CGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGAc                                                                                     >  1:738075/1‑139 (MQ=255)                                                     tCGCTCCCATGTCATTAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATAT‑CGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAAt                                                                                <  2:397036/139‑1 (MQ=255)                                                                  aTTAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATAT‑CGTGCCGGGGATCGAAGCTATTGCCAGCGATATGCTGATCCAACTTGAAGAGGAAGGAc                                                                                     <  1:478503/121‑1 (MQ=255)                                                                  aTTAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATAT‑AGTGCCGGGGATCGAAGCTATTGCCACCGTTATGCTGATCCAACTTGAAGAGGAAGGAc                                                                                     >  2:478503/1‑121 (MQ=255)                                                                      atatGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATAT‑CGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTcg                                                               <  2:660001/139‑1 (MQ=255)                                                                         tGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATAT‑CGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCgcgc                                                            >  1:73000/1‑139 (MQ=255)                                                                           cTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATAT‑CGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCaa                                                          >  1:539608/1‑139 (MQ=255)                                                                            ttGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATAT‑CGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAAc                                                         >  2:536443/1‑139 (MQ=255)                                                                             tGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATAT‑CGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAACg                                                        <  1:790380/139‑1 (MQ=255)                                                                              gATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATAT‑CGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAAtgt                                                                              >  2:251294/1‑116 (MQ=255)                                                                              gATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATAT‑CGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAAtgt                                                                              <  1:251294/116‑1 (MQ=255)                                                                              gATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATAT‑CGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAACGa                                                       >  2:398997/1‑139 (MQ=255)                                                                              gATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATAT‑CGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAACGa                                                       >  1:236013/1‑139 (MQ=255)                                                                                   tGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATAT‑CGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAACGAAATTa                                                  >  2:365097/1‑139 (MQ=255)                                                                                       gCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATAT‑CGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAACGAAATTAACGa                                              <  2:539608/139‑1 (MQ=255)                                                                                       gCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATAT‑CGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAACGAAATTAACGa                                              <  1:798976/139‑1 (MQ=255)                                                                                             cGCCGTGTGGTTGAACTGGAAAAACCACATTATAT‑CGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAACGAAATTAACGATGAATc                                        <  1:406742/139‑1 (MQ=255)                                                                                                cGTGTGGTTGAACTGGAAAAACCACATTATAT‑CGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAACGAAATTAACGATGAATcgcg                                     >  2:318299/1‑139 (MQ=255)                                                                                                  tgtgGTTGAACTGGAAAAACCACATTATAT‑CGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAACg                                                        >  2:94806/1‑118 (MQ=255)                                                                                                  tgtgGTTGAACTGGAAAAACCACATTATAT‑CGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAACg                                                        <  1:94806/118‑1 (MQ=255)                                                                                                  tgtgGTTGAACTGGAAAAACCACATTATAT‑CGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAACGAAATTAACGATGAATCGCGAg                                   >  2:329532/1‑139 (MQ=255)                                                                                                          aaCTGGAAAAACCACATTATAT‑CGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAACGAAATTAACGATGAATcg                                       <  1:389972/127‑1 (MQ=255)                                                                                                          aaCTGGAAAAACCACATTATAT‑CGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAACGAAATTAACGATGAATcg                                       >  2:389972/1‑127 (MQ=255)                                                                                                                aaaaaCCACATTATAT‑CGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAACGAAATTAACGATGAATCGCGAGGGTATCCGTCGCCt                     <  1:398997/139‑1 (MQ=255)                                                                                                                             tat‑CGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAACGAAATTAACGATGAATCGCGAGGGTATCCGTCGCCTGGCGGCAGAagag        <  1:329532/139‑1 (MQ=255)                                                                                                                                   tGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAACGAAATTAACGATGAATCGCGAGGGTATCCGTCGCCTGGCGGCAGAAGagctgca   >  2:437884/1‑139 (MQ=255)                                                                                                                                   tGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAACGAAATTAACGATGAATCGCGAGGGTATCCGTCGCCTGGCGGCAGAAGagctgca   >  1:37191/1‑139 (MQ=255)                                                                                                                                    gCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAACGAAATTAACGATGAATCGCGAGGGTATCCGTCGCCTGGCGGCAGAAGagctgcag  <  2:738075/139‑1 (MQ=255)                                                                                                                                    gCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAACGAAATTAACGATGAATCGCGAGGGTATCCGTCGCCTGGCGGCAGAAGagctgcag  >  2:802663/1‑139 (MQ=255)                                                                                                                                     ccGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAACGAAATTAACGATGAATCGCGAGGGTATCCGTCGCCTggc                  >  1:146007/1‑122 (MQ=255)                                                                                                                                     ccGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAACGAAATTAACGATGAATCGCGAGGGTATCCGTCGCCTggc                  <  2:146007/122‑1 (MQ=255)                                                                                                                                         |                                                                                                                                      AGGTGATTGCCGTCGATCGCTATGCCGACGCACCAGCCATGCATGTCGCGCATCGCTCCCATGTCATTAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATAT‑CGTGCCGGAGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAACGAAATTAACGATGAATCGCGAGGGTATCCGTCGCCTGGCGGCAGAAGAGCTGCAG  >  NC_000913/1930990‑1931258

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