Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,458,617 A→G A118A (GCA→GCG 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,458,6170AG100.0% 187.7 / NA 50A118A (GCA→GCGpaaG1,2‑epoxyphenylacetyl‑CoA isomerase, oxepin‑CoA‑forming
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (25/25);  total (25/25)

GCCCCGCACCGGATTTAGGGATGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCACTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCTTAATACCCGATTGCGGTGGAACCTGGTTACTGCCACGCGTTGCCGGACGAGCGCGCGCCATGGGGC  >  NC_000913/1458487‑1458750
                                                                                                                                  |                                                                                                                                     
gCCCCGCACCGGATTTAGGGATGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCCGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTgggggg                                                                                                                               <  1:508160/139‑1 (MQ=255)
 ccccGCACCGGATTTAGGGATGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCCGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGc                                                                                                                              <  1:1386644/139‑1 (MQ=255)
 ccccGCACCGGATTTAGGGATGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCCGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGc                                                                                                                              <  1:330680/139‑1 (MQ=255)
   ccGCACCGGATTTAGGGATGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCCGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGa                                                                                                                            <  1:645418/139‑1 (MQ=255)
   ccGCACCGGATTTAGGGATGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCCGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGa                                                                                                                            >  1:1404385/1‑139 (MQ=255)
      cACCGGATTTAGGGATGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCCGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACAt                                                                                                                         >  2:826428/1‑139 (MQ=255)
                     tGTCAGTTGAACGTTTCTATAACCAACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCCGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCg                                                                                                          <  2:1345685/139‑1 (MQ=255)
                      gTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCCGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGa                                                                                                                            <  1:1351460/120‑1 (MQ=255)
                      gTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCCGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGa                                                                                                                            >  2:1351460/1‑120 (MQ=255)
                       tCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCCGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGtt                                                                                                        <  1:1414070/139‑1 (MQ=255)
                                     cTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCCGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTc                                                                                          <  1:900183/139‑1 (MQ=255)
                                      tataACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCCGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCa                                                                                         <  1:676189/139‑1 (MQ=255)
                                             cACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCCGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGAc                                                                                                                           >  2:93725/1‑98 (MQ=255)
                                             cACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCCGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGAc                                                                                                                           <  1:93725/98‑1 (MQ=255)
                                              aCTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCCGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGttt                                                                                 >  2:285001/1‑139 (MQ=255)
                                               cTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCCGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTa                                                                                >  2:1450147/1‑139 (MQ=255)
                                               cTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCCGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTa                                                                                >  1:661423/1‑139 (MQ=255)
                                                tGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCCGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCg                                                                                    <  1:996766/134‑1 (MQ=255)
                                                tGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCCGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCg                                                                                    >  2:996766/1‑134 (MQ=255)
                                                tGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCCGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAg                                                                               <  2:638620/139‑1 (MQ=255)
                                                 ggTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCCGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGt                                                                              <  2:1404385/139‑1 (MQ=255)
                                                    aCGTCGCCTGGCAAAACTGCCAAAACCGGTGATCCGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAg                                                                           <  2:661423/139‑1 (MQ=255)
                                                     cGTCGCCTGGCAAAACTGCCAAAACCGGTGATCCGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCg                                                                                    >  2:677489/1‑129 (MQ=255)
                                                     cGTCGCCTGGCAAAACTGCCAAAACCGGTGATCCGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCg                                                                                    <  1:677489/129‑1 (MQ=255)
                                                      gTCGCCTGGCAAAACTGCCAAAACCGGTGATCCGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGtt                                                                         >  1:672697/1‑139 (MQ=255)
                                                       tCGCCTGGCAAAACTGCCAAAACCGGTGATCCGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTa                                                                        <  2:10568/139‑1 (MQ=255)
                                                       tCGCCTGGCAAAACTGCCAAAACCGGTGATCCGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTa                                                                        >  2:942092/1‑139 (MQ=255)
                                                          cgtGGCAAAACTGCCAAAACCGGTGATCCGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTACCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGc                                                                     >  2:904706/3‑139 (MQ=255)
                                                          ccTGGCAAAACTGCCAAAACCGGTGATCCGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGGAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTa                                                                                <  2:1029397/128‑1 (MQ=255)
                                                          ccTGGCAAAACTGCCAAAACCGGTGATCCGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGa                                                                                                                            >  1:991643/1‑84 (MQ=255)
                                                          ccTGGCAAAACTGCCAAAACCGGTGATCCGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGa                                                                                                                            <  2:991643/84‑1 (MQ=255)
                                                          ccTGGCAAAACTGCCAAAACCGGTGATCCGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTa                                                                                >  1:1029397/1‑128 (MQ=255)
                                                          ccTGGCAAAACTGCCAAAACCGGTGATCCGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGt                                                                          >  2:679594/1‑134 (MQ=255)
                                                          ccTGGCAAAACTGCCAAAACCGGTGATCCGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGt                                                                          <  1:679594/134‑1 (MQ=255)
                                                             ggCAAAACTGCCAAAACCGGTGATCCGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCTTa                                                                  >  1:30206/1‑139 (MQ=255)
                                                                 aaaCTGCCAAAACCGGTGATCCGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCTTAATAc                                                              <  2:1444285/139‑1 (MQ=255)
                                                                             ccGGTGATCCGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCTTAATACCCGATTGCGGTg                                                  <  1:871380/139‑1 (MQ=255)
                                                                                        tgCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCaaaa                                                                                                <  2:1306613/82‑1 (MQ=255)
                                                                                        tgCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCaaaa                                                                                                >  1:1306613/1‑82 (MQ=255)
                                                                                        tgCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCTTAATACCCGATTGCGGTGGAACCTGGTTa                                       >  1:348717/1‑139 (MQ=255)
                                                                                             tCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCa                                                                                         <  2:231224/84‑1 (MQ=255)
                                                                                             tCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCa                                                                                         >  1:231224/1‑84 (MQ=255)
                                                                                               aaTGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCAACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCTTAATACCCGATTGCGGTGGAACCTGGTTACTGCCAc                                >  2:929908/1‑139 (MQ=255)
                                                                                                  ggcgtggcgGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGAc                                                                                                                           >  1:1021954/1‑45 (MQ=255)
                                                                                                  ggcgtggcgGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGAc                                                                                                                           <  2:1021954/45‑1 (MQ=255)
                                                                                                       ggcggcAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCTTAATACCCGATTGCGGTGGAACCTGGTTACTGCCACGCGTTGcc                        >  2:1465833/1‑139 (MQ=255)
                                                                                                         cggcAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCTTAATACCCGATTGCGGTGGAACCTGGTTACTGCCACGCGTTGCCgg                      >  2:1184853/1‑139 (MQ=255)
                                                                                                             aggggcaggCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCTTAATACCCGATTGCGGTGGAACCTGGTTACTGCCACGCGTTGCCGGACGa                  >  2:1261774/1‑139 (MQ=255)
                                                                                                               gggcaggCGCAACACTGGCGCTGGGGGGTGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCTTAATACCCGATTGCGGTGGAACCTGGTTACTGCCACGCGTTGCCGGACGAgc                <  1:904706/139‑1 (MQ=255)
                                                                                                                     gcgcAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCTTAATACCCGATTGCGGTGGAACCTGGTTACTGCCACGCGTTGCCGGACGAgcgcgcgc          >  1:904179/1‑139 (MQ=255)
                                                                                                                             ctggcgctggGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCTTAATACCCGATTGCGGTGGAACCTGGTTACTGCCACGCGTTGCCGGACGAGCGCGCGCCATGGGGc  <  2:1157529/139‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                     
GCCCCGCACCGGATTTAGGGATGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCACTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCTTAATACCCGATTGCGGTGGAACCTGGTTACTGCCACGCGTTGCCGGACGAGCGCGCGCCATGGGGC  >  NC_000913/1458487‑1458750

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