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%	116.3 / 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 (19/15);  total (19/15)

ACTGCGGAAAACATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAATTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGGATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTGGATCACTCGCGCACCAACCTTGCGGCGTAAAGCCATTCTGTTGGCCAAAGTGCAGGATGAAGCCGGTCATGGTCTGTATCTCTATAGCGCCGCCGAAACACTGGGCTGCGCCAGAGAAGACATCTACCAAAAAATGCTCGACGGACGGATGAAATAT  >  NC_000913/1453749‑1454238                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                       aCTGCGGCAAACCTTCATGCAATGATAAATAATGACTTATAAACCTAGGGTGCGAAATCAGTCACCGTTAAAACATACACCATTTCTGATTTAACTTAACTATTGTGTCACTTTCCTCAAACAATATGATTCGTGTTTTTAATTAATTACCGAAAACTGGAATCGTAAAGGTGATGCCGTGACAACAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCAAAGGACTGAATGCCCg                                                                                                                                                                                                                                                 >  2:451462/1‑251 (MQ=255)  atGCGGAAAACATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAATTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGa                                                                                                                                                                                                                                                <  2:63649/250‑1 (MQ=255)            cATTATTGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAATTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCg                                                                                                                                                                                                                                       >  1:361336/1‑250 (MQ=255)                   gCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAATTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACAt                                                                                                                                                                                                                               >  1:254059/1‑251 (MQ=255)                        gATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAATTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATc                                                                                                                                                                                                                          <  2:303247/251‑1 (MQ=255)                          taaattaTGACTTATAAAAATAGGGTTCGCAATCCGTCACAGTTAAAACATACCAACATTGAAATTATGCTTAACTATTGCGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGACTCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTATGAGCAACGAATAGCACAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCg                                                                                                                                                                                                                        >  2:437818/1‑251 (MQ=255)                            aataatGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAATTTGTGATTT‑ACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCa                                                                                                                                                                                                                     <  1:109345/251‑1 (MQ=255)                                 tGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAATTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGAtt                                                                                                                                                                                                                 >  2:421749/1‑251 (MQ=255)                                             atcgggTGCGAAATGCATCACAGTTCAACCATACAAAATTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTACCGAAAACTGGAATCGCACAGGCGATGACGTGACCCAAGACGAACGCTTTGAGCAACGGATAGACCAGGAACCGGCTACAGAGACACAGGACTGAATGCCCGATGCTTACGGCAAGCCATTTATACGCCAGATTGCGCAGCATg                                                                                                                                                                                                       >  2:151490/4‑249 (MQ=255)                                                  gTGCGAAATCCGTCACAGTTCAAACATACAAAATTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTc                                                                                                                                                                                                >  1:192913/1‑251 (MQ=255)                                                    gCGAAATCCGTCACAGTTCAAACATACAAAATTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGAGTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGAGGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCTCTCCg                                                                                                                                                                                              <  2:443127/251‑1 (MQ=255)                                                           ccGTCACAGTTCAAACATACAAAATTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGAttg                                                                                                                                                                                        <  1:90726/250‑1 (MQ=255)                                                               cacgtttaaaacaggcAAAATTTGTGATTTTACTTAAGTATTGTGTAACTTGCATAAAACAATTGGAGTCGTGTTTTAAATTAATTCACGAAAACGGGAATCGAAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAACCCAGGAAACGGCTAATGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATGGATCCGCCAGATTGGGCAGCAGGCGCGATCCGCGATTGTTGGt                                                                                                                                                                                   <  2:380554/236‑1 (MQ=255)                                                                  aGTTCAAACATACAAAATTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAACACTGGAATAGTAAAGGTGATGACGTGAACCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATg                                                                                                                                                                                >  1:445962/1‑251 (MQ=255)                                                                               aaaaTTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGt                                                                                                                                                                    >  2:4957/1‑250 (MQ=255)                                                                                     tgtgATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTgg                                                                                                                                                              >  1:409064/1‑250 (MQ=255)                                                                                                 tAACTATTGTGTAATTTTCATAAAATAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCATATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTGGATCACTCGCGCa                                                                                                                                                  <  1:150863/250‑1 (MQ=255)                                                                                                        tgtgTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAAATGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGCCTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTGGATCACTCGCGCACCAACCtt                                                                                                                                          >  1:154838/1‑251 (MQ=255)                                                                                                          tgtAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTGGATCACTCGCGCACCAACCTTgc                                                                                                                                        >  1:337608/1‑251 (MQ=255)                                                                                                          tgtAACTTTCATAAAAAAATGAGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGAAAAGGTGATGACGTGAACCAAGAAGAAAGCTTTGAGCAACGGATAGCCCAGGAAACGGATATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTGGATCACTCGCGCACCAACCTTgc                                                                                                                                        >  2:463578/1‑251 (MQ=255)                                                                                                           gtAACTTTCATAAAACAATGTGAATCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTGGATCACTCGCGCACCAACCTTgc                                                                                                                                        <  1:238160/250‑1 (MQ=255)                                                                                                                      aaaaCAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTGGATCACTCGCGCACCAACCTTgcggcg                                                                                                                                    >  2:221395/1‑243 (MQ=255)                                                                                                                      aaaaCAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTGGATCACTCGCGCACCAACCTTgcggcg                                                                                                                                    <  1:221395/243‑1 (MQ=255)                                                                                                                                      tgtTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGACGGTAACTGGATCACTCGCGCACCAACCTTGCGGCGTAAAGCCATTCTGTTGGCCAAAGt                                                                                                            >  1:154397/1‑251 (MQ=255)                                                                                                                                             aattaattCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTGGATCACTCGCGCACCAACCTTGCGGCGTAAAGCCATTCTGTTGGCCAAAGTGCAGGAt                                                                                                     <  1:463578/251‑1 (MQ=255)                                                                                                                                                       cGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTGGATCACTCGCGCGCCAACCTTGCGGCGTAAAGCCATTCTGTTTGCCAAAGTGCAGGATGAAGCCGGTc                                                                                           <  2:61095/251‑1 (MQ=255)                                                                                                                                                               ggAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTGGATCACTCGCGCACCAACCTTGCGGCGTAAAGCCATTCTGTTGGCCAAAGTGCAGGATGAAGCCGGTCATGGTc                                                                                     <  1:83622/249‑1 (MQ=255)                                                                                                                                                                              atgaCGTTACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTTATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTGGATCACTCGCGCACCAACCTTGCGGCGTAAAGCCATTCTGTTGGCTAAAGTGCAGGATGAAGCCGGTCATGGTCTGTATCTCTATATcgcc                                                                    <  2:89800/251‑1 (MQ=255)                                                                                                                                                                                  cGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTGGATCACTCGCGCACCAACCTTGCGGCGTAAAGCCATTCTGTTGGCCAAAGTGCAGGATGAAGCCGGTCATGGTCTGTATCTCTATAGCGccgccg                                                                <  2:154838/251‑1 (MQ=255)                                                                                                                                                                                                       tGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTGGATCACTCGCGCACCAACCTTGCGGCGTAAAGCCATTCTGTTGGCCAAAGTGCAGGATGAAGCCGGTCATGGTCTGTATCTCTATAGCGCCGCCGAAACACTGGGCTGCGCCag                                             >  1:326613/1‑249 (MQ=255)                                                                                                                                                                                                             ccGGATAGCCCAGGGGACAGCTATCGAGACACAGGACTGACTGGCCGATGCTTACCGAAAGACATTGATCCGCCAGCTTGGGCAGCATGCTCACTCCGAGATTGTTGGTATGTTGCCTGAAGTCAACTGGATCACTCGCGCACCAACCTAGCGGCGTCACGCCATGCTGTTGGCCCACGTGCAGGATGAAGCCGTTCATGGTCTGGATCTCTAAAGCGCCGCCGACACACTGGGCTGCGCCAGAGCAGAcc                                     >  2:127136/2‑250 (MQ=255)                                                                                                                                                                                                                    gCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTCCCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTGGATCACTAGCGCACCAACCTTGCGGCGTAAAGCCATTCTGTTGGCCAAAGTGCAGGATGAAGCCGGTCATGGTCTGTATCTCTATAGCGCCGCCGAAACACTGGGCTGCGCCAGAGAAGACATCTAcc                               >  1:203041/1‑250 (MQ=255)                                                                                                                                                                                                                               ggCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTGGATCACTCGCGCACCAACCTTGCGGCGTAAAGCCATTCTGTTGGCCAAAGTGCAGGATGAAGCCGGTCATGGTCTGTATCTCTATAGCGCCGCCGAAACACTGGGCTGCGCCAGAGAAGACATCTACCAAAAAATGCTCg                   >  1:439566/1‑251 (MQ=255)                                                                                                                                                                                                                                                aCTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTGGATCACTCGCGCACCAACCTTGCGGCGTAAAGCCATTCTGTTGGCCAAAGTGCAGGATGAAGCCGGTCATGGTCTGTATCTCTATAGCGCCGCCGAAACACTGGGCTGCGCCAGAGAAGACATCTACCAAAAAATGCTCGACGGACGGATGAAatat  <  1:6806/251‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                       ACTGCGGAAAACATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAATTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGGATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTGGATCACTCGCGCACCAACCTTGCGGCGTAAAGCCATTCTGTTGGCCAAAGTGCAGGATGAAGCCGGTCATGGTCTGTATCTCTATAGCGCCGCCGAAACACTGGGCTGCGCCAGAGAAGACATCTACCAAAAAATGCTCGACGGACGGATGAAATAT  >  NC_000913/1453749‑1454238

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