Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,787,439	A→G	9.3%	L2149P (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,787,439	0	A	G	9.3%	117.9 / 3.0	45	L2149P (CTG→CCG)	pvdI	non‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base A (20/20);  new base G (2/2);  total (23/22)
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.93e-01

GCGCTGCTCCAGCGGCATTGCTACGAAACGTCTGGCGATGCGTTCTGCGGTTGTCTTGTCCATTTCAAGCCTCAGTCAATTCATTCATCAGTTGTTCGATACCGGCGAGGCCTTCTTCGGTCAG‑GGCCTGGCCCGCGCCTTCCAGCGCCTGGGCAAAGGCCCGCAGAACCGGTGT‑TTCGAAGATGAGACGGATCGGCATATCGATGCCCAGTTGGCTGTTGATACGTGAGAGGGCCTGGGCCGCCAGCAGCGAGTGACCACCCAGTTCAAAGAAGTTGTCGTCCAGCCCGA  >  NC_002947/4787296‑4787586                                                                                                                                                 |                                                                                                                                                     gcgcTGCTCCAGCGGCATTGCTACGAAACGTCTGGCGATGCGTTCTGCGGTTGTCTTGTCCATTTCAAGCCTCAGTCAATTCATTCATCAGTTGTTCGATACCGGCGAGGCCTTCTTCGGTCAG‑GGCCTGGCCCGCGCCTTCCAGCGCCt                                                                                                                                                <  1:594108/150‑1 (MQ=255)           aGCGGCATTGCTACGAAACGTCTGGCGATGCGTTCTGCGGTTGTCTTGTCCATTTCAAGCCTCAGTCAATTCATTCATCAGTTGTTCGATACCGGCGAGGCCTTCTTCGGTCAG‑GGCCTGGCCCGCGCCTTCCAGCGCCTGGGCAAAGGc                                                                                                                                      <  2:1001670/150‑1 (MQ=255)             cGGCATTGCTACGAAACGTCTGGCGATGCGTTCTGCGGTTGTCTTGTCCATTTCAAGCCTCAGTCAATTCATTCATCAGTTGTTCGATACCGGCGAGGCCTTCTTCGGTCAG‑GGCCTGGCCCGCGCCTTCCAGCGCCTGGGCAAAGGccc                                                                                                                                    <  1:829849/150‑1 (MQ=255)                        cGAAACGTCTGGCGATGCGTTCTGCGGTTGTCTTGTCCATTTCAAGCCTCAGTCAATTCATTCATCAGTTGTTCGATACCGGCGAGGCCTTCTTCGGTCAG‑GGCCTGGCCCGCGCCTTCCAGCGCCTGGGCAAAGGCCCGCAGAACCGgt                                                                                                                         <  1:225637/150‑1 (MQ=255)                                   gCGATGCGTTCTGCGGTTGTCTTGTCCATTTCAAGCCTCAGTCAATTCATTCATCAGTTGTTCGATACCGGCGAGGCCTTCTTCGGTCAG‑GGCCTGGCCCGCGCCTTCCAGCGCCTGGGCAAAGGCCCGCAGAACCGGTGT‑TTCGAAGAt                                                                                                             <  2:895434/150‑1 (MQ=255)                                      aTGCGTTCTGCGGTTGTCTTGTCCATTTCAAGCCTCAGTCAATTCATTCATCAGTTGTTCGATACCGGCGAGGCCTTCTTCGGTCAG‑GGCCTGGCCCGCGCCTTCCAGCGCCTGGGCAAAGGCCCGCAGAACCGGTGT‑TTCGAAGATgag                                                                                                          >  1:486255/1‑150 (MQ=255)                                       tGCGTTCTGCGGTTGTCTTGTCCATTTCAAGCCTCAGTCAATTCATTCATCAGTTGTTCGATACCGGCGAGGCCTTCTTCGGTCAG‑GGCCTGGCCCGCGCCTTCCAGCGCCTGGGCAAAGGCCCGCAGAACCGGTGt‑tt                                                                                                                    >  2:365264/1‑139 (MQ=255)                                       tGCGTTCTGCGGTTGTCTTGTCCATTTCAAGCCTCAGTCAATTCATTCATCAGTTGTTCGATACCGGCGAGGCCTTCTTCGGTCAG‑GGCCTGGCCCGCGCCTTCCAGCGCCTGGGCAAAGGCCCGCAGAACCGGTGt‑tt                                                                                                                    <  1:365264/139‑1 (MQ=255)                                          gTTCTGCGGTTGTCTTGTCCATTTCAAGCCTCAGTCAATTCATTCATCAGTTGTTCGATACCGGCGAGGCCTTCTTCGGTCAG‑GGCCTGGCCCGCGCCTTCCAGCGCCTGGGCAAAGGCCCGCAGAACCGGTGT‑TTCGAAGATGAGACgg                                                                                                      >  1:776732/1‑150 (MQ=255)                                                        ttgtcCATTTCAAGCCTCAGTCAATTCATTCATCAGTTGTTCGATACCGGCGAGGCCTTCTTCGGTCAG‑GGCCTGGCCCGCGCCTTCCAGCGCCTGGGCAAAGGCCCGCAGAACCGGTGT‑TTCGAAGATGAGACGGATCGGCATATCGAt                                                                                        >  2:368690/1‑150 (MQ=255)                                                           tcCATTTCAAGCCTCAGTCAATTCATTCATCAGTTGTTCGATACCGGCGAGGCCTTCTTCGGTCAG‑GGCCTGGCCCGCGCCTTCCAGCGCCTGGGCAAAGGCCCGCAGAACCGGTGT‑TTCGAAGATGAGACGGATCGGCATATCGATGcc                                                                                     >  2:371150/1‑150 (MQ=255)                                                            ccATTTCATGCCTCAGTCAATTCATTCATCAGTTGTTCGATACCGGCGAGGCCTTCTTCGGTCAG‑GGCCTGGCCCGCGCCTTCCAGCGCCTGGGCAAAGGCCCGCAGAACCGGTGT‑TTCGAAGATGAGACGGATCGGCATATCGATGccc                                                                                    <  2:608825/150‑1 (MQ=255)                                                             cATTTCAAGCCTCAGTCAATTCATTCATCAGTTGTTCGATACCGGCGAGGCCTTCTTCGGTCAG‑GGCCTGGCCCGCGCCTTCCAGCGCCTGGGCAAAGGCCCGCAGAACCGGTGT‑TTCGAAGATGAGACGGATCGGCatat                                                                                            >  2:208412/1‑141 (MQ=255)                                                             cATTTCAAGCCTCAGTCAATTCATTCATCAGTTGTTCGATACCGGCGAGGCCTTCTTCGGTCAG‑GGCCTGGCCCGCGCCTTCCAGCGCCTGGGCAAAGGCCCGCAGAACCGGTGT‑TTCGAAGATGAGACGGATCGGCatat                                                                                            <  1:208412/141‑1 (MQ=255)                                                             cATTTCAAGCCTCAGTCAATTCATTCATCAGTTGTTCGATACCGGCGAGGCCTTCTTCGGTCAG‑GGCCTGGCCCGCGCCTTCCAGCGCCTGGGCAAAGGCCCGCAGAACCGGTGT‑TTCGAAGATGAGACGGATCGGCATATCGATGCCCa                                                                                   <  1:119693/150‑1 (MQ=255)                                                              aTTTCAAGCCTCAGTCAATTCATTCATCAGTTGTTCGATACCGGCGAGGCCTTCTTCGGTCAG‑GGCCTGGCCCGCGCCTTCCAGCGCCTGGGCAAAGGCCCGCAGAACCGGTGT‑TTCGAAGATGAGACGGATCGGCATATCGATGCCCAg                                                                                  >  1:553422/1‑150 (MQ=255)                                                                       cTCAGTCAATTCATTCATCAGTTGTTCGATACCGGCGAGGCCTTCTTCGGTCAG‑GGCCTGGCCCGCGCCTTCCAGCGCCTGGGCAAAGGCCCGCAGAACCGGTGT‑TTCGAAGATGAGACGGATCGGCATATCGATGCCCAGTTGGCTGtt                                                                         <  2:553422/150‑1 (MQ=255)                                                                            tCAATTCATTCATCAGTTGTTCGATACCGGCGAGGCCTTCTTCGGTCAG‑GGCCTGGCCCGCGCCTTCCAGCGCCTGGGCAAAGGCCCGCAGAACCGGTGT‑TTCGAAGATGAGACGGATCGGCATATCGATGCCCAGTTGGCTGTTGATAc                                                                    <  2:1105908/150‑1 (MQ=255)                                                                              aaTTCATTCATCAGTTGTTCGATACCGGCGAGGCCTTCTTCGGTCAG‑GGCCTGGCCCGCGCCTTCCAGCGCCTGGGCAAAGGCCCGCAGAACCGGTGT‑TTCGAAGATGAGACGGATCGGCATATCGATGCCCAGTTGGCTGTTGATACGt                                                                  <  1:223574/150‑1 (MQ=255)                                                                               attcattcATCAGTTGTTCGATACCGGCGAGGCCTTCTTCGGTCAG‑GGCATGGCCCGCGCCTTCCAGCGCCTGGGCAAAGGCCCGCAGAACCGGTGT‑TTCGAAGATGAGACGGATCGGCATATCGATGCCCAGTTGGCTGTTGATACGTg                                                                 >  1:136813/1‑150 (MQ=255)                                                                                   attcatCAGTTGTTCGATACCGGCGAGGCCTTCTTCGGTCAG‑GGCCTGGCCCGCGCCTTCCAGCGCCTGGGCAACTTCCCGCAGCACCGGTGTCTGGGTTGAACTGA‑TAATCATCGTTTGGATGCCCAGTTGGCTGTTGATACGTGagag                                                             >  1:1097712/1‑150 (MQ=255)                                                                                       atcaGTTGTTCGATACCGGCGAGGCCTTCTTCGGTCAG‑GGCCTGGCCCGCGCCTTCCAGCGCCTGGGCAAAGGCCCGCAGAACCGGTGT‑TTCGAAGATGAGACGGATCGGCATATCGATGCCCAGTTgg                                                                              <  2:985287/129‑1 (MQ=255)                                                                                       atcaGTTGTTCGATACCGGCGAGGCCTTCTTCGGTCAG‑GGCCTGGCCCGCGCCTTCCAGCGCCTGGGCAAAGGCCCGCAGAACCGGTGT‑TTCGAAGATGAGACGGATCGGCATATCGATGCCCAGTTgg                                                                              >  1:985287/1‑129 (MQ=255)                                                                                          aGTTGTTCGATACCGGCGAGGCCTTCTTCGGTCAG‑GGCCTGGCCCGCGCCTTCCAGCGCCTGGGCAAAGGCCCGCAGTTTCGGTGT‑TTCGAAGATGAGACGGATCGGCATATCGATGCCCTGTTGGCTGTTGATACGTGAGAGGGCCTgg                                                      >  1:623795/1‑150 (MQ=255)                                                                                           gttgttCGATACCGGCGAGGCCTTCTTCGATAAGAGGTTAGGAACGC‑AGGTACTCTTCCTGGGCAAAGGCCGACTGAACCGGTGT‑TTCGAAGATGAGACGGATCGGCATATCGATGCCCAGTTGGCTGTTGATACGTGAGAGGGCCTggg                                                     >  1:193069/1‑150 (MQ=255)                                                                                               ttCGATACCGGCGAGGCCTTCTTCGGTCAG‑GGCCTGGCCCGCGCCTTCCAGCGCCTGGGCAAAGGCCCGCAGAACCGGTGT‑TTCGaa                                                                                                                <  2:93223/87‑1 (MQ=255)                                                                                               ttCGATACCGGCGAGGCCTTCTTCGGTCAG‑GGCCTGGCCCGCGCCTTCCAGCGCCTGGGCAAAGGCCCGCAGAACCGGTGT‑TTCGaa                                                                                                                >  1:93223/1‑87 (MQ=255)                                                                                               ttCGATACCGGCGAGGCCTTCTTCGGTCAG‑GGCCTGGCCCGCGCCTTCCAGCGCCTGGGCAAAGGCCCGCAGAACCGGTGT‑TTCGAAGATGAGACGGATCGGCATATCGATGCCCa                                                                                   >  2:877055/1‑116 (MQ=255)                                                                                               ttCGATACCGGCGAGGCCTTCTTCGGTCAG‑GGCCTGGCCCGCGCCTTCCAGCGCCTGGGCAAAGGCCCGCAGAACCGGTGT‑TTCGAAGATGAGACGGATCGGCATATCGATGCCCa                                                                                   <  1:877055/116‑1 (MQ=255)                                                                                                        ggCGAGGCCTTCTTCGGTCAG‑GGCCTGGCCCGCGCCTTCCAGCGCCTGGGCAAAGGCCCGCAGAACCGGTGT‑TTCGAAGATGAGACGGATCGGCATATCGATGCCCAGTTGGCTGTTGATACGTGAGAGGGCCTGGGCCGCCAGCAGCGa                                        <  2:486255/150‑1 (MQ=255)                                                                                                             ggCCTTCTTCGGTCAG‑GGCCTGGCCCGCGCCTTCCAGCGCCTGGGCAAAGGCCCGCAGAACCGGTGT‑TTCGAAGATGAGACGGATCGGCATATCGATGCCCAGTTGGCTGTTGATACGTGAGAGGGCCTGGGCCGCCAGCAGCGAGTGac                                   >  2:1093980/1‑150 (MQ=255)                                                                                                              gCCTTCTTCGGTCAG‑CGCCTGGCCCGCGCCTTCCAGCGCCTGCGCAAAGGCCCGCAGAACCGGTGT‑TTCGAAGATGAGACGGATCGGCATATCGATGCCCAGTTGGCTGTTGATACGTGAGAGGGCCTGGGCCGCCAGCAGCGAGTGacc                                  >  2:739152/1‑150 (MQ=255)                                                                                                                          cAG‑GGCCTGGCCCGCGCCTTCCAGCGCCTGGGCAAAGGCCCGCAGAACCGGTGT‑TTCGAAGATGAGACGGATCGGCATATCGATGCCCAGTTGGCTGTTGATACGTGAGAGGGCCTGGGCCGCCAGCAGCGAGTGACCACCCAGTTCaaa                      <  2:112175/150‑1 (MQ=255)                                                                                                                                ccTGGCCCGCGCCTTCCAGCGCCTGGGCAAAGGCCCGCAGAACCGGTGT‑TTCGAAGATGAGACGGATCGGCATATCGATGCCCAGTTGGCTGTTGATACGTGAGAGGGCCTGGGCCGCCAGCAGCGAGTGACCACCCAGTTCAGAGAAGt                 >  1:294574/1‑150 (MQ=255)                                                                                                                                 cTGGCCCGCGCCTTCCAGCGCCTGGGCAAAGGCCCGCAGAACCGGTGT‑TTCGAAGATGAGACGGATCGGCATATCGATGCCCAGTTGGCTGtt                                                                         <  2:474071/93‑1 (MQ=255)                                                                                                                                 cTGGCCCGCGCCTTCCAGCGCCTGGGCAAAGGCCCGCAGAACCGGTGT‑TTCGAAGATGAGACGGATCGGCATATCGATGCCCAGTTGGCTGtt                                                                         >  1:474071/1‑93 (MQ=255)                                                                                                                                 cTGGCCCGCGCCTTCCAGCGCCTGGGCAAAGGCCCGCAGAACCGGTGT‑TTCGAAGATGAGACGGATCGGCATATCGATGCCCAGTTGGCTGTTGATACGTGAGAGGGCCTGGGCCGCCAGCAGCGAGTGACCACCCAGTTCAAAGAAGtt                >  1:539329/1‑150 (MQ=255)                                                                                                                                        gcgcCTTCCGGCGCCTGGGCAAAGGCCCGCAGAACCGGTGT‑TTCGAAGATGAGACGGATCGGCATATCGATGCCCAGTTGGCTGTTGAt                                                                      >  2:199817/1‑89 (MQ=255)                                                                                                                                        gcgcCTTCCGGCGCCTGGGCAAAGGCCCGCAGAACCGGTGT‑TTCGAAGATGAGACGGATCGGCATATCGATGCCCAGTTGGCTGTTGAt                                                                      >  2:198603/1‑89 (MQ=255)                                                                                                                                        gcgcCTTCCGGCGCCTGGGCAAAGGCCCGCAGAACCGGTGT‑TTCGAAGATGAGACGGATCGGCATATCGATGCCCAGTTGGCTGTTGAt                                                                      <  1:199817/89‑1 (MQ=255)                                                                                                                                        gcgcCTTCCGGCGCCTGGGCAAAGGCCCGCAGAACCGGTGT‑TTCGAAGATGAGACGGATCGGCATATCGATGCCCAGTTGGCTGTTGAt                                                                      <  1:198603/89‑1 (MQ=255)                                                                                                                                             ttCCAGCGCCTGGGCAAAGGCCCGCAGAACCGGTGT‑TTCGAAGATGAGACGGATCGGCATATCGATGCCCAGTTGGCTGTTGATACGTGAGAGGGCCTGGGCCGCCAGCAGCGAGTGACCACCCAGTTCAAAGAAGTTGTCGTCCAGccc    <  1:371150/150‑1 (MQ=255)                                                                                                                                              tCCAGCGCCTGGGCAAAGGCCCGCAGAACCGGTGT‑TTCGAAGATGAGACGGATCGGCATATCGATGCCCAGTTGGCTGTTGATACGTGAGAGGGCCTGGGCCGCCAGCAGCGAGTGa                                    >  1:986337/1‑117 (MQ=255)                                                                                                                                              tCCAGCGCCTGGGCAAAGGCCCGCAGAACCGGTGT‑TTCGAAGATGAGACGGATCGGCATATCGATGCCCAGTTGGCTGTTGATACGTGAGAGGGCCTGGGCCGCCAGCAGCGAGTGa                                    <  2:986337/117‑1 (MQ=255)                                                                                                                                               ccAGCGCCTGGGCAAAGGCCCGCAGAACCGGTGT‑TTCGAAGATGAGACGTATCGGCATATCGATGCCCAGTTGGCTGTTGATACGTGAGAGGGCCTGGGCCGCCAGCAGCGAGTGACCACCCAGTTCAAAGAAGTTGTCGTCCAGCCCGa  >  1:528241/1‑150 (MQ=255)                                                                                                                                                 |                                                                                                                                                     GCGCTGCTCCAGCGGCATTGCTACGAAACGTCTGGCGATGCGTTCTGCGGTTGTCTTGTCCATTTCAAGCCTCAGTCAATTCATTCATCAGTTGTTCGATACCGGCGAGGCCTTCTTCGGTCAG‑GGCCTGGCCCGCGCCTTCCAGCGCCTGGGCAAAGGCCCGCAGAACCGGTGT‑TTCGAAGATGAGACGGATCGGCATATCGATGCCCAGTTGGCTGTTGATACGTGAGAGGGCCTGGGCCGCCAGCAGCGAGTGACCACCCAGTTCAAAGAAGTTGTCGTCCAGCCCGA  >  NC_002947/4787296‑4787586

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑