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% 172.6 / NA 45A118A (GCA→GCG) paaG1,2‑epoxyphenylacetyl‑CoA isomerase, oxepin‑CoA‑forming
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (22/23);  total (22/23)

CAGGCCCCGCACCGGATTTAGGGATGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCACTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCTTAATACCCGATTGCGGTGGAACCTGGTTACTGCCACGCGTTGCCGGACGAGCGCG  >  NC_000913/1458484‑1458739
                                                                                                                                     |                                                                                                                          
caGGCCCCGCACCGGATTTAGGGATGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTggg                                                                                                                       <  2:64861/139‑1 (MQ=255)
     cccGCACCGGATTTAGGGATGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCg                                                                                                                  <  2:685793/139‑1 (MQ=255)
         cACCGGATTTAGGGATGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACAt                                                                                                              >  2:555176/1‑139 (MQ=255)
                                 aaCGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCg                                                                                                                  >  1:857725/1‑111 (MQ=255)
                                 aaCGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCg                                                                                                                  <  2:857725/111‑1 (MQ=255)
                                           taACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATg                                                                            <  1:303772/139‑1 (MQ=255)
                                           taACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATg                                                                            <  1:234378/139‑1 (MQ=255)
                                                  cTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATtgctgc                                                                                                  <  2:835673/110‑1 (MQ=255)
                                                  cTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATtgctgc                                                                                                  >  1:835673/1‑110 (MQ=255)
                                                   tGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGc                                                                                                                   <  1:877580/92‑1 (MQ=255)
                                                   tGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGc                                                                                                                   >  2:877580/1‑92 (MQ=255)
                                                     gTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCGGTGCAGTCGATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCGGGGGGGCGACATCGTGATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTa                                                                  >  1:775911/1‑139 (MQ=255)
                                                          tCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGAc                                                                                                                <  1:611437/88‑1 (MQ=255)
                                                          tCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAATCAATCGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGAc                                                                                                                >  2:611437/1‑88 (MQ=255)
                                                             ccTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGc                                                          >  1:274576/1‑139 (MQ=255)
                                                             ccTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCAGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGc                                                          >  2:610403/1‑139 (MQ=255)
                                                               tGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATgg                                                                           <  1:1779/120‑1 (MQ=255)
                                                               tGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATgg                                                                           >  2:1779/1‑120 (MQ=255)
                                                                   aaaaCTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCTTAATa                                                    <  2:274576/139‑1 (MQ=255)
                                                                           cAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGACAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCTTAATACCCGATTg                                            <  1:473841/139‑1 (MQ=255)
                                                                                                tCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGa                                                                                                                 <  1:89333/49‑1 (MQ=255)
                                                                                                tCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGa                                                                                                                 >  2:89333/1‑49 (MQ=255)
                                                                                                tCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCTTAATACCCGATTGCGGTGGAACCTGGTTACTGcc                       >  1:516126/1‑139 (MQ=255)
                                                                                                tCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCG‑TATTGCTGCCCGTTCAGCAAAATTCGTc                                                                               >  2:312327/1‑82 (MQ=255)
                                                                                                tCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCG‑TATTGCTGCCCGTTCAGCAAAATTCGTc                                                                               <  1:312327/82‑1 (MQ=255)
                                                                                                     ggcgtggcgGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGc                                                                                                                   >  1:889000/1‑42 (MQ=255)
                                                                                                     ggcgtggcgGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGc                                                                                                                   <  2:889000/42‑1 (MQ=255)
                                                                                                     ggcgtggcgGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGa                                                                                                                 <  1:893422/44‑1 (MQ=255)
                                                                                                     ggcgtggcgGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGa                                                                                                                 >  2:893422/1‑44 (MQ=255)
                                                                                                     ggcgtggcgGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCt                                                         <  2:366902/100‑1 (MQ=255)
                                                                                                     ggcgtggcgGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCt                                                         >  1:366902/1‑100 (MQ=255)
                                                                                                     ggcgtggcgGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCTTAATACCCGATTGCGGt                                        <  2:17014/117‑1 (MQ=255)
                                                                                                     ggcgtggcgGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCTTAATACCCGATTGCGGt                                        >  1:17014/1‑117 (MQ=255)
                                                                                                          ggcggcAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATgg                                                                           >  1:142654/1‑77 (MQ=255)
                                                                                                          ggcggcAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATgg                                                                           <  2:142654/77‑1 (MQ=255)
                                                                                                            cggcAGGGGCAGGCGCAACACTGGCGCTGGGGGGc                                                                                                                   <  2:438720/35‑1 (MQ=255)
                                                                                                            cggcAGGGGCAGGCGCAACACTGGCGCTGGGGGGc                                                                                                                   >  1:438720/1‑35 (MQ=255)
                                                                                                            cggcAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCTTAATACCCGATTGCGGTGGAACCTGGTTACTGCCACGCGTTGCCgg           >  2:176698/1‑139 (MQ=255)
                                                                                                                 ggggcaggCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCTTAATACCCGATTGCGGTGGAACCTGGTTACTGCCACGCGTTGCCGGACGAg      >  2:502780/1‑139 (MQ=255)
                                                                                                                   ggcaggCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCTTAATACCCGATTGCGGTGGAACCTGGTTACTGCCACGCGTTGCCGGACGAgcg    <  2:302236/139‑1 (MQ=255)
                                                                                                                     caggCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCTTAATACCCGATTGCGGTGGAACCTGGTTACTGCCACGCGTTGCCGGACGAgcgcg  <  1:176698/139‑1 (MQ=255)
                                                                                                                        gcgcAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCg                                                                                 <  2:134265/57‑1 (MQ=255)
                                                                                                                        gcgcAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCg                                                                                 >  1:134265/1‑57 (MQ=255)
                                                                                                                             acacTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGt                                                               >  2:66248/1‑70 (MQ=255)
                                                                                                                             acacTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGt                                                               <  1:66248/70‑1 (MQ=255)
                                                                                                                                     |                                                                                                                          
CAGGCCCCGCACCGGATTTAGGGATGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCACTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCTTAATACCCGATTGCGGTGGAACCTGGTTACTGCCACGCGTTGCCGGACGAGCGCG  >  NC_000913/1458484‑1458739

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