Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,787,291	C→T	7.8%	R8H (CGC→CAC)	pvdJ ←	non‑ribosomal peptide synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,787,291	0	C	T	7.8%	146.8 / 2.4	51	R8H (CGC→CAC)	pvdJ	non‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base C (26/21);  new base T (2/2);  total (28/23)
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.12e-01

CTGGGCATAGGACAGTGGGATACGCGCCACGTCATGACGGCTGGCGAC‑AATTGGCAGCAGGCGGAAGCTCTGGCCGCCAGCGCTCATTTTTTCGACCACCTGGCGG‑CGCTGCTCCAGCGGCATTGCTACGAAACGTCTGGCGATGCGTTCTGCGGTTGTCTT‑GTCCATTTCAAGCCTCAGTCAATTCATTCATCAGTTGTTCGATACCGGCGAGGCCTTCTTCGGTCAGGGCCTGGCCCGCGCCTTCCAGCGCCTGGGCAAAGGCCCGCAGAACCGGTGTTTCGAAGATGAGA  >  NC_002947/4787146‑4787438                                                                                                                                                    |                                                                                                                                                     ctgGGCATAGGACAGTGGGATACGCGCCACGTCATGACGGCTGGCGACAAATTGCCAGCAGCCGGAACAAATGGCCGCCAGGGGTCATTTTTACGACAACATGCCGG‑CCATGCGCCACCGGCATTGCTACGAAACTAGTGCCGATGCGt                                                                                                                                                    >  1:63899/1‑149 (MQ=255) ctgGGCATAGGACAGTGGGATACGCGCCACGTCATGACGGCTGGCGAC‑AATTGGCAGCAGGCGGAAGCTCTGGCCGCCAGCGCTCATTTTTTCGACCACCTGGCGG‑CGCTGCTCCAGCGGCATTGCTACGAAACGTCTGGCGATGCGtt                                                                                                                                                   <  4:258021/149‑1 (MQ=255)    ggCATAGGACAGTGGGATACGCGCCACGTCATGACGGCTGGCGAC‑AATTGGCAGCAGGCGGAAGCTCTGGCCGCCAGCGCTCATTTTTTCGACCACCTGGCGG‑CGCTGCTCCAGCGGCATTGCTACGAAACGTCTGGCGATGCGTTCt                                                                                                                                                 <  2:239394/148‑1 (MQ=255)           gACAGTGGGATACGCGCCACGTCATGACGGCTGGCGAC‑AATTGGCAGCAGGCGGAAGCTCTGGCCGCCAGCGCTCATTTTTTCGACCACCTGGCGG‑CGCTGCTCCAGCGGCATTGCTACGAAACGTCTGGCGATGCGTTCTGCGGttgt                                                                                                                                         >  2:84763/1‑149 (MQ=255)            aCAGTGGGATACGCGCCACGTCATGACGGCTGGCGAC‑AATTGGCAGCAGGCGGAAGCTCTGGCCGCCAGCGCTCATTTTTTCGACCACCTGGCGG‑CGCTGCTCCAGCGGCATTGCTACGAAACGTCTGGCGATGCGTTCTGCGGttgtc                                                                                                                                        <  2:157819/149‑1 (MQ=255)                tGGGATACGCGCCACGTCATGACGGCTGGCGAC‑AATTGGCAGCAGGCGGAAGCTCTGGCCGCCAGCGCTCATTTTTTCGACCACCTGGCGG‑CGCTGCTCCAGCGGCATTGCTACGAAACGTCTGGCGATGCGTTCTGCGGttgtctt‑gt                                                                                                                                   >  7:225179/1‑149 (MQ=255)                       cgcgCCACTTCATGACGGCTGGCGAC‑AATTGGCAGCAGGCGGAAGCTCTGGCCGCCAGCGCTCATTTTTTCGACCACCTGGCGG‑CGCTGCTCCAGCGGCATTGCTACGAAACGTCTGGCGATGCGTTCTGCGGTTGTCTT‑GTCCAttt                                                                                                                             <  1:92106/148‑1 (MQ=255)                       cgcgCCACGTCATGACGGCTGGCGAC‑AATTGGCAGCAGGCGGAAGCTCTGGCCGCCAGCGCTCATTTTTTCGACCACCTGGCGG‑CGCTGCTCCAGCGGCATTGCTACGAAACGTCTGGCGATGCGTTCTGCGGTTGTCTT‑GTCCAttt                                                                                                                             <  8:225179/148‑1 (MQ=255)                         cgcCACGTCATGACGGCTGGCGAC‑AATTGGCAGCAGGCGGAAGCTCTGGCCGCCAGCGCTCATTTTTTCGACCACCTGGCGG‑CGCTGCTCCAGCGGCATTTCTACGAAACGTCTCGCGATGCGTTCTGCGGTTGTCTT‑GTCCATTTCaa                                                                                                                          <  8:101323/149‑1 (MQ=255)                         cgcCACGTCATGACGGCTGGCGAC‑AATTGGCAGCAGGCGGAAGCTCTGGCCGCCAGCGCTCATTTTTTCGACCACCTGGCGG‑CGCTGCTCCAGCGGCATTGCTACGAAACGTCTGGCGATGCGTTCTGCGGTTGTCTT‑GTCCATTTCaa                                                                                                                          <  2:63899/149‑1 (MQ=255)                               gTCATGACGGCTGGCGAC‑AATTGGCAGCAGGCGGAAGCTCTGGCCGCCAGCGCTCATTTTTTCGACCACCTGGCGG‑CGCTGCTCCAGCGGCATTGCTACGAAACGTCTGGCGATGCGTTCTGCGGTTGTCTT‑GTCCATTTCAAGCCt                                                                                                                      >  7:140030/1‑147 (MQ=255)                               gTCATGACGGCTGGCGAC‑AATTGGCAGCAGGCGGAAGCTCTGGCCGCCAGCGCTCATTTTTTCGACCACCTGGCGG‑CGCTGCTCCAGCGGCATTGCTACGAAACGTCTGGCGATGCGTTCTGCGGTTGTCTT‑GTCCATTTCAAGCCTCa                                                                                                                    <  2:203006/149‑1 (MQ=255)                                tCATGACGGCTGGCGAC‑AATTGGCAGCAGGCGGAAGCTCTGGCCGCCAGCGCTCATTTTTTCGACCACCTGGCGG‑CGCTGCTCCAGCGGCATTGCAACGAAACGTCTGGCGATGCGTTCTGCGGTTGTCTT‑GTCCATTTCAAGCCTCAg                                                                                                                   >  5:44148/1‑149 (MQ=255)                                  aTGACGGCTGGCGAC‑AATTGGCAGCAGGCGGAAGCTCTGGCCGCCAGCGCTCATTTTTTCGACCACCTGGCGG‑CGCTGCTCCAGCGGCATTGCTACGAAACGTCTGGCGATGCGTTCTGCGGTTGTCTT‑GTCCATTTCAAGCCTCAGt                                                                                                                  >  5:120565/1‑148 (MQ=255)                                    gACGGCTGGCGAC‑AATTGGCAGCAGGCGGAAGCTCTGGCCGCCAGCGCTCATTTTTTCGACCACCTGGCGG‑CGCTGCTCCAGCGGCATTGCTACGAAACGTCTGGCGATGCGTTCTGCGGtt                                                                                                                                           <  1:195895/122‑1 (MQ=255)                                    gACGGCTGGCGAC‑AATTGGCAGCAGGCGGAAGCTCTGGCCGCCAGCGCTCATTTTTTCGACCACCTGGCGG‑CGCTGCTCCAGCGGCATTGCTACGAAACGTCTGGCGATGCGTTCTGCGGtt                                                                                                                                           >  2:195895/1‑122 (MQ=255)                                          tGGCGAC‑AATTGGCAGCAGGCGGAAGCTCTGGCCGCCAGCGCTCATTTTTTCGACCACCTGGCGG‑CGCTGCTCCAGCGGCATTGCTACGAAACGTCTGGCGATGCGTTCTGCGGTTGTCTT‑GTCCATTTCAAGCCTCAGTCAattcat                                                                                                          >  3:145216/1‑148 (MQ=255)                                              gAC‑AATTGGCAGCAGGCGGAAGCTCTGGCCGCCAGCGCTCATTTTTTCGACCACCTGGCGG‑CTCTGCTCCAGCGGCATTGCTACGAAACGTCTGGCGATGCGTTCTGCGGTTGTCTT‑GTCCATTTCAAGCCTCAGTCAATTCATtcatc                                                                                                     <  2:236891/149‑1 (MQ=255)                                              gAC‑AATTGGCAGCAGGCGGAAGCTCTGGCCGCCAGCGCTCATTTTTTCGACCACCTGGCGG‑CGCTGCTCCAGCGGCATTGCTACGAAACGTCTGGCGATGTGTTCTGCGGTTGTCTT‑GTCCATTTCAAGCCTCAGTCAATtcattcat                                                                                                      >  1:146375/1‑148 (MQ=255)                                              gAC‑AATTGGCAGCAGGCGGAAGCTCTGGCCGCCAGCGCTCATTTTTTCGACCACCTGGCGG‑CGCTGCTCCAGCGGCATTGCTACGAAACGTCTGGCGATGCGTTCTGGGGTTGTCTT‑GTCCATTTCAAGCCTCAGTCAATtcattcat                                                                                                      >  5:232538/1‑148 (MQ=255)                                               aC‑AATTGGCAGCAGGCGGAAGCTCTGGCCGCCAGCGCTCATTTTTTCGACCACCTGGCGG‑CGCTGCTCCAGCGGCATTGCTACGAAACGTCTGGCGATGTGTTCTGCGGTTGTCTT‑GTCCATTTCAAGCCTCAGTCAATTCATtcatca                                                                                                    <  2:146375/149‑1 (MQ=255)                                               aC‑AATTGGCAGCAGGCGGAAGCTCTGGCCGCCAGCGCTCATTTTTTCGACCACCTGGCGG‑CGCTGCTCCAGCGGCATTGCTACGAAACGTCTGGCGATGCGTTCTGCGGTTGTCTT‑GTCCATTTCAAGCCTCAGTCAATtcattcat                                                                                                      >  5:256036/1‑147 (MQ=255)                                                       gcagcaGGCGGAAGCTCTGGCCGCCAGCGCTCATTTTTTCGACCACCTGGCGG‑CGCTGCTCCAGCGGCATTGCTACGAAACGTCTGGCGATGCGTTCTGCGGTTGTCTT‑GTCCATTTCAAGCCTCAGTCAATTCATTCATCAGTTGTTc                                                                                             <  8:140030/149‑1 (MQ=255)                                                            agGCGGAAGCTCTGGCCGCCAGCGCTCATTTTTTCGACCACCTGGCGG‑CGCTGCTCCAGCGGCATTGCTACGAAACGTCTGGCGATGCGTTCTGCGGTTGTCTT‑GTCCATTTCAAGCCTCAGTCAATTCATTCATCAGTTGTTCGATAc                                                                                        >  8:80003/1‑149 (MQ=255)                                                              gCGGAAGCTCTGGCCGCCAGCGCTCATTTTTTCGACCACCTGGCGG‑CGCTGCTCCAGCGGCATTGCTACGAAACGTCTGGCGATGCGTTCTGCGGTTGTCTT‑GTCCATTTCAAGCCTCAGTCAATTCATTCATCAGTTGTTCGATACCg                                                                                      <  6:167855/149‑1 (MQ=255)                                                              gCGGAAGCTCTGGCCGCCAGCGCTCATTTTTTCGACAACCTGGCGG‑CGCTGCTCCAGCGGCATTGCTACGAAACGTCTGGCGATGCGTTCTGCGGtg                                                                                                                                           >  2:216724/1‑96 (MQ=255)                                                                 gAAGCTCTGGCCGCCAGCGCTCATTTTTTCGACCACCTGGCGG‑CGCTGCTCCAGCGGCATTGCTACGAAACGTCTGGCGATGCGTTCTGCGGTTGTCTT‑GTCCAGTTCAAGCCTCAGTCAATTCATTCATCAGTTGTTCGATACCGGCg                                                                                   >  7:153585/1‑149 (MQ=255)                                                                  aaGCTCTGGCCGCCAGCGCTCATTTTTTCGACCACCTGGCGG‑CGCTGCTCCAGCGGCATTGCTACGAAACGTCTGGCGATGCGTTCTGCGGTTGTCTT‑GTCCATTTCAAGCCTCAGTCAATTCATTCATCAGTTGTTCGATACCGGCGa                                                                                  <  4:145216/149‑1 (MQ=255)                                                                   aGCTCTGGCCGCCAGCGCTCATTTTTTCGACCACCTGGCGG‑CGCTGCTCCAGCGGCATTGCTACGAAACGTGTGGCGATGCGTTCTGCGGTTGTCTT‑GTCCATTTCAAGCCTCAGTCAATTCATTCATCAGTTGTTCGATACCGGCGAg                                                                                 <  1:84763/149‑1 (MQ=255)                                                                           ccgccAGCGCTCATTTTTTCGACCACCTGGCGGACGCTGCTCCAGCGGCATTGCTACGAAACGTCTGGCGATGCGTTCTGCGGTTGTCTT‑GTCCATTTCAAGCCTCAGTCAATTCATTCATCAGTTGTTCGATACCGGCGAGGCcttct                                                                          >  1:37071/1‑149 (MQ=255)                                                                             gccAGCGCTCATTTTTTCGACCACCTGGCGG‑CGCTGCTCCAGCGGCATTGCTACGAAACGTCTGGCGATGCGTTCTGCGGTTGTCTT‑GTCCATTTCAAGCCTCAGTCAATTCATTCATCAGTTGTTCGATACCGGCGAGGCCTTCTTCg                                                                       <  7:80003/149‑1 (MQ=255)                                                                             gccAGCGCTCATTTTTTCGACCACCTGGCGG‑CGCTGCTCCAGCGGCATTGCTACGAAACGTCTGGCGATGCGTTCTGCGGTTGTCTT‑GTCCATTTCAAGCCTCAGTCAATTCATTCATCAGTTGTTCGATACCGGCGAGGCCTTCTTCg                                                                       <  6:256036/149‑1 (MQ=255)                                                                                   gcTCATTTTTTCGACCACCTGGCGG‑CGCTGCTCCAGCGGCATTGCTACGAAACGTCTGGCGATGCGTTCTGCGGTTGTCTT‑GTCCATTTCAAGCCTCAGTCAATTCATTCATCAGTTGTTCGATACCGGCGAGGCCTTCTTCGGTCAgg                                                                 <  2:37071/149‑1 (MQ=255)                                                                                           tttCGACCACCTGGCGG‑CGCTGCTCCAGCGGCATTGCTACGAAACGTCTGGCGATGCGTTCTGCGGTTGTCTT‑GTCCATTTCAAGCCTCAGTCAATTCATTCATCAGTTGTTCGATACCGGCGAGGCCTTCTTCGGTCAGGGCCTGGcc                                                         >  8:58052/1‑149 (MQ=255)                                                                                           tttCGACCACCTGGCGG‑CGCTGCTCCAGCGGCATTGCTACGAAACGTCTGGCGATGCGTTCTGCGGTTGTCTT‑GTCCATTTCAAGCCTCAGTCAATTCATTCATCAGTTGTTCGATACCGGCGAGGCCTTCTTCGGTCAGGGCCTGGcc                                                         >  4:195775/1‑149 (MQ=255)                                                                                             tCGACCACCTGGCGG‑CGCTGCTCCAGCGGCATTGCTACGAAACGTCTGGCGATGCGTTCTGCGGTTGTCTT‑GTCCATTTCAAGCCTCAGTCAATTCATTCATCAGTTGTTCGATACCGGCGAGGCCTTCTTCGGTCAGGGCCTGGCCcg                                                       >  2:266020/1‑149 (MQ=255)                                                                                                  caccTGGCGG‑CGCTGCTCCAGCGGCATTGCTACGAAACGTCTGGCGATGCGTTCTGCGGTTGTCTT‑GTCCATTTCAAGCCTCAGTCAATTCATTCATCAGTTGTTCGATACCGGCGAGGCCTTCTTCGGTCAGGGCCTGGCCCGCGCCt                                                  >  7:107702/1‑149 (MQ=255)                                                                                                    ccTGGCGG‑CGCTGCTCCAGCGGCATTGCTACGAAACGTCTGGCGATGTGTTCTGCGGTTGTCTT‑GTCCATTTCAAGCCTCAGTCAATTCATTCATCAGTTGTTCGATACCGGCGAGGCCTTCTTCGGTCAGGGCCTGGCCCGCGCCTTc                                                >  7:215404/1‑149 (MQ=255)                                                                                                     cTGGCGG‑CGCTGCTCCAGCGGCATTGCTACGAAACGTCTGGCGATGCGTTCTGCGGTTGTCTT‑GTCCATTTCAAGCCTCAGTCAATTCATTCATCAGTTGTTCGATACCGGCGAg                                                                                 <  4:180588/115‑1 (MQ=255)                                                                                                     cTGGCGG‑CGCTGCTCCAGCGGCATTGCTACGAAACGTCTGGCGATGCGTTCTGCGGTTGTCTT‑GTCCATTTCAAGCCTCAGTCAATTCATTCATCAGTTGTTCGATACCGGCGAg                                                                                 >  3:180588/1‑115 (MQ=255)                                                                                                        gcgg‑cgCTGCTCCAGCGGCATTGCTACGAAACGTCTGGCGATGCGTTCTGCGGTTGTCTT‑GTCCATTTCAAGCCTCAGTCAATTCATTCATCAGTTGTTCGATACCGGCGAGGCCTTCTTCGGTCAGGGCCTGGCCCGCGCCTTCCAgc                                            >  2:82153/1‑149 (MQ=255)                                                                                                               ctgctCCAGCGGCATTGCTACGAAACGTCTGGCGATGCGTTCTGCGGTTGTCTT‑GTCCATTTCAAGCCTCAGTCAATTCATTCATCAGTTGTTCGATACCGGCGAGGCCTTCTTCGGTCAGGGCCTGGCCCGCGCCTTCCAGCGCCTgg                                      <  1:266020/149‑1 (MQ=255)                                                                                                                 gctCCAGCGGCATTGCTACGAAACGTCTGGCGATGCGTTCTGCGGTTGTCTT‑GTCCATTTCAAGCCTCAGTCAATTCATTCATCAGTTGTTCGATACCGGCGAGGCCTTCTTCGGTCAGGGCCTGGCCCGCGCCTTCCAGCGCCTggg                                     >  7:52917/1‑148 (MQ=255)                                                                                                                   tCCAGCGGCATTGCTACGAAACGTCTGGCGATGCGTTCTGCGGTTGTCTT‑GTCCATTTCAAGCCTCAGTCAATTCATTCATCAGTTGTTCGATACCGGCGAGGCCTTCTTCGGTCAGGGCCTGGCCCGCGCCTTCCAGCGCCTGGGCa                                   >  7:149527/1‑148 (MQ=255)                                                                                                                         ggCATTGCTACGAAACGTCTGGCGATGCGTTCTGCGGTTGTCTT‑GTCCATTTCAAGCCTCAGTCAATTCATTCATCAGTTGTTCGATACCGGCGAGGCCTTCTTCGGTCAGGGCCTGGCCCGCGCCTTCCAGCGCCTGGGCAAAGGccc                            <  6:232538/149‑1 (MQ=255)                                                                                                                                cTACGAAACGTCTGGCGATGCGTTCTGCGGTTGTCTT‑GTCCATTTCAAGCCTCAGTCAATTCATTCATCAGTTGTTCGATACCGGCGAGGCCTTCTTCGGTCAGGGCCTGGCCCGCGCCTTCCAGCGCCTGGGCAAAGGCCCGCAGAAc                     >  6:121334/1‑149 (MQ=255)                                                                                                                                  aCGAAACGTCTGGCGATGCGTTCTGCGGTTGTCTT‑GTCCATTTCAAGCCTCAGTCAATTCATTCATCAGTTGTTCGATACCGGCGAGGCCTTCTTCGGTCAGGGCCTGGCCCGCGCCTTCCAGCGCCTGGGCAAAGGCCCGCAGAACCg                   >  4:103440/1‑149 (MQ=255)                                                                                                                                     aaaCGTCTGGCGATGTGTTCTGCGGTTGTCTT‑GTCCATTTCAAGCCTCAGTCAATTCATTCATCAGTTGTTCGATACCGGCGAGGCCTTCTTCGGTCAGGGCCTGGCCCGCGCCTTCCAGCGCCTGGGCAAAGGCCCGCAGAACCGgtg                <  8:215404/149‑1 (MQ=255)                                                                                                                                     aaaCGTCTGGCGATGCGTTCTGCGGTTGTCTT‑GTCCATTTCAAGCCTCAGTCAATTCATTCATCAGTTGTTCGATACCGGCGAGGCCTTCTTCGGTCAGGGCCTGGCCCGCGCCTTCCAGCGCCTGGGCAAAGGCCCGCAGAACCGgtg                >  2:102978/1‑149 (MQ=255)                                                                                                                                         gTCTGGCGATGCGTTCTGCGGTTGTCTTAGTCCATTTCAAGCCTCAGTCAATTCATTCATCAGTTGTTCGATACCGGCGAGGCCTTCTTCGGTCAGGGCCTGGCCCGCGCCTTCCAGCGCCTGGGCAAAGGCCCGCAGAACCGGTGttt             >  1:153655/1‑149 (MQ=255)                                                                                                                                               cGATGCGTTCTGCGGTTGTCTT‑GTCCATTTCAAGCCTCAGTCAATTCATTCATCAGTTGTTCGATACCGGCGAGGCCTTCTTCGGTCAGGGCCTGGCCCGCGCCTTCCAGCGCCTGGGCAAAGGCCCGCAGAACCGGTGTTTCGAAGAt      <  1:102978/149‑1 (MQ=255)                                                                                                                                                    cGTTCTGCGGTTGTCTT‑GTCCATTTCAAGCCTCAGTCAATTCATTCATCAGTTGTTCGATACCGGCGAGGCCTTCTTCGGTCAGGGCCTGGCCCGCGCCTTCCAGCGCCTGGGCAAAGGCCCGCAGAACCGGTTTTTCGAAGATgaga  <  1:82153/148‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                     CTGGGCATAGGACAGTGGGATACGCGCCACGTCATGACGGCTGGCGAC‑AATTGGCAGCAGGCGGAAGCTCTGGCCGCCAGCGCTCATTTTTTCGACCACCTGGCGG‑CGCTGCTCCAGCGGCATTGCTACGAAACGTCTGGCGATGCGTTCTGCGGTTGTCTT‑GTCCATTTCAAGCCTCAGTCAATTCATTCATCAGTTGTTCGATACCGGCGAGGCCTTCTTCGGTCAGGGCCTGGCCCGCGCCTTCCAGCGCCTGGGCAAAGGCCCGCAGAACCGGTGTTTCGAAGATGAGA  >  NC_002947/4787146‑4787438

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