Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_pae_pgi 1,453,830 A→G intergenic (‑188/‑97) paaZ ← / → paaA oxepin‑CoA hydrolase and 3‑oxo‑5,6‑dehydrosuberyl‑CoA semialdehyde dehydrogenase/ring 1,2‑phenylacetyl‑CoA epoxidase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_pgi1,453,8300AG98.1% 162.9 / NA 52intergenic (‑188/‑97)paaZ/paaAoxepin‑CoA hydrolase and 3‑oxo‑5,6‑dehydrosuberyl‑CoA semialdehyde dehydrogenase/ring 1,2‑phenylacetyl‑CoA epoxidase subunit
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (20/31);  total (21/31)
Fisher's exact test for biased strand distribution p-value = 4.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATCGCGAATCATAAAAAGCAAAGGATCTTTTAACGAAATGTTAACTATGCGATCTGTATAGCAACTGCGGAAAACATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAATTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCT  >  NC_000913_3_pae_pgi/1453685‑1453974
                                                                                                                                                 |                                                                                                                                                
aaTCGCGAATCATAAAAAGCAAAGGATCTTTTAAAGAAATGTTAACTATGCGATCTGTATAGCAACTGCGGAAAACATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGttt                                                                                                                                               <  2:534998/149‑1 (MQ=255)
  tCGCGAATCATAAAAAGCAAAGGATCTTTTAACGAAATGTTAACTATGCGATCTGTATAGCAACTGCGGAAAACATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTtgt                                                                                                                                             <  2:696331/149‑1 (MQ=255)
     cgAATCATAAAAAGCAAAGGATCTTTTAACGAAATGTTAACTATGCGATCTGTATAGCAACTGCGGAAAACATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGAt                                                                                                                                          >  1:385684/1‑149 (MQ=255)
         tCATAAAAAGCAAAGGATCTTTTAACGAAATGTTAACTATGCGATCTGTATAGCAACTGCGGAAAACATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTa                                                                                                                                      >  1:225123/1‑149 (MQ=255)
           aTAAAAAGCAAAGGATCTTTTAACGAAATGTTAACTATGCGATCTGTATAGCAACTGCGGAAAACATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTACt                                                                                                                                    <  2:1190139/149‑1 (MQ=255)
               aaaGCAAAGGATCTTTTAACGAAATGTTAACTATGCGATCTGTATAGCAACTGCGGAAAACATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTACTTAAc                                                                                                                                >  1:383274/1‑149 (MQ=255)
                   cAAAGGATCTTTTAACGAAATGTTAACTATGCGATCTGTATAGCAACTGCGGAAAACATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTACTTAACTAtt                                                                                                                            >  1:1197421/1‑149 (MQ=255)
                   cAAAGGATCTTGTAACGAAATGTTAAAGATGCGATCTGTATAGCAACTGCGGAAAACATTAATGCACTGATAAAAAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTACTTAACTAtt                                                                                                                            <  2:471775/149‑1 (MQ=255)
                     aaGGATCTTTTAACGAAATGTTAACTATGCGATCTGTATAGCAACTGCGGAAAACATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAACTCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTACTTAACTATtgt                                                                                                                          <  1:1173495/149‑1 (MQ=255)
                     aaGGATCTTTTAACGAAATGTTAACTATGCGATCTGTATAGCAACTGCGGAAAACATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAATTTGTGATTTTACTTAACTATtgt                                                                                                                          >  2:748139/1‑149 (MQ=255)
                       ggATCTTTTAACGAAATGTTAACTATGCGATCTGTATAGCAACTGCGGAAAACATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTACTTAACTATTgtgt                                                                                                                        >  1:1217879/1‑149 (MQ=255)
                               tAACGAAATGTTAACTATGCGATCTGTATAGCAACTGCGGAAAACATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTACTTAACTATTGTGTAACTTTCa                                                                                                                <  1:1105727/149‑1 (MQ=255)
                                aaCGAAATGTTGAGTAGGCGGTCTGTATAGCAACTGCGGAAAACATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTACTTAACTATTGTGTAACTTTCAt                                                                                                               <  2:517005/149‑1 (MQ=255)
                                  cGAAATGTTAACTATGCGATCTGTATAGCAACTGCGGAAAACATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTACTTAACTATTGTGTAACTTTCATaa                                                                                                             <  2:225353/149‑1 (MQ=255)
                                  cGAAAAGTTAAATATGCGATCTGTATAGCAACTGCGAAAAACAGTAATGCACTGATAAATAATGATTTATAAAAATAGTGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTACTTAACTATTGTGTAACTTTCATaa                                                                                                             <  2:134286/149‑1 (MQ=255)
                                   gAAATGTTAACTATGCGATCTGTATAGCAACTGCGGAAAACATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAAACCGTCACAGTTAAAACATCCAAAGATCGTGATTTTACTGAACTATTGTGTAACTTTCATaaa                                                                                                            <  2:779470/149‑1 (MQ=255)
                                        gTTAACTATGCAATCTGTAGAGCAACTGCTGAAAACATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGAGTTCACTTAACTATTGTGTAACTTTCATAAAACAAt                                                                                                       <  1:825496/149‑1 (MQ=255)
                                         ttAACTATGCGATCTGTATAGCAACTGCGGAAAACATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAAtg                                                                                                      >  2:835621/1‑149 (MQ=255)
                                           aaCTATGCGATCTGTATAGCAACTGCGGAAAACATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAAggtt                                                                                                    >  2:588032/1‑148 (MQ=255)
                                             cTATGCGATCTGTATAGCAACTGCGGAAAACATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTACTTTACTATTGTGTAACTTTCATAAAACAATGTGAt                                                                                                  <  1:318330/149‑1 (MQ=255)
                                              tATGCGATCTGTATAGCAACTGCGGAAAACATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGAtg                                                                                                 >  2:427117/1‑148 (MQ=255)
                                                     tCTGTATAGCAACTGCGGAAAACATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGttt                                                                                          >  2:971181/1‑149 (MQ=255)
                                                       tGTATAGCAACTGCGGAAAACATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTACTTaa                                                                                                                                 >  1:168889/1‑108 (MQ=255)
                                                       tGTATAGCAACTGCGGAAAACATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTACTTaa                                                                                                                                 <  2:168889/108‑1 (MQ=255)
                                                         tataGCAACTGCGGAAAACATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTACATAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTttaa                                                                                      >  1:575126/1‑149 (MQ=255)
                                                               aaCTGCGGAAAACATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTaattaatt                                                                                <  2:495532/149‑1 (MQ=255)
                                                                aCTGCGGAAAACATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACGGTTCAAACATACAAAGTTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTc                                                                               <  1:775385/149‑1 (MQ=255)
                                                                  tGCGGAAAACATTAATGCACTTAAAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCAc                                                                             <  2:717552/149‑1 (MQ=255)
                                                                  tGCGGAAAACATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCAc                                                                             <  2:385684/149‑1 (MQ=255)
                                                                    cGGAAAACATGAACGCACTGATAAATAATGATGTATAAAAATAGGGTTCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCTCGa                                                                           <  2:506806/149‑1 (MQ=255)
                                                                          aCATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTg                                                                     <  2:383274/149‑1 (MQ=255)
                                                                             ttAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGaa                                                                  >  2:747205/1‑149 (MQ=255)
                                                                              tAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTACTTAATTATTTTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAAt                                                                 <  2:972089/149‑1 (MQ=255)
                                                                                   aaCTTATAAAAAATGATTTATAAAAATAGGGTGCGAAATCCGTTACAGTTCAAAAATACAAAGTTTGTGATTTTAGTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTaa                                                            <  2:647363/148‑1 (MQ=255)
                                                                                       gATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGt                                                        <  2:351290/149‑1 (MQ=255)
                                                                                          aaataaTGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGtgat                                                     <  2:414774/149‑1 (MQ=255)
                                                                                             taacGATGTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATTCTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGAc                                                  <  2:1168337/149‑1 (MQ=255)
                                                                                                  aTTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGAc                                             >  1:1199976/1‑149 (MQ=255)
                                                                                                   tttATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGCCGTGAcc                                            >  1:128458/1‑149 (MQ=255)
                                                                                                   tttATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGAcc                                            <  2:59465/149‑1 (MQ=255)
                                                                                                    gtaaaaaacTAGTGTCGGAAATCCATCAAATATCAAACATACAAAGTTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAAGGTGATTCCTGTGTTTAATTAAGTCACGAAAACTGGAATCGTAAAGGGGATGACGTGAccc                                           <  2:959583/145‑1 (MQ=255)
                                                                                                      ataAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCaa                                         <  1:747205/149‑1 (MQ=255)
                                                                                                                  tGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTATTTAATTAATTCACGAAAACTGCAATCGTAAAGGTGATGACGTAACCCAAGAAGAACACtct                             >  2:1285506/1‑147 (MQ=255)
                                                                                                                    cGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGa                           >  2:626253/1‑149 (MQ=255)
                                                                                                                     gAAATCCGTCACATTTCAAACATACAAAGTTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAg                          <  2:899232/149‑1 (MQ=255)
                                                                                                                     gAAATCCGTCACAGTTCAAAAATACAAAGTTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAg                          <  2:501671/149‑1 (MQ=255)
                                                                                                                            gTCACAGTTCAAACATACAAAGATTGTGAGTTTACTGAACTATTGTGTAACTTGCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTCAGCAACGGa                   >  2:1134516/1‑149 (MQ=255)
                                                                                                                             tCACAGTTCAAACATACAAAGTTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGGGTTTTTAATTAATTCACGAAAACTGGAAACGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGAt                  >  1:864270/1‑149 (MQ=255)
                                                                                                                                   ttCAAACATACAAAGTTTGTGATTTTACTTAACTATTGTGTAACTTTCATTAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCa            <  1:835621/149‑1 (MQ=255)
                                                                                                                                   ttCAAACATACAAAGTTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGAGGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCa            >  1:245351/1‑149 (MQ=255)
                                                                                                                                         cATACAAAGTTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAAc      <  1:748139/149‑1 (MQ=255)
                                                                                                                                             cAAAGTTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCt  >  2:157140/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
AATCGCGAATCATAAAAAGCAAAGGATCTTTTAACGAAATGTTAACTATGCGATCTGTATAGCAACTGCGGAAAACATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAATTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCT  >  NC_000913_3_pae_pgi/1453685‑1453974

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑