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%	194.8 / NA	56	G92D (GGC→GAC)	nanM	N‑acetylneuraminic acid mutarotase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (27/29);  total (27/29)

GTGTACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATCCAGCTTGTACCATGCCGTACCTGCG  >  NC_000913/4538358‑4538625                                                                                                                                      |                                                                                                                                       gtgtACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGATTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCa                                                                                                                                   <  1:541829/139‑1 (MQ=255) gtgtACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCa                                                                                                                                   <  1:915684/139‑1 (MQ=255) gtgtACAAAAGTAACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCa                                                                                                                                   <  1:452829/139‑1 (MQ=255)            tCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAg                                                                                                                        >  1:295335/1‑139 (MQ=255)             cacaTGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGa                                                                                                                       <  2:64995/139‑1 (MQ=255)                        ccATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAAt                                                                                                            <  2:818780/139‑1 (MQ=255)                                    gcgcGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGt                                                                                                <  2:448988/139‑1 (MQ=255)                                       cgTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGc                                                                                             >  2:412264/1‑139 (MQ=255)                                                       aaCCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCg                                                                             <  1:73250/139‑1 (MQ=255)                                                        aCCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGc                                                                            >  2:275100/1‑139 (MQ=255)                                                          ccAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCa                                                                          <  2:379485/139‑1 (MQ=255)                                                          ccAACTATTGATTTTTGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCa                                                                          <  2:353716/139‑1 (MQ=255)                                                              cTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGcc                                                                                                                  >  2:77655/1‑95 (MQ=255)                                                              cTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGcc                                                                                                                  <  1:77655/95‑1 (MQ=255)                                                                aTTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGCCCGCCAGGGAAt                                                                    <  1:631342/139‑1 (MQ=255)                                                                  tGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGGAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGc                                                                  <  1:223976/139‑1 (MQ=255)                                                                              tGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGt                                                      >  2:292662/1‑139 (MQ=255)                                                                                     gtgtACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATctct                                                                                    >  2:717497/1‑102 (MQ=255)                                                                                     gtgtACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATctct                                                                                    <  1:717497/102‑1 (MQ=255)                                                                                       gtACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGa                                                                                                                       >  1:209838/1‑65 (MQ=255)                                                                                       gtACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGa                                                                                                                       <  2:209838/65‑1 (MQ=255)                                                                                       gtACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCAtttttt                                             >  1:82122/1‑139 (MQ=255)                                                                                        tACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTa                                            >  2:868639/1‑139 (MQ=255)                                                                                          cGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGcc                                                                                                                  <  2:594993/67‑1 (MQ=255)                                                                                          cGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGcc                                                                                                                  >  1:594993/1‑67 (MQ=255)                                                                                          cGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATc                                          >  1:386110/1‑139 (MQ=255)                                                                                           gTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAAt                                                                                                            <  1:242011/72‑1 (MQ=255)                                                                                           gTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAAt                                                                                                            >  2:242011/1‑72 (MQ=255)                                                                                           gTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCt                                         <  2:82122/139‑1 (MQ=255)                                                                                             cATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCaaa                                                                                                                                 <  1:652908/49‑1 (MQ=255)                                                                                             cATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCaaa                                                                                                                                 >  2:652908/1‑49 (MQ=255)                                                                                              aTTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATgtcgc                                                                                                                                     <  2:799634/44‑1 (MQ=255)                                                                                              aTTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATgtcgc                                                                                                                                     >  1:799634/1‑44 (MQ=255)                                                                                                tAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAAt                                                                                                            >  1:289740/1‑67 (MQ=255)                                                                                                tAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAAt                                                                                                            <  2:289740/67‑1 (MQ=255)                                                                                                tAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGc                                    >  2:456245/1‑139 (MQ=255)                                                                                                       cTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAAcaca                                                                                                                             >  2:632924/1‑43 (MQ=255)                                                                                                       cTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAAcaca                                                                                                                             <  1:632924/43‑1 (MQ=255)                                                                                                       cTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAAt                                                                    <  2:548773/100‑1 (MQ=255)                                                                                                       cTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAAt                                                                    >  1:548773/1‑100 (MQ=255)                                                                                                       cTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTa                             >  2:54221/1‑139 (MQ=255)                                                                                                       cTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTa                             >  1:228966/1‑139 (MQ=255)                                                                                                        tGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTAt                            <  2:255596/139‑1 (MQ=255)                                                                                                          aGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATcc                          >  2:25096/1‑139 (MQ=255)                                                                                                            tCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATCCAg                        >  2:915878/1‑139 (MQ=255)                                                                                                               aGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATCCAGCtt                     >  1:193483/1‑139 (MQ=255)                                                                                                               aGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATCCAGCtt                     <  1:456245/139‑1 (MQ=255)                                                                                                                   cTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATCCAGCTTGTAc                 <  1:464041/139‑1 (MQ=255)                                                                                                                     cGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATCCAGCTTGTACCa               >  2:337987/1‑139 (MQ=255)                                                                                                                     cGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATCCAGCTTGTACCa               >  2:25382/1‑139 (MQ=255)                                                                                                                     cGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATCCAGCTTGTACCa               <  2:228966/139‑1 (MQ=255)                                                                                                                     cGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATCCAGCTTGTACCa               <  1:915878/139‑1 (MQ=255)                                                                                                                     cGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATCCAGCTTGTACCa               >  1:370783/1‑139 (MQ=255)                                                                                                                          tttttGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATCCAGCTTGTACCATGCCg          >  2:404745/1‑139 (MQ=255)                                                                                                                            tttGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATCCAGCTTGTACCATGCCGTa        >  1:135791/1‑139 (MQ=255)                                                                                                                             ttGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATCCAGCTTGTACCATGCCGTAc       <  1:868639/139‑1 (MQ=255)                                                                                                                                  aaTGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATCCAGCTTGTACCATGCCGTACCTgcg  <  2:37864/139‑1 (MQ=255)                                                                                                                                  aaTGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATCCAGCTTGTACCATGCCGTACCTgcg  <  1:337987/139‑1 (MQ=255)                                                                                                                                      |                                                                                                                                       GTGTACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATCCAGCTTGTACCATGCCGTACCTGCG  >  NC_000913/4538358‑4538625

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