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	96.8%	91.8 / NA	31	W22* (TGG→TGA)	paaA	phenylacetyl‑CoA 1,2‑epoxidase, monooxygenase subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (19/11);  minor base C (1/0);  total (20/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

AAACATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAATTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGGATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTGGATCACTCGCGCACCAACCTTGCGGCGTAAAGCCATTCTGTTGGCCAAAGTGCAGGATGAAGCCGGTCATGGTCTGTATCTCTATAGCGCCGCCGAAACACTGGGCTGCGCCAGAGAAGACATCTACCAAAAAATGCTCGACGGAC  >  NC_000913/1453757‑1454227                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                            aaaCATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAATTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTAc                                                                                                                                                                                                                              >  1:326081/1‑251 (MQ=255)       tAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAATTTGTGATTTTACTTAACTATTGTGTAACTTTCATACAACAATGTGATTCGTGTTTTTAATTAATTCACGAACACTGGAATAGTAAAGGTGATGACGTGACCCAAGAAGACCGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGACCGATGCTAACCGaaa                                                                                                                                                                                                                         >  2:226620/1‑250 (MQ=255)           gCACTGATAAATAATAATTTATAAAAATAGGGTGCGAAATCCGTCACATTTCAAACATACAAAATTTGTGAGTTTACTTAACTATTGTGTAACTGTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAATCTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGAAACGGACAGCCCAGGAAACGGCTATCGCGCCACAGGACTGAATGCCCGTTGCTTTCCTAAAGACa                                                                                                                                                                                                                     <  2:188158/250‑1 (MQ=255)             aCTGATCAATAATGATTTATAAAAATAGGGTGCGAACTCCGTCACAGTTCAAACATACAAAATTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTg                                                                                                                                                                                                                  >  2:82398/1‑251 (MQ=255)                     ataatGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAATTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGcc                                                                                                                                                                                                           >  1:220556/1‑250 (MQ=255)                                              gAAATCCGTCACAGTTCAAACATACAAAATTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCgag                                                                                                                                                                                 <  2:191801/251‑1 (MQ=255)                                              gAAATCCGTCACAGTTCAAACATACAAAATTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCg                                                                                                                                                                                   >  1:45020/1‑249 (MQ=255)                                                                       aaaaTTTGTGATTTTACTTAACTATTGTGTAACTTTCATAATACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGACGGt                                                                                                                                                         >  2:71951/1‑250 (MQ=255)                                                                           tttGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTg                                                                                                                                                    >  1:127892/1‑251 (MQ=255)                                                                                    ttACTTAACTAATGTGAAACTTTCATAAAACACTGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGCTTGGGCAGCATGCGAACTCCGAGATTGTTGGTATGTTGCCTGGAGGTAACTGGATCAcacg                                                                                                                                           >  2:239278/1‑248 (MQ=255)                                                                                                  tgtAACTTTCATAAAACAATGTGATTCGGGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGAATGGTGGTATGTTGCCTGAAGGTAACTGGATCACTCGCGCACAAACCTTgc                                                                                                                             >  1:296747/1‑251 (MQ=255)                                                                                                  tgtAAATTTCATAAAACAATGAGATTCGTGTTTTTAATTAATTCACGAACACTGGAATCGTAAAGGTGATGCCGTGACCAAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGCGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTGGATCACTCGCGCACCAACCTTgc                                                                                                                             >  2:267304/1‑251 (MQ=255)                                                                                                               aaaCAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGTATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGAAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGGAACTGGATCACTCGCGCACCAACCTTGCGGCGTAAAGccct                                                                                                                >  2:292052/1‑249 (MQ=255)                                                                                                                  cAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTAACTCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTGGATCCCTCGCGCACCAACCTTGCGGCGTAAAGCCATTCt                                                                                                             >  1:39185/1‑251 (MQ=255)                                                                                                                         aTTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTGGATCACTCGCGCACCAACCTTGCGGCGTAAAGACATTCTGTTGGcc                                                                                                      >  2:111187/1‑251 (MQ=255)                                                                                                                                  tttAATTAATTCACGAAAACTGGAATCGTAAAGGGGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACCGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTGGATCACTCGCGCACCAACCTTGCGGCGTAAAGCCATTCTGTTGGCCAAAGTGCa                                                                                              <  1:235522/250‑1 (MQ=255)                                                                                                                                     aattaattCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTGGATCACTCGCGCACCAACCTTGCGGCGTAAAGCCATTCTGTTGGCCAAAGTGCAGGAt                                                                                          >  2:70284/1‑251 (MQ=255)                                                                                                                                        taattCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTGGATCACTCGCGCACCAACCTTGCGGCGTAAAGCCATTCTGTTGGCCAAAGTGCAGGATGa                                                                                        >  2:286369/1‑250 (MQ=255)                                                                                                                                         aattCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGCAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTGGATCACTCGCGCACCAACCTTGCGGCGTAAAGCCATTCTGTTGGCCAAAGTGCAGGATGa                                                                                        >  1:245065/1‑249 (MQ=255)                                                                                                                                             cACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTGGATCACTCGCGCACCAACCTTGCGGCGTAAAGCCATTCTGTTGGCCAAAGTGCAGGATGAAGCCg                                                                                   <  1:284369/250‑1 (MQ=255)                                                                                                                                                 aaaaCTGGAATCGTAAAGGTGATGACTTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTAGCGAGCCACAGTACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTTGGCAGCATTCGCACTCCGATATTGTTGGTATTTTGCCTGAAGGTAACTGGATCACTCGCGCACCAACCTTGCGGCGTAAAGCCATTCTGTTGGCCAAAGTGCAGGATGAAGCCGGTCAt                                                                              <  1:45205/251‑1 (MQ=255)                                                                                                                                                                    tgatgaCGTGACCCAAGCAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACCGGACTGCATGCCCGATGCTTACCGAACGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGCGATTGTTGGTATGTTGCCTGAAGGTAACTGGATCACTCGCGCACCAACCTTGCGGCGTAAAGCCATTCTGTTGGCCAAAGGGCAGGATGAAGCCGGTCATGGTCTGTATCTCTATcgcg                                                           >  2:198626/1‑251 (MQ=255)                                                                                                                                                                                      agaaCGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTGGATCACTCGCGCACCAACCTTGCGGCGTAAAGCCATTCTGTTGGCCAAAGTGCAGGATGAAGCCGGTCATGGTCTGTATCTCTATAGCGCCGCCGAAACACTGGGCt                                         <  1:111187/251‑1 (MQ=255)                                                                                                                                                                                                 aGCAACGGATAGCCCAGGCAACGGCTCTCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTGGATCACTCGCGCACCAACCTTGCGGCGTAAAGCCATTCTGTTGGCCAAAGTGCAGGATGAAGCCGGTCATGGTCTGTATCTCTATAGCGCCGCCGACACACTGGGCTGCGCCAgaga                               >  2:275046/1‑250 (MQ=255)                                                                                                                                                                                                       ggATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTGGATCACTCGCGCACCAACCTTGCGGCGTAAAGCCATTCTGTTGGCCAAAGTGCAGGATGAAGCCGGTCATGGTCTGTATCTCTATAGCGCCGCCGAAACACTGGGCTGCGCCAGAGAAGACATc                        >  1:127128/1‑251 (MQ=255)                                                                                                                                                                                                            gCCCAGGAAACGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTGGATCACTCGCGCACCAACCTTGCGGCGTAAAGCCATTCTGTTGGCCAAAGTGCAGGATGAAGCCGGTCATGGTCTGTATCTCTATAGCGCCGCCGAAACACTGGGCTGCGCCAGAGAAGACATCTACCa                   <  1:70284/251‑1 (MQ=255)                                                                                                                                                                                                                aGGAAACGGCTATAGCGCCACAGGACTGAATGCCCGATGCTTACCGAAAAACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTGGATCACTCGCGAACCAACCTTGCGGCGTAAAGCCATTCTGTTGGCCAAAGTGCAGGATGAACCCGGTCATGGTCTGTATCTCTCTAGCGCCGCCGAAACACTGGGCTGCGCCAGAGAAGACATCTCCCaaaa                >  2:209941/1‑250 (MQ=255)                                                                                                                                                                                                                   aaaCGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTGGATCACTCGCGCACCAACCTTGCGGCGTAAAGCCATTCTGTTGGCCAAAGTGCAGGATGAAGCCTGTCATGGTCTGTATCTCTATAGCTCCGCCGAAACACTGGGCTGCGCCAGAGAAGACATCTACCAAAAAATg            <  1:267304/251‑1 (MQ=255)                                                                                                                                                                                                                     aCGGCTATCGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTGGATCACTCGCGCACCAACCTTGCGGCGTAAAGCCATTCTGTTGGCCAAAGTGCAGGATGAAGCCGGTCATGGTCTGTATCTCTATAGCGCCGCCGAAACACTGGGCTGCGCCAGAGAAGACATCTACCAAAAAATGCt          <  1:54939/251‑1 (MQ=255)                                                                                                                                                                                                                       gtgtATCTAGCCACAGGACTGAATGCCCGATGCTTACCGAATGATAGTGATCCGCCAGATTGGTCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTGGATCACTCGCGCACCAACCTTGCGGCGTAAAGCCATTCTGTTGGCCAAAGTGCAGGATGAAGCCGGTCATGGTCTGTATCTCTATAGCGCCGCCGAAACACTGGGCTGCGCCAGAGAAGACATCTACCAAAAAATGCTCg        <  2:215715/248‑1 (MQ=255)                                                                                                                                                                                                                             cGAGCCACAGGACTGAATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTGGATCACTCGCGCACCAACCGTGCGGCGTAAAGCCATTCTGTTGGCCAAAGTGCAGGATGAAGCCGGTCATGGTCTGTAGCTCTATAGCGCCGCCGAAACACGGGGCTGCGCCAGAGAAGACTGCTACCAAAGAATGCTCgacggac  <  2:245065/251‑1 (MQ=255)                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                            AAACATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAATTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTATCGAGCCACAGGACTGGATGCCCGATGCTTACCGAAAGACATTGATCCGCCAGATTGGGCAGCATGCGCACTCCGAGATTGTTGGTATGTTGCCTGAAGGTAACTGGATCACTCGCGCACCAACCTTGCGGCGTAAAGCCATTCTGTTGGCCAAAGTGCAGGATGAAGCCGGTCATGGTCTGTATCTCTATAGCGCCGCCGAAACACTGGGCTGCGCCAGAGAAGACATCTACCAAAAAATGCTCGACGGAC  >  NC_000913/1453757‑1454227

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