Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,789,896	C→G	12.8%	G1315A (GGA→GCA)	pvdI ←	non‑ribosomal peptide synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,789,896	0	C	G	12.8%	105.8 / 6.9	39	G1315A (GGA→GCA)	pvdI	non‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base C (17/17);  new base G (3/2);  total (20/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

GCCGATGTTGATCTTCTCGTTGCCGCTGTAGCGCGCCAGCAGCAGGCCGAACGAGGCTAGGAAGACGTGGAACAGCGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAATCATCTTCCAGGTGCGCCTGCCAGTAGCCGAGCTGGGCTTGCATCTGGCCTTCGCTCAGGCGCTCGCGTTGCCAGCTGGCATAGTC  >  NC_002947/4789748‑4790044                                                                                                                                                     |                                                                                                                                                     gCCGATGTTGATCTTCTCGTTGCCGCTGTAGCGCGCCAGCAGCAGGCCGAACGAGGCTAGGAAGACGTGGAACAGCGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTc                                                                                                                                                      >  3:184585/1‑149 (MQ=255)   cGATGTTGATCTTCTCGTTGCCGCTGTAGCGCGCCAGCAGCAGGCCGAACGAGGCTAGGAAGACGTGGAACAGCGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCg                                                                                                                                                    >  7:275032/1‑149 (MQ=255)     aTGTTGATCTTCTCGTTGCCGCTGTAGCGCGCCAGCAGCAGGCCGAACGAGGCTAGGAAGACGTGGAACAGCGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGcc                                                                                                                                                  <  4:184585/149‑1 (MQ=255)      tGTTGATCTTCTCGTTGCCGCTGTAGCGCGCCAGCAGCAGGCCGAACGAGGCTAGGAAGACGTGGAACAGCGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGccc                                                                                                                                                 <  1:192096/149‑1 (MQ=255)              ttctCGTTGCCGCTGTAGCGCGCCAGCAGCAGGCCGAACGAGGCTAGGAAGACGTGGAACAGCGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCt                                                                                                                                         >  7:202599/1‑149 (MQ=255)                ctcGTTGCCGCTGTAGCGCGCCAGCAGCAGGCCGAACGAGGCTAGGAAGACGTGGAACAGCGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCtt                                                                                                                                        <  8:275032/148‑1 (MQ=255)                 tcGTTGCCGCTGTAGCGCGCCAGCAGCAGGCCGAACGAGGCTAGGAAGACGTGGAACAGCGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTg                                                                                                                                      >  3:23523/1‑149 (MQ=255)                   gTTGCCGCTGTAGCGCGCCAGCAGCAGGCCGAACGAGGCTAGGAAGACGTGGAACAGCGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCttt                                                                                                                                       <  3:281160/146‑1 (MQ=255)                   gTTGCCGCTGTAGCGCGCCAGCAGCAGGCCGAACGAGGCTAGGAAGACGTGGAACAGCGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTg                                                                                                                                      <  3:24346/147‑1 (MQ=255)                   gTTGCCGCTGTAGCGCGCCAGCAGCAGGCCGAACGAGGCTAGGAAGACGTGGAACAGCGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCa                                                                                                                                    <  4:94682/149‑1 (MQ=255)                    ctgctgcTGGCGCGCTACAGCGGCAGGCCGAACGAGGCTAGGAAGACGTGGAACAGCGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCAg                                                                                                                                   >  5:312467/2‑149 (MQ=255)                                cgcgcCAGCAGCAGGCCGAACGAGGCTAGGAAGACGTGGAACAGCGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTgg                                                                                                                       >  2:156584/1‑149 (MQ=255)                                     cagcggcagGCCGAACGAGGCTAGGAAGACGTGGAACAGCGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCg                                                                                                                    <  6:312467/147‑1 (MQ=255)                                          gcagGCCGAACGAGGCTAGGAAGACGTGGAACAGCGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCa                                                                                                             <  2:258150/149‑1 (MQ=255)                                                          tAGGAAGACGTGGAACAGCGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCg                                                                                                                    <  7:3021/126‑1 (MQ=255)                                                          tAGGAAGACGTGGAACAGCGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCg                                                                                                                    >  8:3021/1‑126 (MQ=255)                                                                  cGTGGAACAGCGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAatca                                                                                     >  6:43733/1‑149 (MQ=255)                                                                    tGGAACAGCGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGc                                                                                                              <  8:205872/122‑1 (MQ=255)                                                                    tGGAACAGCGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGc                                                                                                              >  7:205872/1‑122 (MQ=255)                                                                    tGGAACAGCGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAatcatc                                                                                   <  1:156584/149‑1 (MQ=255)                                                                    tGGAACAGCGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAatcatc                                                                                   <  5:43733/149‑1 (MQ=255)                                                                        aCAGCGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTcc                                                                                                                                                     >  7:61976/1‑79 (MQ=255)                                                                        aCAGCGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTcc                                                                                                                                                     <  8:61976/79‑1 (MQ=255)                                                                                 cGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAATCATCTTCCAGGTGCGcc                                                                      >  8:287682/1‑149 (MQ=255)                                                                                  gTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAATCATCTTCCAGGTGCGCCt                                                                     >  6:236463/1‑149 (MQ=255)                                                                                    tGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAATCATCTTCCAGGTGCGCCTGc                                                                   <  4:23523/149‑1 (MQ=255)                                                                                      cAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAATCATCTTCCAGGTGCGCCTGc                                                                   >  1:192995/1‑147 (MQ=255)                                                                                               gcagtgcgcTGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAATCATCTTCCAGGTGCGCCTGCCAGTAGCCGAg                                                        <  4:18112/144‑1 (MQ=255)                                                                                                              tCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAATCATCTTCCAGGTGCGCCTGCCAGTAGCCGAGCTGGGCTTGCATCTg                                         <  8:202599/149‑1 (MQ=255)                                                                                                               ccgccgTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAATCATCTTCCAGGTGCGCCTGCCAGTAGCCGAGCTGGGCTTGCATCTgg                                        >  4:259463/1‑149 (MQ=255)                                                                                                                  ccgTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAATCATCTTCCAGGTGCGCCTGCCAGTAGCCGAGCTGGGCTTGCATCTGGCCt                                     >  7:89525/1‑149 (MQ=255)                                                                                                                         cTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAATCATCTTCCAGGTGCGCCTGCCAGTAGCCGAGCTGGGCTTGCATCTGGCCTTCGCt                                <  6:179518/147‑1 (MQ=255)                                                                                                                         cTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAATCATCTTCCAGGTGCGCCTGCCAGTAGCCGAGCTGGGCTTGCATCTGGCCTTCGCt                                >  5:179518/1‑147 (MQ=255)                                                                                                                                    ggCGGACGTCTGGACGTGCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAATCATCTTCCAGGTGCGCCTGCCAGTAGCCGAGCTGGGCTTGCATCTGGCCTTCGCTCAGGCGCTCGCGt                   >  1:122565/1‑149 (MQ=255)                                                                                                                                     gcggACGTCTGGCCGTGCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAATCATCTTCCAGGTGCGCCTGCCAGTAGCCGAGCTGGGCTTGCATCTGGCCTTCGCTCAGGCGCTCGCGtt                  <  2:224196/149‑1 (MQ=255)                                                                                                                                     gcggACGTCTGGACGTGCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAATCATCTTCCAGGTGCGCCTGCCAGTAGCCGAGCTGGGCTTGCATCTg                                         <  8:252484/126‑1 (MQ=255)                                                                                                                                     gcggACGTCTGGACGTGCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAATCATCTTCCAGGTGCGCCTGCCAGTAGCCGAGCTGGGCTTGCATCTg                                         >  7:252484/1‑126 (MQ=255)                                                                                                                                         aCGTCTGGACGTGCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAATCATCTTCCAGGTGCGCCTGCCAGTAGCCGAGCTGGGCTTGCATCTGGCCTTCGCTCAGGCGCTCGCGTTGc                >  4:280571/1‑147 (MQ=255)                                                                                                                                           gTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAATCATCTTCCAGGTGCGCCTGCCAGTAGCCGAGCTGGGCTTGCATCTGGCCTTCGCTCAGGCGCTCGCGTTGCCAGc            >  8:285284/1‑149 (MQ=255)                                                                                                                                           gTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAATCATCTTCCAGGTGCGCCTGCCAGTAGCCGAGCTGGGCTTGCATCTGGCCTTCGCTCAGGCGCTCGCGTTGCCAGc            >  1:165435/1‑149 (MQ=255)                                                                                                                                             ccAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAATCATCTTCCAGGTGCGCCTGCCAGTAGCCGAGCTGGGCTTGCATCTGGCCTTCGCTCGGGCGCTc                       >  3:40889/1‑136 (MQ=255)                                                                                                                                             ccAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAATCATCTTCCAGGTGCGCCTGCCAGTAGCCGAGCTGGGCTTGCATCTGGCCTTCGCTCGGGCGCTc                       <  4:40889/136‑1 (MQ=255)                                                                                                                                                     ccGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAATCATCTTCCAGGTGCGCCTGCCAGTAGCCGAGCTGGGCTTGCATCTGGCCTTCGCTCAGGCGCTCGCGTTGCCAGCTGGCATAGTc  >  4:133104/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                     GCCGATGTTGATCTTCTCGTTGCCGCTGTAGCGCGCCAGCAGCAGGCCGAACGAGGCTAGGAAGACGTGGAACAGCGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAATCATCTTCCAGGTGCGCCTGCCAGTAGCCGAGCTGGGCTTGCATCTGGCCTTCGCTCAGGCGCTCGCGTTGCCAGCTGGCATAGTC  >  NC_002947/4789748‑4790044

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