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% 143.1 / NA 39A118A (GCA→GCGpaaG1,2‑epoxyphenylacetyl‑CoA isomerase, oxepin‑CoA‑forming
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (20/19);  total (20/19)

CACAGGCCCCGCACCGGATTTAGGGATGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCACTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCTTAATACCCGATTGCGGTGGAACCTGGTTACTGCCACGCGTTGCCGGACGAGCGCGCGCCAT  >  NC_000913/1458482‑1458745
                                                                                                                                       |                                                                                                                                
cacaGGCCCCGCACCGGATTTAGGGATGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAAcactggcgctg                                                                                                                               <  2:733206/139‑1 (MQ=255)
    ggCCCCGCACCGGATTTAGGGATGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTggggg                                                                                                                           >  1:455777/1‑139 (MQ=255)
      ccccGCACCGGATTTAGGGATGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGc                                                                                                                         >  2:27830/1‑139 (MQ=255)
             ccGGATTTAGGGATGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCg                                                                                                                  <  2:522348/139‑1 (MQ=255)
                       ggATGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATtgctgc                                                                                                        <  1:27830/139‑1 (MQ=255)
                           gTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTgggggg                                                                                                                          >  2:17368/1‑117 (MQ=255)
                           gTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTgggggg                                                                                                                          <  1:17368/117‑1 (MQ=255)
                                          cTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTc                                                                                     <  2:324167/139‑1 (MQ=255)
                                             taACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCGGCCCGTTCAGCAAAATTCGTCATg                                                                                  >  1:449982/1‑139 (MQ=255)
                                                 ccACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGt                                                                              <  2:449982/139‑1 (MQ=255)
                                                    cTGGTACGTCGCCTGGCAAAACTGCTAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGGGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCg                                                                               <  1:544961/135‑1 (MQ=255)
                                                    cTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCg                                                                               >  2:544961/1‑135 (MQ=255)
                                                    cTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTa                                                                           >  2:422152/1‑139 (MQ=255)
                                                    cTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTa                                                                           >  1:78867/1‑139 (MQ=255)
                                                           gTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTggg                                                                                                                             >  1:674770/1‑82 (MQ=255)
                                                           gTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTggg                                                                                                                             <  2:674770/82‑1 (MQ=255)
                                                           gTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGtt                                                                    <  2:457239/139‑1 (MQ=255)
                                                           gTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGtt                                                                    <  1:166168/139‑1 (MQ=255)
                                                            tCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTa                                                                   >  1:508842/1‑139 (MQ=255)
                                                             cGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAg                                                                  >  1:104515/1‑139 (MQ=255)
                                                                                aaCCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCTTAATACCCGATTGCgg                                               <  2:565001/139‑1 (MQ=255)
                                                                                      tGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCTTAATACCCGATTGCGGTGGAAc                                         <  1:784420/139‑1 (MQ=255)
                                                                                             tgCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCTTAATACCCGATTGCGGTGGAACCTGGTTa                                  >  2:767750/1‑139 (MQ=255)
                                                                                             tgCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCTTAATACCCGATTGCGGTGGAACCTGGTTa                                  >  2:10672/1‑139 (MQ=255)
                                                                                               cAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATg                                                                                  <  1:310940/89‑1 (MQ=255)
                                                                                               cAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATg                                                                                  >  2:310940/1‑89 (MQ=255)
                                                                                                aGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCTTAATACCCGATTGCGGTGGAACCTGGTTACTg                               <  2:165443/139‑1 (MQ=255)
                                                                                                    aaTGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCTTAATACCCGATTg                                                  >  2:745732/1‑116 (MQ=255)
                                                                                                    aaTGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCTTAATACCCGATTg                                                  <  1:745732/116‑1 (MQ=255)
                                                                                                    aaTGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCTTAATACCCGATTGCGGTGGAACCTGGTTACTGCCAc                           >  1:454328/1‑139 (MQ=255)
                                                                                                    aaTGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCTTAATACCCGATTGCGGTGGAACCTGGTTACTGCCAc                           >  1:404482/1‑139 (MQ=255)
                                                                                                    aaTGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCTTAATACCCGATTGCGGTGGAACCTGGTTACTGCCAc                           >  2:281202/1‑139 (MQ=255)
                                                                                                     aTGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCTTAATACCCGATTGCGGTGGAACCTGGTTACTGCCAcg                          <  1:431818/139‑1 (MQ=255)
                                                                                                       ggcgtggcgGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCTTAATACCCGATTGCGGTGGAACCTGGTTACTGCCAcgcg                        >  1:382555/1‑139 (MQ=255)
                                                                                                              cggcAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGtt                                                                    <  1:215657/88‑1 (MQ=255)
                                                                                                              cggcAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGtt                                                                    >  2:215657/1‑88 (MQ=255)
                                                                                                                       caggCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCTTAATACCCGATTGCGGTGGAACCTGGTTACTGCCACGCGTTGCCGGACGAgcgcg        <  1:767750/139‑1 (MQ=255)
                                                                                                                          gcgcAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCTTAATACCCGATTGCGGTGGAACCTGGTTACTGCCACGCGTTGCCGGACGAgcgcgcgc     >  1:698907/1‑139 (MQ=255)
                                                                                                                             cAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCTTAATACCCGATTGCGGTGGAACCTGGTTACTGCCACGCGTTGCCGGACGAGCGCGCGCCAt  <  1:281202/139‑1 (MQ=255)
                                                                                                                             cAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCTTAATACCCGATTGCGGTGGAACCTGGTTACTGCCACGCGTTGCCGGACGAGCGCGCGCCAt  <  1:140038/139‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                
CACAGGCCCCGCACCGGATTTAGGGATGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCACTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCTTAATACCCGATTGCGGTGGAACCTGGTTACTGCCACGCGTTGCCGGACGAGCGCGCGCCAT  >  NC_000913/1458482‑1458745

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