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

ATTGCCGTCGATCGCTATGCCGACGCACCAGCCATGCATGTCGCGCATCGCTCCCATGTCATTAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGAGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAACGAAATTAACGATGAATCGCGAGGGTATCCGTCGCCTGGCGGCAGAAGAGCTGCAGC  >  NC_000913/1930995‑1931259                                                                                                                                   |                                                                                                                                       aTTGCCGTCGATCGCTATGCCGACGCACCAGCCATGCATGTCGCGCATCGCTCCCATGTCATTAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAg                                                                                                                                <  2:360215/139‑1 (MQ=255) aTTGCCGTCGATCGCTATGCCGACGCACCAGCCATGCATGTCGCGCATCGCTCCCATGTCATTAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAg                                                                                                                                >  1:668230/1‑139 (MQ=255)     ccGTCGATCGCTATGCCGACGCACCAGCCATGCATGTCGCGCATCGCTCCCATGTCATTAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTAt                                                                                                                            <  2:958751/139‑1 (MQ=255)                       aCGCACCAGCCATGCATGTCGCGCATCGCTCCCATGTCATTAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTAt                                                                                                                            <  1:701384/121‑1 (MQ=255)                       aCGCACCAGCCATGCATGTCGCGCATCGCTCCCATGTCATTAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTAt                                                                                                                            >  2:701384/1‑121 (MQ=255)                       aCGCACCAGCCATGCATGTCGCGCATCGCTCCCATGTCATTAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCg                                                                                                                    >  2:1550774/1‑129 (MQ=255)                       aCGCACCAGCCATGCATGTCGCGCATCGCTCCCATGTCATTAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCg                                                                                                                    <  1:1550774/129‑1 (MQ=255)                            ccAGCCATGCATGTCGCGCATCGCTCCCATGTCATTAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAAc                                                                                                     <  2:672421/139‑1 (MQ=255)                              aGCCATGCATGTCGCGCATCGCTCCCATGTCATTAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACtt                                                                                                   >  2:1608987/1‑139 (MQ=255)                                             gcATCGCTCCCATGTCATTAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGc                                                                                                                         >  1:655626/1‑102 (MQ=255)                                             gcATCGCTCCCATGTCATTAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGc                                                                                                                         <  2:655626/102‑1 (MQ=255)                                                 cGCTCCCATGTCATTAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAAtg                                                                                <  1:675693/139‑1 (MQ=255)                                                    tCCCATGTCATTAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAAtgttg                                                                             <  1:226993/139‑1 (MQ=255)                                                      ccATGTCATTAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTc                                                                           >  1:269698/1‑139 (MQ=255)                                                         tGTCATTAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTcccc                                                                        >  2:827325/1‑139 (MQ=255)                                                                 atatGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTcg                                                                <  1:1608987/139‑1 (MQ=255)                                                                   atGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGc                                                                                                                         >  2:1491508/1‑80 (MQ=255)                                                                   atGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGc                                                                                                                         <  1:1491508/80‑1 (MQ=255)                                                                     gCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCg                                                                                                                                   <  2:1363772/68‑1 (MQ=255)                                                                     gCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCg                                                                                                                                   >  1:1363772/1‑68 (MQ=255)                                                                     gCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCa                                                            >  1:313986/1‑139 (MQ=255)                                                                      cTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTa                                                                                                                             <  2:739542/73‑1 (MQ=255)                                                                      cTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTa                                                                                                                             >  1:739542/1‑73 (MQ=255)                                                                      cTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGc                                                                                                                         >  1:1163620/1‑77 (MQ=255)                                                                      cTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGc                                                                                                                         <  2:1163620/77‑1 (MQ=255)                                                                      cTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCaa                                                           <  2:668230/139‑1 (MQ=255)                                                                           tGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAACGaaa                                                      >  2:1375848/1‑139 (MQ=255)                                                                                       aCGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGagg                                                                                            >  2:253221/1‑89 (MQ=255)                                                                                       aCGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGagg                                                                                            <  1:253221/89‑1 (MQ=255)                                                                                        cGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAg                                                                                                                                >  2:515991/1‑52 (MQ=255)                                                                                        cGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAg                                                                                                                                <  1:515991/52‑1 (MQ=255)                                                                                           cGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAACGAAATTAACGATGAATcgcg                                      <  2:58679/139‑1 (MQ=255)                                                                                             tgtgGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGaa                                                                                                                                 <  1:1216112/46‑1 (MQ=255)                                                                                             tgtgGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGaa                                                                                                                                 >  2:1216112/1‑46 (MQ=255)                                                                                             tgtgGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGa                                                                                                                   <  1:1007583/60‑1 (MQ=255)                                                                                             tgtgGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAAc                                                          <  1:764411/117‑1 (MQ=255)                                                                                             tgtgGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAAc                                                          >  2:764411/1‑117 (MQ=255)                                                                                             tgtgGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAACGAAATTAACGATGAAt                                          <  1:742365/133‑1 (MQ=255)                                                                                             tgtgGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAACGAAATTAACGATGAAt                                          >  2:742365/1‑133 (MQ=255)                                                                                             tgtgGTTGAACTGGAAAAACCACATTAGATCGCGCCGGGGATCGAAGCTATTGCCACCGa                                                                                                                   >  2:1007583/1‑60 (MQ=255)                                                                                                       cTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCg                                                                                                                    >  1:1014304/1‑49 (MQ=255)                                                                                                       cTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCg                                                                                                                    <  2:1014304/49‑1 (MQ=255)                                                                                                       cTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGagg                                                                                            <  2:429862/73‑1 (MQ=255)                                                                                                       cTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGagg                                                                                            >  1:429862/1‑73 (MQ=255)                                                                                                       cTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTcg                                                                <  1:1035648/101‑1 (MQ=255)                                                                                                       cTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTcg                                                                >  2:1035648/1‑101 (MQ=255)                                                                                                         ggAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAACGAAATTAACGATGAATCGCGAGGGTATCCGTCGc                        <  2:313986/139‑1 (MQ=255)                                                                                                           aaaaaCCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCTCTGCGCTCGCGCAACGAAATTAACGATGAATCGCGAGGGTATCCGTCGCCt                      >  2:1253083/1‑139 (MQ=255)                                                                                                                   caTTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAACGAAATTAACGATGAATCGCGAGGGTATCCGTCGCCTGGCGGCag              <  1:1253083/139‑1 (MQ=255)                                                                                                                             tGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAACGAAATTAACGATGAATcgcg                                      <  2:1466989/105‑1 (MQ=255)                                                                                                                             tGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAACGAAATTAACGATGAATcgcg                                      >  1:1466989/1‑105 (MQ=255)                                                                                                                             tGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAACGAAATTAACGATGAATCGCGAGGGTATCCGTCGCCTGGCGGCAGAAGagctgca    >  1:446850/1‑139 (MQ=255)                                                                                                                             tGCCGGGGATCGAAGCTATTGCCACCGATATGCCGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAACGAAATTAACGATGAATCGCGAGGGTATCCGTCGCCTGGCGGCAGAAGagctgca    >  2:1334030/1‑139 (MQ=255)                                                                                                                              gCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAACGAAATTAACGATGAATCGCGAGGGTATCCGTCGCCTGGCGGCAGAAGagctgcag   >  2:835360/1‑139 (MQ=255)                                                                                                                              gCCGGGGATCGAAGCTATTGCCACCGATATGCCGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAACGAAATTAACGATGAATCGCGAGGGTATCCGTCGCCTGGCGGCAGAAGagctgcag   <  1:1334030/139‑1 (MQ=255)                                                                                                                               ccGGGGATCGAAGCTTTTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACt                                                                                     >  2:1130759/1‑56 (MQ=255)                                                                                                                               ccGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAACGAAATTAACGATGAATCGCGAGGGTATCCGTCGCCTGGCGGCAGAAGagctgcagc  >  2:861246/1‑139 (MQ=255)                                                                                                                               ccGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAACGAAATTAACGATGAATCGCGAGGGTATCCGTCGCCTGGCGGCAGAAGagctgcagc  >  2:1513068/1‑139 (MQ=255)                                                                                                                               ccGGGGATCGAAGCTATTCCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACt                                                                                     <  1:1130759/56‑1 (MQ=255)                                                                                                                                   |                                                                                                                                       ATTGCCGTCGATCGCTATGCCGACGCACCAGCCATGCATGTCGCGCATCGCTCCCATGTCATTAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGAGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAACGAAATTAACGATGAATCGCGAGGGTATCCGTCGCCTGGCGGCAGAAGAGCTGCAGC  >  NC_000913/1930995‑1931259

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