Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,789,890	A→G	11.1%	L1317P (CTG→CCG) ‡	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,890	0	A	G	11.1%	116.2 / 5.5	45	L1317P (CTG→CCG) ‡	pvdI	non‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base A (21/19);  new base G (3/2);  total (24/21)
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.15e-01

CGCCGATGTTGATCTTCTCGTTGCCGCTGTAGCGCGCCAGCAGCAGGCCGAACGAGGCTAGGAAGACGTGGAACAGCGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAATCATCTTCCAGGTGCGCCTGCCAGTAGCCGAGCTGGGCTTGCATCTGGCCTTCGCTCAGGCGCTCGCGTTGCCAGC  >  NC_002947/4789747‑4790034                                                                                                                                                |                                                                                                                                                 cGCCGATGTTGATCTTCTCGTTGCCGCTGTAGCGCGCCAGCAGCAGGCCGAACGAGGCTAGGAAGACGTGGAACAGCGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGt                                                                                                                                             >  4:24346/1‑149 (MQ=255)  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)                    gTTGCCGCTGTAGCGCGCCAGCAGCAGGCCGAACGAGGCTAGGAAGACGTGGAACAGCGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAg                                                                                                                                                 >  4:279642/1‑126 (MQ=255)                    gTTGCCGCTGTAGCGCGCCAGCAGCAGGCCGAACGAGGCTAGGAAGACGTGGAACAGCGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAg                                                                                                                                                 <  3:279642/126‑1 (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                                                                                                    >  7:205872/1‑122 (MQ=255)                                                                     tGGAACAGCGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGc                                                                                                    <  8:205872/122‑1 (MQ=255)                                                                     tGGAACAGCGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAatcatc                                                                         <  5:43733/149‑1 (MQ=255)                                                                     tGGAACAGCGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAatcatc                                                                         <  1:156584/149‑1 (MQ=255)                                                                         aCAGCGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTcc                                                                                                                                           <  8:61976/79‑1 (MQ=255)                                                                         aCAGCGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTcc                                                                                                                                           >  7:61976/1‑79 (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  >  1:165435/1‑149 (MQ=255)                                                                                                                                            gTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAATCATCTTCCAGGTGCGCCTGCCAGTAGCCGAGCTGGGCTTGCATCTGGCCTTCGCTCAGGCGCTCGCGTTGCCAGc  >  8:285284/1‑149 (MQ=255)                                                                                                                                              ccAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAATCATCTTCCAGGTGCGCCTGCCAGTAGCCGAGCTGGGCTTGCATCTGGCCTTCGCTCGGGCGCTc             >  3:40889/1‑136 (MQ=255)                                                                                                                                              ccAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAATCATCTTCCAGGTGCGCCTGCCAGTAGCCGAGCTGGGCTTGCATCTGGCCTTCGCTCGGGCGCTc             <  4:40889/136‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                 CGCCGATGTTGATCTTCTCGTTGCCGCTGTAGCGCGCCAGCAGCAGGCCGAACGAGGCTAGGAAGACGTGGAACAGCGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAATCATCTTCCAGGTGCGCCTGCCAGTAGCCGAGCTGGGCTTGCATCTGGCCTTCGCTCAGGCGCTCGCGTTGCCAGC  >  NC_002947/4789747‑4790034

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