Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,297,782	G→A	P529S (CCT→TCT)	fdhF ←	formate dehydrogenase H

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,297,782	0	G	A	93.9%	87.4 / ‑3.6	34	P529S (CCT→TCT)	fdhF	formate dehydrogenase H
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (19/12);  total (21/13)
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.24e-01

GATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTCAACCGCCTTGAAGAAACGCTGGAAGCCACGGTCAGCCGCAGTAAACACGCCTTCATGCTCGCCCCACG  >  NC_000913/4297537‑4298027                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      gATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCAc                                                                                                                                                                                                                                                   <  1:371481/250‑1 (MQ=255)            aGCCAGGGTCATCAGCCAGCGCCGCCAGAGCCGCACAGTTACCGGTCATCGAAAGGCAAGAGTAGTGACCAAATTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGCTCGGGGTATCAAACTTCTCTTTACACCGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGAGGCAAGCCCCCTGAATGAAggcc                                                                                                                                                                                                                                       >  1:230690/1‑248 (MQ=255)                  gtgcatcaGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCg                                                                                                                                                                                                                                 <  1:13574/248‑1 (MQ=255)                                       gTGCCGCAAAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTATCTCGTAAGTCGCa                                                                                                                                                                                                            <  1:357337/251‑1 (MQ=255)                                                  ttACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGACATc                                                                                                                                                                                                 >  1:225875/1‑251 (MQ=255)                                                                aCGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGAc                                                                                                                                                                                   <  1:84255/251‑1 (MQ=255)                                                                          tagtGACCCACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGGGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGCAGAATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTTAGTCGCACCGTAGAAATCCGGGCACAGGTGACGCAACCCATc                                                                                                                                                                         >  2:203315/1‑251 (MQ=255)                                                                                       cACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTg                                                                                                                                                             <  1:379862/250‑1 (MQ=255)                                                                                                   aCAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGATCCGTAGAAATCCGGGCACAGATGACGCATCTCATCCCAGATCTACTGGGcgttgttgt                                                                                                                                                  >  2:205262/1‑249 (MQ=255)                                                                                                           aTCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCgg                                                                                                                                        >  1:105444/1‑251 (MQ=255)                                                                                                              gggTGCTCGTCGGTGAGTTTGTCGATGAGCGCTCCCCAGGCGCAGGTTAAGAACTGCGCCAGACCTTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCTGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCTGGCACAGATGACGCAACTCATCCCAGATCTTCTGGGTGTTGTTGTAGTGCATCGGATa                                                                                                                                     <  2:261629/251‑1 (MQ=255)                                                                                                                gTACTCGGCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTGCTCGGAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTTTTGTTGTAGTGCATCGGATAAc                                                                                                                                   <  1:378634/251‑1 (MQ=255)                                                                                                                     cgtcgGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGCACTGCACCAGACCGTTCGGGGTATCAACCTACTATTTAAACAAATAAAAAGTCCCCTGATCGGCATCTGAAGCATCGCCGCAAGACCACTGACAGAAGCCCACTTCGCCCATTTTCTCGAACGCCGCACCGAAGAAATCCGGCAACAGATGACGACACTCATCCCATATATCCTGAGCGATGTTGTCGTGCATCAGACAAccccg                                                                                                                               >  2:32261/1‑248 (MQ=255)                                                                                                                        cgGTGAGTTCGTCGATTGGCGCTACCCAGTCTCAGTTGCAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTACCCTGATCGGCATCTGAAGTCTCGCGGCAAGACCACTGAATGAAGCCCCGTTCGCCCATTTTCTCGTCAGTCGAACCGTAGAAATCCGGTCACCGATGACGCAACTCATCCCAGATCACCTGGGTGTTGTTGAAGTTCATCGGATACTCCATACgg                                                                                                                           >  2:343414/1‑251 (MQ=255)                                                                                                                            gAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTgg                                                                                                                       >  2:109107/1‑251 (MQ=255)                                                                                                                                     gATTGGCGCTACCCAGTCGCAGGAGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCCCAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGAtttaac                                                                                                              >  2:382259/1‑248 (MQ=255)                                                                                                                                      aTTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTAGCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCGCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACt                                                                                                             <  1:72540/251‑1 (MQ=255)                                                                                                                                         ggCGCTACCCAGTCGCAGGTGAAGAACTTCGCCATACCTGTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCTGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCTCCGTAGAAATCCGTGCACAGATGACTCAACTCTTCCCAGAGCTCCTGTGGGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACtga                                                                                                           <  2:187397/250‑1 (MQ=255)                                                                                                                                                cccAGACGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGGTGTAGTGCATCGGATAACCCATGCGGGTGGCGATTTCACTGATGATTTGc                                                                                                   <  1:104899/251‑1 (MQ=255)                                                                                                                                                  cAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCa                                                                                                 <  1:203315/251‑1 (MQ=255)                                                                                                                                                  cAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCa                                                                                                 <  1:109107/251‑1 (MQ=255)                                                                                                                                                         aGGTGAAGAACTGCGCCAAACCGTTCGGGGAATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGCGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGtt                                                                                          >  1:265298/1‑251 (MQ=255)                                                                                                                                                            tGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTc                                                                                       >  1:194359/1‑251 (MQ=255)                                                                                                                                                                              cGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTAGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGCTCCCATTTCGGTTc                                                                      >  1:365346/1‑250 (MQ=255)                                                                                                                                                                                ttCGGGGTATCAAACTTCTCTTTAAACAGATACGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCACGACCAATGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCCAATCATCCAAGATCTCCTGGGTGTTGTTGTAGTGCATCGGCTAACCCATACGGGTGGCGATTTCACTGGTGCTTTGCCAGTCCGTTTTCAGATCCACTTTCGGTTCacc                                                                   >  2:261962/1‑249 (MQ=255)                                                                                                                                                                                  cGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTCAAcc                                                                  <  1:31145/250‑1 (MQ=255)                                                                                                                                                                                         tCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGTTCCCATTTCGGTTCAACCGCCTTgc                                                           >  1:59502/1‑249 (MQ=255)                                                                                                                                                                                                                       cccTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGGGCATCGGCTAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTCAACCGCCTTGAAGCAACGCTGGAAGCCACGGTAAGCCGCAGt                            >  2:135652/1‑251 (MQ=255)                                                                                                                                                                                                                        ccTGATCGTAATCTGAAGTATCGCGGACAGACCACTGAATGAAGCCCAGTTCGCCCATTTTATAGTACGTCGCACCGTCGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATCCGGATGGCGGTTTCCCTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTAAAACGCCTTAGAGAAACGATGGAAGCAACgacacgccacagt                            >  2:353499/1‑238 (MQ=255)                                                                                                                                                                                                                         cTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCAAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCCCAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTACGTTTTCAGATCCCATTTCGGTTCAACCGCCTTGAAGAAACGCTGGAAGCCGCGGTCAGCCGCCGTa                           >  1:209575/1‑250 (MQ=255)                                                                                                                                                                                                                            aTCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTCAACCGCCTTGAAGAAACGCTGGAAGCCACGGTCAGCCGCAGTAAac                        >  1:327197/1‑250 (MQ=255)                                                                                                                                                                                                                                      gcaGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCCTTTCGGTTCAACCGCCTTGCAGAAACGCTGGAAGCCACGGTCAGCCGCAGTAAACACGCCTTCATg             >  1:123702/3‑251 (MQ=255)                                                                                                                                                                                                                                      gAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTCAACCGCCTTGAAGAAACGCTGGAAGCCACGGTCAGCCGCAGTAAACACGCCTTCATg             >  1:123698/1‑251 (MQ=255)                                                                                                                                                                                                                                                 gCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTCAACCGCCTTGAAGAAACGCTGGAAGCCACGGTCAGCCGCAGTAAACACGCCTTCATGCTCGCCCCacg  >  1:374467/1‑251 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      GATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTCAACCGCCTTGAAGAAACGCTGGAAGCCACGGTCAGCCGCAGTAAACACGCCTTCATGCTCGCCCCACG  >  NC_000913/4297537‑4298027

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