Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,551,531	C→T	7.8%	G300G (GGC→GGT)	pydA →	NADP‑dependent dihydropyrimidine dehydrogenase subunit PreA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,551,531	0	C	T	7.8%	147.1 / 2.3	51	G300G (GGC→GGT)	pydA	NADP‑dependent dihydropyrimidine dehydrogenase subunit PreA
Reads supporting (aligned to +/- strand):  ref base C (21/26);  new base T (2/2);  total (23/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCGGGGCTTGCCGATTTGTGGCATTGGCGGCATCGGCAACTGGCGTGACGCCGCCGAGTTCATGGCCTTGGGCAGTGGCGCCGTGCAGGTGTGCACGGCGGCGATGCTGCATGGCTTCCGTATC‑GTCGAGGACATGCAGGACGGCCTGGCACGCTGGATGGACCAGCATGGGCATGCCACTGTCGAAGCGTTCCGCGGGCAGGCGGTGGGGCATACCACCGACTGGAAATACCTGGACATCAACTACAAGTCGGTCGCGCACATTGACCAGGAGGCGTGCATTGGCTGCGGGC  >  NC_002947/4551386‑4551679                                                                                                                                                   |                                                                                                                                                     cgctGGGCTTGCCGATTTGTGGCATTGGCGGCATCGGCAACTGGCGTGACGCCGCCGAGTTCATGGCCTTGGGCAGTGGCGCCGTGCAGGTGTGCACGGCGGCGATGCTGCATGGCTTCCGTATC‑GTCGAGGACATGCAGGACGGCCTg                                                                                                                                                   <  3:246579/149‑1 (MQ=255)  gcgGGGCTTGCCGATTTGTGGCATTGGCGGCATCGGCAACTGGCGTGACGCCGCCGAGTTCATGGCCTTGGGCAGTGGCGCCGTGCAGGTGTGCACGGCGGCGATGCTGCATGGCTTCCGTATC‑GTCGAGGACATGCAGGACGGCCTgg                                                                                                                                                  >  5:125311/1‑149 (MQ=255)  gcgGGGCTTGCCGATTTGTGGCATTGGCGGCATCGGCAACTGGCGTGACGCCGCCGAGTTCATGGCCTTGGGCAGTGGCGCCGTGCAGGTGTGCACGGCGGCGATGCTGCATGGCTTCCGTATC‑GTCGAGGACATGCAGGACGGCCTgg                                                                                                                                                  >  2:52515/1‑149 (MQ=255)            ccGATTTGTGGCATTGGCGGCATCGGCAACTGGCGTGACGCCGCCGAGTTCATGGCCTTGGGCAGTGGCGCCGTGCAGGTGTGCACGGCGGCGATGCTGCATGGCTTCCGTATC‑GTCGAGGACATGCAGGACGGCCTGGCACGCtggat                                                                                                                                        <  8:202902/149‑1 (MQ=255)             cGATTTGTGGCATTGGCGGCATCGGCAACTGGCGTGACGCCGCCGAGTTCATGGCCTTGGGCAGTGGCGCCGTGCAGGTGTGCACGGCGGCGATGCTGCATGGCTTCCGTATC‑GTCGAGGACATGCAGGACGGCCTGGCACGCtggatg                                                                                                                                       >  5:254499/1‑149 (MQ=255)             cGATTTGTGGCATTGGCGGCATCGGCAACTGGCGTGACGCCGCCGAGTTCATGGCCTTGGGCAGTGGCGCCGTGCAGGTGTGCACGGCGGCGATGCTGCATGGCTTCCGTATC‑GTCGAGGACATGCAGGACGGCCTGGCACGCtggatg                                                                                                                                       >  7:48439/1‑149 (MQ=255)             cGATTTGGGGCATTGGCGGCATCGGCAATTGGCGTGACGCCGCCGAGTTCATGGCCTTGGGCAGTGGCGCCGTGCAGGTGTGCACGGCGGCGATGCTGCATGGCTTCCGTATC‑GTCGAGGACATGCAGGACGGCCTGGCACGCtggatg                                                                                                                                       <  7:98709/149‑1 (MQ=255)                  tgtgGCATTGGCGGCATCGGCAACTGGCGTGACGCCGCCGAGTTCATGGCCTTGGGCAGTGGCGCCGTGCAGGTGTGCACGGCGGCGATGCTGCATGGCTTCCGTATC‑GTCGAGGACATCCAGGACGGCCTGGCACGCTGGATGGAcc                                                                                                                                   <  6:88673/148‑1 (MQ=255)                               gcATCGGCAACTGGCGTGACGCCGCCGAGTTCATGGCCTTGGGCAGTGGCGCCGTGCAGGTGTGCACGGCGGCGATGCTGCATGGCTTCCGTATC‑GTCGAGGACATGCAGGACGGCCTGGCACGCTGGATGGACCAGCATGGGCATGc                                                                                                                      >  4:223034/1‑148 (MQ=255)                                       aaCTGGCGTGACGCCGCCGAGTTCATGGCCTTGGGCAGTGGCGCCGTGCAGGTGTGCACGGCGGCGATGCTGCATGGCTTCCGTATC‑GTCGAGGACATGCAGGACGGTCTGGCACGCTGGATGGACCAGCATGGGCATGCCACTGTCg                                                                                                              >  3:38078/1‑148 (MQ=255)                                       aaCTGGCGTGACGCCGCCGAGTTCATGGCCTTGGGCAGTGGCGCCGTGCAGGTGTGCACGGCGGCGATGCTGCATGGCTTCCGTATC‑GTCGAGGACATGCAGGACGGTCTGGCACGCTGGATGGACCAGCATGGGCATGCCACTGTCg                                                                                                              >  3:34616/1‑148 (MQ=255)                                            gCGTGACGCCGCCGAGTTCATGGCCTTGGGCAGTGGCGCCGTGCAGGTGTGCACGGTGGCGATGCTGCATGGCTTCCGTATC‑GTCGAGGACATGCAGGACGGCCTGGCACGCTGGATGGACCAGCATGGGCATGCCACTGTCGAAGc                                                                                                          <  8:141148/147‑1 (MQ=255)                                            gCGTGACGCCGCCGAGTTCATGGCCTTGGGCAGTGGCGCCGTGCAGGTGTGCACGGCGGCGATGCTGCATGGCTTCCGTATC‑GTCGAGGACATGCAGGACGGCCTGGCACGCTGGATGGACCAGCATGGGCATGCCACTGTCGAAGc                                                                                                          <  3:223034/147‑1 (MQ=255)                                                    ccgccgAGTTCATGGCCTTGGGCAGTGGCACCGTGCAGGTGTGCACGGCGGCGATGCTGCATGGCTTCCGTATC‑GTCGAGGACATGCAGGACGGCCTGGCACGCTGGATGGACCAGCATGGGCATGc                                                                                                                      <  6:233035/127‑1 (MQ=255)                                                    ccgccgAGTTCATGGCCTTGGGCAGTGGCACCGTGCAGGTGTGCACGGCGGCGATGCTGCATGGCTTCCGTATC‑GTCGAGGACATGCAGGACGGCCTGGCACGCTGGATGGACCAGCATGGGCATGc                                                                                                                      >  5:233035/1‑127 (MQ=255)                                                      gccgAGTTCATGGCCTTGGGCAGTGGCGCCGTGCAGGTGTGCACGGCGGCGATGCTGCATGGCTTCCGTATC‑GTCGAGGACATGCAGGACGGCCTGGCACGCtggatgga                                                                                                                                     <  7:123162/110‑1 (MQ=255)                                                      gccgAGTTCATGGCCTTGGGCAGTGGCGCCGTGCAGGTGTGCACGGCGGCGATGCTGCATGGCTTCCGTATC‑GTCGAGGACATGCAGGACGGCCTGGCACGCtggatgga                                                                                                                                     >  8:123162/1‑110 (MQ=255)                                                         gAGTTCATGGCCTTGGGCAGTGGCGCCGTGCAGGTGTGCACGGCGGCGATGCTGCATGGCTTCCGTATC‑GTCGAGGACATGCAGGACGGCCTGGCACGCTGGATGGACCAGCATGGGCATGCCACTGTCGAAGCGTTCCGCGGGCAGGc                                                                                           <  8:48439/149‑1 (MQ=255)                                                             tCATGGCCTTGGGCAGTGGCGCCGTGCAGGTGTGCACGGCGGCGATGCTGCATGGCTTCCGTATC‑GTCGAGGACATGCAGGACGGCCTGGCACGCTGGATGGACCAGCATGGGCATGCCACTGt                                                                                                                >  3:155125/1‑124 (MQ=255)                                                             tCATGGCCTTGGGCAGTGGCGCCGTGCAGGTGTGCACGGCGGCGATGCTGCATGGCTTCCGTATC‑GTCGAGGACATGCAGGACGGCCTGGCACGCTGGATGGACCAGCATGGGCATGCCACTGt                                                                                                                <  4:155125/124‑1 (MQ=255)                                                              cATGGCCTTGGGCAGTGGCGCCGTGCAGGTGTGCACGGCGGCGATGCTGCATGGCTTCCGTATC‑GTCGAGGACATGCAGGACGGCCTGGCACGCTGGATGGACCAGCATGGGCATGCCACTGTCGAAGCGTTCCGCGGGCAGGCGGTgg                                                                                      >  5:136009/1‑149 (MQ=255)                                                               aTGGCCTTGGGCAGTGGCGCCGTGCAGGTGTGCACGGCGGCGATGCTGCATGGCTTCCGTATC‑GTCGAGGACATGCAGGACGGTCTGGCACGCTGGATGGACCAGCATGGGCATGCCACTGTCGAAGCGTTCCGCGGGCAGGCGGTgg                                                                                      <  4:34616/148‑1 (MQ=255)                                                               aTGGCCTTGGGCAGTGGCGCCGTGCAGGTGTGCACGGCGGCGATGCTGCATGGCTTCCGTATC‑GTCGAGGACATGCAGGACGGTCTGGCACGCTGGATGGACCAGCATGGGCATGCCACTGTCGAAGCGTTCCGCGGGCAGGCGGTgg                                                                                      <  4:38078/148‑1 (MQ=255)                                                                               gcgcCGTGCAGGTGTGCACGGCGGCGATGCTGCATGGCTTCCGTATC‑GTCGAGGACATGCAGGACGGCCTGGCACGCTGGATGGACCAGCATGGGCATGCCACTGTCGAAGCGTTCCGCGGGCAGGCGGTGGGGCATACCACCGACTgg                                                                     <  1:52515/149‑1 (MQ=255)                                                                               gcgcCGTGCAGGTGTGCACGGCGGCGATGCTGCATGGCTTCCGTATC‑GTCGAGGACATGCAGGACGGCCTGGCACGCTGGATGGACCAGCATGGGCATGCCACTGTCGAAGCGTTCCGCGGGCAGGCGGTGGGGCATACCACCGACTg                                                                      >  2:90094/1‑148 (MQ=255)                                                                               gcgcCGTGCAGGTGTGCACGGCGGCGATGCTGCATGGCTTCCGTATC‑GTCGAGGACATGCAGGACGGCCTGGCACGCTGGATGGACCAGCATGGGCATGCCACTGTCGAAGCGTTCCGCGGGCAGGCGGTGGGGCATACCACCGACTg                                                                      >  5:85635/1‑148 (MQ=255)                                                                                cgcCGTGCAGGTGTGCACGGCGGCGATGCTGCATGGCTTCCGTATCAGTCGAGGACATGCAGGACGGCCTGGCACGCTGGATGGACCAGCATGGGCATGCCACTGTCGAAGCGTTCCGCGGGCAGGCGGTggg                                                                                     >  1:196102/1‑133 (MQ=255)                                                                                cgcCGTGCAGGTGTGCACGGCGGCGATGCTGCATGGCTTCCGTATC‑GTCGAGGACATGCAGGACGGCCTGGCACGCTGGATGGACCAGCATGGGCATGCCACTGTCGAAGCGTTCCGCGGGCAGGCGGTggg                                                                                     <  2:196102/132‑1 (MQ=255)                                                                                    gTGCAGGTGTGCACGGCGGCGATGCTGCATGGCTTCCGTATC‑GTCGAGGACATGCAGGACGGCCTGGCACGCTGGATGGACCAGCATGGGCATGCCACTGTCGAAGCGTTCCGCGGGCAGGCGGTGGGGCATACCACCGACTGGAAATa                                                                >  1:108374/1‑149 (MQ=255)                                                                                     tGCAGGTGTGCACGGCGGCGATGCTGCATGGCTTCCGTATC‑GTCGAGGACATGCAGGACGGCCTGGCACGCTGGATGGACCAGCATGGGCATGCCACTGTCGAAGCGTTCCGCGGGCAGGCGGTGGGGCATACCACCGACTGGAAATAc                                                               <  2:221998/149‑1 (MQ=255)                                                                                     tGCAGGTGTGCACGGCGGCGATGCTGCATGGCTTCCGTATC‑GTCGAGGACATGCAGGACGGCCTGGCACGCTGGATGGACCAGCATGGGCATGCCACTGTCGAAGCGTTCCGCGGGCAGGCGGTGGGGCATACCACCGACTGGAAATAc                                                               <  6:200724/149‑1 (MQ=255)                                                                                               cACGGCGGCGATGCTGCATGGCTTCCGTATC‑GTCGAGGACATGCAGGACGGCCTGGCACGCTGGATGGACCAGCATGGGCATGCCACTGTCGAAGCGTTCCGCGGGCAGGCGGTGGGGCATACCACCGACTGGAAATACCTGGACATc                                                      <  1:103803/148‑1 (MQ=255)                                                                                                          tgctgcATGGCTTCCGTATC‑GTCGAGGACATGCAGGACGGCCTGGCACGCTGGATGGACCAGCATGGGCATGCCACTGTCGAAGCGTTCCGCGGGCAGGCGGTGGGGCATACCACCGACTGGAAATACCTGGACATCAACTACAAgtcg                                          <  6:136009/149‑1 (MQ=255)                                                                                                             tgcATGGCTTCCGTATC‑GTCGAGGACATGCAGGACGGCCTGGCACGCTGGATGGACCAGCATGGGCATGCCACTGTCGAAGCGTTCCGCGGGCAGGCGGTGGGGCATACCACCGACTGGAAATACCTGGACATCAACTACaa                                              <  7:164928/142‑1 (MQ=255)                                                                                                             tgcATGGCTTCCGTATC‑GTCGAGGACATGCAGGACGGCCTGGCACGCTGGATGGACCAGCATGGGCATGCCACTGTCGAAGCGTTCCGCGGGCAGGCGGTGGGGCATACCACCGACTGGAAATACCTGGACATCAACTACaa                                              >  8:164928/1‑142 (MQ=255)                                                                                                               cATGGCTTCCGTATC‑GTCGAGGACATGCAGGACGGCCTGGCACGCTGGATGGACCAGCATGGGCATGCCACTGTCGAAGCGTTCCGCGGGCAGGCGGTGGGGCATACCACCGACTGGAAATACCTGGACATCAACTACAAGTCGGTcgc                                     >  3:47155/1‑149 (MQ=255)                                                                                                                       ccGTATC‑GTCGAGGACATGCAGGACGGCCTGGCACGCTGGATGGACCAGCATGGGCATGCCACTGTCGAAGCGTTCCGCGGGCAGGCGGTGGGGCATACCACCGACTGGAAATACCTGGACATCAACTACAAGTCGGTCGCGCACAtt                              <  6:254499/148‑1 (MQ=255)                                                                                                                       ccGTATC‑GTCGAGGACATGCAGGACGGCCTGGCACGCTGGATGGACCAGCATGGGCATGCCACTGTCGAAGCGTTCCGCGGGCAGGCGGTGGGGCATACCACCGACTGGAAATACCTGGACATCAACTACAAGTCGGTCGCGCACAtt                              >  1:133150/1‑148 (MQ=255)                                                                                                                               gtcgAGGACATGCAGGACGGCCTGGCACGCTGGATGGACCAGCATGGGCATGCCACTGTCGAAGCGTTCCGCGGGCAGGCGGTGGGGCGTACCACCGACTGGAAATACCTGGACATCAACTACAAGTCGGTCGCGCACATTGACCagga                      <  1:35239/149‑1 (MQ=255)                                                                                                                                tcgAGGACATGCAGGACGGCCTGGCACGCTGGATGGACCAGCATGGGCATGCCACTGTCGAAGCGTTCCGCGGGCAGGCGGTGGGGCATACCACCGACTGGAAATACCTGGACATCAACTACAAGTCGGTCGCGCACATTGACCaggag                     <  6:85635/149‑1 (MQ=255)                                                                                                                                tcgAGGACATGCAGGACGGCCTGGCACGCTGGATGGACCAGCATGGGCATGCCACTGTCGAAGCGTTCCGCGGGCAGGCGGTGGGGCATACCACCGACTGGAAATACCTGGACATCAACTACAAGTCGGTCGCGCACATTGACCaggag                     >  2:243098/1‑149 (MQ=255)                                                                                                                                    ggACATGCAGGACGGCCTGGCACGCTGGATGGACCAGCATGGGCATGCCACTGTCGAAGCGTTCCGCGGGCAGGCGGTGGGGCATACCACCGACTGGAAATACCTGGACATCAACTACAAGTCGGTCGCGCACATTGACCAGGAGGCGt                 <  6:125311/149‑1 (MQ=255)                                                                                                                                    ggACATGCAGGACGGCCTGGCACGCTGGATGGACCAGCATGGGCATGCCACTGTCGAAGCGTTCCGCGGGCAGGCGGTGGGGCATACCACCGACTGGAAATACCTGGACATCAACTACAAGTCGGTCGCGCACATTGACCAGGAGGCGt                 <  2:133150/149‑1 (MQ=255)                                                                                                                                       cATGCAGGACGGCCTGGCACGCTGGATGGACCAGCATGGGCATGCCACTGTCGAAGCGTTCCGCGGGCAGGCGGTGGGGCATACCACCGACTGGAAATACCTGGACATCAACTACAAg                                             >  6:108533/1‑118 (MQ=255)                                                                                                                                       cATGCAGGACGGCCTGGCACGCTGGATGGACCAGCATGGGCATGCCACTGTCGAAGCGTTCCGCGGGCAGGCGGTGGGGCATACCACCGACTGGAAATACCTGGACATCAACTACAAg                                             <  5:108533/118‑1 (MQ=255)                                                                                                                                           cAGGACGGCCTGGCACGCTGGATGGACCAGCATGGGCATGCCACTGTCGAAGCGTTCCGCGGGCAGGCGGTGGGGCATACCACCGACTGGAAATACCTGGACATCAACTACAAGTCGGTCGCGCACATTGACCAGGAGGCGTGCATTgg          <  1:243098/149‑1 (MQ=255)                                                                                                                                           cAGGACGGCCTGGCACGCTGGATGGACCAGCATGGGCATGCCACTGTCGAAGCGTTCCGCGGGCAGGCGGTGGGGCATACCACCGACTGGAAATACCTGGACATCAACTACAAGTCGGTCGCGCACATTGACCAGGAGGCGTGCATTgg          >  4:130678/1‑149 (MQ=255)                                                                                                                                                cGGCCTGGCACGCTGGATGGACCAGCATGGGCATGCCACTGTCGAAGCGTTCCGCGGGCAGGCGGTGGGGCATACCACCGACTGGAAATACCTGGACATCAACTACAAGTCGGTCGCGCACATTGAcc                          >  3:10335/1‑128 (MQ=255)                                                                                                                                                cGGCCTGGCACGCTGGATGGACCAGCATGGGCATGCCACTGTCGAAGCGTTCCGCGGGCAGGCGGTGGGGCATACCACCGACTGGAAATACCTGGACATCAACTACAAGTCGGTCGCGCACATTGAcc                          <  4:10335/128‑1 (MQ=255)                                                                                                                                                 ggCCTGGCACGCTGGATGGACCAGCATGGGCATGCCACTGTCGAAGCGTTCCGCGGGCAGGCGGTGGGGCATACCACCGACTGGAAATACCTGGACATCAACTACAAGTCGGTCGCGCACATTGACCAGGAGGCGTGCATTGGCTGCgg    >  6:181334/1‑149 (MQ=255)                                                                                                                                                  gCCTGGCACGCTGGATGGACCAGCATGGGCATGCCACTGTCGAAGCGTTCCGCGGGCAGGCGGTGGGGCATACCACCGACTGGAAATACCTGGACATCAACTACAAGTCGGTCGCGCACATTGACCAGGAGGCGTGCATTGGCTGCggg   <  2:108374/149‑1 (MQ=255)                                                                                                                                                   ccTGGCACGCTGGATGGACCAGCATGGGCATGCCACTGTCGAAGCGTTCCGCGGGCAGGCGGTGGGGCATACCACCGACTGGAAATACCTGGACATCAACTACAAGTCGGTCGCGCACATTGACCAGGAGGCGTGCATTGGCTACGGgc  >  8:93466/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                     CGCGGGGCTTGCCGATTTGTGGCATTGGCGGCATCGGCAACTGGCGTGACGCCGCCGAGTTCATGGCCTTGGGCAGTGGCGCCGTGCAGGTGTGCACGGCGGCGATGCTGCATGGCTTCCGTATC‑GTCGAGGACATGCAGGACGGCCTGGCACGCTGGATGGACCAGCATGGGCATGCCACTGTCGAAGCGTTCCGCGGGCAGGCGGTGGGGCATACCACCGACTGGAAATACCTGGACATCAACTACAAGTCGGTCGCGCACATTGACCAGGAGGCGTGCATTGGCTGCGGGC  >  NC_002947/4551386‑4551679

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