Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,797,855	C→G	11.1%	V47L (GTG→CTG)	vdh ←	vanillin dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,797,855	0	C	G	11.1%	130.3 / 7.7	45	V47L (GTG→CTG)	vdh	vanillin dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (14/26);  new base G (3/2);  total (17/28)
Fisher's exact test for biased strand distribution p-value = 3.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.70e-01

AGTTGGCCATGGCCCCGGTTTCACCGGCGGCGGCGATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCA‑GGGCGGCCCAGGCCGGAAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGGCGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCCATCGCTGGCGGGGCGCGACTGCCCGCCAATCAGCAAAGGCACCTGCAAC  >  NC_002947/3797708‑3797991                                                                                                                                                     |                                                                                                                                         aGTTGGCCATGGCCCCGGTTTCACCGGCGGCGGCGATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCA‑GGGCGGCCCAGGCCGGAAACGCCGCGCTGGCAGCAGCCAcc                                                                                                                                         >  2:252774/1‑149 (MQ=255)  gtgggCCATGGCCCCGGTTGCACCGGCGGCGGCGATGAACTCGGCGGCCCTCGCGTGCAACAGATCATCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGCTGCCA‑GGGCGGCCCAGGCCGGAAACGCCGCGCTGGCAGCAGCCACcg                                                                                                                                        <  4:139913/146‑1 (MQ=255)   ttGGCCATGGCCCCGGTTTCACCGGCGGCGGCGATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCA‑GGGCGGCCCAGGCCGGAAACGCCGCGCTGGCAGCAGCCACcgc                                                                                                                                       <  3:107478/149‑1 (MQ=255)        cATGGCCCCGGTTTCACCGGCGGCGGCGATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCA‑GGGCGGCCCAGGCCGGAAACGCCGCGCTGGCAGCAGCCACCGcggcat                                                                                                                                  >  3:72146/1‑149 (MQ=255)          tGGCCCCGGTTTCACCGGCGGCGGCGATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCA‑GGGCGGCCCAGGCCGGAAACGCCGCGCTGGCAGCAGCCACCGcggcatcg                                                                                                                                <  7:10594/149‑1 (MQ=255)               ccGGTTTCACCGGCGGCGGCGATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCA‑GGGCGGCCCAGGCCGGAAACGCCGCGCTGGCAGCAGCCACCGCggcatcggcatc                                                                                                                           <  3:115991/149‑1 (MQ=255)                cGGTTTCACCGGCGGCGGCGATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCA‑GGGCGGCCCAGGCCGGAAACGCCGCGCTGGCAGCAGCCACCGCGgcatcggcatcg                                                                                                                          <  7:14245/149‑1 (MQ=255)                  gTTTCACCGGCGGCGGCGATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCA‑GGGCGGCCCAGGCCGGAAACGCCGCGCTGGCAGCAGCCACCGCGGCatcggcatcggc                                                                                                                        <  1:251359/149‑1 (MQ=255)                          ggcggcggcGATTAACTCGGCGGCCCTCGCCGGCAACACTTCAGCGCCTGCAAGCAAGGGGCTGCGCCGCTCGCCCGGTGCCATCGGCGGCCCAGGCCGGAAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGtgtg                                                                                                                 <  1:252774/149‑1 (MQ=255)                                ggcgATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCA‑GGGCGGCCCAGGCCGGAAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTggcggcgg                                                                                                          <  7:52386/149‑1 (MQ=255)                                ggcgATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCA‑GGGCGGCCCAGGCCGGAAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTggcggcgg                                                                                                          <  3:186182/149‑1 (MQ=255)                                  cgATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCA‑GGGCGGCCCAGGCCGGAAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGGCGGCa                                                                                                        <  8:164701/149‑1 (MQ=255)                                                gcCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCA‑GGGCGGCCCAGGCCGGAAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGGCGGCAGCCTGCGAcaccac                                                                                          <  4:72146/149‑1 (MQ=255)                                                  ccTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCA‑GGGCGGCCCAGGCCGGAAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGGCGGCAGCCTGCGACACCACCt                                                                                        >  8:205782/1‑149 (MQ=255)                                                       ccTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCA‑GGGCGGCCCAGGCCGGAAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGGCGGCAGCCTGCGACACCACCTCGCCa                                                                                   <  3:59901/149‑1 (MQ=255)                                                               aGATCAGTGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCA‑GGGCGGCCCAGGCCGGAAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGGCGGCAGCCTGCGACACCACCTCGCCAGTCACCgg                                                                           <  5:195425/149‑1 (MQ=255)                                                                    aGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCG‑GGGCGGCCCAGGCCGGAAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGGCGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTa                                                                       >  7:186420/1‑148 (MQ=255)                                                                       gcCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCA‑GGGCGGCCCAGGCCGGAAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGGCGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAgc                                                                   >  6:52007/1‑149 (MQ=255)                                                                         cTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCA‑GGGCGGCCCAGGCCGGAAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGGCGGTAGCCTGCGACACCACCTCGCCAGTCACCg                                                                            <  8:56583/138‑1 (MQ=255)                                                                         cTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCA‑GGGCGGCCCAGGCCGGAAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGGCGGTAGCCTGCGACACCACCTCGCCAGTCACCg                                                                            >  7:56583/1‑138 (MQ=255)                                                                            ccAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCA‑GGGCGGCCCAGGCCGGAAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGGCGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTc                                                               >  4:224045/1‑148 (MQ=255)                                                                            ccAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCA‑GGGCGGCCCAGGCCGGAAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGGCGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTc                                                               >  4:156332/1‑148 (MQ=255)                                                                                   gCGGCTGCGCCGCTCGCCCGGTGCCG‑GGGCGGCCCAGGCCGGAAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGGCGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTc                                                       <  8:186420/149‑1 (MQ=255)                                                                                   gCGGCTGCGCCGCTCGCCCGGTGCCA‑GGGCGGCCCAGGCCGGAAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGGCGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTc                                                       <  3:224045/149‑1 (MQ=255)                                                                                       cTGCGCCGCTCGCCCGGTGCCA‑GGGCGGCCCAGGCCGGAAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGGCGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTcc                                                   <  2:280784/149‑1 (MQ=255)                                                                                             cgcTCGCCCGGTGCCA‑GGGCGGCCCAGGCCGGAAACGCCGCGCTGGCAGCCGCCACCGCGGCCTCGGCATCGGCCCGTGTGGCGGCGGCCGCCTGCGACCCCACCTCGCCCGGGCCCGGGGTACAGCGCGCGACGGTTCGTCCCTCGCt                                             >  5:49688/1‑149 (MQ=255)                                                                                              gcTCGCCCGGTGCCA‑GGGCGGCCCAGGCCGGAAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGGCGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCCATCGCTg                                            >  1:258310/1‑149 (MQ=255)                                                                                              gcTCGCCCGGTGCCA‑GGGCGGCCCAGGCCGGAAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGGCGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCCATCGCTg                                            >  1:76217/1‑149 (MQ=255)                                                                                                  gCCCGGTGCCA‑GGGCGGCCCAGGCCGGAAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGGCGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCCATCGCTGGCgg                                        <  2:258310/149‑1 (MQ=255)                                                                                                   cccGGTGCCA‑GGGCGGCCCAGGCCGGAAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGGCGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCCATCGCTGGCggg                                       <  8:261270/149‑1 (MQ=255)                                                                                                    ccGGTGCCA‑GGGCGGCCCAGGCCGGAAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGGCGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGaa                                                            >  4:304660/1‑127 (MQ=255)                                                                                                    ccGGTGCCA‑GGGCGGCCCAGGCCGGAAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGGCGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGaa                                                            <  3:304660/127‑1 (MQ=255)                                                                                                          ccA‑GGGCGGCCCAGGCCGGAAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGGCGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCCATCGCTGGCGGGGCGCGa                                 <  2:76217/148‑1 (MQ=255)                                                                                                                      cAGGCCGGAAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGGCGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCCATCGCTGGCGGGGCGCGACTGCCCGCCAAt                     <  3:156332/149‑1 (MQ=255)                                                                                                                       aGGCCGGAAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGGCGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCCATCGCTGGCGGGGCGCGACTGCCCGCCAAt                     >  3:216938/1‑148 (MQ=255)                                                                                                                        ggCCGGAAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGGCGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCCATCGCTGGCGGGGCGCGACTGCCCGCCAATCa                   <  4:216938/149‑1 (MQ=255)                                                                                                                        ggCCGGAAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGGCGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCCATCGCTGGCGGGGCGCGACTGCCCGCCAATCa                   <  5:52007/149‑1 (MQ=255)                                                                                                                             gAAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGGCGGCAGCCTGCTACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCCATCGCTGGCGGGGCGCGACTGCCCGCCAATCAGCaaa              <  7:205782/149‑1 (MQ=255)                                                                                                                             gAAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGGCGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCCATCGCTGGCGGGGCGCGACTGCCCGCCAATCAGCa                >  2:313108/1‑147 (MQ=255)                                                                                                                               gatgccgcGGTGGCTGCTGCCAGCGCGGCATCGGCATCGGCCAGTGTGGCGGCGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCCATCGCTGGCGGGGCGCGACTGCCCGCCAATCAGCAAAgg            >  6:243804/4‑149 (MQ=255)                                                                                                                                    cgcgCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGGCGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCCATCGCTGGCGGGGCGCGACTGCCCGCCAATCAGCAAAGGCACCt       >  8:263903/1‑149 (MQ=255)                                                                                                                                     gcggTGGCTGCTGCCAGCGCGGCATCGGCATCGGCCAGTGTGGCGGCGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCCATCGCTGGCGGGGCGCGACTGccc                           >  7:260531/1‑128 (MQ=255)                                                                                                                                     gcggTGGCTGCTGCCAGCGCGGCATCGGCATCGGCCAGTGTGGCGGCGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCCATCGCTGGCGGGGCGCGACTGccc                           <  8:260531/128‑1 (MQ=255)                                                                                                                                           gcagcaGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGGCGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCCATCGCTGGCGGGGCGCGACTGCCCGCCAATCAGCAAAGGCACCTGCAAc  <  6:49688/147‑1 (MQ=255)                                                                                                                                                 gcCAGCGCGGCATCGGCATCGGCCAGTGTGGCGGCGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCCATCGCt                                             <  7:105750/98‑1 (MQ=255)                                                                                                                                                 gcCAGCGCGGCATCGGCATCGGCCAGTGTGGCGGCGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCCATCGCt                                             >  8:105750/1‑98 (MQ=255)                                                                                                                                                     |                                                                                                                                         AGTTGGCCATGGCCCCGGTTTCACCGGCGGCGGCGATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCA‑GGGCGGCCCAGGCCGGAAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGGCGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCCATCGCTGGCGGGGCGCGACTGCCCGCCAATCAGCAAAGGCACCTGCAAC  >  NC_002947/3797708‑3797991

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