Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,455,774	A→G	N204S (AAC→AGC)	paaC →	phenylacetyl‑CoA 1,2‑epoxidase, structural subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,455,774	0	A	G	100.0%	97.6 / NA	27	N204S (AAC→AGC)	paaC	phenylacetyl‑CoA 1,2‑epoxidase, structural subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (13/14);  total (13/14)

TTGTGGCGTTTTACCGCCGAACTGTTCGATGCCGACGAGATTGATATTGCACTGAGTGAAGAGGGTATTGCGGTTGATCCACGCACTTTACGCGCAGCGTGGGAAGCCGAAGTTTTTGCCGGGATCAACGAAGCCACATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAAGGACTGCATACCGAACACCTTGGACCCATGCTGGCAGAAATGCAGTATCTCCAGCGTGTCTTGCCCGGTCAGCAAT  >  NC_000913/1455647‑1455905                                                                                                                                |                                                                                                                                    ttGTGGCGTTTTACCGCCGAACTGTTCGATGCCGACGAGATTGATATTGCACTGAGTGAAGAGGGTATTGCGGTTGATCCACGCACTTTACGCGCAGCGTGGGAAGCCGAAGTTTTTGCCGGGATCAGCGAAGCCACAtt                                                                                                                         <  2:15943/140‑1 (MQ=255)    tgGCGTTTTACCGCCGAACTGTTCGATGCCGACGAGATTGATATTGCACTGAGTGAAGAGGGTATTGCGGTTGATCCACGCACTTTACGCGCAGCGTGGGAAGCCGAAGTTTTTGCCGGGATCAGCGAAGCCACATTGaa                                                                                                                      >  2:336242/1‑140 (MQ=255)           ttACCGCCGAACTGTTCGATGCCGACGAGATTGATATTGCACTGAGTGAAGAGGGTATTGCGGTTGATCCACGCACTTTACGCGCAGCGTGGGAAGCCGAAGTTTTTGCCGGGATCAGCGAAGCCACATTGAACGTACCg                                                                                                               >  1:246489/1‑140 (MQ=255)                gccgAACTGTTCGATGCCGACGAGATTGATATTGCACTGAGTGAAGAGGGTATTGCGGTTGATCCACGCACTTTACGCGCAGCGTGGGAAGCCGAAGTTTTTGCCGGGATCAGCGAAGCCACATTGAACGTACCGCAaga                                                                                                          <  2:53488/140‑1 (MQ=255)                    aaCTGTTCGATGCCGACGAGATTGATATTGCACTGAGTGAAGAGGGTATTGCGGTTGATCCACGCACTTTACGCGCAGCGTGGGAAGCCGAAGTTTTTGTCGGGATCAGCGAAGCCACATTGAACGTACCGCAAGAGCAg                                                                                                      <  1:189530/140‑1 (MQ=255)                    aaCTGTTCGATGCCGACGAGATTGATATTGCACTGAGTGAAGAGGGTATTGCGGTTGATCCACGCACTTTACGCGCAGCGTGGGAAGCCGAAGTTTTTGCCGGGATCAGCGAAGCCACATTGAACGTACCGCAAGAGCAg                                                                                                      <  1:283199/140‑1 (MQ=255)                                     gagaTTGATATTGCACTGAGTGAAGAGGGTATTGCGGTTGATCCACGCACTTTACGCGCAGCGTGGGAAGCCGAAGTTTTTGCCGGGATCAGCGAAGCCACATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCgg                                                                                     <  1:224725/140‑1 (MQ=255)                                                       aGTGAAGAGGGTATTGCGGTTGATCCACGCACTTTACGCGCAGCGTGGGAAGCCGAAGTTTTTGCCGGGATCAGCGAAGCCACATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAAGGACTGCATAc                                                                   >  2:320036/1‑140 (MQ=255)                                                       aGTGAAGAGGGTATTGCGGTTGATCCACGCACTTTACGCGCAGCGTGGGAAGCCGAAGTTTTTGCCGGGATCAGCGAAGCCACATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAAGGACTGCATAc                                                                   >  2:532703/1‑140 (MQ=255)                                                                ggTATTGCGGTTGATCCACGCACTTTACGCGCAGCGTGGGAAGCCGAAGTTTTTGCCGGGATCAGCGAAGCCACATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAAGGACTGCATACCGAACACCt                                                          >  1:525959/1‑140 (MQ=255)                                                                  tATTGCGGTTGATCCACGCACTTTACGCGCAGCGTGGGAAGCCGAAGTTTTTGCCGGGATCAGCGAAGCCACATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAAGGACTGCATACCGAACACCTTg                                                        >  1:279357/1‑140 (MQ=255)                                                                      gCGGTTGATCCACGCACTTTACGCGCAGCGTGGGAAGCCGAAGTTTTTGCCGGGATCAGCGAAGCCACATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAAGGACTGCATACCGAACACCTTGGAcc                                                    <  2:279357/140‑1 (MQ=255)                                                                       cGGTTGATCCACGCACTTTACGCGCAGCGTGGGAAGCCGAAGTTTTTGCCGGGATCAGCGAAGCCACATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAAGGACTGCATACCGAACACCTTGGAccc                                                   <  1:336242/140‑1 (MQ=255)                                                                                       tttACGCGCAGCGTGGGAAGCCGAAGTTTTTGCCGGGATCAGCGAAGCCACATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAAGGACTGCATACCGAACACCTTGGACCCATGCTGGCAGAAATGc                                   <  1:320036/140‑1 (MQ=255)                                                                                         tacacgcAGCGTGGGAAGCCGAAGTTTTTGCCGGGATCAGCGAAGCCACATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAAGGACTGCATACCGAACACCTTGGACCCATGCTGGCAGAAATGCAg                                 >  1:382689/5‑140 (MQ=255)                                                                                         tACGCGCAGCGTGGGAAGCCGAAGTTTTTGCCGGGATCAGCGAAGCCACATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAAGGACTGCATACCGAACACCTTGGACCCATGCTGGCAGAAATGCAg                                 >  2:478843/1‑140 (MQ=255)                                                                                             cgcAGCGTGGGAGGCCGAAGTTTTTGCCGGGATCAGCGAAGCCACATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAAGGACTGCATACCGAACACCTTTGACCCATGCTAGCAGAAATGCacgatt                             >  2:579268/1‑135 (MQ=255)                                                                                             cgcAGCGTGGGAAGCCGAAGTTTTTGCCGGGATCAGCGAAGCCACATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAAGGACTGCATACCGAACACCTTGGACCCATGCTGGCAGAAATGCAGTAtc                             >  2:308656/1‑140 (MQ=255)                                                                                                    tGGGAAGCCGAAGTTTTTGCCGGGATCAGCGAAGCCACATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAAGGACTGCATACCGAACACCTTGGACCCATGCTGGCAGAAATGCAGTATCTCCAGCg                      >  1:49936/1‑140 (MQ=255)                                                                                                            cGAAGTTTTTGCCGGGATCAGCGAAGCCACATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAAGGACTGCATACCGAACACCTTGGACCCATGCTGGCAGAAATGCAGTATCTCCAGCGTGTCTTGc              >  2:333864/1‑140 (MQ=255)                                                                                                             gAAGTTTTTGCCGGGATCAGCGAAGCCACATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAAGGACTGCATACCGAACACCTTGGACCCATGCTGGCAGAAATGCAGTATCTCCAGCGTGTCTTGcc             <  2:525959/140‑1 (MQ=255)                                                                                                              aaGTTTTTGCCGGGATCAGCGAAGCCACATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAAGGACTGCATACCGAACACCTTGGACCCATGCTGGCAGAAATGCAGTATCTCCAGCGTGTCTTGccc            <  2:147338/140‑1 (MQ=255)                                                                                                               aGTTTTTGCCGGGATCAGCGAAGCCACATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAAGGACTGCATACCGAACACCTTGGACCCATGCTGGCAGAAATGCAGTATCTCCTGCGTGTCTTGCCCg           >  2:588046/1‑140 (MQ=255)                                                                                                                   tttGCCGGGATCAGCGAAGCCACATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAAGGACTGCATACCGAACACCTTGGACCCATGCTGGCAGAAATGCAGTATCTCCAGCGTGTCTTGCCCGGTCa       <  1:532703/140‑1 (MQ=255)                                                                                                                     tGCCGGGATCAGCGAAGCCACATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAAGGACTGCATACCGAACACCTTGGACCCATGCTGGCAGAAATGCAGTATCTCCAGCGTGTCTTGCCCGGTCAGc     >  2:562807/1‑140 (MQ=255)                                                                                                                      gCCGGGATCAGCGAAGCCACATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAAGGACTGCATACCGAACACCTTGGACCCATGCTGGCAGAAATGCAGTATCTCCAGCGTGTCTTGCCCGGTCAGCa    <  2:382689/140‑1 (MQ=255)                                                                                                                       ccGGGATCAGCGAAGCCACATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAAGGACTGCATACCGAACACCTTGGACCCATGCTGGCAGAAATGCAGTATCTCCAGCGTGTCTTGCCCGGTCAGCaa   <  1:478843/140‑1 (MQ=255)                                                                                                                        cGGGATCAGCGAAGCCACATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAAGGACTGCATACCGAACACCTTGGACCCATGCTGGCAGAAATGCAGTATCTCCAGCGTGTCTTGCCCGGTCAGCAAt  <  2:11869/140‑1 (MQ=255)                                                                                                                                |                                                                                                                                    TTGTGGCGTTTTACCGCCGAACTGTTCGATGCCGACGAGATTGATATTGCACTGAGTGAAGAGGGTATTGCGGTTGATCCACGCACTTTACGCGCAGCGTGGGAAGCCGAAGTTTTTGCCGGGATCAACGAAGCCACATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAAGGACTGCATACCGAACACCTTGGACCCATGCTGGCAGAAATGCAGTATCTCCAGCGTGTCTTGCCCGGTCAGCAAT  >  NC_000913/1455647‑1455905

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