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,6170AG98.2% 197.7 / NA 56A118A (GCA→GCG) paaG1,2‑epoxyphenylacetyl‑CoA isomerase, oxepin‑CoA‑forming
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (26/29);  total (27/29)
Fisher's exact test for biased strand distribution p-value = 4.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

CCGCACCGGATTTAGGGATGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCACTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCTTAATACCCGATTGCGGTGGAACCTGGTTACTGCCACGCGTTGCCGGACGAGCGCGCG  >  NC_000913/1458490‑1458741
                                                                                                                               |                                                                                                                            
ccGCACCGGATTTAGGGATGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGa                                                                                                                   >  2:478585/1‑139 (MQ=255)
       ggATTTAGGGATGTCAGTTGAACGTTTCTATAACCCACTGGTACCTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGtt                                                                                                            >  2:1067946/1‑139 (MQ=255)
                       gTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAg                                                                                            <  1:1165170/139‑1 (MQ=255)
                                 tCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGt                                                                                  >  2:759037/1‑139 (MQ=255)
                                     taACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATg                                                                              <  2:256038/139‑1 (MQ=255)
                                      aaCCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATgg                                                                             <  2:533427/139‑1 (MQ=255)
                                      aaCCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATgg                                                                             <  1:759037/139‑1 (MQ=255)
                                        cccACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCg                                                                           >  2:269724/1‑139 (MQ=255)
                                        cccACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCg                                                                           <  1:269724/139‑1 (MQ=255)
                                           aCTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGACCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCAGGGGGGCGAAATCGTTATTGCTGCCCGTTCAGCAAAATTAGTCATGGCGtct                                                                        >  2:119300/1‑137 (MQ=255)
                                                       ccTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTgggggg                                                                                                                      <  2:501090/81‑1 (MQ=255)
                                                       ccTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTgggggg                                                                                                                      >  1:501090/1‑81 (MQ=255)
                                                       ccTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGt                                                                          >  2:37334/1‑125 (MQ=255)
                                                       ccTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGt                                                                          <  1:37334/125‑1 (MQ=255)
                                                       ccTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGc                                                            >  1:1136004/1‑139 (MQ=255)
                                                       ccTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGc                                                            >  1:1418995/1‑139 (MQ=255)
                                                       ccTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGc                                                            >  2:1106827/1‑139 (MQ=255)
                                                        cTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCt                                                           <  1:1067946/139‑1 (MQ=255)
                                                               aaCTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATtgct                                                                                                      <  2:874126/89‑1 (MQ=255)
                                                               aaCTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATtgct                                                                                                      >  1:874126/1‑89 (MQ=255)
                                                                         aCCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCTTAATACCCGATTGCGGt                                          <  2:810479/139‑1 (MQ=255)
                                                                             gTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCTTAATACCCGATTg                                              <  1:585458/131‑1 (MQ=255)
                                                                             gTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCTTAATACCCGATTg                                              >  2:585458/1‑131 (MQ=255)
                                                                              tGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACa                                                                                                                 >  1:663795/1‑63 (MQ=255)
                                                                              tGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACa                                                                                                                 <  2:663795/63‑1 (MQ=255)
                                                                                     tgCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCa                                                                                <  1:415086/89‑1 (MQ=255)
                                                                                     tgCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCa                                                                                >  2:415086/1‑89 (MQ=255)
                                                                                     tgCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCTTAATACCCGATTGCGGTGGAACCTGGTTa                              <  1:1106827/139‑1 (MQ=255)
                                                                                         gacAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCTTAATACCCGATTGCGGTGGAACCTGGTTACTGc                          <  2:1418995/137‑1 (MQ=255)
                                                                                         gTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCTTAATACCCGATTGCGGTGGAACCTGGTTACTGc                          <  2:1136004/139‑1 (MQ=255)
                                                                                          tCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAg                                                                                            <  1:112327/72‑1 (MQ=255)
                                                                                          tCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAg                                                                                            >  2:112327/1‑72 (MQ=255)
                                                                                            actggcgtggcgGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGt                                                                 <  2:979255/95‑1 (MQ=255)
                                                                                            aaTGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGAc                                                                                                                  >  1:409013/1‑48 (MQ=255)
                                                                                            aaTGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGAc                                                                                                                  <  2:409013/48‑1 (MQ=255)
                                                                                            aaTGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGt                                                                 >  1:979255/1‑97 (MQ=255)
                                                                                            aaTGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCTTAATa                                                      <  1:859962/108‑1 (MQ=255)
                                                                                            aaTGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCTTAATa                                                      >  2:859962/1‑108 (MQ=255)
                                                                                               ggcgtggcgGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATtgctgc                                                                                                    >  1:1311358/1‑59 (MQ=255)
                                                                                               ggcgtggcgGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATtgctgc                                                                                                    <  2:1311358/59‑1 (MQ=255)
                                                                                               ggcgtggcgGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCa                                                                                <  1:575980/79‑1 (MQ=255)
                                                                                               ggcgtggcgGCAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCa                                                                                >  2:575980/1‑79 (MQ=255)
                                                                                                    ggcggcAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGtt                                                                <  2:1283188/90‑1 (MQ=255)
                                                                                                    ggcggcAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGtt                                                                >  1:1283188/1‑90 (MQ=255)
                                                                                                    ggcggcAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCTTAATACCCGATTGCGGTGGAACCTGGTTACTGCCACGCGTTGcc               >  1:134991/1‑139 (MQ=255)
                                                                                                      cggcAGGGGCAGGCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCTTAATACCCGATTGCGGTGGAACCTGGTTACTGCCACGCGTTGCCgg             >  2:239774/1‑139 (MQ=255)
                                                                                                       ggcaggggcaggCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCTTAATACCCGATTGCGGTGGAACCTGGTTACTGCCACGCGTTGCCGGa            <  2:1359783/139‑1 (MQ=255)
                                                                                                            gggcaggCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGt                                                                 <  2:1007019/81‑1 (MQ=255)
                                                                                                            gggcaggCGCAACACTGGCACTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTCGTAAGt                                                                 >  1:1007019/1‑81 (MQ=255)
                                                                                                               caggCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCTTAATACCCGATTGCGGTGGAACCTGGTTACTGCCACGCGTTGCCGGACGAgcgcg    >  1:31600/1‑139 (MQ=255)
                                                                                                               caggCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCTTAATACCCGATTGCGGTGGAACCTGGTTACTGCCACGCGTTGCCGGACGAgcgcg    <  1:701088/139‑1 (MQ=255)
                                                                                                                 ggCGCAACACTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCTTAATACCCGATTGCGGTGGAACCTGGTTACTGCCACGCGTTGCCGGACGAgcgcgcg  <  1:478585/139‑1 (MQ=255)
                                                                                                                       acacTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGt                                                                                                >  1:139031/1‑39 (MQ=255)
                                                                                                                       acacTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGt                                                                                                <  2:139031/39‑1 (MQ=255)
                                                                                                                       acacTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCTTAATACCCGATTGCGGTgg                                        >  2:230804/1‑95 (MQ=255)
                                                                                                                       acacTGGCGCTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCTTAATACCCGATTGCGGTgg                                        <  1:230804/95‑1 (MQ=255)
                                                                                                                               |                                                                                                                            
CCGCACCGGATTTAGGGATGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGCAGGCGCAACACTGGCACTGGGGGGCGACATCGTTATTGCTGCCCGTTCAGCAAAATTCGTCATGGCGTTTAGTAAGTTAGGCTTAATACCCGATTGCGGTGGAACCTGGTTACTGCCACGCGTTGCCGGACGAGCGCGCG  >  NC_000913/1458490‑1458741

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: ‑