Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,789,889	C→T	11.6%	L1317L (CTG→CTA) ‡	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,889	0	C	T	11.6%	111.7 / 4.2	43	L1317L (CTG→CTA) ‡	pvdI	non‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base C (20/18);  new base T (3/2);  total (23/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.78e-01

TGGGCACGCCGATGTTGATCTTCTCGTTGCCGCTGTAGCGCGCCAGCAGCAGGCCGAACGAGGCTAGGAAGACGTGGAACAGCGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAATCATCTTCCAGGTGCGCCTGCCAGTAGCCGAGCTGGGCTTGCATCTGGCCTTCGCTCAGGCGCTCGCGTTGCCAGC  >  NC_002947/4789741‑4790034                                                                                                                                                     |                                                                                                                                                  tGGGCACGCCGATGTTGATCTTCTCGTTGCCGCTGTAGCGCGCCAGCAGCAGGCCGAACGAGGCTAGGAAGACGTGGAACAGCGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTcc                                                                                                                                                   <  1:133590/149‑1 (MQ=255)       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                                                                                                          >  8:3021/1‑126 (MQ=255)                                                                 tAGGAAGACGTGGAACAGCGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCg                                                                                                          <  7:3021/126‑1 (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                                                                                                                                           <  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                               >  7:252484/1‑126 (MQ=255)                                                                                                                                            gcggACGTCTGGACGTGCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAATCATCTTCCAGGTGCGCCTGCCAGTAGCCGAGCTGGGCTTGCATCTg                               <  8:252484/126‑1 (MQ=255)                                                                                                                                                aCGTCTGGACGTGCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAATCATCTTCCAGGTGCGCCTGCCAGTAGCCGAGCTGGGCTTGCATCTGGCCTTCGCTCAGGCGCTCGCGTTGc      >  4:280571/1‑147 (MQ=255)                                                                                                                                                  gTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAATCATCTTCCAGGTGCGCCTGCCAGTAGCCGAGCTGGGCTTGCATCTGGCCTTCGCTCAGGCGCTCGCGTTGCCAGc  >  1:165435/1‑149 (MQ=255)                                                                                                                                                  gTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAATCATCTTCCAGGTGCGCCTGCCAGTAGCCGAGCTGGGCTTGCATCTGGCCTTCGCTCAGGCGCTCGCGTTGCCAGc  >  8:285284/1‑149 (MQ=255)                                                                                                                                                    ccAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAATCATCTTCCAGGTGCGCCTGCCAGTAGCCGAGCTGGGCTTGCATCTGGCCTTCGCTCGGGCGCTc             <  4:40889/136‑1 (MQ=255)                                                                                                                                                    ccAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAATCATCTTCCAGGTGCGCCTGCCAGTAGCCGAGCTGGGCTTGCATCTGGCCTTCGCTCGGGCGCTc             >  3:40889/1‑136 (MQ=255)                                                                                                                                                     |                                                                                                                                                  TGGGCACGCCGATGTTGATCTTCTCGTTGCCGCTGTAGCGCGCCAGCAGCAGGCCGAACGAGGCTAGGAAGACGTGGAACAGCGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAATCATCTTCCAGGTGCGCCTGCCAGTAGCCGAGCTGGGCTTGCATCTGGCCTTCGCTCAGGCGCTCGCGTTGCCAGC  >  NC_002947/4789741‑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: ‑