Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,459,052 G→A *263* (TAG→TAA paaG → 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_0009131,459,0520GA100.0% 101.1 / NA 32*263* (TAG→TAApaaG1,2‑epoxyphenylacetyl‑CoA isomerase, oxepin‑CoA‑forming
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (13/19);  total (13/19)

AAACCAATACGCTCGATACGCAACTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAGCATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAGTTGCTGCCAGTCATGGACACCAGGTTTTACTGTATGACATTTCTGCTGAAGCGC  >  NC_000913/1458922‑1459174
                                                                                                                                  |                                                                                                                          
aaaCCAATACGCTCGATACGCAACTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACAtgatga                                                                                                                    >  1:952475/1‑139 (MQ=255)
        aCGCTCGATACGCAACTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAAtgtg                                                                                                            <  1:842018/139‑1 (MQ=255)
              gATACGCAACTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACt                                                                                                      >  1:158857/1‑139 (MQ=255)
                tACGCAACTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACtgt                                                                                                    <  1:71545/139‑1 (MQ=255)
                    cAACTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCa                                                                                                <  2:952475/139‑1 (MQ=255)
                            tCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTgg                                                                                        <  2:942898/139‑1 (MQ=255)
                                  aCGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCgg                                                                                  <  2:767057/139‑1 (MQ=255)
                                    gTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACt                                                                                                      <  2:523489/117‑1 (MQ=255)
                                    gTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACt                                                                                                      >  1:523489/1‑117 (MQ=255)
                                               cGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAAtgtg                                                                                                            >  1:369500/1‑100 (MQ=255)
                                               cGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAAtgtg                                                                                                            <  2:369500/100‑1 (MQ=255)
                                               cGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCa                                                                             <  2:1038772/131‑1 (MQ=255)
                                               cGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCa                                                                             >  1:1038772/1‑131 (MQ=255)
                                                      ccGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTg                                                              <  1:1008817/139‑1 (MQ=255)
                                                             gaGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAAtgtg                                                                                                            >  1:1098873/1‑86 (MQ=255)
                                                             gaGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAAtgtg                                                                                                            <  2:1098873/86‑1 (MQ=255)
                                                                           gTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAGTTGCTGCCAGTCATg                                         <  1:642928/139‑1 (MQ=255)
                                                                               aGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAGTTGCTGCCAGTCATGGaca                                     >  2:427700/1‑139 (MQ=255)
                                                                               aGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAGTTGCTGCCAGTCATGGaca                                     <  1:427700/139‑1 (MQ=255)
                                                                                  tgtCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAGTTGCTGCCAGTCATGGACACCa                                  <  1:1093390/139‑1 (MQ=255)
                                                                                   gtCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACat                                                                                                                         >  2:247684/1‑51 (MQ=255)
                                                                                   gtCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACat                                                                                                                         <  1:247684/51‑1 (MQ=255)
                                                                                   gtCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTgg                                                                                                  >  2:162633/1‑74 (MQ=255)
                                                                                   gtCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTgg                                                                                                  <  1:162633/74‑1 (MQ=255)
                                                                                   gtCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTggg                                                                                       <  2:659639/85‑1 (MQ=255)
                                                                                   gtCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTggg                                                                                       >  1:659639/1‑85 (MQ=255)
                                                                                                  gCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCa                                                                                                <  2:1055188/61‑1 (MQ=255)
                                                                                                  gCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCa                                                                                                >  1:1055188/1‑61 (MQ=255)
                                                                                                  gCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAGTTGCTGCCAGTCATGGACACCAGGTTTTACTGTATGAc                  <  1:289935/139‑1 (MQ=255)
                                                                                                            cACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAGTTGCTGCCAGTCATGGACACCAGGTTTTACTGTATGACATTTctgctg        <  2:158857/139‑1 (MQ=255)
                                                                                                                  aGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAGTTGCTGCCAGTCATGGACACCAGGTTTTACTGTATGACATTTCTGCTGAAgcgc  >  1:4955/1‑139 (MQ=255)
                                                                                                                  aGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAGTTGCTGCCAGTCATGGACACCAGGCTTTACTGTATGACATTTCTGCTGAAgcgc  >  2:886362/1‑139 (MQ=255)
                                                                                                                                  |                                                                                                                          
AAACCAATACGCTCGATACGCAACTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAGCATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAGTTGCTGCCAGTCATGGACACCAGGTTTTACTGTATGACATTTCTGCTGAAGCGC  >  NC_000913/1458922‑1459174

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: