Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,780,121	A→G	10.0%	V2398A (GTG→GCG) ‡	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,121	0	A	G	10.0%	133.2 / 5.9	50	V2398A (GTG→GCG) ‡	pvdJ	non‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base A (23/22);  new base G (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 = 8.04e-01

CGCATCCACCACCACGGCGCTGTGCACACCTTCATGGGCCAGCAGGCGCGATTCGATTTCGCCCAGTTCTATGCGGTAGCCGCGGATCTTCACCTGGAAGTCGGCGCGCCCCACATACTCCAGGGTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGT‑ACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCGGGGACAAAGCTCAATGCCGTGCGCACCGGATCCCCAAGGTAACCGCGCCCGACACCTACGCCGGCCACATGCAACTCGCCGGTGGCCCG  >  NC_002947/4779979‑4780269                                                                                                                                               |                                                                                                                                                      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                                                   <  3:46126/131‑1 (MQ=255)                                                                                                                cacaTACTCCAGGGTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGT‑ACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCGGGGACAAAGCTCAATGCCGTGCGCACCGGATCCCCAAGGt                                                   >  4:46126/1‑131 (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                                                                         >  5:275424/1‑107 (MQ=255)                                                                                                                  caTACTCCAGGGTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGT‑ACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCGGGGACAAAGCTCAATGcc                                                                         <  6:275424/107‑1 (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)                                                                                                                                               |                                                                                                                                                      CGCATCCACCACCACGGCGCTGTGCACACCTTCATGGGCCAGCAGGCGCGATTCGATTTCGCCCAGTTCTATGCGGTAGCCGCGGATCTTCACCTGGAAGTCGGCGCGCCCCACATACTCCAGGGTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGT‑ACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCGGGGACAAAGCTCAATGCCGTGCGCACCGGATCCCCAAGGTAACCGCGCCCGACACCTACGCCGGCCACATGCAACTCGCCGGTGGCCCG  >  NC_002947/4779979‑4780269

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