Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,718,098	G→T	16.6%	intergenic (‑111/+88)	paaA ← / ← paaK	ring 1,2‑phenylacetyl‑CoA epoxidase subunit alpha/phenylacetate‑CoA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,718,098	0	G	T	16.6%	79.5 / 2.6	29	intergenic (‑111/+88)	paaA/paaK	ring 1,2‑phenylacetyl‑CoA epoxidase subunit alpha/phenylacetate‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base G (16/8);  new base T (2/3);  total (18/11)
Fisher's exact test for biased strand distribution p-value = 3.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

TTTCCACTAGCTGTGCGTACATGTGCTGCGACTCCTGCTAAGGCATGAGTCATTTGTAAGCGATACAAAACATAACGTCAAACACATTTATGTGTATCAAGTAAGATTTGTGTATCACGTTTACAGCCAAGAGGGCTGTCGGGCTGCAGCGTATGCCCGGCAGAACGCAGGCAAAAAAAAGCCCCGGACAATGCCGGGGCCGAAATGGCAGCCGAATGGGTCAGCTGGCCAACCGTCTGTCGAACACGTGGCACGCCTTGCCCTCGG  >  NC_002947/3717966‑3718232                                                                                                                                     |                                                                                                                                       tttCCACTAGCTGTGCGTACATGTGCTGCGACTCCTGCTAAGGCATGAGTCATTTGTAAGCGATACAAAACATAACGTCAAACACATTTATGTGTATCAAGTAAGATTTGTGTATCACGTTTACAGCCAAGAGGGc                                                                                                                                     <  1:73266/136‑1 (MQ=255)         aGCTGTGCGTACATGTGCTGCGACTCCTGCTAAGGCATGAGTCATTTGTAAGCGATACAAAACATAACGTCAAACACATTTATGTGTATCAAGTATGATTTGTGTATCACGTTTACAGACAAGAGTGCTGTCGGGc                                                                                                                             >  2:61177/1‑136 (MQ=255)         aGCTGTGCGTACATGTGCTGCGACTCCTGCTAAGGCATGAGTCATTTGTAAGCGATACAAAACATAACGTCAAACACATTTATGTGTATCAAGTAAGATTTGTGTATCACGTTTACAGCCAAGAGGGCTGTCGGGc                                                                                                                             >  2:324558/1‑136 (MQ=255)         aGCTGTGCGTACATGTGCTGCGACTCCTGCTAAGGCATGAGTCATTTGTAAGCGATACAAAACATAACGTCAAACACATTTATGTGTATCAAGTAAGATTTGTGTATCACGTTTACAGCCAAGAGGGCTGTCGGGc                                                                                                                             >  2:59968/1‑136 (MQ=255)                  tACATGTGCTGCGACTCCTGCTAAGGCATGAGTCATTTGTAAGCGATACAAAACATAACGTCAAACACATTTATGTGTATCAAGTAAGATTTGTGTATCACGTTTACAGCCAAGAGGGCTGTCGGGCTGCAGCGTa                                                                                                                    >  2:97365/1‑136 (MQ=255)                   aCATGTGCTGCGACTCCTGCTAAGGCATGAGTCATTTGTAAGCGATACAAAACATAACGTCAAACACATTTATGTGTATCAAGTAAGATTTGTGTATCACGTTTACAGCCAAGAGGGCTGTCGGGCTGCAGCGTAt                                                                                                                   >  2:20490/1‑136 (MQ=255)                             cGACTCCTGCTAAGGCATGAGTCATTTGTAAGCGATACAAAACATAACGTCAAACACATTTATGTGTATCAAGTAAGATTTGTGTATCACGTTTACAGCCAACAGGGCTGTCGGGCTGCAGCGTATGCCCGGCAGa                                                                                                         >  2:219527/1‑136 (MQ=255)                                       tAAGGCATGAGTCATTTGTAAGCGATACAAAACATAACGTCAAACACATTTATGTGTATCAAGTAAGATTTGTGTATCACGTTTACAGCCAAGAGGGCTGTCGGGCTGCAGCGTATGCCCGGCAGAACGCAGGCaa                                                                                               >  2:189814/1‑136 (MQ=255)                                                        gTAAGCGATACAAAACATAACGTCAAACACATTTATGTGTATCAAGTAAGATTTGTGTATCACGTTTACAGCCAAGAGGGCTGTCGGGCTGCAGCGTATGCCCGGCAGAACGCAGGCAAAAAAAAGCCCCGGACaa                                                                              >  2:361538/1‑136 (MQ=255)                                                                   aaaaCATAACGTCAAACACATTTATGTGTATCAAGTAAGATTTGTGTATCACGTTTACAGCCAAGATGGCTGTCGGGCTGCAGCGTATGCCCGGCAGAACGCAGGCAAAAAAAAGCCCCGGACAATGCCGGGGCCg                                                                   >  2:55705/1‑136 (MQ=255)                                                                   aaaaCATAACGTCAAACACATTTATGTGTATCAAGTAAGATTTGTGTATCACGTTTACAGCCAAGATGGCTGTCGGGCTGCAGCGTATGCCCGGCAGAACGCAGGCAAAAAAAAGCCCCGGACAATGCCGGGGCCg                                                                   >  2:347815/1‑136 (MQ=255)                                                                               cAAACACATTTATGTGTATCAAGTAAGATTTGTGTATCACGTTTACAGCCAAGAGGGCTGTCGGGCTGCAGCGTATGCCCGGCAGAACGCAGGCAAAAAAAAGCCCCGGACAATGCCGGGGCCGAAATGGCAGCCg                                                       >  2:130655/1‑136 (MQ=255)                                                                                  acacATTTATGTGTATCAAGTAAGATTTGTGTATCACGTTTACAGCCAAGAGGGCTGTCGGGCTGCAGCGTATGCCCGGCAGAACGCAGGCAAAAAAAAGCCCCGGACAATGCCGGGGCCGAAATGGCAGCCGAAt                                                    <  1:324558/136‑1 (MQ=255)                                                                                  acacATTTATGTGTATCAAGTAAGATTTGTGTATCACGTTTACAGCCAAGAGGGCTGTCGGGCTGCAGCGTATGCCCGGCAGAACGCAGGCAAAAAAAAGCCCCGGACAATGCCGGGGCCGAAATGGCAGCCGAAt                                                    <  1:59968/136‑1 (MQ=255)                                                                                  acacATTTATGTGTATCAAGTAAGATTTGTGTATCACGTTTACAGCCAAGAGGGCTGTCGGGCTGCAGCGTATGCCCGGCAGAACGCAGGCAAAAAAAAGCCCCGGACAATGCCGGGGCCGAAATGGCAGCCGAAt                                                    <  1:61177/136‑1 (MQ=255)                                                                                    acaTTTATGTGTATCAAGTAAGATTTGTGTATCACGTTTACAGCCAAGAGGGCTGTCGGGCTGCAGCGTATGCCCGGCAGAACGCAGGCAAAAAAAAGCCCCGGACAATGCCGGGGCCGAAATGGCAGCCGAATgg                                                  >  2:497264/1‑136 (MQ=255)                                                                                          aTGTGTATCAAGTAAGATTTGTGTATCACGTTTACAGCCAAGAGGGCTGTCGGGCTGCAGCGTATGCCCGGCAGAACGCAGGCAAAAAAAAGCCCCGGACAATGCCGGGGCCGAAATGGCAGCCGAATGGGTCAGc                                            >  1:173692/1‑136 (MQ=255)                                                                                                 tCAAGTAAGATTTGTGTATCACGTTTACAGCCAAGAGGGCTGTCGGGCTGCAGCGTATGCCCGGCAGAACGCAGGCAAAAAAAAGCCCCGGACAATGCCGGGGCCGAAATGGCAGCCGAATGGGTCAGCTGGCCaa                                     >  1:365923/1‑136 (MQ=255)                                                                                                   aaGTAAGATTTGTGTATCACGTTTACAGCCAAGAGGGCTGTCGGGCTGCAGCGTATGCCCGGCAGAACGCAGGCAAAAAAAAGCCCCGGACAATGCCGGGGCCGAAATGGCAGCCGAATGGGTCAGCTGGCCAAcc                                   >  1:95881/1‑136 (MQ=255)                                                                                                                 tATCACGTTTACAGCCAAGATGGCTGTCGGGCTGCAGCGTATGCCCGGCAGAACGCAGGCAAAAAAAAGCCCCGGACAATGCCGGGGCCGAAATGGCAGCCGAATGGGTCAGCTGGCCAACCGTCTGTCGAACACg                     <  1:55705/136‑1 (MQ=255)                                                                                                                 tATCACGTTTACAGCCAAGATGGCTGTCGGGCTGCAGCGTATGCCCGGCAGAACGCAGGCAAAAAAAAGCCCCGGACAATGCCGGGGCCGAAATGGCAGCCGAATGGGTCAGCTGGCCAACCGTCTGTCGAACACg                     <  1:347815/136‑1 (MQ=255)                                                                                                                  aTCACGTTTACAGCCAAGAGGGCTGTCGGGCTGCAGCGTATGCCCGGCAGAACGCAGGCAAAAAAAAGCCCCGGACAATGCCGGGGCCGAAATGGCAGCCGAATGGGTCAGCTGGCCAACCGTCTGTCGAACACGt                    >  1:340107/1‑136 (MQ=255)                                                                                                                   tCACGTTTACAGCCAAGAGGGCTGTCGGGCTGCAGCGTATGCCCGGCAGAACGCAGGCAAAAAAAAGCCCCGGACAATGCCGGGGCCGAAATGGCAGCCGAATGGGTCAGCTGGCCAACCGTCTGTCGAACACGTg                   <  1:189814/136‑1 (MQ=255)                                                                                                                   tCACGTTTACAGCCAAGAGGGCTGTCGGGCTGCAGCGTATGCCCGGCAGAACGCAGGCAAAAAAAAGCCCCGGACAATGCCGGGGCCGAAATGGCAGCCGAATGGGTCAGCTGGCCAACCGTCTGTCGAACACGTg                   <  1:20490/136‑1 (MQ=255)                                                                                                                       gTTTACAGCCAAGAGGGCTGTCGGGCTGCAGCGTATGCCCGGCAGAACGCAGGCAAAAAAAAGCCCCGGACAATGCCGGGGCCGAAATGGCAGCCGAATGGGTGAGCTGGCCAACCGTCTGTCGAACACGTGGCAc               >  1:108484/1‑136 (MQ=255)                                                                                                                       gTTTACAGCCAAGAGGGCTGTCGGGCTGCAGCGTATGCCCGGCAGAACGCAGGCAAAAAAAAGCCCCGGACAATGCCGGGGCCGAAATGGCAGCCGAATGGGTCAGCTGGCCAACCGTCTGTCGAACACGTGGCAc               >  2:3982/1‑136 (MQ=255)                                                                                                                        tttACATCCAAGATGGCTGGCGGGCTGCAGCGTATGCCCGGCAGAACGCAGGCAAAAAAAAGCCCCGGACAATGCCGGTGCCGAAATGGCAGCCGAATGGGTCAGCTGGCCAACCGTCTGTCGAACACGTGGCACg              <  2:185042/136‑1 (MQ=255)                                                                                                                                cAAGAGGGCTGTCGGGCTGCAGCGTATGCCCGGCAGAACGCAGGCAAAAAAAAGCCCCGGACAATGCCGGGGCCGAAATGGCAGCCGAATGGGTGAGCTGGCCAACCGTCTGTCGAACACGTGGCACGCCTTGccc      <  2:108484/136‑1 (MQ=255)                                                                                                                                aaaGAGGGCTGTCGGGCTGCAGCGTATGCCCGGCAGAACGCAGGCAAAAAAAAGCCCCGGACAATGCCGGGGCCGAAATGGCAGCCGAATGGGTCAGCTGGCCAACCGTCTGTCGAACACGTGGCACGCCTTGccc      <  1:221831/135‑1 (MQ=255)                                                                                                                                  agagGGCTGTCGGGCTGCAGCGTATGCCCGGCAGAACGCAGGCAAAAAAAAGCCCCGGACAATGCCGGGGCCGAAATGGCAGCCGAATGGGTCAGCTGGCCAACCGTCTGTCGAACACGTGGCACGCCTTGCCCTc    >  2:303796/1‑136 (MQ=255)                                                                                                                                  agagGGCTGTCGGGCTGCAGCGTATGCCCGGCAGAACGCAGGCAAAAAAAAGCCCCGGACAATGCCGGGGCCGAAATGGCAGCCGAATGGGTCAGCTGGCAAACCGTCTGTCGAACACGTGGCACGCCTTGCCCTc    <  1:130655/136‑1 (MQ=255)                                                                                                                                    agGGCTGTCGGGCTGCAGCGTATGCCCGGCAGAACGCAGGCAAAAAAAAGCCCCGGACAATGCCGGGGCCGAAATGGCAGCCGAATGGGTCAGCTGGCCAACCGTCTGTCGAACACGTGGCACGCCTTGCCCTCgg  >  1:146668/1‑136 (MQ=255)                                                                                                                                     |                                                                                                                                       TTTCCACTAGCTGTGCGTACATGTGCTGCGACTCCTGCTAAGGCATGAGTCATTTGTAAGCGATACAAAACATAACGTCAAACACATTTATGTGTATCAAGTAAGATTTGTGTATCACGTTTACAGCCAAGAGGGCTGTCGGGCTGCAGCGTATGCCCGGCAGAACGCAGGCAAAAAAAAGCCCCGGACAATGCCGGGGCCGAAATGGCAGCCGAATGGGTCAGCTGGCCAACCGTCTGTCGAACACGTGGCACGCCTTGCCCTCGG  >  NC_002947/3717966‑3718232

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