Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,538,474	C→T	E98K (GAG→AAG)	nanM ←	N‑acetylneuraminate mutarotase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,538,474	0	C	T	100.0%	205.3 / NA	57	E98K (GAG→AAG)	nanM	N‑acetylneuraminate mutarotase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (36/21);  total (36/21)

ACATAAGCCTTGCCGTTGTGTACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATCCAGCTTGTAC  >  NC_000913/4538341‑4538610                                                                                                                                      |                                                                                                                                         aCATAAGCCTTGCCGTTGTGTACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTTGCTGtt                                                                                                                                    <  2:230439/140‑1 (MQ=255)     aaGCCTTGCCGTTGTGTACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTTGCTGTTTTTg                                                                                                                                >  1:205298/1‑140 (MQ=255)      aGCCTTGCCGTTGTGTACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTTGCTGTTTTTGc                                                                                                                               >  1:222255/1‑140 (MQ=255)           tGCCGTTGTGTACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTTGCTGTTTTTGCCAATg                                                                                                                          <  1:668702/140‑1 (MQ=255)           tGCCGTTGTGTACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTTGCTGTTTTTGCCAATg                                                                                                                          >  2:210789/1‑140 (MQ=255)           tGCCGTTGTGTACAAAAGTCACATGGCCCGCCATGCCCATCGGCACGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTTGCTGTTTTTGCCAATg                                                                                                                          <  1:749382/140‑1 (MQ=255)            gCCGTTGTGTACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTTGCTGTTTTTGCCAATgc                                                                                                                         >  2:433109/1‑140 (MQ=255)             ccGTTGTGTACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTTGCTGTTTTTGCCAATgcc                                                                                                                        >  2:44598/1‑140 (MQ=255)                 tgtgTACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTTGCTGTTTTTGCCAATGCCGCCa                                                                                                                    <  2:57244/140‑1 (MQ=255)                  gtgtACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTTGCTGTTTTTGCCAATGCCGCCa                                                                                                                    >  2:336357/1‑139 (MQ=255)                  gtgtACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTTGCTGTTTTTGCCAATGCCGCCa                                                                                                                    <  1:336357/139‑1 (MQ=255)                   tgtACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTTGCTGTTTTTGCCAATGCCGCCaaa                                                                                                                  <  2:222255/140‑1 (MQ=255)                             tCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTTGCTGTTTTTGCCAATGCCGCCAAAcaca                                                                                                              <  2:453824/134‑1 (MQ=255)                             tCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTTGCTGTTTTTGCCAATGCCGCCAAAcaca                                                                                                              >  1:453824/1‑134 (MQ=255)                               acaTGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTTGCTGTTTTTGCCAATGCCGCCAAACACATACAGAtt                                                                                                      >  1:675177/1‑140 (MQ=255)                                       cGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTTGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCaa                                                                                              >  1:152538/1‑140 (MQ=255)                                                        cgTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTTGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCt                                                                             >  2:397772/1‑140 (MQ=255)                                                        cgTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTTGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCt                                                                             <  1:210789/140‑1 (MQ=255)                                                        cgTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTTGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCt                                                                             <  1:623922/140‑1 (MQ=255)                                                         gTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTTGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCtt                                                                            >  2:169063/1‑140 (MQ=255)                                                               cATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTTGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTg                                                                           <  1:512418/135‑1 (MQ=255)                                                               cATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTTGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTg                                                                           >  2:512418/1‑135 (MQ=255)                                                                         aCCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTTGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTg                                                                   >  1:15322/1‑133 (MQ=255)                                                                         aCCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTTGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTg                                                                   <  2:15322/133‑1 (MQ=255)                                                                          cccAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTTGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCa                                                           >  2:568396/1‑140 (MQ=255)                                                                            cAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTTGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAgg                                                         <  1:169063/140‑1 (MQ=255)                                                                               cTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTTGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGaa                                                      >  1:386118/1‑140 (MQ=255)                                                                                         tGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTTGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAAc                                            >  1:153388/1‑140 (MQ=255)                                                                                         tGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTTGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAAc                                            >  2:319732/1‑140 (MQ=255)                                                                                                          aCGTCATTAAATACCTGAGTCAAGCCCTTGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATc                           >  2:779896/1‑140 (MQ=255)                                                                                                           cGTCATTAAATACCTGAGTCAAGCCCTTGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCt                          >  2:725240/1‑140 (MQ=255)                                                                                                           cGTCATTAAATACCTGAGTCAAGCCCTTGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCt                          >  2:677308/1‑140 (MQ=255)                                                                                                            gTCATTAAATACCTGAGTCAAGCCCTTGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCtt                         >  2:21953/1‑140 (MQ=255)                                                                                                             tCATTAAATACCTGAGTCAAGCCCTTGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCttt                        <  1:319732/140‑1 (MQ=255)                                                                                                              cATTAAATACCTGAGTCAAGCCCTTGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTg                       >  1:40680/1‑140 (MQ=255)                                                                                                               aTTAAATACCTGAGTCAAGCCCTTGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTgg                      <  1:433109/140‑1 (MQ=255)                                                                                                               aTTAAATACCTGAGTCAAGCCCTTGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTgg                      >  2:326398/1‑140 (MQ=255)                                                                                                                ttAAATACCTGAGTCAAGCCCTTGCTG‑TTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTc                                      >  1:415428/1‑122 (MQ=255)                                                                                                                ttAAATACCTGAGTCAAGCCCTTGCTG‑TTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTc                                      <  2:415428/122‑1 (MQ=255)                                                                                                                 tAAATACCTGAGTCAAGCCCTTGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGcc                    <  2:539276/140‑1 (MQ=255)                                                                                                                  aaaTACCTGAGTCAAGCCCTTGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCt                   <  1:677308/140‑1 (MQ=255)                                                                                                                     tACCTGAGTCAAGCCCTTGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCtgtg                >  1:340704/1‑140 (MQ=255)                                                                                                                      aCCTGAGTCAAGCCCTTGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTgtgt               >  1:625749/1‑140 (MQ=255)                                                                                                                        cTGAGTCAAGCCCTTGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTAt             >  1:264925/1‑140 (MQ=255)                                                                                                                        cTGAGTCAAGCCCTTGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTAt             >  1:93525/1‑140 (MQ=255)                                                                                                                        cTGAGTCAAGCCCTTGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTAt             >  2:666868/1‑140 (MQ=255)                                                                                                                        cTGAGTCAAGCCCTTGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTACCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCAGTCCATTTTTTATCTTTGGCCTGGGTAt             >  1:534696/1‑140 (MQ=255)                                                                                                                          gAGTCAAGCCCTTGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATcc           >  1:99397/1‑140 (MQ=255)                                                                                                                           aGTCAAGCCCTTGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATCCa          >  1:278864/1‑140 (MQ=255)                                                                                                                           aGTCAAGCCCTTGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATCCa          >  2:447187/1‑140 (MQ=255)                                                                                                                           aGTCAAGCCCTTGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATCCa          >  1:74434/1‑140 (MQ=255)                                                                                                                             tCAAGCCCTTGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATCCAGc        <  1:477918/140‑1 (MQ=255)                                                                                                                             tCAAGCCCTTGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATCCAGc        >  2:477918/1‑140 (MQ=255)                                                                                                                             tCAAGCCCTTGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATCCAGc        <  1:725240/140‑1 (MQ=255)                                                                                                                                aGCCCTTGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATCCAGCTTg     >  1:556646/1‑140 (MQ=255)                                                                                                                                  cccTTGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATCCAGCTTGTa   <  2:340704/140‑1 (MQ=255)                                                                                                                                   ccTTGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATCCAGCTTGTAc  <  2:625749/140‑1 (MQ=255)                                                                                                                                      |                                                                                                                                         ACATAAGCCTTGCCGTTGTGTACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATCCAGCTTGTAC  >  NC_000913/4538341‑4538610

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