Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,538,491	C→T	G92D (GGC→GAC)	nanM ←	N‑acetylneuraminic acid mutarotase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,538,491	0	C	T	100.0%	153.1 / NA	44	G92D (GGC→GAC)	nanM	N‑acetylneuraminic acid mutarotase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (24/20);  total (24/20)

GTACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATCCAGCTTGTACCATGCCGTACCTGCG  >  NC_000913/4538360‑4538625                                                                                                                                    |                                                                                                                                       gtACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCaaa                                                                                                                                 <  1:132824/139‑1 (MQ=255) gtACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCaaa                                                                                                                                 <  2:573313/139‑1 (MQ=255) gtACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCaaa                                                                                                                                 <  2:635903/139‑1 (MQ=255)          tCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATgtcgc                                                                                                                                     <  1:457944/126‑1 (MQ=255)          tCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAg                                                                                                                        >  2:381589/1‑139 (MQ=255)          tCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTATATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATgtcgc                                                                                                                                     >  2:457944/1‑126 (MQ=255)                   ccGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATc                                                                                                               >  1:38199/1‑139 (MQ=255)                    cGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATgtcgc                                                                                                                                     >  1:541673/1‑116 (MQ=255)                    cGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATgtcgc                                                                                                                                     <  2:541673/116‑1 (MQ=255)                               tCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCaga                                                                                                   >  2:551426/1‑139 (MQ=255)                                  gcgcGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGt                                                                                                <  1:121408/139‑1 (MQ=255)                                     cgTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGc                                                                                             >  2:217383/1‑139 (MQ=255)                                        gCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATgtcgc                                                                                                                                     >  1:184089/1‑96 (MQ=255)                                        gCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATgtcgc                                                                                                                                     <  2:184089/96‑1 (MQ=255)                                                            cTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCa                                                                                                                                   >  1:515001/1‑78 (MQ=255)                                                            cTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCa                                                                                                                                   <  2:515001/78‑1 (MQ=255)                                                                      tGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTaa                                                            >  1:472252/1‑139 (MQ=255)                                                                                        cGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTg                                                                                  <  2:723501/99‑1 (MQ=255)                                                                                        cGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTg                                                                                  >  1:723501/1‑99 (MQ=255)                                                                                        cGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATc                                          >  1:604321/1‑139 (MQ=255)                                                                                        cGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATc                                          >  2:231935/1‑139 (MQ=255)                                                                                        cGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATc                                          >  1:11096/1‑139 (MQ=255)                                                                                           cATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCttt                                       >  2:386642/1‑139 (MQ=255)                                                                                           cATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATACAGCTAACGCTGTCCATTTTTTATCttt                                       >  1:223514/1‑139 (MQ=255)                                                                                             ttAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTgg                                     <  2:38199/139‑1 (MQ=255)                                                                                              tAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCaaa                                                                                                                                 <  2:113637/46‑1 (MQ=255)                                                                                              tAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCaaa                                                                                                                                 >  1:113637/1‑46 (MQ=255)                                                                                              tAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGc                                    <  1:231935/139‑1 (MQ=255)                                                                                               aaaTACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGcc                                   <  1:217383/139‑1 (MQ=255)                                                                                                     cTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGa                                                                                                                       <  1:353258/49‑1 (MQ=255)                                                                                                     cTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGa                                                                                                                       >  2:353258/1‑49 (MQ=255)                                                                                                     cTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCagag                                                                                                  >  2:746483/1‑70 (MQ=255)                                                                                                     cTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCagag                                                                                                  <  1:746483/70‑1 (MQ=255)                                                                                                     cTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCa                                                                                                     >  2:176461/1‑67 (MQ=255)                                                                                                     cTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCa                                                                                                     <  1:176461/67‑1 (MQ=255)                                                                                                     cTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTa                             >  1:5862/1‑139 (MQ=255)                                                                                                     cTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCACTTTTTATCTTTGGCCTGTGTa                             >  1:460281/1‑139 (MQ=255)                                                                                                      tGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTAt                            <  1:575831/139‑1 (MQ=255)                                                                                                      tGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTAt                            <  2:604321/139‑1 (MQ=255)                                                                                                       gAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAAt                                                                    >  2:432210/1‑98 (MQ=255)                                                                                                       gAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAAt                                                                    <  1:432210/98‑1 (MQ=255)                                                                                                        aGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTg                                      >  2:412693/1‑127 (MQ=255)                                                                                                        aGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTg                                      <  1:412693/127‑1 (MQ=255)                                                                                                            aaGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATCCAGCt                      >  1:234179/1‑139 (MQ=255)                                                                                                               cccTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATCCAGCTTGt                   >  1:384713/1‑139 (MQ=255)                                                                                                                 cTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATCCAGCTTGTAc                 <  2:384713/139‑1 (MQ=255)                                                                                                                  tCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATCCAGCTTGTAcc                <  2:5862/139‑1 (MQ=255)                                                                                                                   cGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATCCAGCTTGTACCa               >  2:110809/1‑139 (MQ=255)                                                                                                                   cGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATCCAGCTTGTACCa               >  2:231108/1‑139 (MQ=255)                                                                                                                                aaTGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATCCAGCTTGTACCATGCCGTACCTgcg  <  2:472252/139‑1 (MQ=255)                                                                                                                                    |                                                                                                                                       GTACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATCCAGCTTGTACCATGCCGTACCTGCG  >  NC_000913/4538360‑4538625

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