Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,797,420	C→A	9.8%	G192C (GGC→TGC)	vdh ←	vanillin dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,797,420	0	C	A	9.8%	154.0 / 6.3	51	G192C (GGC→TGC)	vdh	vanillin dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (25/21);  new base A (2/3);  total (27/24)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

CCGCCAGTTCGCCGACGATGCGCCCGACGTGCGTCGAACCGGTGAAGTTGACCCGGCGTACCGCAGGGTTGGCGATCAGCC‑‑G‑CTCGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGC  >  NC_002947/3797273‑3797563                                                                                                                                                       |                                                                                                                                                ccgccAGTTCGCCGACGATGCGCCCGACGTGCGTCGAACCGGTGAAGTTGACCCGGCGTACCGCAGGGTTGGCGATCAGCC‑‑G‑CTCGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCt                                                                                                                                                <  1:236314/149‑1 (MQ=255) ccgccAGTTCGCCGACGATGCGCCCGACGTGCGTCGAACCGGTGAAGTTGACCCGGCGTACCGCAGGGTTGGCGATCAGCC‑‑G‑CTCGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCt                                                                                                                                                >  2:236314/1‑149 (MQ=255)           gCCGACGATGCGCCCGACGTGCGTCGAACCGGTGAAGTTGACCCGGCGTACCGCAGGGTTGGCGATCAGCC‑‑G‑CTCGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTg                                                                                                                                        >  2:294053/1‑147 (MQ=255)             cgacgaTGCGCCCGACGTGCGTCGAACCGGTGAAGTTGACCCGGCGTACCGCAGGGTTGGCGATCAGCC‑‑G‑CTCGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCAGGCATCggcgacg                                                                                                                                    <  8:244471/149‑7 (MQ=255)              gacgaTGCGCCCGACGTGCGTCGAACCGGTGAAGTTGACCCGGCGTACCGCAGGGTTGGCGATCAGCCAGGAGTAGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTg                                                                                                                                        >  1:290301/1‑147 (MQ=255)                  aTGCGCCCGACGTGCGTCGAACCGGTGAAGTTGACCCGGCGTACCGCAGGGTTGGCGATCAGCC‑‑G‑CTCGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTg                                                                                                                               <  8:28713/149‑1 (MQ=255)                       cccGACGTGCGTCGAACCGGTGAAGTTGACCCGGCGTACCGCAGGGTTGGCGATCAGCC‑‑G‑CTCGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCg                                                                                                                           >  3:103768/1‑148 (MQ=255)                                  tCGAACCGGTGAAGTTGACCCGGCGTACCGCAGGGTTGGCGATCAGCC‑‑G‑CTCGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATg                                                                                                                >  1:4885/1‑148 (MQ=255)                                     aaCCGGTGAAGTTGACCCGGCGTACCGCAGGGTTGGCGATCAGCC‑‑G‑CTCGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGAcc                                                                                                            >  7:334742/1‑149 (MQ=255)                                     aaCCGGTGAAGTTGACCCGGCGTACCGCAGGGTTGGCGATCAGCC‑‑G‑CTCGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGAcc                                                                                                            >  4:154688/1‑149 (MQ=255)                                        cGGTGAAGTTGACCCGGCGTACCGCAGGGTTGGCGATCAGCC‑‑G‑CTCGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGAccgcc                                                                                                         <  5:144235/149‑1 (MQ=255)                                           tGAAGTTGACCCGGCGTACCGCAGGGTTGGCGATCAGCC‑‑G‑CTCGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCggg                                                                                                      <  2:290301/149‑1 (MQ=255)                                            gAAGTTGACCCGGCGTACCGCAGGGTTGGCGATCAGCC‑‑G‑CTCGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCC‑GCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCt                                                                                                    >  7:203837/1‑149 (MQ=255)                                              aGTTGACCCGGCGTACCGCAGGGTTGGCGATCAGCC‑‑G‑CTCGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCt                                                                                                    >  1:174732/1‑148 (MQ=255)                                                ttGACCCGGCGTACCGCAGGGTTGGCGATCAGCC‑‑G‑CTCGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGAc                                                                                                             >  1:185630/1‑137 (MQ=255)                                                ttGACCCGGCGTACCGCAGGGTTGGCGATCAGCC‑‑G‑CTCGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGAc                                                                                                             <  2:185630/137‑1 (MQ=255)                                                      cGGCGTACCGCAGGGTTGGCGATCAGCC‑‑G‑CTCGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGAc                                                                                                             <  3:216543/131‑1 (MQ=255)                                                      cGGCGTACCGCAGGGTTGGCGATCAGCC‑‑G‑CTCGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGAc                                                                                                             >  4:216543/1‑131 (MQ=255)                                                          gTACCGCAGGGTTGGCGATCAGCC‑‑G‑CTCGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGcc                                                                                       <  2:174732/149‑1 (MQ=255)                                                              cGCAGGGTTGGCGATCAGCC‑‑G‑CTCGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGccat                                                                                     >  6:27662/1‑145 (MQ=255)                                                              cGCAGGGTTGGCGATCAGCC‑‑G‑CTCGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCtt                                                                                     >  5:15024/1‑147 (MQ=255)                                                              cGCAGGGTTGGCGATCAGCC‑‑G‑CTCGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCtt                                                                                     <  6:15024/147‑1 (MQ=255)                                                              cGCAGGGTTGGCGATCAGCC‑‑G‑CTCGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGa                                                                                   <  8:334742/149‑1 (MQ=255)                                                              cGCAGGGTTGGCGATCAGCC‑‑G‑CTCGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGa                                                                                   <  1:230647/149‑1 (MQ=255)                                                              cGCAGGGTTGGCGATCAGCC‑‑G‑CTCGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCAGGCATCgg                                                                                                                                         <  6:282308/95‑2 (MQ=255)                                                              cGCAGGGTTGGCGATCAGCC‑‑G‑CTCGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCAGGCATCgg                                                                                                                                         >  5:282308/1‑94 (MQ=255)                                                                cAGGGTTGGCGATCAGCG‑‑G‑CTCGACGATGGCGGGGGCATCCTGCGGCGCTTTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGc                                                                                 <  4:122740/149‑1 (MQ=255)                                                                 aGGGTTGGCGATCAGCC‑‑G‑CTCGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAg                                                                                          >  7:87880/1‑139 (MQ=255)                                                                 aGGGTTGGCGATCAGCC‑‑G‑CTCGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAg                                                                                          <  8:87880/139‑1 (MQ=255)                                                                   ggTTGGCGATCAGCC‑‑G‑CTCGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGcacc                                                                              <  3:154688/149‑1 (MQ=255)                                                                              aGCC‑‑G‑CTCGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGcc                                                                   <  2:4885/149‑1 (MQ=255)                                                                                 c‑‑G‑CTCGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCC‑GCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAg                                                               <  8:203837/149‑1 (MQ=255)                                                                                                  gggggCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGG                                                  <  7:310220‑M1/149‑1 (MQ=255)                                                                                                         tCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGC                                           >  5:13782‑M1/1‑149 (MQ=255)                                                                                                           cTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCAc                                         <  1:288567/149‑1 (MQ=255)                                                                                                               ggcgCATTGCTGATGACATTGACCACGCCGTCGCAGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGt                                     <  4:103768/149‑1 (MQ=255)                                                                                                                         tgatgaCATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGc                                                                    >  1:25624/1‑108 (MQ=255)                                                                                                                         tgatgaCATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGc                                                                    <  2:25624/108‑1 (MQ=255)                                                                                                                           atgaCATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCAc                         >  2:323/1‑149 (MQ=255)                                                                                                                            tgaCATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCAcc                        <  1:323/149‑1 (MQ=255)                                                                                                                             gaCATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCg                       >  7:60687/1‑149 (MQ=255)                                                                                                                               cATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCggg                     >  5:331619/1‑149 (MQ=255)                                                                                                                                      cACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCcctgcgcgcgccctgcggcgtggtgttgg              >  6:267208‑M1/1‑120 (MQ=255)                                                                                                                                      cACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTcc              >  8:338934/1‑149 (MQ=255)                                                                                                                                      cACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCA                                                >  1:82174‑M1/1‑115 (MQ=255)                                                                                                                                      cACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCA                                                <  2:82174‑M1/115‑1 (MQ=255)                                                                                                                                              gcacGATGCAGGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGtt                <  4:168175/136‑1 (MQ=255)                                                                                                                                              gcacGATGCAGGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGtt                >  3:168175/4‑139 (MQ=255)                                                                                                                                                ccGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCcctgcgcgcgccctgcggcgtggtgttgggcatcgcacc    >  7:132058‑M1/1‑110 (MQ=255)                                                                                                                                                 cGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATg   >  5:332649/1‑149 (MQ=255)                                                                                                                                                  gATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGc  >  7:107131/1‑149 (MQ=255)                                                                                                                                                       |                                                                                                                                                CCGCCAGTTCGCCGACGATGCGCCCGACGTGCGTCGAACCGGTGAAGTTGACCCGGCGTACCGCAGGGTTGGCGATCAGCC‑‑G‑CTCGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGC  >  NC_002947/3797273‑3797563

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: ‑