Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,418,862	C→T	G255D (GGT→GAT)	ulaG ←	L‑ascorbate‑6‑phosphate lactonase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,418,862	0	C	T	96.9%	101.8 / NA	32	G255D (GGT→GAT)	ulaG	L‑ascorbate‑6‑phosphate lactonase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (20/11);  total (20/12)
Fisher's exact test for biased strand distribution p-value = 3.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTGGTACTCGAAGTTGTCTTTATCCAGCGGCCAGGTAAATTTGCCACCCACCTGCCAGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGATCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATACCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCACCGCGCGATCGTCCATGCCATCTGGCAGTACGCCAGCCGCTTTTTGATCGGCAGGCAGGGTGATCAGTGCAGTACGGTCGAAAGCATCAAGCGCATGAATTTCAATGTCT  >  NC_000913/4418626‑4419107                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                      aGTGGTACTCGAAGTTGTCTTTATCCAGCGGCCAGGTAAATTTGCCACCCACCTGCCAGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGATCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGAGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCg                                                                                                                                                                                                                                         >  2:37915/1‑251 (MQ=255)                         ccAGCGGCCAGGTAAATTTGCCACCCACCTGCCAGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGATCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCg                                                                                                                                                                                                                 >  1:301489/1‑251 (MQ=255)                               gCCAGGTAAATTTGCCACCCACCTGCCAGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGATCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTg                                                                                                                                                                                                            >  2:15381/1‑250 (MQ=255)                                 cAGGTAAATTTGCCACCCACCTGCCAGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGATCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGGGGTGGAACGGGATCACTACTTTCGCATTCATCGCTGCACCCATACGCAGCATATCGGCGCTGGTCATGGTGTCGGTGATATCGCGCGGGTTTGCGCCGTACGATCCGGACGCCACGTCGATCTGATg                                                                                                                                                                                                         <  2:191197/251‑1 (MQ=255)                                    gTAAATTTGCCACCCACCTGCCAGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGATCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTc                                                                                                                                                                                                      <  1:22063/251‑1 (MQ=255)                                               acccaccTGCCAGATCAACGGCTTGAAGCCATACTTCTGGCGATCTTTTTTCATCTCCCACAGCACGCGGATCTCTTGCGGATCGGCCTGGACGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCCTACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGCTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCt                                                                                                                                                                                           >  2:154549/1‑251 (MQ=255)                                                      tGCCAGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGATCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCat                                                                                                                                                                                     >  1:23172/1‑250 (MQ=255)                                                           gATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGATCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTc                                                                                                                                                                                                      <  2:119157/228‑1 (MQ=255)                                                           gATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGATCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTc                                                                                                                                                                                                      >  1:119157/1‑228 (MQ=255)                                                                            cATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGATCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATAAGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTcg                                                                                                                                                              >  1:303215/1‑251 (MQ=255)                                                                                       cGATCTTTTTTCATCTCCCACAGCACGCGGATCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGAt                                                                                                                                                    >  2:55731/1‑250 (MQ=255)                                                                                                  cATCTCCCACAGCACGCGGATCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAgccgc                                                                                                                                        >  1:189950/1‑251 (MQ=255)                                                                                                        ccACAGCACGCGGATCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATCACTTCTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGATGAGTGGGAGTCTCCTCTGTGATTCAGTGAGTCTCCAGGCg                                                                                                                                  <  2:303215/251‑1 (MQ=255)                                                                                                                    gATCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGt                                                                                                                      >  2:183204/1‑251 (MQ=255)                                                                                                                      tctcTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGt                                                                                                                      >  1:134158/1‑249 (MQ=255)                                                                                                                                    gCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCACCGCGCGATc                                                                                                      >  1:145728/1‑251 (MQ=255)                                                                                                                                           aggttGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCACCGCGCGATCGTCCATg                                                                                               <  2:234217/248‑1 (MQ=255)                                                                                                                                                    aGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGTCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCACCGCGCGATCGTCCATGCCATCTGGc                                                                                      >  2:267941/1‑251 (MQ=255)                                                                                                                                                                      ggATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCACCGCGCGATCGTCCATGCCATCTGGCAGTACGCCAGCCGCtttt                                                                    >  1:89873/1‑251 (MQ=255)                                                                                                                                                                      ggATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCACCGCGCGATCGTCCATGCCATCTGGCAGTACGCCAGCCGCtttt                                                                    <  1:109860/251‑1 (MQ=255)                                                                                                                                                                        aTCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCACCGCGCGATCGTCCATGCCATCTGGCAGTACGCCAGCCGCTTTTTg                                                                  <  2:249990/251‑1 (MQ=255)                                                                                                                                                                                    gCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCACCGCGCGATCGTCCATGCCATCTGGCAGTACGCCAGCCGCTTTTTGATCggcaggc                                                        >  1:233236/1‑249 (MQ=255)                                                                                                                                                                                    gCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCACCGCGCGATCGTCCATGCCATCTGGCAGTACGCCAGCCGCTTTTTGATCGgcaggcag                                                      <  1:277264/251‑1 (MQ=255)                                                                                                                                                                                        tCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGACACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCACCGCGCGATCGTCCATGCCATCTGGCAGTACGCCAGCCGCTTTTTGATCGGCAGGCAggggg                                                  >  1:156526/1‑249 (MQ=255)                                                                                                                                                                                             gcTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCACCGCGCGATCGTCCATGCCATCTGGCAGTACGCCAGCCGCTTTTTGATCGGCAGGCAGGGTGATc                                               >  1:42094/1‑249 (MQ=255)                                                                                                                                                                                                 cACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCACCGCGCGATCGTCCATGCCATCTGGCAGTACGCCAGCCGCTCTTTGATCGGCAGGCAGGGTGATCAGTGc                                          >  1:90684/1‑250 (MQ=255)                                                                                                                                                                                                     cATACGCAGCATATCGGCGCTGGTCATGTTGTCGGTGATACCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCACCGCGCGATCGTCCATGCCATCTGGCAGTACGCCAGCCGCTCTTTGATCGGCAGGCAGGGTGATCAGTGCAGTAc                                     <  2:90684/251‑1 (MQ=255)                                                                                                                                                                                                        aCGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCACCGCGCGATCGTCCATGCCATCTGGCAGTACGCCAGCCGCTTTTTGATCGGCAGGCAGGGTGATCAGTGCAGTACGGt                                  <  1:277409/251‑1 (MQ=255)                                                                                                                                                                                                         cGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCACCGCGCGATCGTCCATGCCATCTGGCAGTACGCCAGCCGCTTTTTGATCGGCAGGCAGGGTGATCAGTGCAGTACGGTc                                 <  1:135828/251‑1 (MQ=255)                                                                                                                                                                                                                         tGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCACCGCGCGATCGTCCATGCCATCTGGCAGTACGCCAGCCGCTTTTTGATCGGCAGGCAGGGGGATCAGTGCAGTACGGTCGAAAGCATCAAgcgc                  >  2:50489/1‑250 (MQ=255)                                                                                                                                                                                                                          ggTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCACCGCGCGATCGTCCATGCCATCTGGCAGTACGCCAGCCGCTTTTTGATCGGCAGGCAGGGTGATCAGTGCAGTACGGTCGAAAGCATCAAGCGCAt                >  1:152250/1‑251 (MQ=255)                                                                                                                                                                                                                                         gATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCACCGCGCGATCGTCCATGCCATCTGGCAGTACGCCAGCCGCTTTTTGATCGGCAGGCAGGGTGATCAGTGCAGTACGGTCGAAAGCATCAAGCGCATGAATTTCAATGTCt  <  1:95942/250‑1 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                      AGTGGTACTCGAAGTTGTCTTTATCCAGCGGCCAGGTAAATTTGCCACCCACCTGCCAGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGATCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATACCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCACCGCGCGATCGTCCATGCCATCTGGCAGTACGCCAGCCGCTTTTTGATCGGCAGGCAGGGTGATCAGTGCAGTACGGTCGAAAGCATCAAGCGCATGAATTTCAATGTCT  >  NC_000913/4418626‑4419107

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