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% 166.7 / NA 52*263* (TAG→TAA) paaG1,2‑epoxyphenylacetyl‑CoA isomerase, oxepin‑CoA‑forming
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (25/27);  total (25/27)

CTGAAACCAATACGCTCGATACGCAACTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAGCATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAGTTGCTGCCAGTCATGGACACCAGGTTTTACTGTATGACATTTCTGCTGAAGCGCTGACCCGCGCAA  >  NC_000913/1458919‑1459186
                                                                                                                                     |                                                                                                                                      
cTGAAACCAATACGCTCGATACGCAACTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACatga                                                                                                                                   >  1:1244601/1‑139 (MQ=255)
   aaaCCAATACGCTCGATACGCAACTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACAtgatga                                                                                                                                <  2:1518260/139‑1 (MQ=255)
     aCCAATACGCTCGATACGCAACTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATa                                                                                                                              <  1:1425397/139‑1 (MQ=255)
        aaTACGCTCGATACGCAACTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATa                                                                                                                              <  1:1105442/136‑1 (MQ=255)
        aaTACGCTCGATACGCAACTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATa                                                                                                                              >  2:1105442/1‑136 (MQ=255)
             gCTCGATACGCAACTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCa                                                                                                                      <  1:1132616/139‑1 (MQ=255)
             gCTCGATACGCAACTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCa                                                                                                                      <  1:1319791/139‑1 (MQ=255)
             gCTCGATACGCAACTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCa                                                                                                                      <  1:918398/139‑1 (MQ=255)
                 gATACGCAACTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACt                                                                                                                  >  1:1552420/1‑139 (MQ=255)
                        aaCTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATaa                                                                                                                             <  2:1333815/121‑1 (MQ=255)
                        aaCTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATaa                                                                                                                             >  1:1333815/1‑121 (MQ=255)
                                  ggAGCGTGACTATCAGCGACCTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATCAATGTGCAAACTGTGGCAGTGATTGGGAg                                                                                                 >  2:1218977/1‑139 (MQ=255)
                                  ggAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACat                                                                                                                                     <  1:1052604/103‑1 (MQ=255)
                                  ggAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACat                                                                                                                                     >  2:1052604/1‑103 (MQ=255)
                                  ggAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAg                                                                                                 >  1:1234140/1‑139 (MQ=255)
                                  ggAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAg                                                                                                 <  2:1244601/139‑1 (MQ=255)
                                     aCGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCgg                                                                                              <  2:1437233/139‑1 (MQ=255)
                                       gTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAAtgtg                                                                                                                        >  2:1072506/1‑111 (MQ=255)
                                       gTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAAtgtg                                                                                                                        <  1:1072506/111‑1 (MQ=255)
                                           cTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGATCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCAt                                                                                        >  1:1039569/1‑139 (MQ=255)
                                                aGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATggggg                                                                                   <  2:1234140/139‑1 (MQ=255)
                                                     cTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGggcaggc                                                                              <  1:1380885/139‑1 (MQ=255)
                                                       tGCCGGCCGGACCGCAGCTTATCGGCAAGGTGTCAGTGCGCTCCTGCCTAAACGCTCACCGCAGTACACGGGGAAATAACATGATGATAAATGTGCAAACAGTGGCAGTGATAGGGAGCGCCACCATGGGGGCAGGCAt                                                                            >  1:347468/1‑139 (MQ=255)
                                                        gCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCAtt                                                                           >  1:1652133/1‑139 (MQ=255)
                                                               ggaGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAg                                                                    <  2:222736/139‑1 (MQ=255)
                                                                 aGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTGACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGTGAGCGGCACCATGGGggc                                                                                  >  1:403322/1‑123 (MQ=255)
                                                                 aGCGCAGAGTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGggc                                                                                  <  2:403322/123‑1 (MQ=255)
                                                                           aTCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAGTTGCTGCCAGTc                                                        >  1:503582/1‑139 (MQ=255)
                                                                                  aGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAGTTGCTGCCAGTCATGGaca                                                 <  2:1039569/139‑1 (MQ=255)
                                                                                      gtCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAGTTGCTGCCAGTCATGGACACCAg                                             >  1:601884/1‑139 (MQ=255)
                                                                                           tGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAGTTGCTGCCAGTCATg                                                     >  1:809893/1‑126 (MQ=255)
                                                                                           tGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAGTTGCTGCCAGTCATg                                                     <  2:809893/126‑1 (MQ=255)
                                                                                           tGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAGTTGCTGCCAGTCATGGACACCAGGtttt                                        >  2:938107/1‑139 (MQ=255)
                                                                                               ttCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAGTTGCTGCCAGTCATGGACACCAGGTTTTACTg                                    <  2:601884/139‑1 (MQ=255)
                                                                                                          aCGCTCACCGCAGTTCACGGGGAAATAACATGATGATaaa                                                                                                                            <  2:558182/40‑1 (MQ=255)
                                                                                                          aCGCTCACCGCAGTTCACGGGGAAATAACATGATGATaaa                                                                                                                            >  1:558182/1‑40 (MQ=255)
                                                                                                          aCGCTCACCGCAGTTCACGGGGAAATAACATGATGATa                                                                                                                              >  1:359916/1‑38 (MQ=255)
                                                                                                          aCGCTCACCGCAGTTCACGGGGAAATAACATGATGATa                                                                                                                              <  2:359916/38‑1 (MQ=255)
                                                                                                              tCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCgg                                                                                              >  1:1297172/1‑66 (MQ=255)
                                                                                                              tCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCgg                                                                                              <  2:1297172/66‑1 (MQ=255)
                                                                                                               cACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAGTTGCTGCCAGTCATGGACACCAGGttt                                         <  1:36587/118‑1 (MQ=255)
                                                                                                               cACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAGTTGCTGCCAGTCATGGACACCAGGttt                                         >  2:36587/1‑118 (MQ=255)
                                                                                                                        tCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAGTTGCTGCCAGTCATGGACACCAGGTTTTACTGTATGACATTTCTGCTGAAGCGCTGa           <  1:938107/139‑1 (MQ=255)
                                                                                                                        tCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAGTTGCTGCCAGTCATGGACACCAGGTTTTACTGTATGACATTTCTGCTGAAGCGCTGa           >  2:1340566/1‑139 (MQ=255)
                                                                                                                              ggAAATAACATGATGATAAATGTGCAAACTGTGGCGGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAGTtgctgc                                                             <  2:1410502/83‑1 (MQ=255)
                                                                                                                              ggAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAGTtgctgc                                                             >  1:1410502/1‑83 (MQ=255)
                                                                                                                              ggAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAGTTGCTGCCAGTCATGGACACCAGGTTTTACTGTATGACATTTCTGCTGAAGCGCTGACCcgcg     <  2:1552420/139‑1 (MQ=255)
                                                                                                                               gAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCAt                                                                            <  1:695340/67‑1 (MQ=255)
                                                                                                                               gAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCAt                                                                            >  2:695340/1‑67 (MQ=255)
                                                                                                                                 aaTAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAGTTGCTGCCAGTCATGGACAcc                                               >  1:746702/1‑94 (MQ=255)
                                                                                                                                 aaTAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAGTTGCTGCCAGTCATGGACAcc                                               <  2:746702/94‑1 (MQ=255)
                                                                                                                                 aaTAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAGTTGCTGCCAGTCATGGACACCAGGTTTTACTGTATGACATTTCTGCTGAAGCGCTGACCCGCGCaa  >  1:217582/1‑139 (MQ=255)
                                                                                                                                     |                                                                                                                                      
CTGAAACCAATACGCTCGATACGCAACTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAGCATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAGTTGCTGCCAGTCATGGACACCAGGTTTTACTGTATGACATTTCTGCTGAAGCGCTGACCCGCGCAA  >  NC_000913/1458919‑1459186

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