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%	112.6 / NA	34	W22* (TGG→TGA)	paaA	phenylacetyl‑CoA 1,2‑epoxidase, monooxygenase subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (21/13);  total (21/13)

CTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAATTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGGATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTGGATCACTCGCGCACCAACCTTGCGGCGTAAAGCCATTCTGTTGGCCAAAGTGCAGGATGAAGCCGGTCATGGTCTGTATCTCTATAGCGCCGCCGAAACACTGGGCTGCGCCAGAGAAGACATCTACCAAAAAATGCTCGACGGACGGATG  >  NC_000913/1453770‑1454232                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                 cTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAATTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTg                                                                                                                                                                                                                       >  1:373875/1‑250 (MQ=255)                                      ccGTCACAGTTCAAACATACAAAATTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGAttgt                                                                                                                                                                                 >  1:470941/1‑251 (MQ=255)                                       cGTCACAGTTCAAACATACAAAATTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTcc                                                                                                                                                                                         >  2:453156/1‑242 (MQ=255)                                       cGTCACAGTTCAAACATACAAAATTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTcc                                                                                                                                                                                         <  1:453156/242‑1 (MQ=255)                                              gTTCAAACATACAAAATTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATGCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGt                                                                                                                                                                         <  2:368283/251‑1 (MQ=255)                                                                    aTTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGCGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTGGATCAc                                                                                                                                                   >  1:522725/1‑251 (MQ=255)                                                                       ttACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTGGATCACTc                                                                                                                                                 <  1:456570/250‑1 (MQ=255)                                                                                        aaCTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGCATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCGTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCGCAGGACTGAATGCCCGATGCTTACCGAAAGACATGGATCCGCCAGATTGGGCTGCATGCGCACTCCGAGATTGTTGGTATGTGGCCTGAAGGGAACTGGATCACTCGCGCGCCAACCTTgcggc                                                                                                                               <  2:236622/251‑1 (MQ=255)                                                                                         aCTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTGGATCACTCGCGCACCAACCTTgcggcg                                                                                                                              <  1:127391/251‑1 (MQ=255)                                                                                                  aaaCAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTGGATCACTCGCGCACCAACCTTGCGGCGTAAAGCCAt                                                                                                                     <  2:137641/251‑1 (MQ=255)                                                                                                         gtgATTCGTGGTTTTAATTAATTCACGAAAACTGGCATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGCAACTGGATCACTCGCGCACCAACCTTGCGGCGTAAAGCCATTCTGTTg                                                                                                              >  1:536681/1‑251 (MQ=255)                                                                                                          tgatgcgggtTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGCCGTTACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCAATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGAAATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTGGATCACTCGCGCACCAACCTTGCGGCGTAAAGCCATAATGTTgg                                                                                                             >  1:257691/9‑251 (MQ=255)                                                                                                                gtgtTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGCTAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTGGATCACTCGCGCACCAACCTTGCGGCGTAAAGCCATTCTGTTGGCCac                                                                                                         >  2:557548/1‑248 (MQ=255)                                                                                                                         attaattCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTGGATCACTCGCGCACCAACCTTGCGGCGTAAAGCCATTCTGTTGGCCAAAGTGCAGGATg                                                                                              <  1:557548/251‑1 (MQ=255)                                                                                                                             attCACGAAAACTGGAAGCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTGGATCACTCGCGCACCAACCTTGCGGCGTAAAGCCATTCTGTTGGCCAAAGTGCAGGATGAAGc                                                                                          <  2:470941/251‑1 (MQ=255)                                                                                                                                cACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACCCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTGGATCACTCGCGCACCAACCTTGCGGCGTAAAGCCATTCTGTTGGCCAGAGAGCAGGATGAAGCCgg                                                                                       >  2:188509/1‑251 (MQ=255)                                                                                                                                 aCGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTGGATCACTCGCGCACCAACCTTGCGGCGTAAAGCCATTCTGTTGGCCAAAGTGCAGGATGAAGCCggc                                                                                      >  1:70923/1‑250 (MQ=255)                                                                                                                                  cGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTGGATCACTCGCGCACCAACCTTGCGGCGTAAAGCCATTCTGTTGGCCAAAGTGCAGGATGAAGCCGGTc                                                                                     >  1:127594/1‑251 (MQ=255)                                                                                                                                     aaaCTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTGGATCACTCGCGCACCAACCTTGCGGCGTAAAGCCATTCTGTTGGCCAAAGTGCAGGATGAAGCCGGTCATg                                                                                  >  1:250407/1‑251 (MQ=255)                                                                                                                                     aaaCTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCACAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTGGATCACTCGCGCACCAACCTTGCGGCGTAAAGCCATTCTGTTGGCCAAAGTGCAGGATGACGCCGGTCATg                                                                                  >  2:484937/1‑251 (MQ=255)                                                                                                                                       aCTGGAATCGGAAAGTTGATGCCGTGACCCAAGAAAAACGCTTCGAGCAGCGGGAAGCCCAGGAAATGGCTATCCAGCCACAGGACTGAATGCCCGAGGTTGACCTAAAGACATTGAGCCGCCAGAGTGGGGATCATGCGCAGTCCGGGAAGGTTGGTATGTTGCCTGAAGGTAACTGGATCACGGGCGCACCAACCTTGCTGCGTAAAGCCATTCTGTTGGCCAAAGTGCAGGATGAAGCCGGTCATGGt                                                                                <  1:584620/251‑1 (MQ=255)                                                                                                                                         tGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTGGATCACTCGCGCACCAACCTTGCGGCGTAAAGCCATTCTGTTGGCCAAAGTGCAGGATGAAGCCGGTCATGGTCt                                                                              >  2:582265/1‑251 (MQ=255)                                                                                                                                                     ggTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACCTTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTGGATCACTCGCGCACCAACCTTGCGGCGTAAAGCCATTCTGTTGGCCAAAGTGCAGGATGAAGCCGGTCATGGTCTGTATCTCtat                                                                    >  1:280082/1‑249 (MQ=255)                                                                                                                                                      gTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTGGATCACTCGCGCACCAACCTTGCGGCGTAAAGCCATTCTGTTGGCCAAAGTGCAGGATGAAGCCGGTCATGGTCTGTATCTCTATAgc                                                                 >  2:137490/1‑251 (MQ=255)                                                                                                                                                               tGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTGGATCACTCGCGCACCAACCTTGCGGCGTAAAGCCATTCTGTTGGCCAAAGTGCAGGATGAAGCCGGTCATGGTCTGTATCTCTATAGCGCCGCCGa                                                         >  1:503375/1‑250 (MQ=255)                                                                                                                                                                     aagaagAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGCAGGTAACTGGATCACTCGCGCACCAACCTTGCGGCGCAAAGCCATTCTGTTGGCCAAAGTGCAGGATGAAGCCGGTCATGGTCTGTATCTCTATAGCGCCGCCGAAACACTg                                                  >  1:257015/1‑251 (MQ=255)                                                                                                                                                                           aaCGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTGGATCCCTCGCGCACCAACCTTGCGGCGTACAGCCATTCTGTTGGCCAAAGTGCAGGATGAAGCCGGTCATGGTCTGTATCTCTATAGCGCCGCCGAAACACTGGGCTgc                                            >  2:529067/1‑251 (MQ=255)                                                                                                                                                                                     gCAACGGCGAGCCCAGGCAACGGCTCTCGAGCCACAGGACTGAATGCCCGTTGCTTACCGAAAGACATGGATCGGCCAGATTGGGCAGCATGCGCACTCCGAGATGGTTGGTATGTGGCCTGAAGTTAAATGGATCACTCGCGCACCACCCTTGCGGCGTAAAGCCATTCTGTTGGCCAAAGTGCAGGATGAAGCCGGTCATGGTCTGTATCGCTATAGCGCCGCCGAGTCACTGGGCTGAGCCAGagaag                                  <  2:257691/251‑1 (MQ=255)                                                                                                                                                                                          ggATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTGGATCACTCGCGCACCAACCTTGCGGCGTAAAGCCATTCTGTTGGCCAAAGTGCAGGATGAAGCCGGTCATGGTCTGTATCTCTATAGCGCCGCCGAAACACTGGGCTGCGCCAGAGAAGACATc                             <  1:57262/251‑1 (MQ=255)                                                                                                                                                                                              aGCCCAGGAACCGGATATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAACGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTGGATCACTCGCGCACCAACCTTGCGGCGTAAAGCCATTCTGTTGGCCAAAGTGCAGGATGAAGCCGGTCATGGTCTGTATCTCTATAGCGCCGCCGAAACACTGGGCTGCGCCAGAGAAGACATCTAcc                         >  2:75847/1‑251 (MQ=255)                                                                                                                                                                                              aGCCCAGGAAACGGCTATCGAGCCACAGGACGGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTGGATCACTCGCGCACCAACCTTGCGGCGTAAAGCCATTCTGTTGGCCAAAGTGCAGGATGAAGCCGGTCATGGTCTGTATCTCTATAGCGCCGCCGAAACACTGGGCTGCGCCAGAGAAGACATCTAc                          <  2:339255/250‑1 (MQ=255)                                                                                                                                                                                                 ccAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTGGATCACTCGCGCACCAACCTTGCGGCGTAAAGCCATTCTGTTGGCCAAAGTGCAGGATGAAGCCGGTCATGGTCTGTATCTCTATAGCGCCGCCGAAACACTGGGCTGCGCCAGAGAAGACATCTACCaaa                      <  1:582265/251‑1 (MQ=255)                                                                                                                                                                                                 ccAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTGGATCACTCGCGCACCAACCTTGCGGCGTAAAGCCATTCTGTTGGCCAAAGTGCAGGATGAAGCCGGTCATGGTCTGTATCTCTATAGCGCCGCCGAAACACTGGGCTGCGCCAGAGAAGACATCTACCaa                       >  1:470823/1‑250 (MQ=255)                                                                                                                                                                                                                      acaGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTGGATCACTCGCGCACCAACCTTGCGGCGTAAAGCCATTCTGTTGGCCAAAGTGCAGGATGAAGCCGGTCATGGTCTGTATCTCTATAGCGCCGCCGAAACACTGGGCTGCGCCAGAGAAGACATCTACCAAAAAATGCTCGACGGACGGATg  >  1:311125/1‑250 (MQ=255)                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                 CTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAATTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGGATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTGGATCACTCGCGCACCAACCTTGCGGCGTAAAGCCATTCTGTTGGCCAAAGTGCAGGATGAAGCCGGTCATGGTCTGTATCTCTATAGCGCCGCCGAAACACTGGGCTGCGCCAGAGAAGACATCTACCAAAAAATGCTCGACGGACGGATG  >  NC_000913/1453770‑1454232

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