Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,537,169	A→G	D142D (GAT→GAC)	nanS ←	N‑acetyl‑9‑O‑acetylneuraminate esterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,537,169	0	A	G	100.0%	190.9 / NA	53	D142D (GAT→GAC)	nanS	N‑acetyl‑9‑O‑acetylneuraminate esterase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (22/31);  total (22/31)

GTCATTAAGTCAAATTCGCCTTGCATCCAGCATGCGCCGAGGAATTTGTTCTGCGGATTTTTTGCCAGTGCGGCTCGCGTTCTGCTGACTAAATCCTGGTACAGCGGAGTATCCGTTCCCCAACGACAAGCATCATGGCTGGCTCCGTGCCGTTCTGAATATGTCCCTTCGCTGCCCGCGGTAAAAGCCGATCCGCCACGGCAACACGGAACGATAAGAATCCCTGCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCAATA  >  NC_000913/4537038‑4537301                                                                                                                                    |                                                                                                                                     gTCATTAAGTCAAATTCGCCTTGCATCCAGCATGCGCCGAGGAATTTGTTCTGCGGATTTTTTGCCAGTGCGGCTCGCGTTCTGCTGACTAAATCCTGGTACAGCGGAGTATCCGTTCCCCAACGACAAGCGTCAtggct                                                                                                                              >  2:391981/1‑140 (MQ=255)   cATTAAGTCAAATTCGCCTTGCATCCAGCATGCGCCGAGGAATTTGTTCTGCGGATTTTTTGCCAGTGCGGCTCGCGTTCTGCTGACTAAATCCTGGTACAGCGGAGTATCCGTTCCCCAACGACAAGcgtc                                                                                                                                    >  1:136552/1‑132 (MQ=255)   cATTAAGTCAAATTCGCCTTGCATCCAGCATGCGCCGAGGAATTTGTTCTGCGGATTTTTTGCCAGTGCGGCTCGCGTTCTGCTGACTAAATCCTGGTACAGCGGAGTATCCGTTCCCCAACGACAAGcgtc                                                                                                                                    <  2:136552/132‑1 (MQ=255)         gTCAAATTCGCCTTGCATCCAGCATGCGCCGAGGAATTTGTTCTGCGGATTTTTTGCCAGTGCGGCTCGCGTTCTGCTGACTAAATCCTGGTACAGCGGAGTATCCGTTCCCCAACGACAAGCGTCATGGCTGGCTCCGt                                                                                                                      >  2:646482/1‑140 (MQ=255)          tCAAATTCGCCTTGCATCCAGCCTGAGCCGAGGAATTCGTTCTGCGGATTTTTTGCCAGTGCGGCTCGCGTTCTGCTGACTAAATCCTGGTACAGCGGAGTATCCGTTCCCCAACGACAAGCGTCATGGCTGGCTCCGTg                                                                                                                     <  1:272840/140‑1 (MQ=255)          tCAAATTCGCCTTGCATCCAGCATGCGCCGAGGAATTTGTTCTGCGGATTTTTTGCCAGTGCGGCTCGCGTTCTGCTGACTAAATCCTGGTACAGCGGAGTATCCGTTCCCCAACGACAAGCGTCATGGCTGGCTCCGTg                                                                                                                     <  2:331897/140‑1 (MQ=255)           cAAATTCGCCTTGCATCCAGCATGCGCCGAGGAATTTGTTCTGCGGATTTTTTGCCAGTGCGGCTCGCGTTCTGCTGACTAAATCCTGGTACAGCGGAGTATCCGTTCCCAAACTACAAGCGTCATGGCTGGCTCCGTGc                                                                                                                    >  1:279048/1‑140 (MQ=255)            aaaTTCGCCTTGCATCCAGCATGCGCCGAGGAATTTGTTCTGCGGATTTTTTGCCAGTGCGGCTCGCGTTCTGCTGACTAAATCCTGGTACAGCGGAGTATCCGTTCCCCAACGACAAGCGTCATGGCTGGCTCCGTGcc                                                                                                                   >  1:582733/1‑140 (MQ=255)               ttCGCCTTGCATCCAGCATGCGCCGAGGAATTTGTTCTGCGGATTTTTTGCCAGTGCGGCTCGCGTTCTGCTGACTAAATCCTGGTACAGCGGAGTATCCGTTCCCCAACGACAAGCGTCATGGCTGGCTCCGTGCCGtt                                                                                                                >  2:395616/1‑140 (MQ=255)               ttCGCCTTGCATCCAGCATGCGCCGAGGAATTTGTTCTGCGGATTTTTTGCCAGTGCGGCTCGCGTTCTGCTGACTAAATCCTGGTACAGCGGAGTATCCGTTCCCCAACGACAAGCGTCATGGCTGGCTCCGTGCCGtt                                                                                                                <  2:106943/140‑1 (MQ=255)                tCGCCTTGCATCCAGCATGCGCCGAGGAATTTGTTCTGCGGATTTTTTGCCAGTGCGGCTCGCGTTCTGCTGACTAAATCCTGGTACAGCGGAGTATCCGTTCCCCAACGACAAGCGTCATGGCTGGCTCCGTGCCGTTc                                                                                                               <  2:82722/140‑1 (MQ=255)                 cGCCTTGCATCCAGCATGCGCCGAGGAATTTGTTCTGCGGATTTTTTGCCAGTGCGGCTCGCGTTCTGCTGACTAAATCCTGGTACAGCGGAGTATCCGTTCCCCAACGACAAGCGTCATGGCTGGCTCCGTGCCGTTCt                                                                                                              <  1:540708/140‑1 (MQ=255)                      tGCATCCATCATGCGCACAGGATTTTGTTCTCCCGATTTATTGCCAGTTCGGCTCGCGTTCTGCTGACTAAATCCTGGTACAGCGGAGTATCCGTTCCCCAACGACAAGCGTCATGGCTGGCTCCGTGCCGTTCTGAata                                                                                                         <  2:320732/140‑1 (MQ=255)                        cATCCAGCATGCGCCGAGGAATTTGTTCTGCGGATTTTTTGCCAGTGCGGCTCGCGTTCTGCTGACTAAATCCTGGTACAGCGGAGTATCCGTTCCCCAACGACAAGCGTCATGGCTGGCTCCGTGCCGTTCTGAATATg                                                                                                       <  1:44778/140‑1 (MQ=255)                            cAGCATGCGCCGAGGAATTTGTTCTGCGGATTTTTTGCCAGTGCGGCTCGCGTTCTGCTGACTAAATCCTGGTACAGCGGAGTATCCGTTCCCCAACGACAAGCGTCATGGCTGGCTCCGTGCCGTTCTGAATATGTccc                                                                                                   <  1:547166/140‑1 (MQ=255)                                ttgcaccGAGGAATTTGTTCTGCGGATTTTTTGCCAGTGCGGCTCGCGTTCTGCTGACTAAATCCTGGTACAGCGGAGTATCCGTTCCCCAACGACAAGCGTCATGGCTGGCTCCGTGCCGTTCTGAATATGTCCCTTCg                                                                                               >  1:358546/6‑140 (MQ=255)                                   cgcCGAGGAATTTGTTCTGCGGATTTTTTGCCAGTGCGGCTCGCGTTCTGCTGACTAAATCCTGGTACAGCGGAGTATCCGTTCCCCAACGACAAGCGTCATGGCTGGCTCCGTGCCGTTCTGAATATGTCCCTTCGCTg                                                                                            >  1:130918/1‑140 (MQ=255)                                   cgcCGAGGAATTTGTTCTGCGGATTTTTTGCCAGTGCGGCTCGCGTTCTGCTGACTAAATCCTGGTACAGCGGAGTATCCGTTCCCCAACGACAAGCGTCATGGCTGGCTCCGTGCCGTTCTGAATATGTCCCTTCGCTg                                                                                            >  1:91/1‑140 (MQ=255)                                      cGAGGAATTTGTTCTGCGGATTTTTTGCCAGTGCGGCTCGCGTTCTGCTGACTAAATCCTGGTACAGCGGAGTATCCGTTCCCCAACGACAAGCGTCATGGCTGGCTCCGTGCCGTTCTGAATATGTCCCTTCGCTGccc                                                                                         <  1:596820/140‑1 (MQ=255)                                        aGGAATTTGTTCTGCGGATTTTTTGCCAGTGCGGCTCGCGTTCTGCTGACTAAATCCTGGTACAGCGGAGTATCCGTTCCCCAACGACAAGCGTCATGGCTGGCTCCGTGCCGTTCTGAATATGTCCCTTCGCTGCCcgc                                                                                       >  2:61616/1‑140 (MQ=255)                                        aGGAATTTGTTCTGCGGATTTTTTGCCAGTGCGGCTCGCGTTCTGCTGACTAAATCCTGGTACAGCGGAGTATCCGTTCCCCAACGACAAGCGTCATGGCTGGCTCCGTGCCGTTCTGAATATGTCCCTTCGCTGCCcgc                                                                                       >  2:89703/1‑140 (MQ=255)                                         ggAATTTGTTCTGCGGATTTTTTGCCAGTGCGGCTCGCGTTCTGCTGACTAAATCCTGGTACAGCGGAGTATCCGTTCCCCAACGACAAGCGTCATGGCTGGCTCCGTGCCGTTCTGAATATGTCCCTTCGCTGCCcgcg                                                                                      >  1:550350/1‑140 (MQ=255)                                           aaTTTGTTCTGCGGATTTTTTGTCAGTGCGGTTCGCGTTTTGCTGACTAAATCCTGGTACAGCGGAGTATCCGTTCCCCAACGACAAGCGTCATGGCTGGCTCCGTGCCGTTCTGAATATGTCCCTTCGCTGCCCGCGGt                                                                                    <  2:251466/140‑1 (MQ=255)                                                    tGCGGATTTTTTGCCAGTGCGGCTCGCGTTCTGCTGACTAAATCCTGGTACAGCGGAGTATCCGTTCCCCAACGACAAGCGTCATGGCTGGCTCCGTGCCGTTCTGAATATGTCCCTTCGCTGCCCGCGGTAAAAGCCGa                                                                           <  1:646482/140‑1 (MQ=255)                                                     gCGGATTTTTTGCCAGTGCGGCTCGCGTTCTGCTGACTAAATCCTGGTACAGCGGAGTATCCGTTCCCCAACGACAAGCGTCATGGCTGGCTCCGTGCCGTTCTGAATATGTCCCTTCGCTGCCCGCGGTAAAAGCCGAt                                                                          >  1:623073/1‑140 (MQ=255)                                                         aTTTTTTGCCAGTGCGGCTCGCGTTCTGCTGACTAAATCCTGGTACAGCGGAGTATCCGTTCCCCAACGACAAGCGTCATGGCTGGCTCCGTGCCGTTCTGAATATGTCCCTTCGCTGCCCGCGGTAAAAGCCGATCCGc                                                                      <  1:146491/140‑1 (MQ=255)                                                         aTTTTTTGCCAGTGCGGCTCGCGTTCTGCTGACTAAATCCTGGTACAGCGGAGTATCCGTTCCCCAACGACAAGCGTCATGGCTGGCTCCGTGCCGTTCTGAATATGTCCCTTCGCTGCCCGCGGTAAAAGCCGATCCGc                                                                      <  2:130918/140‑1 (MQ=255)                                                                      gCGGCTCGCGTTCTGCTGACTAAATCCTGGTACAGCGGAGTATCCGTTCCCCAACGACAAGCGTCATGGCTGGCTCCGTGCCGTTCTGAATATGTCCCTTCGCTGCCCGCGGTAAAAGCCGATCCGCCACGGCAACACgg                                                         <  2:647553/140‑1 (MQ=255)                                                                        ggCTCGCGTTCTGCTGAGTAAATCCTGGTACGGCGTAGTCTCCGTTCCCCAACGACAAGCGTCATGGCTGGCTCCGTGCCGTTCTGAATATGTCCCTTCGCTGCCCGCGGTAAAAGCCGATCCGCCACGGCAACACGGaa                                                       <  2:258558/140‑1 (MQ=255)                                                                        ggCTCGCGTTCTGCTGACTAAATCCTGGTACAGCGGAGTATCCGTTCCCCAACGACAAGCGTCATGGCTGGCTCCGTGCCGTTCTGAATATGTCCCTTCGCTGCCCGCGGTAAAAGCCGATCCGCCACGGCAACACGGaa                                                       <  1:635962/140‑1 (MQ=255)                                                                        ggCTCGCGTTCTGCTGACTAAATCCTGGTACAGCGGAGTATCCGTTCCCCAACGACAAGCGTCATGGCTGGCTCCGTGCCGTTCTGAATATGTCCCTTCGCTGCCCGCGGTAAAAGCCGATCCGCCACGGCAACACGGaa                                                       <  1:61616/140‑1 (MQ=255)                                                                         gCTCGCGTTCTGCTGACTAAATCCTGGTACAGCGGAGTATCCGTTCCCCAACGACAAGCGTCATGGCTGGCTCCGTGCCGTTCTGAATATGTCCCTTCGCTGCCCGCGGTAAAAGCCGATCCGCCACGGCAACACGGAAc                                                      <  2:664503/140‑1 (MQ=255)                                                                                   tgctgACTAAATCCTGGTACAGCGGAGTATCCGTTCCCCAACGACAAGCGTCATGGCTGGCTCCGTGCCGTTCTGAATATGTCCCTTCGCTGCCCGCGGTAAAAGCCGATCCGCCACGGCAACACGGAACGATAAGAATc                                            <  1:22955/140‑1 (MQ=255)                                                                                    gctgACTAAATCCTGGTACAGCGGAGTATCCGTTCCCCAACGACAAGCGTCATGGCTGGCTCCGTGCCGTTCTGAATATGTCCCTTCGCTGCCCGCGGTAAAAGCCGATCCGCCACGGCAACACGGAACGATAAGAATcc                                           <  2:38819/140‑1 (MQ=255)                                                                                    gctgACTAAATCCTGGTACAGCGGAGTATCCGTTCCCCAACGACAAGCGTCATGGCTGGCTCCGTGCCGTTCTGAATATGTCCCTTCGCTGCCCGCGGTAAAAGCCGATCCGCCACGGCAACACGGAACGATAAGAATcc                                           <  2:544890/140‑1 (MQ=255)                                                                                           aaaTCCTGGTACAGCGGAGTATCCGTTCCCCAACGACAAGCGTCATGGCTGGCTCCGTGCCGTTCTGAATATGTCCCTTCGCTGCCCGCGGTAAAAGCCGATCCGCCACGGCAACACGGAACGATAAGAATCCCTGCAtt                                    <  2:652011/140‑1 (MQ=255)                                                                                                cTGGTACAGCGGAGTATCCGTTCCCCAACGACAAGCGTCATGGCTGGCTCCGTGCCGTTCTGAATATGTCCCTTCGCTGCCCGCGGTAAAAGCCGATCCGCCACGGCAACACGGAACGATAAGAATCCCTGCATTATCAg                               <  1:89703/140‑1 (MQ=255)                                                                                                cTGGTACAGCGGAGTATCCGTTCCCCAACGACAAGCGTCATGGCTGGCTCCGTGCCGTTCTGAATATGTCCCTTCGCTGCCCGCGGTAAAAGCCGATCCGCCACGGCAACACGGAACGATAAGAATCCCTGCATTATCAg                               <  1:523937/140‑1 (MQ=255)                                                                                                    tACAGCGGAGTATCCGTTCCCCAACGACAAGCGTCATGGCTGGCTCCGTGCCGTTCTGAATATGTCCCTTCGCTGCCCGCGGTAAAAGCCGATCCGCCACGGCAACACGGAACGATAAGAATCCCTGCATTATCAGGAAt                           >  1:606853/1‑140 (MQ=255)                                                                                                       aGCGGAGTATCCGTTCCCCAACGACAAGCGTCATGGCTGGCTCCGTGCCGTTCTGAATATGTCCCTTCGCTGCCCGCGGTAAAAGCCGATCCGCCACGGCAACACGGAACGATAAGAATCCCTGCATTATCAGGAATGaa                        <  1:603177/140‑1 (MQ=255)                                                                                                        gCGGAGTATCCGTTCCCCAACGACAAGCGTCATGGCTGGCTCCGTGCCGTTCTGAATATGTCCCTTCGCTGCCCGCGGTAAAAGCCGATCCGCCACGGCAACACGGAACGATAAGAATCCCTGCATTATCAGGAATGAAg                       >  2:500088/1‑140 (MQ=255)                                                                                                        gCGGAGTATCCGTTCCCCAACGACAAGCGTCATGGCTGGCTCCGTGCCGTTCTGAATATGTCCCTTCGCTGCCCGCGGTAAAAGCCGATCCGCCACGGCAACACGGAACGATAAAAATCCCTGCATTATCAGGAATGAAg                       >  2:519239/1‑140 (MQ=255)                                                                                                         cGGAGTATCCGTTCCCCAACGACAAGCGTCATGGCTGGCTCCGTGCCGTTCTGAATATGTCCCTTCGCTGCCCGCGGTAAAAGCCGATCCGCCACGGCAACACGGAACGATAAGAATCCCTGCATTATCAGGAATGAAgg                      >  1:485543/1‑140 (MQ=255)                                                                                                         cGGAGTATCCGTTCCCCAACGACAAGCGTCATGGCTGGCTCCGTGCCGTTCTGAATATGTCCCTTCGCTGCCCGCGGTAAAAGCCGATCCGCCACGGCAACAAGGAACGATAAGAATCCCTGCATTATCAGGAATGAAgg                      >  1:276039/1‑140 (MQ=255)                                                                                                          ggAGTATCCGTTCCCCAACGACAAGCGTCATGGCTGGCTCCGTGCCGTTCTGAATATGTCCCTTCGCTGCCCGCGGTAAAAGCCGATCCGCCACGGCAACACGGAACGATGAGAGTCTCTGCATGATCAGGAATGAAggg                     >  1:173882/1‑140 (MQ=255)                                                                                                          ggAGTATCCGTTCCCCAACGACAAGCGTCATGGCTGGCTCCGTGCCGTTCTGAATATGTCCCTTCGCTGCCCGCGGTAAAAGCCGATCCGCCACGGCAACACGGAACGATAAGAATCCCTGCATTATCAGGAATGAAggg                     >  1:485665/1‑140 (MQ=255)                                                                                                              tATCCGTTCCCCAACGACAAGCGTCATGGCTGGCTCCGTGCCGTTCTGAATATGTCCCTTCGCTGCCCGCGGTAAAAGCCGATCCGCCACGGCAACACGGAACGATAAGAATCCCTGCATTATCAGGAATGaa                        >  2:572195/1‑133 (MQ=255)                                                                                                              tATCCGTTCCCCAACGACAAGCGTCATGGCTGGCTCCGTGCCGTTCTGAATATGTCCCTTCGCTGCCCGCGGTAAAAGCCGATCCGCCACGGCAACACGGAACGATAAGAATCCCTGCATTATCAGGAATGaa                        <  1:572195/133‑1 (MQ=255)                                                                                                              tATCCGTTCCCCAACGACAAGCGTCATGGCTGGCTCCGTGCCGTTCTGAATATGTCCCTTCGCTGCCCGCGGTAAAAGCCGATCCGCCACGGCAACACGGAACGATAAGAATCCCTGCATTATCAGGAATGAAGGGCAGt                 >  1:558294/1‑140 (MQ=255)                                                                                                                tCCGTTCCCCAACGACAAGCGTCATGGCTGGCTCCGTGCCGTTCTGAATATGTCCCTTCGCTGCCCGCGGTAAAAGCCGATCCGCCACGGCAACACGGAACGATAAGAATCCCTGCATTATCAGGAATGAAGGGCAGTaa               >  1:505724/1‑140 (MQ=255)                                                                                                                   gTTCCCCAACGACAAGCGTCATGGCTGGCTCCGTGCCGTTCTGAATATGTCCCTTCGCTGCCCGCGGTAAAAGCCGATCCGCCACGGCAACACGGAACGATAAGAATCCCTGCATTATCAGGAATGAAGGGCAGTAAttt            <  1:395616/140‑1 (MQ=255)                                                                                                                     tCCCCAACGACAAGCATCATGGCTGACTCCGTGCCGTCCTGAATATGTCCCTTCGCTGCCCGCGGTAAAAGCCGATCCGCCACGGCAACACGGAACGATAAGAATCCCTGCATTATCAGGAATGAAGGGCAGTAATTTcc          <  1:313617/140‑1 (MQ=255)                                                                                                                      ccccAACGACAAGCGTCATGGCTGGCTCCGTGCCGTTCTGAATATGTCCCTTCGCTGCCCGCGGTAAAAGCCGATCCGCCACGGCAACACGGAACGATAATAATCCCTGGATTATCAGGAATGCAGGGCAGTAATTTCCg         >  2:318046/1‑140 (MQ=255)                                                                                                                      ccccAACGACAAGCGTCATGGCTGGCTCCGTGCCGTTCTGAATATGTCCCTTCGCTGCCCGCGGTAAAAGCCGATCCGCCACGGCAACACGGAACGATAAGAATCCCTGCATTATCAGGAATGAAGGGCAGTAATTTCCg         <  2:550350/140‑1 (MQ=255)                                                                                                                          aaCGACAAGCGTCATGGCTGGCTCCGTGCCGTTCTGAATATGTCCCTTCGCTGCCCGCGGTAAAAGCCGATCCGCCACGGCAACACGGAACGATAAGAATCCCTGCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCa     <  2:582733/140‑1 (MQ=255)                                                                                                                           aCGACAAGCGTCATGGCTGGCTCCGTGCCGTTCTGAATATGTCCCTTCGCTGCCCGCGGTAAAAGCCGATCCGCCACGGCAACACGGAACGATAAGAATCCCTGCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCaa    <  2:606853/140‑1 (MQ=255)                                                                                                                           aCGACAAGCGTCATGGCTGGCTCCGTGCCGTTCTGAATATGTCCCTTCGCTGCCCGCGGTAAAAGCCGATCCGCCACGGCAACACGGAACGATAAGAATCCCTGCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCaa    <  2:276039/140‑1 (MQ=255)                                                                                                                             gACAAGCGTCATGGCTGGCTCCGTGCCGTTCTGAATATGTCCCTTCGCTGCCCGCGGTAAAAGCCGATCCGCCACGGCAACACGGAACGATAAGAATCCCTGCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCAata  <  1:620018/140‑1 (MQ=255)                                                                                                                                    |                                                                                                                                     GTCATTAAGTCAAATTCGCCTTGCATCCAGCATGCGCCGAGGAATTTGTTCTGCGGATTTTTTGCCAGTGCGGCTCGCGTTCTGCTGACTAAATCCTGGTACAGCGGAGTATCCGTTCCCCAACGACAAGCATCATGGCTGGCTCCGTGCCGTTCTGAATATGTCCCTTCGCTGCCCGCGGTAAAAGCCGATCCGCCACGGCAACACGGAACGATAAGAATCCCTGCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCAATA  >  NC_000913/4537038‑4537301

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