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% 148.2 / NA 46*263* (TAG→TAA) paaG1,2‑epoxyphenylacetyl‑CoA isomerase, oxepin‑CoA‑forming
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (21/25);  total (21/25)

CGCTGAAACCAATACGCTCGATACGCAACTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAGCATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAGTTGCTGCCAGTCATGGACACCAGGTTTTACTGTATGACATTTCTGCTGAAGCGCTGACCCGCG  >  NC_000913/1458917‑1459183
                                                                                                                                       |                                                                                                                                   
cgcTGAAACCAATACGCTCGATACGCAACTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACat                                                                                                                                  >  1:1305315/1‑139 (MQ=255)
cgcTGAAACCAATACGCTCGATACGCAACTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACat                                                                                                                                  <  2:1434205/139‑1 (MQ=255)
cgcTGAAACCAATACGCTCGATACGCAACTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACat                                                                                                                                  >  2:609082/1‑139 (MQ=255)
     aaaCCAATACGCTCGATACGCAACTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACAtgatga                                                                                                                             >  1:721726/1‑139 (MQ=255)
       aCCAATACGCTCGATACGCAACTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATa                                                                                                                           <  2:800382/139‑1 (MQ=255)
                   gATACGCAACTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACt                                                                                                               <  1:609082/139‑1 (MQ=255)
                           aCTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGt                                                                                                       >  1:131058/1‑139 (MQ=255)
                               gATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATaa                                                                                                                          <  1:471373/116‑1 (MQ=255)
                               gATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATaa                                                                                                                          >  2:471373/1‑116 (MQ=255)
                                  cTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTggg                                                                                                <  1:238725/139‑1 (MQ=255)
                                       aCGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCgg                                                                                           <  1:798501/139‑1 (MQ=255)
                                                                     cgcAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATggggg                                                                                <  2:715082/120‑1 (MQ=255)
                                                                     cgcAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATggggg                                                                                >  1:715082/1‑120 (MQ=255)
                                                                        aGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAGTtgctgc                                                          <  2:721726/139‑1 (MQ=255)
                                                                          aTTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAGTTGCTGCCa                                                        >  1:57603/1‑139 (MQ=255)
                                                                            tATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAGTTGCTGCCAGt                                                      <  1:1128308/139‑1 (MQ=255)
                                                                                    aGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAAtgtg                                                                                                                     >  1:266871/1‑68 (MQ=255)
                                                                                    aGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAAtgtg                                                                                                                     <  2:638419/68‑1 (MQ=255)
                                                                                    aGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAAtgtg                                                                                                                     <  2:266871/68‑1 (MQ=255)
                                                                                    aGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAAtgtg                                                                                                                     >  1:638419/1‑68 (MQ=255)
                                                                                    aGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAGTTGCTGCCAGTCATGGaca                                              <  2:131058/139‑1 (MQ=255)
                                                                                    aGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCAGCACCATGGGGGCAGGCATTGCTGAAGTTGCTGCCAGTCATGGaca                                              <  1:922161/139‑1 (MQ=255)
                                                                                      gtgtCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAGTTGCTGCCAGTCATGGACAcc                                            <  2:389495/139‑1 (MQ=255)
                                                                                        gtCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATg                                                                                    <  2:690222/97‑1 (MQ=255)
                                                                                        gtCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATg                                                                                    >  1:690222/1‑97 (MQ=255)
                                                                                        gtCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAGTTGCTGCCAGTCATGGACACCAg                                          >  1:103777/1‑139 (MQ=255)
                                                                                         tCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAGTTGCTGCCAGTCATGGaca                                              >  2:1514112/1‑134 (MQ=255)
                                                                                         tCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAGTTGCTGCCAGTCATGGaca                                              <  1:1514112/134‑1 (MQ=255)
                                                                                                   ccTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCAtt                                                                        >  2:149359/1‑98 (MQ=255)
                                                                                                   ccTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCAtt                                                                        <  1:149359/98‑1 (MQ=255)
                                                                                                       gCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAGTTGCTGCCAGTCATGGaca                                              >  2:987294/1‑120 (MQ=255)
                                                                                                       gCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAGTTGCTGCCAGTCATGGaca                                              <  1:987294/120‑1 (MQ=255)
                                                                                                            aCGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGggcagg                                                                            <  2:221048/85‑1 (MQ=255)
                                                                                                            aCGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGggcagg                                                                            >  1:221048/1‑85 (MQ=255)
                                                                                                            aCGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAGTTGCTGCCAGTCATGGACACCAGGTTTTACTGTATGACATTTc                      >  2:755526/1‑139 (MQ=255)
                                                                                                                 cACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAGTTGCTGCCAGTCATGGACACCAGGTTTTACTGTATGACATTTctgctg                 <  1:755526/139‑1 (MQ=255)
                                                                                                                  aCCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCAt                                                                         >  2:592486/1‑82 (MQ=255)
                                                                                                                  aCCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCAt                                                                         <  1:592486/82‑1 (MQ=255)
                                                                                                                         ttCACGGGAAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAGTTGCTGCCAGTCATGGACACCAGGTTTTACTGTATGACATTTCTGCTGAAGCGCTg         >  1:622138/1‑139 (MQ=255)
                                                                                                                              ggggAAATAACATGATGATAAATGTGCAAACTGTGGCa                                                                                                         <  1:1528902/38‑1 (MQ=255)
                                                                                                                              ggggAAATAACATGATGATAAATGTGCAAACTGTGGCa                                                                                                         >  2:1528902/1‑38 (MQ=255)
                                                                                                                              ggggAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGc                                                                                          <  2:113891/53‑1 (MQ=255)
                                                                                                                              ggggAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGc                                                                                          >  1:113891/1‑53 (MQ=255)
                                                                                                                              ggggAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAGTTGCTGCCAGTCATGGACACCAgg                                         <  1:274154/102‑1 (MQ=255)
                                                                                                                              ggggAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAGTTGCTGCCAGTCATGGACACCAgg                                         >  2:274154/1‑102 (MQ=255)
                                                                                                                                gaaaaTAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAGTTGCTGCCAGTCATGGACACCAGGTTTTACTGTATGACATTTCTGCTGAAGCGCTGACCcgcg  <  2:622138/137‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                   
CGCTGAAACCAATACGCTCGATACGCAACTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAGCATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAGTTGCTGCCAGTCATGGACACCAGGTTTTACTGTATGACATTTCTGCTGAAGCGCTGACCCGCG  >  NC_000913/1458917‑1459183

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