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%	239.4 / NA	67	E82G (GAG→GGG)	purT	phosphoribosylglycinamide formyltransferase 2
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (37/30);  total (37/30)

AGGTGATTGCCGTCGATCGCTATGCCGACGCACCAGCCATGCATGTCGCGCATCGCTCCCATGTCATTAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGAGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAACGAAATTAACGATGAATCGCGAGGGTATCCGTCGCCTGGCGGCAGAAGAGCTGCAG  >  NC_000913/1930990‑1931258                                                                                                                                        |                                                                                                                                      agGTGATTGCCGTCGATCGCTATGCCGACGCACCAGCCATGCATGTCGCGCATCGCTCCCATGTCATTAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGAt                                                                                                                                    <  1:375534/139‑1 (MQ=255)      aTTGCCGTCGATCGCTATGCCGACGCACCAGCCATGCATGTCGCGCATCGCTCCCATGTCATTAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAg                                                                                                                               >  1:848814/1‑139 (MQ=255)          ccGTCGATCGCTATGCCGACGCACCAGCCATGCATGTCGCGCATCGCTCCCATGTCATTAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTAt                                                                                                                           <  1:1241518/139‑1 (MQ=255)             tCGATCGCTATGCCGACGCACCAGCCATGCATGTCGCGCATCGCTCCCATGTCATTAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGc                                                                                                                        <  1:449595/139‑1 (MQ=255)             tCGATCGCTATGCCGACGCACCAGCCATGCATGTCGCGCATCGCTCCCATGTCATTAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGc                                                                                                                        <  1:1417218/139‑1 (MQ=255)             tCGATCGCTATGCCGACGCACCAGCCATGCATGTCGCGCATCGCTCCCATGTCATTAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGc                                                                                                                        <  1:993048/139‑1 (MQ=255)                              gCACCAGCCATGCATGTCGCGCATCGCTCCCATGTCATTAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGa                                                                                                                  <  2:281980/128‑1 (MQ=255)                              gCACCAGCCATGCATGTCGCGCATCGCTCCCATGTCATTAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGa                                                                                                                  >  1:281980/1‑128 (MQ=255)                               cACCAGCCATGCATGTCGCGCATCGCTCCCATGTCATTAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCa                                                                                                      >  2:111513/1‑139 (MQ=255)                                   aGCCATGCATGTCGCGCATCGCTCCCATGTCATTAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACtt                                                                                                  >  1:2837/1‑139 (MQ=255)                                               cgcgCATCGCTCCCATGTCATTAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGaggaagga                                                                                      >  1:857748/1‑139 (MQ=255)                                               cgcgCATCGCTCCCATGTCATTAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGaggaagga                                                                                      >  1:78753/1‑139 (MQ=255)                                                gcgcATCGCTCCCATGTCATTAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGAc                                                                                     <  2:78753/139‑1 (MQ=255)                                                gcgcATCGCTCCCATGTCATTAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGAc                                                                                     <  2:857748/139‑1 (MQ=255)                                                             aTGTCATTAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTccc                                                                        <  2:1459710/139‑1 (MQ=255)                                                              tGTCATTAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTcccc                                                                       >  1:460989/1‑139 (MQ=255)                                                                       tatGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTcgc                                                              >  2:895949/1‑139 (MQ=255)                                                                        atGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATTTTGTCCCCTGCGCTcgcg                                                             <  1:895949/139‑1 (MQ=255)                                                                         tGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCgcgc                                                            >  2:654567/1‑139 (MQ=255)                                                                          gCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGACGCTATTGCTACCAATATGCTGATCAAACTTGAAGAGGAAGTACTGAATGTTGTGCCGTGCGCTCGCGCa                                                           >  1:402910/1‑139 (MQ=255)                                                                          gCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCa                                                           >  2:971850/1‑139 (MQ=255)                                                                          gCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCa                                                           >  1:1561828/1‑139 (MQ=255)                                                                          gCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCa                                                           >  1:1050384/1‑139 (MQ=255)                                                                           cTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAAc                                                                                                    >  2:429928/1‑97 (MQ=255)                                                                           cTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAAc                                                                                                    <  1:429928/97‑1 (MQ=255)                                                                           cTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCaa                                                          >  2:637188/1‑139 (MQ=255)                                                                            ttGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAAc                                                         <  2:769369/139‑1 (MQ=255)                                                                              gATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCa                                                                                                                      >  2:1354271/1‑76 (MQ=255)                                                                              gATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCa                                                                                                                      <  1:1354271/76‑1 (MQ=255)                                                                              gATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGagga                                                                                          >  2:219531/1‑104 (MQ=255)                                                                              gATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGagga                                                                                          <  1:219531/104‑1 (MQ=255)                                                                              gATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAACGa                                                       >  1:564252/1‑139 (MQ=255)                                                                              gATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAACGa                                                       >  2:1102373/1‑139 (MQ=255)                                                                                tGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAACGaaa                                                     >  2:565138/1‑139 (MQ=255)                                                                                 ggTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAACGAAAt                                                    <  2:662474/139‑1 (MQ=255)                                                                                   tGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAAc                                                         <  1:738235/132‑1 (MQ=255)                                                                                   tGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAAc                                                         >  2:738235/1‑132 (MQ=255)                                                                                     aTGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATg                                                                                                              >  2:1369364/1‑77 (MQ=255)                                                                                     aTGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATg                                                                                                              <  1:1369364/77‑1 (MQ=255)                                                                                      tGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAACGAAATTAACg                                               <  1:1611277/139‑1 (MQ=255)                                                                                         aTTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAACGAAATTAACGATg                                            <  1:1558492/139‑1 (MQ=255)                                                                                           tACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAACGAAATTAACGATGaa                                          <  1:1363455/139‑1 (MQ=255)                                                                                            aCGCCGTGTGGCTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAACGAAATTAACGATGAAt                                         <  2:1050384/139‑1 (MQ=255)                                                                                             cGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCg                                                                                                                                  >  2:251181/1‑49 (MQ=255)                                                                                             cGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCg                                                                                                                                  <  1:251181/49‑1 (MQ=255)                                                                                                cGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCGCCTGCGCTCGCGCAACGAAATTAACGATGAATcgcg                                     <  1:111513/139‑1 (MQ=255)                                                                                                cGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAACGAAATTAACGATGAATcgcg                                     >  1:360103/1‑139 (MQ=255)                                                                                                     ggTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAACGAAATTAACGATGAATCGCGAGGGt                                <  1:1102373/139‑1 (MQ=255)                                                                                                            cTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAAc                                                         >  1:130902/1‑107 (MQ=255)                                                                                                            cTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAAc                                                         <  2:130902/107‑1 (MQ=255)                                                                                                               gAAAAACGACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAACGAAATTAACGATGAATCGCGAGGGTATCCGTCGcc                      >  2:984433/1‑139 (MQ=255)                                                                                                                    aCCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAACGAAATTAACGATGAATCGCGAGGGTATCCGTCGCCTggcg                 >  1:1538837/1‑139 (MQ=255)                                                                                                                      cacaTTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAACGAAATTAACGATGAATCGCGAGGGTATCCGTCGCCTggcggc               >  2:379779/1‑139 (MQ=255)                                                                                                                         aTTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATg                                                                                                              >  1:882808/1‑41 (MQ=255)                                                                                                                         aTTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATg                                                                                                              <  2:882808/41‑1 (MQ=255)                                                                                                                         aTTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGa                                                                                                <  1:864274/55‑1 (MQ=255)                                                                                                                         aTTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGa                                                                                                >  2:864274/1‑55 (MQ=255)                                                                                                                                  tGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAACGAAATTAACGATGAATCGCGAGGGTATCCGTCGCCTGGCGGCAGAAGagctgca   >  2:389165/1‑139 (MQ=255)                                                                                                                                  tGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAACGAAATTAACGATGAATCGCGAGGGTATCCGTCGCCTGGCGGCAGAAGagctgca   >  2:952153/1‑139 (MQ=255)                                                                                                                                  tGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAACGAAATTAACGAGGAATCGCGAGGGTATCCGTCGCCTGGCGGCAGAAGagctgca   >  1:771159/1‑139 (MQ=255)                                                                                                                                   gCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGa                                                                                  >  2:238061/1‑59 (MQ=255)                                                                                                                                   gCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGa                                                                                  <  1:238061/59‑1 (MQ=255)                                                                                                                                   gCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAACGAAATTAACGATGAATc                                        >  2:1216500/1‑101 (MQ=255)                                                                                                                                   gCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAACGAAATTAACGATGAATc                                        <  1:1216500/101‑1 (MQ=255)                                                                                                                                   gCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAACGAAATTAACGATGAATCGCGAGGGTATCCGTCGCCTGGCGGCAGAAGagctgcag  >  2:1440387/1‑139 (MQ=255)                                                                                                                                    ccGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAACGAAATTAACGATGAATcgcg                                     <  2:355197/103‑1 (MQ=255)                                                                                                                                    ccGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAACGAAATTAACGATGAATcgcg                                     >  1:355197/1‑103 (MQ=255)                                                                                                                                        |                                                                                                                                      AGGTGATTGCCGTCGATCGCTATGCCGACGCACCAGCCATGCATGTCGCGCATCGCTCCCATGTCATTAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGAGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAACGAAATTAACGATGAATCGCGAGGGTATCCGTCGCCTGGCGGCAGAAGAGCTGCAG  >  NC_000913/1930990‑1931258

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: ‑