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

CAGGCCCCGCACCGGATTTAGGGATGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCACTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCTTAATACCCGATTGCGGTGGAACCTGGTTACTGCCACGCGT  >  NC_000913/1458484‑1458724
                                                                                                                                     |                                                                                                            
caGGCCCCGCACCGGATTTAGGGATGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTggg                                                                                                        >  2:426401/1‑139 (MQ=255)
 aGGCCCCGCACCGGATTTAGGGATGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTgggg                                                                                                       >  2:159803/1‑139 (MQ=255)
      ccGCACCGGATTTAGGGATGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTgggggggcg                                                                                                  >  2:292695/1‑136 (MQ=255)
            cGGATTTAGGGATGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGGCGACATCg                                                                                             <  1:426401/139‑1 (MQ=255)
                    gggATGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGGCGACATCGTTATtgct                                                                                     <  1:234436/139‑1 (MQ=255)
                              ttGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGGCGACATCGTTATTGCTGCCCGTTCAg                                                                           <  1:159803/139‑1 (MQ=255)
                                   cGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCaaaa                                                                      >  1:92648/1‑139 (MQ=255)
                                                cACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGt                                                         >  1:119448/1‑139 (MQ=255)
                                                cACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGt                                                         >  1:112886/1‑139 (MQ=255)
                                                 aCTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGt                                                         <  2:381660/138‑1 (MQ=255)
                                                 aCTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTggggggg                                                                                                    <  1:400652/94‑2 (MQ=255)
                                                 aCTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTggggggg                                                                                                    >  2:400652/1‑93 (MQ=255)
                                                 aCTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGAAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGt                                                         >  1:381660/1‑138 (MQ=255)
                                                   tGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTa                                                      <  2:92648/139‑1 (MQ=255)
                                                             ccTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGGCGACATCGTTAtt                                                                                        >  2:455224/1‑95 (MQ=255)
                                                             ccTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGGCGACATCGTTAtt                                                                                        <  1:455224/95‑1 (MQ=255)
                                                             ccTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAgg                                            >  1:344932/1‑139 (MQ=255)
                                                             ccTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAgg                                            >  1:467917/1‑139 (MQ=255)
                                                             ccTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAgg                                            <  2:485381/139‑1 (MQ=255)
                                                             ccTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTa                                                      <  1:461471/128‑1 (MQ=255)
                                                             ccTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTa                                                      >  2:461471/1‑128 (MQ=255)
                                                                                         tgtgCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCTTAATACCCGATTGCGGTGGAACCTgg                >  2:261146/1‑139 (MQ=255)
                                                                                                   aTGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCTTAATACCCGATTGCGGTGGAACCTGGTTACTGCCAc      <  2:467917/139‑1 (MQ=255)
                                                                                                     ggcgtggcgGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCTTAATACCCGATTGCGGTGGAACCTGGTTACTGCCAcgc    >  2:412251/1‑139 (MQ=255)
                                                                                                       cgtggcgGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCTTAATACCCGATTGCGGTGGAACCTGGTTACTGCCACGCGt  <  1:261146/139‑1 (MQ=255)
                                                                                                                                     |                                                                                                            
CAGGCCCCGCACCGGATTTAGGGATGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCACTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCTTAATACCCGATTGCGGTGGAACCTGGTTACTGCCACGCGT  >  NC_000913/1458484‑1458724

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