Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,316,707	T→A	11.3%	V249D (GTC→GAC)	PP_3788 →	non‑ribosomal peptide synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,316,707	0	T	A	11.3%	115.6 / 7.8	45	V249D (GTC→GAC)	PP_3788	non‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base T (18/21);  new base A (2/3);  total (21/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.81e-01

TCACGGCCCTGAAAACCCAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCG‑AGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTGGGGAGTTGAGTTGATCAGTGCGCTTCATCAGC  >  NC_002947/4316560‑4316844                                                                                                                                                     |                                                                                                                                          tCACGGCCCTGAAAACCCAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCG‑AGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCa                                                                                                                                         >  1:997232/1‑150 (MQ=255) tCACGGCCCTGAAAACCCAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCG‑AGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCa                                                                                                                                         <  2:355107/150‑1 (MQ=255)     ggCCCTGAAAACCCAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCG‑AGTCAACTTTCCCAAACACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCAt                                                                                                                                     <  1:473458/150‑1 (MQ=255)       cccTGAAACCCCAGGCCTCTTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCG‑AGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCATGGTGACACCGTGc                                                                                                                                   <  2:616512/150‑1 (MQ=255)         cTGAAAACCCAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCG‑AGTCAACTTTCCCAATTACGTGGCGCCCTGCATCGCTGATCACGGTGTCACCATGCtc                                                                                                                                 >  1:454004/1‑150 (MQ=255)               aCCCAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCG‑AGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATg                                                                                                                           <  1:33260/150‑1 (MQ=255)               aCCCAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCG‑AGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATg                                                                                                                           <  1:148617/150‑1 (MQ=255)                   aGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCG‑AGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTAc                                                                                                                       >  1:408774/1‑150 (MQ=255)                      ccTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCG‑AGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCa                                                                                                                    <  1:1324011/150‑1 (MQ=255)                           ttCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCG‑AGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCAc                                                                                                               <  2:997232/150‑1 (MQ=255)                                 cccTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCG‑AGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTgcag                                                                                                         >  2:428027/1‑150 (MQ=255)                                  ccTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCG‑AGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCATGGTGACACCGTGCtcagc                                                                                                                              <  1:1370203/128‑3 (MQ=255)                                  ccTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCG‑AGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCATGGTGACACCGTGCtcagc                                                                                                                              >  2:1370203/1‑126 (MQ=255)                                            ggaggCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTTGTGCCTGACCACAACTG‑AGTCAGCTTTTCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGTTACCCAGCCCCCTGCAGCATCTTCTCGc                                                                                              >  2:746205/1‑150 (MQ=255)                                               ggCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCG‑AGTTAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAg                                                                                           >  1:263354/1‑150 (MQ=255)                                                         ccTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGAACAAACTCGAACTCGA‑AAACCCAATTACCTGTCGGTCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCAACCTGCAGCATCTTCTCGCCAGCCCGGcaatc                                                                                 >  1:603047/1‑150 (MQ=255)                                                           ttttttttCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCG‑AGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGcaatc                                                                                 >  1:191396/1‑148 (MQ=255)                                                           ttttttttCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCG‑AGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGcaatc                                                                                 <  2:191396/148‑1 (MQ=255)                                                                   ccgccgGCGCCGTGTTCGTGCCTGACCACAACCG‑AGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAAtgcct                                                                          <  2:48441/147‑1 (MQ=255)                                                                   ccgccgGCGCCGTGTTCGTGCCTGACCACAACCG‑AGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAAtgcct                                                                          >  1:48441/1‑147 (MQ=255)                                                                                   cGTGCCTGACCACAACCG‑AGTCAACTTTCCCAATTACCTGTCTGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCatat                                                       <  1:970863/150‑1 (MQ=255)                                                                                    gTGCCTGACCACAACCG‑AGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATc                                                      >  1:438772/1‑150 (MQ=255)                                                                                       ccTGACCACAACCG‑AGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGtt                                                   <  2:740742/150‑1 (MQ=255)                                                                                          gACCACAACCG‑AGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTtgc                                                >  2:1060856/1‑150 (MQ=255)                                                                                             cacaACCG‑AGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCAc                                                                                                               >  2:1378904/1‑84 (MQ=255)                                                                                             cacaACCG‑AGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCAc                                                                                                               <  1:1378904/84‑1 (MQ=255)                                                                                              acaACCG‑AGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCg                                            <  1:428027/150‑1 (MQ=255)                                                                                              acaACCG‑AGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCg                                            <  2:408774/150‑1 (MQ=255)                                                                                                   cG‑AGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGaa                                       <  1:796929/150‑1 (MQ=255)                                                                                                        tCAACTTGCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCCGGCATATCGTTTGCTGCGGAGAACCtt                                   <  2:229406/150‑1 (MQ=255)                                                                                                          aaCTTTCCCAATTACCTGTCGGCCTGCATAGCTGCGCAGAGTGGCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTgg                                 >  1:983624/1‑150 (MQ=255)                                                                                                          aaCTGTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTgg                                 <  2:603047/150‑1 (MQ=255)                                                                                                           aCTTTCCCAATTGCCTGCCGGCCTGCACCGCCGAGCACGGTGTCCCCATGCTCATCATGGTACCCAGCCACCTGCAGCCTCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTggg                                <  1:622231/150‑1 (MQ=255)                                                                                                           aCTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTggg                                <  1:1060856/150‑1 (MQ=255)                                                                                                           aCTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTggg                                >  2:787452/1‑150 (MQ=255)                                                                                                           aCTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTggg                                >  1:532344/1‑150 (MQ=255)                                                                                                             tttCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTGGGga                              >  2:1425098/1‑150 (MQ=255)                                                                                                                   aaTTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGAAATCCAAGCATATCGTTTGCTGCGGAGAACCTTGGGgagttgag                        >  1:1085439/1‑150 (MQ=255)                                                                                                                          tGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTGGGGAGTTGAGTTGATCa                 <  2:983624/150‑1 (MQ=255)                                                                                                                            tCGGCCTGCACCGCTGACAACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTGGGGAGTTGAGTTGATCAGt               <  1:17159/150‑1 (MQ=255)                                                                                                                               gCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCtt                                   >  2:290161/1‑127 (MQ=255)                                                                                                                               gCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCtt                                   <  1:290161/127‑1 (MQ=255)                                                                                                                                     catgatGAGCATGGTGACACCGTGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTgg                                 <  2:585104/118‑1 (MQ=255)                                                                                                                                     catgatGAGCATGGTGACACCGTGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTgg                                 >  1:585104/6‑123 (MQ=255)                                                                                                                                         cTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTGGGGAGTTGAGTTGATCAGTGCGCTTCATCAgc  >  2:22577/1‑150 (MQ=255)                                                                                                                                                     |                                                                                                                                          TCACGGCCCTGAAAACCCAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCG‑AGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTGGGGAGTTGAGTTGATCAGTGCGCTTCATCAGC  >  NC_002947/4316560‑4316844

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