Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	113,831	G→A	V130M (GTG→ATG)	guaC →	GMP reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	113,831	0	G	A	100.0%	179.3 / NA	52	V130M (GTG→ATG)	guaC	GMP reductase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (20/32);  total (20/32)

GCAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACGTGGCGAATGGTTATTCCGAACACTTCGTGCAGTTCGTTGCGAAAGCGCGTGAAGCGTGGCCGACCAAAACCATTTGTGCTGGTAACGTAGTGACTGGTGAAATGTGTGAGGAGCTTATCCTCTCAGGTGCCGA  >  NC_000913/113695‑113964                                                                                                                                         |                                                                                                                                      gCAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACATGAACTCGGCATTAAACTTCGTTTGTATTGACATg                                                                                                                                     >  1:1070496/1‑139 (MQ=255)    aGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACATGGCg                                                                                                                                  <  2:283732/139‑1 (MQ=255)    aGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTCCGTTTGTATTGACATGGCg                                                                                                                                  >  2:1095382/1‑139 (MQ=255)      cGTTTATCAGCAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACATGGCGaa                                                                                                                                <  1:660622/139‑1 (MQ=255)       gTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACATGGCGAAt                                                                                                                               >  1:359191/1‑139 (MQ=255)         ttATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACATGGCGAATgg                                                                                                                             <  2:16351/139‑1 (MQ=255)                caaTTCTTCCGCTGATGTGCTGAAACATTTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACATGGCGAATGTTTATTcc                                                                                                                      >  2:253841/1‑139 (MQ=255)                  aTTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACATGGCGAATGGTTATTCCGa                                                                                                                    <  2:1013580/139‑1 (MQ=255)                   ttcttcCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACATGGCGAATGGTTATTCCGaa                                                                                                                   >  1:1060406/1‑139 (MQ=255)                   ttcttcCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACATGGCGAATGGTTATTCCGaa                                                                                                                   >  2:767228/1‑139 (MQ=255)                   ttcttcCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACATGGCGAATGGTTATTCCGaa                                                                                                                   >  2:382160/1‑139 (MQ=255)                       tcCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACATGGCGAATGGTTATTCCGAACACt                                                                                                               <  1:767228/139‑1 (MQ=255)                       tcCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACATGGCGAATGGTTATTCCGAACACt                                                                                                               <  1:382160/139‑1 (MQ=255)                             gATGTGCTGAAACATGAGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACATGGCGAATGGTTATTCCGAACACTTCGTGc                                                                                                         <  1:253841/139‑1 (MQ=255)                              aTGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACATGGCGAATGGTTATTCCGAACACTTCGTGCa                                                                                                        >  1:983254/1‑139 (MQ=255)                                      aaaCATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACATGGCGAATGGTTATTCCGAACACTTCGTGCAGTTCGTTg                                                                                                >  1:560646/1‑139 (MQ=255)                                       aaCATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCG‑TTGTATTGACATGGCGAATGGTTATTCCGAACACTTCGTGCAGTTCGTTGCg                                                                                              <  2:600566/139‑1 (MQ=255)                                                          cGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACATGGCGAATGGTTATTCCGAACACTTCGTGCAGTTCGTTGCGAAAGCGCGTGAAGCGTgg                                                                            <  2:1060406/139‑1 (MQ=255)                                                          cGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACATGGCGAATGGTTATTCCGAACACTTCGTGCAGTTCGTTGCGAAAGCGAGTGAAGCGTgg                                                                            <  2:707608/139‑1 (MQ=255)                                                            gTACGTCTGATGCGGATTTCGAAAAAACTAGACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACATGGCGaa                                                                                                                                <  1:1081326/85‑1 (MQ=255)                                                            gTACGTCTGATGCGGATTTCGAAAAAACTAGACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACATGGCGaa                                                                                                                                >  2:1081326/1‑85 (MQ=255)                                                            gTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACATGGCGAATGGTTATTCCGAACACTTCGTGCAGTTCGTTGCGAAAGCGCGTGAAGCGTGGcc                                                                          >  1:460524/1‑139 (MQ=255)                                                            gTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACATGGCGAATGGTTATTCCGAACACTTCGTGCAGTTCGTTGCGAAAGCGCGTGAAGCGTGGcc                                                                          >  1:1074163/1‑139 (MQ=255)                                                             tACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACATGGCGAATGGTTATTCCGAac                                                                                                                  <  1:590709/98‑1 (MQ=255)                                                             tACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACATGGCGAATGGTTATTCCGAac                                                                                                                  >  2:590709/1‑98 (MQ=255)                                                                gTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACATGGCGAATGGTTATTCCGAACACTTCGTGCAGTTCGTTGCGAAAGCGCGTGAAGCGTGGCCGAcc                                                                      <  2:213516/139‑1 (MQ=255)                                                                 tCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACATGGCGAATGGTTATTCCGAACACTTCGTGCAGTTCGTTGCGAAAGCGCGTGAAGCGTgg                                                                            <  1:839880/132‑1 (MQ=255)                                                                 tCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACATGGCGAATGGTTATTCCGAACACTTCGTGCAGTTCGTTGCGAAAGCGCGTGAAGCGTgg                                                                            >  2:839880/1‑132 (MQ=255)                                                                     aTGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACATGGCGAATGGTTATTCCGAACACTTCGTGCAGTTCGTTGCGAAAGCGCGTGAAGCGTGGCCGACCAAAAc                                                                 >  1:616859/1‑139 (MQ=255)                                                                          gATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACATGGCGAATGGTTATTCCGAACACTTCGTGCAGTTCGTTGCGAAAGCGCGTGAAGCGTGGCCGACCAAAACCAttt                                                            >  2:127534/1‑139 (MQ=255)                                                                           aTTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACATGGCGAATGGTTATTCCGAACACTTCGTGCAGTTCGTTGCGAAAGCGCGTGAAGCGTGGCCGACCAAAACCATTtg                                                           >  1:79970/1‑139 (MQ=255)                                                                           aTTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACATGGCGAATGGTTATTCCGAACACTTCGTGCAGTTCGTTGCGAAAGCGCGTGAAGCGTGGCCGACCAAAACCATTtg                                                           <  2:709619/139‑1 (MQ=255)                                                                                  aaaaaCTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACATGGCGAATGGTTATTCCGAACACTTCGTGCAGTTCGTTGCGAAAGCGCGTGAAGCGTGGCCGACCAAAACCATTTGTGCTGGt                                                    <  2:924228/139‑1 (MQ=255)                                                                                        tAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACATGGCGAATGGTTATTCCGAACACTTCGTGCAGTTCGTTGCGAAAGCGCGTGAAGCGTGGCCGACCAAAACCATTTGTGCTGGTAACGTa                                              <  2:983254/139‑1 (MQ=255)                                                                                          aaCAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACATGGCGAATGGTTATTCCGAACACTTCGTGCAGTTCGTTGCGAAAGCGCGTGAAGCGTGGCCGACCAAAACCATTTGTGCTGGTAACGTAGt                                            <  2:616859/139‑1 (MQ=255)                                                                                          aaCAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACATGGCGAATGGTTATTCCGAACACTTCGTGCAGTTCGTTGCGAAAGCGCGTGAAGCGTGGCCGACCAAAACCATTTGTGCTGGTAACGTAGt                                            <  1:47962/139‑1 (MQ=255)                                                                                                    cGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACATGGCGAATGGTTATTCCGAACACTTCGTGCAGTTCGTTGCGAAAGCGCGTGAAGCGTGGCCGACCAAAACCATTTGTGCTGGTAACGTAGTGACTGGTGaa                                  <  1:127534/139‑1 (MQ=255)                                                                                                      aCCTGAACCCGGCATTAAACTTCGTTTGTATTGACATGGCGAATGGTTATTCCGAACACTTCGTGCAGTTCGTTGCGAAAGCGCGTGAAGCGTGGCCGACCAAAACCATTTGTGCTGGTAACGTAGTGACTGGTGAAAt                                <  1:679856/139‑1 (MQ=255)                                                                                                          gAACCCGGCATTAAACTTCGTTTGTATTGACATGGCGAATGGTTATTCCGAACACTTCGTGCAGTTCGTTGCGAAAGCGCGTGAAGCGTGGCCGACCAAAACCATTTGTGCTGGTAACGTa                                              >  2:433247/1‑121 (MQ=255)                                                                                                          gAACCCGGCATTAAACTTCGTTTGTATTGACATGGCGAATGGTTATTCCGAACACTTCGTGCAGTTCGTTGCGAAAGCGCGTGAAGCGTGGCCGACCAAAACCATTTGTGCTGGTAACGTa                                              <  1:433247/121‑1 (MQ=255)                                                                                                           aaCCCGGCATTAAACTTCGTTTGTATTGACATGGCGAATGGTTATTCCGAACACTTCGTGCAGTTCGTTGCGAAAGCGCGTGAAGCGTGGCCGACCAAAACCATTTGTGCTGGTAACGTAGTGACTGGTGAAAtgtgtg                           <  2:638193/139‑1 (MQ=255)                                                                                                           aaCCCGGCATTAAACTTCGTTTGTATTGACATGGCGAATGGTTATTCCGAACACTTCGTGCAGTTCGTTGCGAAAGCGCGTGAAGCGTGGCCGACCAAAACCATTTGTGCTGGTAACGTAGTGACTGGTGAAAtgtgtg                           <  1:769888/139‑1 (MQ=255)                                                                                                           aaCCCGGCATTAAACTTCGTTTGTATTGACATGGCGAATGGTTATTCCGAACACTTCGTGCAGTTCGTTGCGAAAGCGCGTGAAGCGTGGCCGACCAAAACCATTTGTGCTGGTAACGTAGTGACTGGTGAAAtgtgtg                           <  1:1041548/139‑1 (MQ=255)                                                                                                              ccGGCATTAGACTTCGTTTGTATTGACATGGCGAATGGTTATTCCGAACACTTCGTGCAGTTCGTTGCGAAAGCGCGTGAAGCGTGGCCGACCAAAACCATTTGTGCTGGTAACGTAGTGACTGGTGAAATGTGTgagg                        <  1:857850/139‑1 (MQ=255)                                                                                                              ccGGCATTAAACTTCGTTTGTATTGACATGGCGAATGGTTATTCCGAACACTTCGTGCAGTTCGTTGCGAAAGCGCGTGAAGCGTGGCCGACCAAAACCATTTGTGCTGGTAACGTAGTGACTGGTGAAATGTGTgagg                        <  2:1074163/139‑1 (MQ=255)                                                                                                                   aTTAAACTTCGTTTGTATTGACATGGCGAATGGTTATTCCGAACACTTCGTGCAGTTCGTTGCGAAAGCGCGTGAAGCGTGGCCGACCAAAACCATTTGTGCTGGTAACGTAGTGACTGGTGAAATGTGTGAGGAGCtt                   <  1:903125/139‑1 (MQ=255)                                                                                                                          ttCGTTTGTATTGACATGGCGAATGGTTATTCCGAACACTTCGTGCAGTTCGTTGCGAAAGCGCGTGAAGCGTGGCCGACCAAAACCATTTGTGCTGGTAACGTAGTGACTGGTGAAATGTGTGAGGAGCTTATCctct            <  2:560646/139‑1 (MQ=255)                                                                                                                          ttCGTTTGTATTGACATGGCGAATGGTTATTCCGAACACTTCGTGCAGTTCGTTGCGAAAGCGCGTGAAGCGTGGCCGACCAAAACCATTTGTGCTGGTAACGTAGTGACTGGTGAAATGTGTGAGGAGCTTATCctct            <  2:613474/139‑1 (MQ=255)                                                                                                                            cGTTTGTATTGACATGGCGAATGGTTATTCCGAACACTTCGTGCAGTTCGTTGCGAAAGCGCGTGAAGCGTGGCCGACCAAAACCATTtgt                                                          >  1:675030/1‑91 (MQ=255)                                                                                                                            cGTTTGTATTGACATAGCGAATGGTTATTCCGAACACTTCGTGCAGTTCGTTGCGAAAGCGCGTGAAGCGTGGCCGACCAAAACCATTtgt                                                          <  2:675030/91‑1 (MQ=255)                                                                                                                              tttGTATTGACATGGCGAATGGTTATTCCGAACACTTCGTGCAGTTCGTTGCGAAAGCGCGTGAAGCGTGGCCGACCAAAACCATTTGTGCTGGTAACGTAGTGACTGGTGAAATGTGTGAGGAGCTTATCCTCTCAgg        >  2:404546/1‑139 (MQ=255)                                                                                                                                    ttGACATGGCGAATGGTTATTCCGAACACTTCGTGCAGTTCGTTGCGAAAGCGCGTGAAGCGTGGCCGACCAAAACCATTTGTGCTGGTAACGTAGTGACTGGTGAAATGTGTGAGGAGCTTATCCTCTCAGGTGCCGa  <  1:1095382/139‑1 (MQ=255)                                                                                                                                         |                                                                                                                                      GCAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACGTGGCGAATGGTTATTCCGAACACTTCGTGCAGTTCGTTGCGAAAGCGCGTGAAGCGTGGCCGACCAAAACCATTTGTGCTGGTAACGTAGTGACTGGTGAAATGTGTGAGGAGCTTATCCTCTCAGGTGCCGA  >  NC_000913/113695‑113964

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