Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,780,122	C→T	9.9%	V2398M (GTG→ATG) ‡	pvdJ ←	non‑ribosomal peptide synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,780,122	0	C	T	9.9%	135.7 / 3.8	50	V2398M (GTG→ATG) ‡	pvdJ	non‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base C (23/22);  new base T (2/3);  total (25/25)
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

CGCATCCACCACCACGGCGCTGTGCACACCTTCATGGGCCAGCAGGCGCGATTCGATTTCGCCCAGTTCTATGCGGTAGCCGCGGATCTTCACCTGGAAGTCGGCGCGCCCCACATACTCCAGGGTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGT‑ACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCGGGGACAAAGCTCAATGCCGTGCGCACCGGATCCCCAAGGTAACCGCGCCCGACACCTACGCCGGCCACATGCAACTCGCCGGTGGCCCGG  >  NC_002947/4779979‑4780270                                                                                                                                                |                                                                                                                                                      cgcATCCACCACCACGGCGCTGTGCACACCTTCATGGGCCAGCAGGCGCGATTCGATTTCGCCCAGTTCTATGCGGTAGCCGCGGATCTTCACCTGGAAGTCGGCGCGCCCCACATACTCCAGGGTCTGCCCATCGGCGCACACGTgc                                                                                                                                                   <  1:284763/148‑1 (MQ=255) cgcATCCACCACCACGGCGCTGTGCACACCTTCATGGGCCAGCAGGCGCGATTCGATTTCGCCCAGTTCTATGCGGTAGCCGCGGATCTTCACCTGGAAGTCGGCGCGCCCCACATACTCCAGGGTCTGCCCATCGGCGCACACGTgc                                                                                                                                                   <  3:319669/148‑1 (MQ=255)  gcATCCACCACCACGGCGCTGTGCACACCTTCATGGGCCAGCAGGCGCGATTCGATTTCGCCCAGTTCTATGCGGTAGCCGCGGATCTTCACCTGGAAGTCGGCGCGCCCCACATACTCCAGGGTCTGCCCATCGGCGCACACGTgcgc                                                                                                                                                 >  6:213689/1‑149 (MQ=255)       caccaccacGGCGCTGTGCACACCTTCATGGGCCAGCAGGCGCGATTCGATTTCGCCCAGTTCTATGCGGTAGCCGCGGATCTTCACCTGGAAGTCGGCGCGCCCCACATACTCCAGGGTCTGCCCATCGGCGCACGTGTGCGccgat                                                                                                                                             <  1:151396/148‑4 (MQ=255)        accaccacGGCGCTGTGCACACCTTCATGGGCCAGCAGGCGCGATTCGATTTCGCCCAGTTCTATGCGGTAGCCGCGGATCTTCACCTGGAAGTCGGCGCGCCCCACATACTCCAGGGTCTGCCCATCGGCGCACGTGTGCGccgatgg                                                                                                                                           <  3:99047/149‑6 (MQ=255)        accaccacGGCGCTGTGCACACCTTCATGGGCCAGCAGGCGCGATTCGATTTCGCCCAGTTCTATGCGGTAGCCGCGGATCTTCACCTGGAAGTCGGCGCGCCCCACATACTCCAGGGTCTGCCCATCGGCGCACGTGTGCGccgatgg                                                                                                                                           >  6:160616/1‑144 (MQ=255)        accaccacGGCGCTGTGCACACCTTCATGGGCCAGCAGGCGCGATTCGATTTCGCCCAGTTCTATGCGGTAGCCGCGGATCTTCACCTGGAAGTCGGCGCGCCCCACATACTCCAGGGTCTGCCCATCGGCGCACACGTGCGCCAGGTc                                                                                                                                           >  6:242874/1‑149 (MQ=255)         ccaccacGGCGCTGTGCACACCTTCATGGGCCAGCAGGCGCGATTCGATTTCGCCCAGTTCTATGCGGTAGCCGCGGATCTTCACCTGGAAGTCGGCGCGCCCCACATACTCCAGGGTCTGCCCATCGGCGCACGTGTGCGccgatggg                                                                                                                                          <  5:160616/149‑7 (MQ=255)                         cacaCCTTCATGGGCCAGCAGGCGCGATTCGATTTCGCCCAGTTCTATGCGGTAGCCGCGGATCTTCACCTGGAAGTCGGCGCGCCCCACATACTCCAGGGTCTGCCCATCGGCGCCCACGTGCGCCAGGTCGCCACTGCGGT‑ACAGgc                                                                                                                         >  4:122936/1‑149 (MQ=255)                            acCTTCATGGGCCAGCAGGCGCGATTCGATTTCGCCCAGTTCTATGCGGTAGCCGCGGATCTTCACCTGGAAGTCGGCGCGCCCCACATACTCCAGGGTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGT‑ACAGGCGCt                                                                                                                      >  4:25317/1‑149 (MQ=255)                              cTTCATGGGCCAGCAGGCGCGATTCGATTTCGCCCAGTTCTATGCGGTAGCCGCGGATCTTCACCTGGAAGTCGGCGCGCCCCACATACTCCAGGGTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGT‑ACAGGCGCTCa                                                                                                                    >  5:39774/1‑149 (MQ=255)                              attCATGGGCCAGCAGGCGCGATTCGATTTCGCCCAGTTCTATGCGGTAGCCGCGGATCTTCACCTGGAAGTCGGCGCGCCCCACATACTCCAGGGTCTGCCCATCGGAGCACACGTGCGCCAGGTCGCCACTGCGGT‑ACAGCCCCTCa                                                                                                                    >  8:46865/2‑149 (MQ=255)                                         agcagGCGCGATTCGATTTCGCCCAGTTCTATGCGGTAGCCGCGGATCTTCACCTGGAAGTCGGCGCGCCCCACATACTCCAGGGTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGT‑ACAGGCGCTCACCCGGGCGGCt                                                                                                         >  8:179372/1‑149 (MQ=255)                                         agcagGCGCGATTCGATTTCGCCCAGTTCTATGCGGTAGCCGCGGATCTTCACCTGGAAGTCGGCGCGCCCCACATACTCCAGGGTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCCGT‑ACAGGCGCTCACCCGGGCGGCt                                                                                                         <  7:46865/149‑1 (MQ=255)                                                gcgATTCGATTTCGCCCAGTTCTATGCGGTAGCCGCGGATCTTCACCTGGAAGTCGGCGCGCCCCACATACTCCAGGGTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGT‑ACAGGCGCTCACCCGGGCGGCTGGCGAAc                                                                                                  >  3:267190/1‑149 (MQ=255)                                                        attTCGCCCAGTTCTATGCGGTAGCCGCGGATCTTCACCTGGAAGTCGGCGCGCCCCACATACTCCAGGGTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGT‑ACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCggg                                                                                          <  6:126954/149‑1 (MQ=255)                                                                     cTATGCGGTAGCCGCGGATCTTCACCTGGAAGTCGGCGCGCCCCACATACTCCAGGGTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGT‑ACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCGGGGACAAAGCTCaa                                                                              <  5:242874/148‑1 (MQ=255)                                                                        tGCGGTAGCCGCGGATCTTCACCTGGAAGTCGGCGCGCCCCACATACTCCAGGGTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGT‑ACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCGGGGACAAAGCTCAATGcc                                                                          <  7:179372/149‑1 (MQ=255)                                                                             tAGCCGCGGATCTTCACCTGGAAGTCGGCGCGCCCCACATACTCCAGGGTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGT‑ACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCgg                                                                                           >  6:282210/1‑127 (MQ=255)                                                                             tAGCCGCGGATCTTCACCTGGAAGTCGGCGCGCCCCACATACTCCAGGGTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGT‑ACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCgg                                                                                           <  5:282210/127‑1 (MQ=255)                                                                                ccGCGGATCTTCACCTGGAAGTCGGCGCGCCCCACATACTCCAGGGTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGT‑ACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCGGGGACAAAGCTCAATGCCGTGCGCAc                                                                  >  1:304412/1‑149 (MQ=255)                                                                                 cgcgGATCTTCACCTGGAAGTCGGCGCGCCCCACATACTCCAGGGTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGT‑ACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCGGGGACAAAGCTCAATGCCGTGCGCAcc                                                                 >  1:139122/1‑149 (MQ=255)                                                                                   cgGATCTTCACCTGGAAGTCGGCGCGCCCCACATACTCCAGGGTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGT‑ACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCGGGGACAAAGCTCAATGCCGTGCGCACCgg                                                               <  3:25317/149‑1 (MQ=255)                                                                                              cTGGAAGTCGGCGCGCCCCACATACTCCAGGGTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGT‑ACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCGGGGACAAAGCTCAATGCCGTGCGCACCGGATcc                                                           >  4:123306/1‑142 (MQ=255)                                                                                              cTGGAAGTCGGCGCGCCCCACATACTCCAGGGTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGT‑ACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCGGGGACAAAGCTCAATGCCGTGCGCACCGGATcc                                                           <  3:123306/142‑1 (MQ=255)                                                                                                          gcgcCCCACATACTCCAGGGTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGT‑ACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCGGGGACAAAGCTCAATGCCGTGCGCACCGGATCCCCAAGGTAACCGCGCCCGa                                        <  4:267190/149‑1 (MQ=255)                                                                                                            gcCCCACATACTCCAGGGTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGT‑ACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCGGGGACAAAGCTCAATGCCGTGCGCACCGGATCCCCAAGGTAACCGCGCCCGaca                                      <  8:91023/149‑1 (MQ=255)                                                                                                            gcCCCACATACTCCAGGGTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGT‑ACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCGGGGACAAAGCTCAATGCCGTGCGCACCGGATCCCCAAGGTAACCGCGCCCGaca                                      <  3:328940/149‑1 (MQ=255)                                                                                                                cacaTACTCCAGGGTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGT‑ACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCgg                                                                                           <  6:184553/92‑1 (MQ=255)                                                                                                                cacaTACTCCAGGGTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGT‑ACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCGGGGACAAAGCTCAATGCCGTGCGCACCGGATCCCCAAGGt                                                    >  4:46126/1‑131 (MQ=255)                                                                                                                cacaTACTCCAGGGTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGT‑ACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCGGGGACAAAGCTCAATGCCGTGCGCACCGGATCCCCAAGGt                                                    <  3:46126/131‑1 (MQ=255)                                                                                                                cacaTACTCCAGGGTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGAGGT‑ACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCgg                                                                                           >  5:184553/1‑92 (MQ=255)                                                                                                                 acaTACTCCAGGGTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGT‑ACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCGGGGACAAAGCTCAATGCCGTGCGCACCGGATCCCCAAGGTAACCGCGCCCGACACCTAc                                 <  6:39774/149‑1 (MQ=255)                                                                                                                 acaTACTCCAGGGTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGT‑ACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCGGGGACAAAGCTCAATGCCGTGCGCACCGGATCCCCAAGGTAACCGCGCCCGACACCTAc                                 <  5:213689/149‑1 (MQ=255)                                                                                                                  caTACTCCAGGGTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGT‑ACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCGGGGACAAAGCTCAATGcc                                                                          <  6:275424/107‑1 (MQ=255)                                                                                                                  caTACTCCAGGGTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGT‑ACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCGGGGACAAAGCTCAATGcc                                                                          >  5:275424/1‑107 (MQ=255)                                                                                                                      cTCCAGGGTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGT‑ACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCGGGGACAAAGCTCAATGCCGTGCGCACCGGATCCCCAAGGTAACCGCGCCCGACACCTACGCCg                             <  2:304412/148‑1 (MQ=255)                                                                                                                       tCCAGGGTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGT‑ACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCGGGGACAAAGCTCAATGCCGTGCGCACCGGATCCCCAAGGTAACCGCGCCCGACACCTACGCCGGc                           >  4:111040/1‑149 (MQ=255)                                                                                                                        ccAGGGTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGT‑ACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCGGGGACAAAGCTCAATGCCGTGCGCACCGGATCCCCAAGGTAACCGCGCCCGACACCTACGCCGGc                           >  5:229721/1‑148 (MQ=255)                                                                                                                          aGGGTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGT‑ACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCGGGGACAAAGCTCAATGCCGTGCGCACCGGATCCCCAAGGTAACCGCGCCCGACACCTACGCCGGCca                         >  1:264646/1‑148 (MQ=255)                                                                                                                             gTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGTAAGAGGCGCTCACCCGGGCGGCTGGCGAACGGATCGGGGACAAAGCTCAATGCCGTGCGCACCGGATCCCCaa                                                       >  1:245271/1‑116 (MQ=255)                                                                                                                             gTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGT‑ACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCGGGGACAAAGCTCAATGCCGTGCGCACCGGATCCCCaa                                                       <  2:245271/115‑1 (MQ=255)                                                                                                                                  cccATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGT‑ACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCGGGGACAAAGCTCAATGCCGTGCGCACCGGATCCCCAAGGTAACCGCGCCCGACACCTACGCCGGCCACATGCAACt                >  4:134520/1‑149 (MQ=255)                                                                                                                                   ccATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGT‑ACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCGGGGACAAAGCTCAATGCCGTGCGCACCGGATCCCCAAGGTAACCGCGCCCGACACCTACGCCGGCCACATGCAACTc               <  2:264646/149‑1 (MQ=255)                                                                                                                                    cATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGT‑ACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCGGGGACAAAGCTCAATGCCGTGCGCACCGGATCCCCAAGGTAACCGCGCCCGACACCTACGCCGGCCACATGCAACTCg              >  4:80919/1‑149 (MQ=255)                                                                                                                                    cATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGT‑ACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCGGGGACAAAGCTCAATGCCGTGCGCACCGGATCCCCAAGGTAACCGCGCCCGACACCTACGCCGGCCACATGCAACTCg              >  6:51427/1‑149 (MQ=255)                                                                                                                                    acctggCGCACGTGTGCGCCGGGTCGCCACTGCGGT‑ACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCGGGGACAAAGCTCAATGCCGTGCGCACCGGATCCCCAAGGTAACCGCGCCCGACACCTACGCCGGCCACATGCAACTCg              >  5:221978/5‑149 (MQ=255)                                                                                                                                        ggCGCACACGTGCGCCAGGTCGCCACTGCGGT‑ACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCGGGGACAAAGCTCAATGCCGTGCGCACCGGATCCCCAAGGTAACCGCGCCCGACACCTACGCCGGCCACATGCAACTCGCCgg          <  3:80919/149‑1 (MQ=255)                                                                                                                                        ggCGCACACGTGCGCCAGGTCGCCACTGCGGT‑ACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCGGGGACAAAGCTCAATGCCGTGCGCACCGGATCCCCAAGGTAACCGCGCCCGACACCTACGCCGGCCACATGCAACTCGCCg           <  2:139122/148‑1 (MQ=255)                                                                                                                                           gcACACGTGCGCCAGGTCGCCACTGCGGT‑ACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCGGGGACAAAGCTCAATGCCGTGCGCACCGGATCCCCAAGGTAACCGCGCCCGACACCTACGCCGGCCACATGCAACTCGCCGGTgg       >  1:60762/1‑149 (MQ=255)                                                                                                                                               acGTGCGCCAGGTCGCCACTGCGGT‑ACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCGGGGACAAAGCTCAATGCCGTGCGCACCGGATCCCCAAGGTAACCGCGCCCGACACCTACGCCGGCCACATGCAACTCGCCGGTGGCCCg   >  7:22699/1‑149 (MQ=255)                                                                                                                                                cGTGCGCCAGGTCGCCACTGCGGT‑ACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCGGGGACAAAGCTCAATGCCGTGCGCACCGGATCCCCAAGGTAACCGCGCCCGACACCTACGCCGGCCACATGCAACTCGCCGGTGGCCCgg  >  3:292183/1‑149 (MQ=255)                                                                                                                                                |                                                                                                                                                      CGCATCCACCACCACGGCGCTGTGCACACCTTCATGGGCCAGCAGGCGCGATTCGATTTCGCCCAGTTCTATGCGGTAGCCGCGGATCTTCACCTGGAAGTCGGCGCGCCCCACATACTCCAGGGTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGT‑ACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCGGGGACAAAGCTCAATGCCGTGCGCACCGGATCCCCAAGGTAACCGCGCCCGACACCTACGCCGGCCACATGCAACTCGCCGGTGGCCCGG  >  NC_002947/4779979‑4780270

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