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% 134.5 / NA 42*263* (TAG→TAApaaG1,2‑epoxyphenylacetyl‑CoA isomerase, oxepin‑CoA‑forming
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (19/23);  total (19/23)

AAACCAATACGCTCGATACGCAACTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAGCATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAGTTGCTGCCAGTCATGGACACCAGGTTTTACTGTATGACATTTCTGCTG  >  NC_000913/1458922‑1459168
                                                                                                                                  |                                                                                                                    
aaaCCAATACGCTCGATACGCAACTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACAtgatga                                                                                                              >  2:1063928/1‑139 (MQ=255)
              gATACGCAACTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACat                                                                                                                   >  2:1364318/1‑120 (MQ=255)
              gATACGCAACTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACat                                                                                                                   <  1:1364318/120‑1 (MQ=255)
              gATACGCAACTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACt                                                                                                >  2:1442406/1‑139 (MQ=255)
              gATACGCAACTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACt                                                                                                <  1:483398/139‑1 (MQ=255)
                tACGCAACTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACtgt                                                                                              >  2:30029/1‑139 (MQ=255)
                  cGCAACTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTgg                                                                                            <  1:1442406/139‑1 (MQ=255)
                    cAACTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCa                                                                                          <  2:1143020/139‑1 (MQ=255)
                          gATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGAtt                                                                                    >  1:144369/1‑139 (MQ=255)
                               ggAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAg                                                                               <  1:1063928/139‑1 (MQ=255)
                                         tATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATg                                                                     <  2:895960/139‑1 (MQ=255)
                                             aGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATggggg                                                                 <  1:30029/139‑1 (MQ=255)
                                                 aCTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGggcagg                                                             >  1:1130259/1‑139 (MQ=255)
                                                                        aTCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCa                                                                                          <  1:716607/87‑1 (MQ=255)
                                                                        aTCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCa                                                                                          >  2:716607/1‑87 (MQ=255)
                                                                        aTCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAGTTGCTGCCAGTc                                      >  1:759057/1‑139 (MQ=255)
                                                                            tGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTgg                                                                                            >  1:1184555/1‑81 (MQ=255)
                                                                            tGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTgg                                                                                            <  2:1184555/81‑1 (MQ=255)
                                                                               aGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAGTTGCTGCCAGTCATGGaca                               <  1:94661/139‑1 (MQ=255)
                                                                                ggTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGc                                                                           <  1:308334/94‑1 (MQ=255)
                                                                                ggTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGc                                                                           >  2:308334/1‑94 (MQ=255)
                                                                                   gtCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGc                                                                              <  2:999153/88‑1 (MQ=255)
                                                                                   gtCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGc                                                                              >  1:999153/1‑88 (MQ=255)
                                                                                    tCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAGTTGCTGCCAGTCATGGACACCAgg                          <  2:1130259/139‑1 (MQ=255)
                                                                                        tGCGTTCCTGGCCAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATa                                                                                                            <  1:1307497/53‑1 (MQ=255)
                                                                                        tGCGTTCCTGGCCAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATa                                                                                                            >  2:1307497/1‑53 (MQ=255)
                                                                                         gCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAGTTGCTGCCAGTCATGGACACCAGGTTTTa                     <  2:1326551/139‑1 (MQ=255)
                                                                                          cGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACat                                                                                                                   <  1:865976/44‑1 (MQ=255)
                                                                                          cGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACat                                                                                                                   >  2:865976/1‑44 (MQ=255)
                                                                                                tGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACtgt                                                                                              >  2:1015741/1‑59 (MQ=255)
                                                                                                tGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACtgt                                                                                              <  1:1015741/59‑1 (MQ=255)
                                                                                                  gCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCgg                                                                            <  2:1159878/75‑1 (MQ=255)
                                                                                                  gCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCgg                                                                            <  2:1093424/75‑1 (MQ=255)
                                                                                                  gCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCgg                                                                            >  1:1159878/1‑75 (MQ=255)
                                                                                                  gCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCgg                                                                            >  1:1093424/1‑75 (MQ=255)
                                                                                                       aCGCTCACCGCAGTTCACGGGGAAATAACATGATGATaa                                                                                                           >  1:754045/1‑39 (MQ=255)
                                                                                                       aCGCTCACCGCAGTTCACGGGGAAATAACATGATGATaa                                                                                                           <  2:754045/39‑1 (MQ=255)
                                                                                                            cACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAGTTGCTGCCAGTCATGGACACCAGGTTTTACTGTATGACATTTctgctg  <  2:582555/139‑1 (MQ=255)
                                                                                                                 cAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAGTTGCTGCCAGTCATGGACAcc                             <  1:840099/107‑1 (MQ=255)
                                                                                                                 cAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAGTTGCTGCCAGTCATGGACAcc                             >  2:840099/1‑107 (MQ=255)
                                                                                                                     tCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAg                                                  <  2:1047190/82‑1 (MQ=255)
                                                                                                                     tCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAg                                                  >  1:1047190/1‑82 (MQ=255)
                                                                                                                                  |                                                                                                                    
AAACCAATACGCTCGATACGCAACTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAGCATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAGTTGCTGCCAGTCATGGACACCAGGTTTTACTGTATGACATTTCTGCTG  >  NC_000913/1458922‑1459168

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