Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,538,491	C→T	G92D (GGC→GAC)	nanM ←	N‑acetylneuraminate 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	97.1%	111.4 / NA	34	G92D (GGC→GAC)	nanM	N‑acetylneuraminate mutarotase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (18/15);  minor base G (0/1);  total (18/16)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATAGCGGTTGAATCTTTTCCAGCCTCGTTGAGATCTTCAAAATAGCCATTGAAGATATTCTGGTTA‑ACACCACCAGTAACATAAGCCTTGCCGTTGTGTACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATCCAGCTTGTACCATGCCGTACCTGCGCTACCTAAACCAATGTAGACAGTGTCGTTATCAATTGCTCCGGTACCACTTTTAAATGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGC  >  NC_000913/4538263‑4538729                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                               ataGCGTTTGAATCTTTTCCAGCCTCGTTGAGATCTTCAAAATAGCCATTGAAGATATTCTGGTTA‑ACCCCACCAGTAACATAAGCCTTGCCGTTGTGTACAAAATTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGcc                                                                                                                                                                                                                          <  1:119472/251‑1 (MQ=255)                          cGTTGAGATCTTCAAAATAGCCATTGAAGATATTCTGGTTA‑ACACCACCAGTAACATAAGCCTTGCCGTTGTGTACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTg                                                                                                                                                                                                  >  2:101680/1‑250 (MQ=255)                                         aagaGCCATGGAAGATATTCGGGTTA‑ACCCCACCAGTAACATAAGCCGTGCCGTTGTGTACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACGATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCa                                                                                                                                                                                  <  2:222304/248‑1 (MQ=255)                                            aGCCATTGAAGATATTCTGGTTA‑ACACCACCAGTAACATAAGCCTTGCCGTTGTGTACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAggg                                                                                                                                                                               >  1:135849/1‑251 (MQ=255)                                            aGCCATTGAAGATATT‑TTGTTACATCCCACCAGTAACTTAAGCCGTGCCGGTGTGTACAAAAGGCACAGGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGTTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTTCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTTCTTGATCTCTTGGTCCGCCAggg                                                                                                                                                                               <  2:377033/251‑1 (MQ=255)                                                         atTCTGGTTA‑ACACCACCAGAAACATAAGCCTTGCCGTTGTGTACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACg                                                                                                                                                                  <  1:365959/251‑1 (MQ=255)                                                          attcttggtaaacccaccgGTAAAATAAGCCGTGCTGTTGTTTACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTGGCTTGATCTCTTGGGCCGCCAGGGAATGCAGCTAACg                                                                                                                                                                  <  2:365598/238‑1 (MQ=255)                                                                  a‑aCACCACCAGTAACATAAGCCTTGCCGTTGTGTACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCAtt                                                                                                                                                         >  2:57161/1‑251 (MQ=255)                                                                    ccccccccgttAACATAAGCCTGTCCGTTTTGTCCAAAGTTCACATGCCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTTTTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTTATCGCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCAttt                                                                                                                                                        <  2:294748/241‑1 (MQ=255)                                                                        caccaGTAACATAAGCCTTGCCGTTGTGTACACAAGACACATCGCCAGCCATACCCATCGGCGCGTGCCACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGAATGCCATCAATAAATGCAGAGGTTTCTTGATCTCTTGGTCCGCCAGGGAATGCAGCAAACGCTGTCCAtttttt                                                                                                                                                     >  1:231720/1‑250 (MQ=255)                                                                                             cGTTGTGTACAAAAGTAACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCGGTCCATTTTTTATCTTTGGCCTGTGTATCCAGc                                                                                                                               >  2:214196/1‑251 (MQ=255)                                                                                                tgtgTACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATCCAGCTTg                                                                                                                            >  1:283657/1‑251 (MQ=255)                                                                                                   gtaaaaaaggcaaggggCCCGCCATGCCCAGCGGCGCGTGCGACATCAAATTAACCCAACTAGTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCGCGCTGTTTTTGCCAATGGCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCTGGGAATGCAGCTAACGCTGTCAATTTTTTATCTTTGGCCTGTGTATCCAGCTTGTa                                                                                                                          <  2:228325/235‑1 (MQ=255)                                                                                                         aaGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATCCAGCTTGTACCATGcc                                                                                                                   <  2:135849/251‑1 (MQ=255)                                                                                                                aTGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATCCAGCTTGTACCATGCCGTACCTg                                                                                                            <  1:82017/251‑1 (MQ=255)                                                                                                                aTGGCCCGCCATGCCCATCGGCGCGGGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTAGCTTTGGCCTGTGTATCCAGCTTGTACCATGCCGTACCTg                                                                                                            >  2:293800/1‑251 (MQ=255)                                                                                                                                            gACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGTAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATCCAGCTTGTACCATGCCGTACCTGCGCTACCTAAACCAATGTAGACAGTGTc                                                                                >  2:265278/1‑251 (MQ=255)                                                                                                                                               aTCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATCCAGCTTGTACCATGCCGTACCTGCGCTACCTAAACCAATGTAGACAGTGTCGtt                                                                             <  2:209516/251‑1 (MQ=255)                                                                                                                                                 cAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTCCGTCATTAAATACCCGCGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATCCAGCTTGTACCATGCCGTACCTGCGCTACCTAAACCAATGTAGACAGTGTCGTTAt                                                                           >  2:128295/1‑251 (MQ=255)                                                                                                                                                     ttAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTT‑ATCTTTGGCCTGTGTATCCAGCTTGTACCATGCCGTACCTGCGCTACCTAAACCCATTTAGACAGGGTCGTTATCCAtt                                                                      >  1:301223/1‑251 (MQ=255)                                                                                                                                                      tAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGAAATTACATCCCTGAGTCAAGACCTAGCTGTTTTTGCCAATGTCACCACACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATCCAGCTTGTACCATGCCGTACCTGCGCCACCTAAACCAATGTAGACAGTGTCGTTATCAAtt                                                                      >  1:43372/1‑251 (MQ=255)                                                                                                                                                      tAACCCAACTATTGGGTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGACTGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATCCAGCTTGTACCATGCCGTACCTGCGCTACCTAAACCAATGTAGACAGTGTCGTTATCAAtt                                                                      >  1:144877/1‑251 (MQ=255)                                                                                                                                                        aCCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATCCAGCTTGTACCATGCCGTACCTGCGCTACCTAAACCAATGTAGACAGTGTCGTTATCAATTGc                                                                    <  1:101680/251‑1 (MQ=255)                                                                                                                                                             aCTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAACTACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATCCAGCTTGTACCATGCCGTACCTGCGCTACCTAAACCAATGTAGACAGTGTCGTTATCAATTGCTCCgg                                                               >  1:254511/1‑251 (MQ=255)                                                                                                                                                                aTTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCATGCCATCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATCCAGCTTGTACCATGCCGTACCTGCTCTACCTCAACCAATGTAGACAGTGTCGTTATCAATTGCTCCGGTAc                                                            >  2:317697/1‑251 (MQ=255)                                                                                                                                                                 ttGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATCCAGCTTGTACCATGCCGTACCTGCGCTACCTAAACCAATGTAGACAGTGTCGTTATCAATTGCTCCGGTAcc                                                           >  1:263416/1‑251 (MQ=255)                                                                                                                                                                         ggggTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATCCAGCTTGTACCATGCCGTACCTGCGCTACCTAAACCAATGTAGACAGTGTCGTTATCAATTGCTCCGGTACCACTTTTa                                                    <  1:316919/250‑1 (MQ=255)                                                                                                                                                                                 aTTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATCCAGCTTGTACCATGCCGTACCTGCGCTACCTAAACCAATGTAGACAGTGTCGTTATCAATTGCTCCGGTACCACTTTTAAATGGcac                                            >  1:35850/1‑250 (MQ=255)                                                                                                                                                                                 aTTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATCCAGCTTGTACCATGCCGTACCTGCGCTACCTAAACCAATGTAGACAGGGTCGTTATCAATTGCTCCGGTACCACTTTTAAATGGcaca                                           <  2:119950/251‑1 (MQ=255)                                                                                                                                                                                     gtgtACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATCAAGCTTGTACCATGCCGTACCTGCGCTACCTAAACCAATGTAGACAGTGTCGTTATCAATTGCTCCGGTACCACTTTTAAATGGCACAgg                                         >  2:165577/1‑249 (MQ=255)                                                                                                                                                                                                 aaaTACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATCCAGCTTGTACCATGCCGTACCTGCGCTACCTAAACCAATGTAGACAGTGTCGTTATCAATTGCTCCGGTACCACTTTTAAATGGCACAGGAGTTTCCGGTAAtt                           <  2:278988/251‑2 (MQ=255)                                                                                                                                                                                                              aaGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATCCAGCTTGTACCATGCCGTACCTGCGCTACCTAAACCAATGTAGACAGTGTCGTTATCAATTGCTCCGGTACCACTTTTAAATGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGc              <  1:219234/251‑1 (MQ=255)                                                                                                                                                                                                                      gCTGTTTTTGCAAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATCCAGCTTGTACCATGCCGTACCTGCGCTACCTAAACCAATGTAGACAGTGTCGTTATCAATTGCTCCGGTACCACTTTTAAATGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATg      >  2:371162/1‑251 (MQ=255)                                                                                                                                                                                                                           ttttGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATCCAGCTTGTACCATGCCGTACCTGCGCTACCTAAACCAATGTAGACAGTGTCGTTATCAATTGCTCCGGTACCACTTTTAAATGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGc  <  1:57161/250‑1 (MQ=255)                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                               ATAGCGGTTGAATCTTTTCCAGCCTCGTTGAGATCTTCAAAATAGCCATTGAAGATATTCTGGTTA‑ACACCACCAGTAACATAAGCCTTGCCGTTGTGTACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATCCAGCTTGTACCATGCCGTACCTGCGCTACCTAAACCAATGTAGACAGTGTCGTTATCAATTGCTCCGGTACCACTTTTAAATGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGC  >  NC_000913/4538263‑4538729

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