Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,454,292	C→A	G122G (GGC→GGA)	paaA →	ring 1,2‑phenylacetyl‑CoA epoxidase subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,454,292	0	C	A	100.0%	116.4 / NA	37	G122G (GGC→GGA)	paaA	ring 1,2‑phenylacetyl‑CoA epoxidase subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (20/17);  total (20/17)

TCTCTATAGCGCCGCCGAAACACTGGGCTGCGCCAGAGAAGACATCTACCAAAAAATGCTCGACGGACGGATGAAATATTCCTCCATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGCTGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGCCCGTATGCCAGAGCAATGGTGAAAATCTGCAAAGAAGAGAGTTTTCACCAGCGTCAGGGGTT  >  NC_000913/1454160‑1454414                                                                                                                                     |                                                                                                                           tctcTATAGCGCCGCCGAAACACTGGGCTGCGCCAGAGAAGACATCTACCAAAAAATGCTCGACGGACGGATGAAATATTCCTCCATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGgatggctg                                                                                                                      <  2:386120/139‑1 (MQ=255)       taGCGCCGCCGAAACACTGGGCTGCGCCAGAGAAGACGTCTACCAAAAAATGCTCGACGGACGGATGAAATATTCCTCCATCTTTAATTATCCGACATAGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGAt                                                                                                                >  1:758748/1‑139 (MQ=255)       taGCGCCGCCGAAACACTGGGCTGCGCCAAAGAAGACATCTACCAAAAAATGCTCGACGGACGGATGAAATATTTCTCCATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGAt                                                                                                                >  1:758053/1‑139 (MQ=255)         gcgcCGCCGAAACACTGGGCTGCGCCAGAGAAGACATCTACCAAAAAATGCTCGACGGACGGATGAAATATTCCTCCATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATgg                                                                                                              >  1:37699/1‑139 (MQ=255)           gccgccGAAACACTGGGCTGCGCCAGAGAAGACATCTACCAAAAAATGCTCGACGGACGGATGAAATATTCCTCCATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGgcg                                                                                                            <  2:224875/139‑1 (MQ=255)                cgAAACACTGGGCTGCGCCAGAGAAGACATCTACCAAAAAATGCTCGACGGACGGATGAAATATTCCTCCATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTTGTGGATGGCGCAGCg                                                                                                       <  2:758748/139‑1 (MQ=255)                         gggCTGCGCCAGAGAAGACATCTACCAAAAAATGCTCGACGGACGGATGAAATATTCCTCCATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAAt                                                                                              >  2:1071127/1‑139 (MQ=255)                          ggCTGCGCCAGAGAAGACATCTACCAAAAAATGCTCGACGGACGGATGAAATATTCCTCCATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATc                                                                                             >  2:263646/1‑139 (MQ=255)                                        agaCATCTACCAAAAAATGCTCGACGGACGGATGAAATATTCCTCCATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGATGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGc                                                                               <  2:37699/139‑1 (MQ=255)                                                 ccAAAAAATGCTCGACGGACGGATGAAATATTCCTCCATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCAc                                                                          >  1:504252/1‑135 (MQ=255)                                                 ccAAAAAATGCTCGACGGACGGATGAAATATTCCTCCATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCAc                                                                          <  2:504252/135‑1 (MQ=255)                                                         tGCTCGACGGACGGATGAAATATTCCTCCATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGccc                                                              <  1:774637/139‑1 (MQ=255)                                                               acggacggATGAAATATTCCTCCATCTTTAATTATCCGACATTGAGTTGGGTTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGCCCGTATGc                                                        <  2:758053/139‑1 (MQ=255)                                                                                 cctccATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCAc                                                                          <  2:839088/103‑1 (MQ=255)                                                                                 cctccATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCAc                                                                          >  1:839088/1‑103 (MQ=255)                                                                                 cctccATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGCCCGTATGCCAGAGCAATCGTGAAAAt                                      >  2:689402/1‑139 (MQ=255)                                                                                       tCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTa                                                                    >  2:939569/1‑103 (MQ=255)                                                                                       tCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTa                                                                    <  1:939569/103‑1 (MQ=255)                                                                                                 tCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCAc                                                                          >  2:379103/1‑87 (MQ=255)                                                                                                 tCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCAc                                                                          <  1:379103/87‑1 (MQ=255)                                                                                                 tCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCAc                                                                          <  2:904564/87‑1 (MQ=255)                                                                                                 tCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCAc                                                                          >  1:904564/1‑87 (MQ=255)                                                                                                   cGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGCCCGTATGCCAGAGCAATGGTGAAAAt                                      >  2:187508/1‑121 (MQ=255)                                                                                                   cGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGCCCGTATGCCAGAGCAATGGTGAAAAt                                      <  1:187508/121‑1 (MQ=255)                                                                                                           aGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGCCCGTATGCCAGAGCAATGGTGAAAATCTGCAAAGAAGAGAGTTTTCACCAGc            >  2:250769/1‑139 (MQ=255)                                                                                                                ggCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGCCCGTATGCCAGAGCAATg                                              <  2:266013/100‑1 (MQ=255)                                                                                                                ggCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGCCCGTATGCCAGAGCAATg                                              >  1:266013/1‑100 (MQ=255)                                                                                                                     atatCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCActtct                                                                      >  2:1042157/1‑71 (MQ=255)                                                                                                                     atatCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCActtct                                                                      <  1:1042157/71‑1 (MQ=255)                                                                                                                     atatCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGc                                                                >  2:409837/1‑77 (MQ=255)                                                                                                                     atatCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGc                                                                <  1:409837/77‑1 (MQ=255)                                                                                                                     atatCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGCCCGTATGCCAGAGCAATGGTGaaa                                        <  1:538547/101‑1 (MQ=255)                                                                                                                     atatCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGCCCGTATGCCAGAGCAATGGTGaaa                                        >  2:538547/1‑101 (MQ=255)                                                                                                                     atatCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGCCCGTATGCCAGAGCAATGGTGAAAATCTGCAAAGAAGAGAGTTTTCACCAGCGTCAGGGGtt  >  1:631483/1‑139 (MQ=255)                                                                                                                      tatCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTa                                                                                  >  2:1108962/1‑58 (MQ=255)                                                                                                                      tatCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTa                                                                                  <  1:1108962/58‑1 (MQ=255)                                                                                                                          ggTGTTATCGGATGGCTGGTGGATGGCGCAGCGAt                                                                                                     >  2:425302/1‑35 (MQ=255)                                                                                                                          ggTGTTATCGGATGGCTGGTGGATGGCGCAGCGAt                                                                                                     <  1:425302/35‑1 (MQ=37)                                                                                                                                     |                                                                                                                           TCTCTATAGCGCCGCCGAAACACTGGGCTGCGCCAGAGAAGACATCTACCAAAAAATGCTCGACGGACGGATGAAATATTCCTCCATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGCTGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGCCCGTATGCCAGAGCAATGGTGAAAATCTGCAAAGAAGAGAGTTTTCACCAGCGTCAGGGGTT  >  NC_000913/1454160‑1454414

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