Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,248,515	C→T	8.9%	V504V (GTG→GTA)	PP_3724 ←	acyl‑CoA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,248,515	0	C	T	8.9%	124.0 / 2.1	45	V504V (GTG→GTA)	PP_3724	acyl‑CoA synthetase
Reads supporting (aligned to +/- strand):  ref base C (26/15);  new base T (2/2);  total (28/17)
Fisher's exact test for biased strand distribution p-value = 6.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

AGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTGGTCCATCAACCGCTGGCGGTTCTTTT  >  NC_002947/4248371‑4248663                                                                                                                                                 |                                                                                                                                                     aGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGTACGc                                                                                                                                                  >  1:263687/1‑149 (MQ=255)     tCCAGGGCGCGGGCATCGCGCAGGTTCAAAGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACAGCCTTCAGTACAAAGGGGTAACCGACGAGGGCGGCGTCTGCAAAGGCGCGTTCGTGCGTACTCATCACGccag                                                                                                                                               >  3:290155/1‑146 (MQ=255)       cAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGTACGCCCGCa                                                                                                                                             <  2:263687/148‑1 (MQ=255)       cAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAg                                                                                                                                            <  6:312160/149‑1 (MQ=255)                        gcAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAg                                                                                                                           >  8:280479/1‑149 (MQ=255)                                        aGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCt                                                                                                           <  2:317499/149‑1 (MQ=255)                                             cAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGAt                                                                                                      >  8:140806/1‑149 (MQ=255)                                                   cGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGcc                                                                                                  <  8:110924/147‑1 (MQ=255)                                                   cGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCa                                                                                                 <  8:205660/148‑1 (MQ=255)                                                      tCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCg                                                                                             >  7:259930/1‑149 (MQ=255)                                                          gtgCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTcc                                                                                          >  6:100772/1‑148 (MQ=255)                                                             cGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCCTTCGCGCGTACTCACCACGCCCGCAGGCACCGGCAGACCGAACGCGCGCACGGCCTGCATGGATTGCCAATCGTCCAACa                                                                                      >  8:118041/1‑149 (MQ=255)                                                                         gcTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGc                                                                          <  2:245526/149‑1 (MQ=255)                                                                             aCCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGt                                                                      >  2:94173/1‑149 (MQ=255)                                                                             aCCGCCTTCAGTACCAAGGGGTAACCGACGCGGACGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGt                                                                      >  5:148088/1‑149 (MQ=255)                                                                              ccGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAATCGGCTTTCGCCCGt                                                                      >  5:33075/1‑148 (MQ=255)                                                                               cGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGt                                                                      >  4:222699/1‑147 (MQ=255)                                                                                 ccTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAgg                                                                  >  8:104726/1‑149 (MQ=255)                                                                                        tACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCg                                                           >  5:286629/1‑149 (MQ=255)                                                                                          ccAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAATCGGCTTTCGCCCGTGAGGGCCTGCgg                                                          <  6:33075/148‑1 (MQ=255)                                                                                           cAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGACGGCGCGCAGTGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCggg                                                         >  8:80829/1‑148 (MQ=255)                                                                                           cAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCggg                                                         >  2:316143/1‑148 (MQ=255)                                                                                            aaGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGca                                                       >  2:22615/1‑149 (MQ=255)                                                                                            aaGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGca                                                       >  7:14777/1‑149 (MQ=255)                                                                                             aGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGcac                                                      >  1:309853/1‑149 (MQ=255)                                                                                                 gTAACCGACGCGGGCGGCGTCTTCCACGGCGCGTTCGGGCGTTCTCAGCACGCCCGCAGGCACCGGCCTACCTAGGGCGCGCAGTGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCTCCCGTGAGGGCCTTCGGGCACAGCg                                                  <  4:290155/149‑1 (MQ=255)                                                                                                        aCGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCg                                           >  3:21528/1‑149 (MQ=255)                                                                                                            gggcggCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGca                                         >  7:157203/1‑147 (MQ=255)                                                                                                                     cTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGccc                              >  3:208901/1‑149 (MQ=255)                                                                                                                     cTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGccc                              <  7:104726/149‑1 (MQ=255)                                                                                                                     cTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGccc                              <  7:45861/149‑1 (MQ=255)                                                                                                                        ccACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCg                                                  >  6:304090/1‑126 (MQ=255)                                                                                                                        ccACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCg                                                  <  5:304090/126‑1 (MQ=255)                                                                                                                        ccACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTgg                           >  4:216768/1‑149 (MQ=255)                                                                                                                           cGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTGGTcc                        <  8:141183/149‑1 (MQ=255)                                                                                                                                  ttCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTTGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTGGTCCATCAAcc                 <  1:94173/149‑1 (MQ=255)                                                                                                                                  ttCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTGGTCCATCAAcc                 >  1:317799/1‑149 (MQ=255)                                                                                                                                  tgcGGGCGTGCTGAGTACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGcaca                                       >  7:272175/3‑127 (MQ=255)                                                                                                                                  tgcGGGCGTGCTGAGTACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGcaca                                       <  8:272175/125‑1 (MQ=255)                                                                                                                                          tACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTGGTCCATCAACCGCTGGCgg         <  7:280479/149‑1 (MQ=255)                                                                                                                                              cAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTGGTCCATCAACCGCTGGCg          >  2:110995/1‑144 (MQ=255)                                                                                                                                              cAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTGGTCCATCAACCGCTGGCg          <  1:110995/144‑1 (MQ=255)                                                                                                                                              cAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTGGTCCATCAACCGCTGGCGGTTCt     >  7:301291/1‑149 (MQ=255)                                                                                                                                              cAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCCCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTGGTCCATCAACCGCTGGCGGTTCt     >  8:280134/1‑149 (MQ=255)                                                                                                                                                gCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTGGTCCATCAACCGCTGGCGGTTCttt   <  3:216768/149‑1 (MQ=255)                                                                                                                                                 cACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACATCGGTACCCGCACAGCCTCGCCCTGGTCCATCAACCGCTGGCGGTTCtttt  >  4:68637/1‑149 (MQ=255)                                                                                                                                                 cACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTGGTCCATCAACCGCTGGCGGTTCtttt  >  2:298166/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                     AGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTGGTCCATCAACCGCTGGCGGTTCTTTT  >  NC_002947/4248371‑4248663

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