Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,453,992	G→A	W22* (TGG→TGA)	paaA →	phenylacetyl‑CoA 1,2‑epoxidase, monooxygenase subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,453,992	0	G	A	100.0%	91.0 / NA	27	W22* (TGG→TGA)	paaA	phenylacetyl‑CoA 1,2‑epoxidase, monooxygenase subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (14/13);  total (14/13)

GCAACTGCGGAAAACATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAATTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGGATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTGGATCACTCGCGCACCAACCTTGCGGCGTAAAGCCATTCTGTTGGCCAAAGTGCAGGATGAAGCCGGTCATGGTCTGTATCTCTATAGCGCCGCCGAAACACTGGGCTGCGCCAGAGAAGACATCTACCAAAAAATGCTCGACGGACGGATGAA  >  NC_000913/1453746‑1454234                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   gCAACTGCGGAAAACATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAATTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGc                                                                                                                                                                                                                                                >  1:33956/1‑251 (MQ=255)    aCTGCGGAAAACATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAATTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGTCCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCg                                                                                                                                                                                                                                             <  1:449034/251‑1 (MQ=255)         ggAAAACATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAATTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCt                                                                                                                                                                                                                                        >  1:55336/1‑251 (MQ=255)                    aTGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAATTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCACAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGCAAGAc                                                                                                                                                                                                                             >  2:194577/1‑251 (MQ=255)                      gCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAATTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACAt                                                                                                                                                                                                                           <  1:106615/251‑1 (MQ=255)                                    tGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAATTTGTGATTTTACTTAACTATTGTATAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGAtt                                                                                                                                                                                                             >  1:376793/1‑251 (MQ=255)                                                                 tCACAGTTCAAACATACAAAATTTGTGATTTTACTTAACTATTGCGTAACTTTCATAAAACAATGTGCTTCGTGTCTTTAATTAATTCACGCAAACTGGAATCGTAACGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTgg                                                                                                                                                                                >  1:35025/1‑251 (MQ=255)                                                                  cacaGTTCCAACATACAAAATTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGCTTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGACGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGt                                                                                                                                                                               >  1:516017/1‑251 (MQ=255)                                                                               tACAAAATTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCCGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGa                                                                                                                                                                    >  2:415952/1‑249 (MQ=255)                                                                                          tgATTTTACTTAACTCTTGTGTAACTTTCTTAACACAATGTGATTCGTGTTTTTAATTAATTCACGACAACTGGAATCGTAAAGGCGATGACGTGACCCAAGAAGAACGCTTTGAGCACCGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGCCATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGACGGTAACTGGATc                                                                                                                                                       >  2:94438/1‑251 (MQ=255)                                                                                           gATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTGGATCa                                                                                                                                                      <  1:460744/251‑1 (MQ=255)                                                                                                     aaCTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTGGATCACTCGCGCAcc                                                                                                                                            <  2:146574/251‑1 (MQ=255)                                                                                                           tgtgTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTGGATCACTCGCGCACCAACCtt                                                                                                                                      >  2:510842/1‑251 (MQ=255)                                                                                                              gtaccTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTGGTTCACTCGCGCACCAACCTTgcg                                                                                                                                   >  2:419254/5‑251 (MQ=255)                                                                                                               tAACTTTCATAAAACAATGTGATTCGTGTTTTTACTTAATTCACGAACACTGGACTCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACCGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGACGGTAACTgg                                                                                                                                                          >  2:377109/1‑227 (MQ=255)                                                                                                               tAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTgg                                                                                                                                                          <  1:377109/227‑1 (MQ=255)                                                                                                                         aaaaCAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTGGATCACTCGCGCACCAACCTTGCGGCGTAAAGcc                                                                                                                         <  1:36948/250‑1 (MQ=255)                                                                                                                          aaaCAATGTGAGTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTGTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTGGATCACTCGCGCACCAACCTTGCGGCGTAAAGCCAt                                                                                                                       <  2:294003/251‑1 (MQ=255)                                                                                                                                             tttAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGATCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGTCAGCATGCGCACTCCGAGATTGGTGGTATGTTGCCTGAAGGTAAATGGATCACTCGCGCACCAACCTTGCGGCGCAAAGCCATTCTGTTGGCCAAAGTgacg                                                                                                    >  2:415408/1‑248 (MQ=255)                                                                                                                                                             aaaCTGGAATCGTAAAGGCGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGACCAGGAAGCGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGACAGACATTGATCCGCCAGCTTGGGCAGCATGCGCACTCCGCGATTGTTGGTATGTTGCCTGAAGGTAACTGGATCACTCGCGCACCAACCTTGCGGCGTACAGCCATTCTGTTGGCCAAAGTGCAGGATGAAGCCGGTAATg                                                                                    >  2:253381/1‑251 (MQ=255)                                                                                                                                                                                                 agaaCGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTGGATCACTCGCGCACCAACCTTGCGGCGTAAAGCCATTCTGTTGGCCAAAGTGCAGGATGAAGCCGGTCATGGTCTGTATCTCTATAGCGCCGCCGAAACACTGGGCt                                                <  1:419254/251‑1 (MQ=255)                                                                                                                                                                                                   aaCGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTGGATCACTCGCGCACCAACCTTGCGGCGTAAAGCCATTCTGTTGGCCAAAGTGCAGGATGAAGCCGGTCATGGTCTGTATCTCTATAGCGCCGCCGAAACACTGGGCTgc                                              <  1:473509/251‑1 (MQ=255)                                                                                                                                                                                                          tGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTGGATCACTCGCGCACCAACCTTGCGGCGTAAAGCCATTCTGTTGGCCAAAGTGCAGGATGAAGCCGGTCATGGTCTGTATCTCTATAGCGCCGCCGAAACACTGGGCTGCGCCaga                                        <  2:179498/250‑1 (MQ=255)                                                                                                                                                                                                                               aaCGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTGGATCACTCGCGCACCAACCTTGCGGCGTAAAGCCATTCTGTTGGCCAAAGTGCAGGATGAAGCCGGTCATGGTCTGTATCTCTATAGCGCCGCCGAAACACTGGGCTGCGCCAGAGAAGACATCTACCAAAAAATGc                  <  1:372086/251‑1 (MQ=255)                                                                                                                                                                                                                                          aGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCTCACTCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTGGATCACTCGCGCACCAACCTTGCGGCGTAAAGCCATTCTGTTGGCCAAAGTGCAGGATGAAGCCGGTCATGGTCTGTATCTCTATAGCGCCGCCGAAACACTGGGCTGCGCCAGAGAAGACATCTTCCAAAAAATGCTCGacggacgg       <  1:94438/251‑1 (MQ=255)                                                                                                                                                                                                                                               caGGACTGAATGCCCGATTCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTGGATCACTCGCGCACCAACCTTGCGGCGTAAAGCCATTCTGTTGGCCAAAGTGCAGGATGAAGCCGGTCATGGTCTGTATCTCTATAGCGCCGCCGAAACACTGGGCTGCGCCAGAGAAGACATCTACCAAAAAATGCTCGACGGACGGATGaa  <  1:291072/251‑1 (MQ=255)                                                                                                                                                                                                                                               caGGACTGAATGAGCGATGCTCACGGAAAGACGTTGATCCGCCAGACTGGGCAGCATGCGCACTACGAGATTGTTGGTATGTTGCCTGAAGGTAACTGGATCACTCGCGCACCACCATTGCGGCGTAAAGCCATTCTGTTAGCCAAAGCGCAGGATGAAGCCGGTCATGGTCTGTATCTCTATAGCGCAGCCGGACCACTGTGCCACGACAGAGAAAACATATACCAAAAAATGCTCgacggacgactga   >  2:273574/1‑245 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   GCAACTGCGGAAAACATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAATTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGGATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTGGATCACTCGCGCACCAACCTTGCGGCGTAAAGCCATTCTGTTGGCCAAAGTGCAGGATGAAGCCGGTCATGGTCTGTATCTCTATAGCGCCGCCGAAACACTGGGCTGCGCCAGAGAAGACATCTACCAAAAAATGCTCGACGGACGGATGAA  >  NC_000913/1453746‑1454234

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