Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,458,944	A→G	Q227Q (CAA→CAG)	paaG →	putative ring 1,2‑epoxyphenylacetyl‑CoA isomerase (oxepin‑CoA forming)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,458,944	0	A	G	100.0%	177.6 / NA	48	Q227Q (CAA→CAG)	paaG	putative ring 1,2‑epoxyphenylacetyl‑CoA isomerase (oxepin‑CoA forming)
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (29/19);  total (29/19)

GGTTGTTGATGATGAAACGCTGGCAGATACCGCGCAACAGCTGGCACGGCATCTGGCGACACAACCGACATTTGGTCTTGGACTTATCAAGCAAGCGATAAATAGCGCTGAAACCAATACGCTCGATACGCAACTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAGCATGATGATAAATGTGCAAACTGT  >  NC_000913/1458812‑1459076                                                                                                                                     |                                                                                                                                     ggTTGTTGATGATGAAACGCTGGCAGATACCGCGCAACAGCTGGCACGGCATCTGGCGACACAACCGACATTTGGTCTTGGACTTATCAAGCAAGCGATAAATAGCGCTGAAACCAATACGCTCGATACGCAGCTGGATc                                                                                                                               <  2:393850/140‑1 (MQ=255)     gttgATGATGAAACGCTGGCAGATACCGCGCAACAGCTGGCACGGCATCTGGCGACACAACCGACATTTGGTCTTGGACTTATCAAGCAAGCGATAAATAGCGCTGAAACCAATACGCTCGATACGCAGCTGGATc                                                                                                                               >  2:1090000/1‑136 (MQ=255)     gttgATGATGAAACGCTGGCAGATACCGCGCAACAGCTGGCACGGCATCTGGCGACACAACCGACATTTGGTCTTGGACTTATCAAGCAAGCGATAAATAGCGCTGAAACCAATACGCTCGATACGCAGCTGGATc                                                                                                                               <  1:1090000/136‑1 (MQ=255)     gttgATGATGAAACGCTGGCAGATACCGCGCAACAGCTGGCACGGCATCTGGCGACACAACCGACATTTGGTCTTGGACTTATCAAGCAAGCGATAAATAGCGCTGAAACCAATACGCTCGATACGCAGCTGGATCTGGa                                                                                                                           <  2:514163/140‑1 (MQ=255)           gatgaAACGCTGGCAGATACCGCGCAACAGCTGGCACGGCATCTGGCGACACAACCGACATTTGGTCTTGGACTTATCAAGCAAGCGATAAATAGCGCTGAAACCAATACGCTCGATACGCAGCTGGATCTGGAACGTGa                                                                                                                     <  1:1086257/140‑1 (MQ=255)           gatgaAACGCTGGCAGATACCGCGCAACAGCTGGCACGGCATCTGGCGACACAACCGACATTTGGTCTTGGACTTATCAAGCAAGCGATAAATAGCGCTGAAACCAATACGCTCGATACGCAGCTGGATCTGGAACGTGa                                                                                                                     <  2:384301/140‑1 (MQ=255)                aaCGCTGGCAGATACCGCGCAACAGCTGGCACGGCATCTGGCGACACAACCGACATTTGGTCTTGGACTTATCAAGCAAGCGATAAATAGCGCTGAAACCAATACGCTCGATACGCAGCTGGATCTGGAACGTGACTATc                                                                                                                >  2:241646/1‑140 (MQ=255)                    cTGGCAGATACCGCGCAACAGCTGGCACGGCATCTGGCGACACAACCGACATTTGGTCTTGGACTTATCAAGCAAGCGATAAATAGCGCTGAAACCAATACGCTCGATACGCAGCTGGATCTGGAACGTGACTATCAGCg                                                                                                            <  2:913061/140‑1 (MQ=255)                     tGGCAGATACCGCGCAACAGCTGGCACGGCATCTGGCGACACAACCGACATTTGGTCTTGGACTTATCAAGCAAGCGATAAATAGCGCTGAAACCAATACGCTCGATACGCAGCTGGATCTGGAACGTGACTATCAGCGa                                                                                                           >  2:4659/1‑140 (MQ=255)                       gCAGATACCGCGCAACAGCTGGCACGGCATCTGGCGACACAACCGACATTTGGTCTTGGACTTATCAAGCAAGCGATAAATAGCGCTGAAACCAATACGCTCGAGACGCAGCTGGATCTGGAACGTGACTATCAGCGACt                                                                                                         >  2:794286/1‑140 (MQ=255)                              ccGCGCAACAGCTGGCACGGCATCTGGCGACACAACCGACATTTGGTCTTGGACTTATCAAGCAAGCGATAAATAGCGCTGAAACCAATACGCTCGATACGCAGCTGGATCTGGAACGTGACTATCAGCGACTTGCcgga                                                                                                  >  1:492159/1‑140 (MQ=255)                                 cgcAACAGCTGGCGCGGCATCTGGCGACACAACCGACATTTGGTCTTGGACTTATCAAGCAAGCGATAAATAGCGCTGAAACCAATACGCTCGATACGCAGCTGGATCTGGAACGTGACTATCAGCGACTTGCcggacgg                                                                                               >  1:613591/1‑140 (MQ=255)                                 cgcAACAGCTGGCACGGCATCTGGCGACACAACCGACATTTGGTCTTGGACTTATCAAGCAAGCGATAAATAGCGCTGAAACCAATACGCTCGATACGCAGCTGGATCTGGAACGTGACTATCAGCGACTTGCcggacgg                                                                                               >  1:240066/1‑140 (MQ=255)                                 cgcAACAGCTGGCACGGCATCTGGCGACACAACCGACATTTGGTCTTGGACTTATCAAGCAAGCGATAAATAGCGCTGAAACCAATACGCTCGATACGCAGCTGGATCTGGAACGTGACTATCAGCGACTTGCcggacgg                                                                                               >  1:160681/1‑140 (MQ=255)                                 cgcAACAGCTGGCACGGCATCTGGCGACACAACCGACATTTGGTCTTGGACTTATCAAGCAAGCGATAAATAGCGCTGAAACCAATACGCTCGATACGCAGCTGGATCTGGAACGTGACTATCAGCGACTTGCcggacgg                                                                                               >  2:348391/1‑140 (MQ=255)                                 cgcAACAGCTGGCACGGCATCTGGCGACACAACCGACATTTGGTCTTGGACTTATCAAGCAAGCGATAAATAGCGCTGAAACCAATACGCTCGATACGCAGCTGGATCTGGAACGTGACTATCAGCGACTTGCcggacgg                                                                                               >  1:71166/1‑140 (MQ=255)                                 cgcAACAGCTGGCACGGCATCTGGCGACACAACCGACATTTGGTCTTGGACTTATCAAGCAAGCGATAAATAGCGCTGAAACCAATACGCTCGATACGCAGCTGGATCTGGAACGTGACTATCAGCGACTTGCcggacgg                                                                                               >  2:1045523/1‑140 (MQ=255)                                 cgcAACAGCTGGCACGGCATCTGGCGACACAACCGACATTTGGTCTTGGACTTATCAAGCAAGCGATAAATAGCGCTGAAACCAATACGCTCGATACGCAGCTGGATCTGGAACGTGACTATCAGCGACTTGCcggacgg                                                                                               >  1:902626/1‑140 (MQ=255)                                 cgcAACAGCTGGCACGGCATCTGGCGACACAACCGACATTTGGTCTTGGACTTATCAAGCAAGCGATAAATAGCGCTGAAACCAATACGCTCGATACGCAGCTGGATCTGGAACGTGACTATCAGCGACTTGCcgga                                                                                                  >  1:194876/1‑137 (MQ=255)                                 cgcAACAGCTGGCACGGCATCTGGCGACACAACCGACATTTGGTCTTGGACTTATCAAGCAAGCGATAAATAGCGCTGAAACCAATACGCTCGATACGCAGCTGGATCTGGAACGTGACTATCAGCGACTTGCcgga                                                                                                  >  2:828434/1‑137 (MQ=255)                                 cgcAACAGCTGGCACGGCATCTGGCGACACAACCGACATTTGGTCTTGGACTTATCAAGCAAGCGATAAATAGCGCTGAAACCAATACGCTCGATACGCAGCTGGATCTGGAACGTGACTATCAGCGACTTGCcgga                                                                                                  <  2:194876/137‑1 (MQ=255)                                 cgcAACAGCTGGCACGGCATCTGGCGACACAACCGACATTTGGTCTTGGACTTATCAAGCAAGCGATAAATAGCGCTGAAACCAATACGCTCGATACGCAGCTGGATCTGGAACGTGACTATCAGCGACTTGCcgga                                                                                                  <  1:828434/137‑1 (MQ=255)                                     aCAGCTGGCACGGCATCTGGCGACACAACCGACATTTGGTCTTGGACTTATCAAGCAAGCGATAAATAGCGCTGAAACCAATACGCTCGATACGCAGCTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAgcg                                                                                           <  1:878166/140‑1 (MQ=255)                                     aCAGCTGGCACGGCATCTGGCGACACAACCGACATTTGGTCTTGGACTTATCAAGCAAGCGATAAATAGCGCTGAAACCAATACGCTCGATACGCAGCTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAgcg                                                                                           <  1:241646/140‑1 (MQ=255)                                                  cATCTGGCGACACAACCGACATTTGGTCTTGGACTTATCAAGCAAGCGATAAATAGCGCTGAAACCAATACGCTCGATACGCAGCTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGa                                                                              <  1:794286/140‑1 (MQ=255)                                                  cATCTGGCGACACAACCGACATTTGGTCTTGGACTTATCAAGCAAGCGATAAATAGCGCTGAAACCAATACGCTCGATACGCAGCTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGa                                                                              <  2:786729/140‑1 (MQ=255)                                                            cacaACCGACATTTGGTCTTGGACTTATCAAGCAAGCGATAAATAGCGCTGAAACCAATACGCTCGATACGCAGCTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGt                                                                    >  2:1151551/1‑140 (MQ=255)                                                               aaCCGACATTTGGTCTTGGACTTATCAAGCAAGCGATAAATAGCGCTGAAACCAATACGCTCGATACGCAGCTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCg                                                                 >  1:152696/1‑140 (MQ=255)                                                                           gTCTTGGACTTATCAAGCAAGCGATAAATAGCGCTGAAACCAATACGCTCGATACGCAGCTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTaaa                                                     >  1:821711/1‑140 (MQ=255)                                                                               tGGACTTATCAAGCAAGCGATAAATAGCGCTGAAACCAATACGCTCGATACGCAGCTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCt                                                 >  1:279983/1‑140 (MQ=255)                                                                                    ttATCAAGCAAGCGATAAATAGCGCTGAAACCAATACGCTCGATACGCAGCTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCg                                            <  2:503324/140‑1 (MQ=255)                                                                                    ttATCAAGCAAGCGATAAATAGCGCTGAAACCAATACGCTCGATACGCAGCTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCg                                            <  1:348391/140‑1 (MQ=255)                                                                                      aTCAAGCAAGCGATAAATAGCGCTGAAACCAATACGCTCGATACGCAGCTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCa                                          >  2:360042/1‑140 (MQ=255)                                                                                         aagcaagcGATAAATAGCGCTGAAACCAATACGCTCGATACGCAGCTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGtt                                       <  2:160681/140‑1 (MQ=255)                                                                                           gcaagcGATAAATAGCGCTGAAACCAATACGCTCGATACGCAGCTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCa                                     <  2:902626/140‑1 (MQ=255)                                                                                                cGATAAATAGCGCTGAAACCAATACGCTCGATACGCAGCTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACgggg                                >  1:63363/1‑140 (MQ=255)                                                                                                cGATAAATAGCGCTGAAACCAATACGCTCGATACGCAGCTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACgggg                                >  1:258429/1‑140 (MQ=255)                                                                                                         gcgcTGAAACCAATACGCTCGATACGCAGCTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAGCat                       <  2:492159/140‑1 (MQ=255)                                                                                                          cgcTGAAACCAATACGCTCGATACGCAGCTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAGCatg                      >  1:62815/1‑140 (MQ=255)                                                                                                          cgcTGAAACCAATACGCTCGATACGCAGCTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAGCatg                      >  1:655947/1‑140 (MQ=255)                                                                                                            cTGAAACCAATACGCTCGATACGCAGCTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAGCatgat                    >  1:485297/1‑140 (MQ=255)                                                                                                               aaaCCAATACGCTCGATACGCAGCTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACg                                   >  1:1032842/1‑122 (MQ=255)                                                                                                               aaaCCAATACGCTCGATACGCAGCTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACg                                   <  2:1032842/122‑1 (MQ=255)                                                                                                               aaaCCAATACGCTCGATACGCAGCTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAGCATgatgat                 >  2:463412/1‑140 (MQ=255)                                                                                                                  ccAATACGCTCGATACGCAGCTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAGCATGATGATaaa              <  1:1045523/140‑1 (MQ=255)                                                                                                                       aCGCTCGATACGCAGCTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAGCATGATGATAAATGTGc         >  2:42246/1‑140 (MQ=255)                                                                                                                           tCGATACGCAGCTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAGCATGATGATAAATGTGCAAAc     >  1:512151/1‑140 (MQ=255)                                                                                                                              aTACGCAGCTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAGCATGATGATAAATGTGCAAACtgt  >  1:838125/1‑140 (MQ=255)                                                                                                                                     |                                                                                                                                     GGTTGTTGATGATGAAACGCTGGCAGATACCGCGCAACAGCTGGCACGGCATCTGGCGACACAACCGACATTTGGTCTTGGACTTATCAAGCAAGCGATAAATAGCGCTGAAACCAATACGCTCGATACGCAACTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAGCATGATGATAAATGTGCAAACTGT  >  NC_000913/1458812‑1459076

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