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,0520GA96.6% 177.0 / ‑1.9 59*263* (TAG→TAApaaG1,2‑epoxyphenylacetyl‑CoA isomerase, oxepin‑CoA‑forming
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (28/29);  total (29/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.66e-01

TACGCTCGATACGCAACTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAGCATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAGTTGCTGCCAGTCATGGACACCAGGTTTTACTGTATGACATTTCTGCTGAAGCGCTGACCCGCGC  >  NC_000913/1458929‑1459184
                                                                                                                           |                                                                                                                                    
tACGCTCGATACGCAACTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAAtgt                                                                                                                       <  1:959964/139‑1 (MQ=255)
   gCTCGATACGCAACTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCa                                                                                                                    <  2:1741270/139‑1 (MQ=255)
     tCGATACGCAACTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCaaa                                                                                                                  >  2:4714/1‑139 (MQ=255)
     tCGATACGCAACTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGATCACGGGGAAATAACATGATGATAAATGTGCaaa                                                                                                                  >  2:4667/1‑139 (MQ=255)
        aTACGCAACTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACtg                                                                                                               >  1:1491133/1‑139 (MQ=255)
           cGCAACTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTgg                                                                                                            <  2:1491133/139‑1 (MQ=255)
                     tCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTgg                                                                                                  >  1:13492/1‑139 (MQ=255)
                      cTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCCGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGa                                                                                                <  2:1050196/139‑1 (MQ=255)
                        ggAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAg                                                                                               >  2:352819/1‑139 (MQ=255)
                           aCGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCgg                                                                                            <  2:1688105/139‑1 (MQ=255)
                             gTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACTGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGAtt                                                                                                    >  2:1239673/1‑129 (MQ=255)
                             gTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGAtt                                                                                                    <  1:1239673/129‑1 (MQ=255)
                                aCTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCa                                                                                       <  2:437500/139‑1 (MQ=255)
                                 cTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCAt                                                                                      <  1:1479198/139‑1 (MQ=255)
                                       gCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGggc                                                                                <  2:13492/139‑1 (MQ=255)
                                        cGACTTACCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGggca                                                                               >  1:582502/1‑139 (MQ=255)
                                              gCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGc                                                                                              <  2:294123/118‑1 (MQ=255)
                                              gCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGc                                                                                              >  1:294123/1‑118 (MQ=255)
                                              gCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGggcagg                                                                             >  2:1518690/1‑135 (MQ=255)
                                              gCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGggcagg                                                                             <  1:1518690/135‑1 (MQ=255)
                                              gCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCAtt                                                                         >  1:212454/1‑139 (MQ=255)
                                                cggacggaGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACat                                                                                                                                   <  2:1696154/79‑1 (MQ=255)
                                                cggacggaGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACat                                                                                                                                   >  1:1696154/1‑79 (MQ=255)
                                                       aGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAGtt                                                                >  2:1243997/1‑139 (MQ=255)
                                                        gcgcAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCa                                                                                                          <  2:1414800/96‑1 (MQ=255)
                                                        gcgcAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCa                                                                                                          >  1:1414800/1‑96 (MQ=255)
                                                         cgcAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGAtt                                                                                                    <  2:1467885/101‑1 (MQ=255)
                                                         cgcAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGAtt                                                                                                    >  1:1467885/1‑101 (MQ=255)
                                                                      gAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAGTTGCTGCCAGTCATGGa                                                 >  2:1060989/1‑139 (MQ=255)
                                                                        aGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACat                                                                                                                                   >  2:1023279/1‑55 (MQ=255)
                                                                        aGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACat                                                                                                                                   <  1:1023279/55‑1 (MQ=255)
                                                                                 tGCGTTCGTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGATCGGCACCATGGGGGCAGGCATTGCTGAAGTTGCTGCCAGTCATGGACACCAGGtttt                                      <  1:455434/139‑1 (MQ=255)
                                                                                   cGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAGTTGCTGCCAGTCATGGACACCAGGttt                                       >  1:92262/1‑136 (MQ=255)
                                                                                   cGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAGTTGCTGCCAGTCATGGACACCAGGttt                                       <  2:92262/136‑1 (MQ=255)
                                                                                    gTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAGCATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCgg                                                                                            <  1:817177/82‑1 (MQ=255)
                                                                                    gTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAGCATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCgg                                                                                            >  2:817177/1‑82 (MQ=255)
                                                                                     ttCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAGTTGCTGCCAGTCATGGACACCAGGTTTTACTg                                  <  2:212454/139‑1 (MQ=255)
                                                                                     ttCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAGTTGCTGCCAGTCATGGACACCAGGTTTTACTg                                  >  2:241383/1‑139 (MQ=255)
                                                                                           gCTAAACGCTCACCGCAGTTCACGGGGAAATAACat                                                                                                                                   >  2:1750360/1‑36 (MQ=255)
                                                                                           gCTAAACGCTCACCGCAGTTCACGGGGAAATAACat                                                                                                                                   <  1:1750360/36‑1 (MQ=255)
                                                                                           gCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAGTTGCTGCCAGTCATGGACACCAGGttt                                       >  2:987962/1‑128 (MQ=255)
                                                                                           gCTAAACGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAGTTGCTGCCAGTCATGGACACCAGGttt                                       <  1:987962/128‑1 (MQ=255)
                                                                                                aCGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGggcaggc                                                                            >  1:1079372/1‑86 (MQ=255)
                                                                                                aCGCTCACCGCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGggcaggc                                                                            <  2:1079372/86‑1 (MQ=255)
                                                                                                         gCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAGTTGCTGCCAGTCATGGACACCAGGTTTTACTg                                  <  2:719954/119‑1 (MQ=255)
                                                                                                         gCAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAGTTGCTGCCAGTCATGGACACCAGGTTTTACTg                                  >  1:719954/1‑119 (MQ=255)
                                                                                                          cAGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAGTTGCTGCCAGTCATGGACACCAGGTTTTACTGTATGACATTTCTGCTGAAgcg             <  1:4714/139‑1 (MQ=255)
                                                                                                          cAGATCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAGTTGCTGCCAGTCATGGACACCAGGTTTTACTGTATGACATTTCTGCTGAAgcg             <  1:4667/139‑1 (MQ=255)
                                                                                                           aGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAGTTGCTGCCAGTCATGGaca                                               >  2:1552427/1‑104 (MQ=255)
                                                                                                           aGTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAGTTGCTGCCAGTCATGGaca                                               <  1:1552427/104‑1 (MQ=255)
                                                                                                            gTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAg                                                                                               >  1:700597/1‑55 (MQ=255)
                                                                                                            gTTCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAg                                                                                               <  2:700597/55‑1 (MQ=255)
                                                                                                             ttCACGGGGAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAGTTGCTGCCAGTCATGGACACCAGGTTTTACTGTATGACATTTCTGCTGAAGCGCTg          >  2:1510950/1‑139 (MQ=255)
                                                                                                                  ggggAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAGTTGCTGCCa                                                         <  2:1247506/87‑1 (MQ=255)
                                                                                                                  ggggAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAGTGGCTGCCa                                                         >  1:1247506/1‑87 (MQ=255)
                                                                                                                   gggAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAGTTGCTGCCAGTCATGGACACCAGGTTTTACTGTATGACATTTCTGCTGAAGCGCTGACCcgc    <  1:1510950/139‑1 (MQ=255)
                                                                                                                    ggAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAGTTGCTGCCAGTCATGGACAcc                                             >  2:1866276/1‑97 (MQ=255)
                                                                                                                    ggAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAGTTGCTGCCAGTCATGGACAcc                                             <  1:1866276/97‑1 (MQ=255)
                                                                                                                     gAAATAACATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAGTTGCTGCCAGTCATGGACACCAGGTTTTACTGTATGACATTTCTGCTGAAGCGCTGACCCgcgc  >  1:1105615/1‑139 (MQ=255)
                                                                                                                           |                                                                                                                                    
TACGCTCGATACGCAACTGGATCTGGAACGTGACTATCAGCGACTTGCCGGACGGAGCGCAGATTATCGTGAAGGTGTCAGTGCGTTCCTGGCTAAACGCTCACCGCAGTTCACGGGGAAATAGCATGATGATAAATGTGCAAACTGTGGCAGTGATTGGGAGCGGCACCATGGGGGCAGGCATTGCTGAAGTTGCTGCCAGTCATGGACACCAGGTTTTACTGTATGACATTTCTGCTGAAGCGCTGACCCGCGC  >  NC_000913/1458929‑1459184

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