Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	536,173	T→C	L155P (CTC→CCC)	hyi →	hydroxypyruvate isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	536,173	0	T	C	97.1%	100.6 / NA	34	L155P (CTC→CCC)	hyi	hydroxypyruvate isomerase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (16/17);  minor base A (1/0);  total (17/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

CCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCTCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAACCAACACCATGACTCAGTGGGCTGATAAAATTGGTCACCTGCAAATTGCCGATAATCCGCATCGCGGCGAACCGGGAACCGGAGAAATTAATTATGATTATCTCTTTAAGGTAATCGAAAATTCTGACTACAAC  >  NC_000913/535930‑536411                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                               ccTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCTGCTGGTTTCAGCAGTGAACATATTCACGCAACGCTTGTAGAAATCCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCgg                                                                                                                                                                                                                                          <  1:479804/250‑1 (MQ=255)                      aagaGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGCAGATATTTTTTTACTGATTGAACCTATTAACCATTTTGATATTCCTGGGTTCCATCCCACCGGAACTCGGCAGGCGCTGAAattg                                                                                                                                                                                                                    >  2:27827/1‑251 (MQ=255)                        gagTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACAATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAattgat                                                                                                                                                                                                                  >  2:88989/1‑251 (MQ=255)                                           cagcagcGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTtgctgc                                                                                                                                                                                               <  1:522419/251‑1 (MQ=255)                                           cagcagcGATTCGTTATGCCCGCGCGCTGGGAAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTACATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTtgctgc                                                                                                                                                                                               >  1:426279/1‑251 (MQ=255)                                            agcagcGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACACGGCCGGCGCTGAAATTGATTGATGATGTTGGGtgcggc                                                                                                                                                                                               >  2:389830/1‑247 (MQ=255)                                                   aTTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGATCTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGaa                                                                                                                                                                                       <  1:303375/251‑1 (MQ=255)                                                     tCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGaaa                                                                                                                                                                                      <  1:745/250‑1 (MQ=255)                                                             cccGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTAt                                                                                                                                                                             <  1:433591/251‑1 (MQ=255)                                                                       gggTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGGTTCAGCAGTGAACAGATTCACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTTGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATc                                                                                                                                                                   <  2:70848/251‑1 (MQ=255)                                                                                  aaaTTAACTGTCTGGTCGGTAAAACGCCGGCTGGGTTCAGAAGTGAACAGATTCACGCAACGCTTGAAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGCAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTGCTGGTTTCCATCACCCCGGAACTCGGCAAGCGCTGAATATGATTGTTGATGTTGGTTGCTGCAATTTGCAAATTCAGAATGACATTAATCATATGCCGCgg                                                                                                                                                        >  2:406137/1‑251 (MQ=255)                                                                                           gTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGACAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGt                                                                                                                                               >  2:60917/1‑251 (MQ=255)                                                                                              tGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGACAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATgttgggt                                                                                                                                                                                                    >  2:431571/1‑193 (MQ=255)                                                                                              tGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATgttggtt                                                                                                                                                                                                    <  1:431571/195‑1 (MQ=255)                                                                                                gTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAAtt                                                                                                                                          >  1:265064/1‑251 (MQ=255)                                                                                                       aaaCGCCGGCTGGTTTCAGCAGTGAACAGATTCACGCCACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAACCAac                                                                                                                                   >  1:439492/1‑251 (MQ=255)                                                                                                                       cagcagTGAACAGATTCACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTACTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAACCAACACCATGACTCAGTggg                                                                                                                   <  1:48801/251‑1 (MQ=255)                                                                                                                          cagTGAACAGATTCACGACACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGACTATATTTTATTACTGATTGAACCTAGTAACCAGTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCGGCAATTTGAAAATTCAGTATGACATTTATCATAGTCAGCGGATGGAAGGTGAAGTAACCACCTCCATGACTCAGTGGGCTg                                                                                                                <  2:133923/251‑1 (MQ=255)                                                                                                                                       cACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGATAATTCAGTATGACATTTATACTATGCAGCGGATGGAAGGTGAATTAACAAACACCATGACTCAGTGGGCTGATAAAATTGGta                                                                                                    >  2:334224/1‑249 (MQ=255)                                                                                                                                                                gTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAACCAACACCATGACTCAGTGGGCTGATAAAATTGGTCACCTGCAAATTGCCGATAATCCGCAt                                                                          >  1:519250/1‑251 (MQ=255)                                                                                                                                                                   aTGCCGCGAATATGCTGATGAAAGAAGAGATTTTATTACTGATTGATCCTATTAACCATTTTGGTATTCCTGGTTTCCATCCGTCCGTAACTCGGCAGGCGCTGAAATGGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAACCAACACCATGACTCAGTGGTCTGATAAAATTGGTCATCTGCAAATTGCCTATAATCCGCATcgt                                                                       <  2:44514/251‑2 (MQ=255)                                                                                                                                                                          gAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAACCAACACCATGACTCAGTGGGCTGATAAAATTGGTCACCTGCAAATTGCCGATAATCCGCATCGCGGCgaac                                                                >  1:2698/1‑251 (MQ=255)                                                                                                                                                                                 cTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGCAGGTGAATTAACCAACACCATGACTCAGTGGGCTGATAAAATTGGTCACCTGCAAATTGCCGATAATCCGCATCGCGGCgaaccgggaac                                                         >  2:310589/1‑251 (MQ=255)                                                                                                                                                                                   gatgaAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAACCAACACCATGACTCAGTGGGCTGATAAAATTGGTCACCTGCAAATTGCCGATAATCCGCATCGCGGCgaaccgggaaccg                                                       >  2:120890/1‑251 (MQ=255)                                                                                                                                                                                      gaAAGAAGATATTTTATTACTGATTGCACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCt                                                                                                                                                                                                                           >  1:131384/1‑84 (MQ=255)                                                                                                                                                                                      gaAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCt                                                                                                                                                                                                                           <  2:131384/84‑1 (MQ=255)                                                                                                                                                                                         agaagaTATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCAACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAACCAACACCATGACTCAGTGGGCTGATAAAATTGGTCACCTGCAAATTGCCGATAATCCGCATCGCGGCGAACCGGGAACCGGAGaaa                                                 <  1:288131/251‑1 (MQ=255)                                                                                                                                                                                         agaacaTATTTTAGTACTGATTGAAACTCTTAACCATTTTGATATTACTGGTTTCCATCCCACCGGAACTCGTCAGGCGCTGAAATTGATTGATGTTGTTGGTCGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGCTGTAAGGTGAATTAACCAACACCATGACTCAGTGGGCAGATAAAATTGGTCACCTGCAAATTGCCGATAATCTGCATCGAGGCGAACCGGGAACCGgagac                                                  >  2:162874/1‑249 (MQ=255)                                                                                                                                                                                                                   cTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAACCAACACCATGACTCAGTGGGCTGATAAAATTGGTCACCTGCAAATTGCCGATAATCCGCATCGCGGCGCACCGGGAACCGGAAAAATTAATTATGATTATCTCTTTAAGGt                        >  2:62201/1‑250 (MQ=255)                                                                                                                                                                                                                     aTTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAACCAACACCATGACTCAGTGGGCTGATAAAATTGGTCACCTGCAAATTGCCGATAATCCGCATCGCGGCGAACCGGGAACCGGAGAAATTAATTATGATTATCTCTTTAAGGTAAt                     <  1:497552/251‑1 (MQ=255)                                                                                                                                                                                                                       tAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATTCAGCGGATGGAAGGTGAATTAACCAACACCATGACTCAGTGGGCTGATAAAATTTGTCACCTGCAAATTGCCGATAATCCGCATCGCGGCGAACCGGGAACCGGAGAAATTAATTATGATTATCTCTTTAAGGTAATc                    <  2:212906/250‑1 (MQ=255)                                                                                                                                                                                                                       tAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAACCAACACCATGACTCAGTGGGCTGATAAAATTGGTCACCTGCAAATTGCCGATAATCCGCATCGCGGCGAACCGGGAACCGGAGAAATTAATTATGATTATCTCTTTAAGGTAATc                    <  1:106083/250‑1 (MQ=255)                                                                                                                                                                                                                                 gATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAACCAACACCATGACTCAGTGGGCTGATAAAATTGGGCACCTGCAAATTGCCGATAATCCGCATCGCGGCGAACCGGGAACCGGAGAAATTAATTATGATTATCTCTTTAAGGTAATCGAAAATTCTGa         <  1:407664/251‑1 (MQ=255)                                                                                                                                                                                                                                          ggTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAACCAACACCATGACTCAGTGGGCTGATAAAATTGGTCACCTGCAAATTGCCGATAATCCGCATCGCGGCGAACCGGGAACCGGAGAAATTAATTATGATTATCTCTTTAAGGTAATCGAAAATTCTGACTACAAc  <  1:297396/249‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                               CCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCTCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAACCAACACCATGACTCAGTGGGCTGATAAAATTGGTCACCTGCAAATTGCCGATAATCCGCATCGCGGCGAACCGGGAACCGGAGAAATTAATTATGATTATCTCTTTAAGGTAATCGAAAATTCTGACTACAAC  >  NC_000913/535930‑536411

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