Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,454,774 A→G D16G (GAT→GGT)  paaI → hydroxyphenylacetyl‑CoA thioesterase PaaI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,454,7740AG98.1% 169.3 / ‑3.7 53D16G (GAT→GGT) paaIhydroxyphenylacetyl‑CoA thioesterase PaaI
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (18/34);  total (18/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.32e-01

GTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGTCACTGCACAAA  >  NZ_CP009273/1454689‑1454857
                                                                                     |                                                                                   
gTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGcc                                                                                 <  1:449282/90‑1 (MQ=255)
     ttGCTGCGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGAATGGGAACGGGGCCTgcgc                                                                            >  2:1315714/1‑90 (MQ=255)
          gcgcCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAg                                                                       >  2:37714/1‑90 (MQ=255)
          gcgcCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAg                                                                       >  1:1171616/1‑90 (MQ=255)
          gcgcCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAg                                                                       >  1:1499258/1‑90 (MQ=255)
                    gcgcTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCAt                                                             <  2:363313/90‑1 (MQ=255)
                      gcTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCg                                                           <  2:510293/90‑1 (MQ=255)
                          cTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACAt                                                       <  2:126803/90‑1 (MQ=255)
                           tGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACAtt                                                      <  1:1128129/90‑1 (MQ=255)
                                 gCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCa                                                >  1:1289177/1‑90 (MQ=255)
                                 gCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCa                                                >  2:244003/1‑90 (MQ=255)
                                 gCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCa                                                >  2:321234/1‑90 (MQ=255)
                                 gCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCa                                                >  1:185576/1‑90 (MQ=255)
                                   ggTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAAt                                              <  2:394777/90‑1 (MQ=255)
                                      tATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCg                                                           >  2:385235/1‑74 (MQ=255)
                                      tATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCg                                                           <  1:385235/74‑1 (MQ=255)
                                          aGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGaa                                       <  1:1074347/90‑1 (MQ=255)
                                             cATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGc                                    <  2:426682/90‑1 (MQ=255)
                                             cATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGc                                    <  2:655354/90‑1 (MQ=255)
                                             cATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGc                                    <  1:1002127/90‑1 (MQ=255)
                                             cATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGc                                    <  2:948552/90‑1 (MQ=255)
                                             cATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGc                                    <  2:1210453/90‑1 (MQ=255)
                                             cATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGc                                    <  1:686778/90‑1 (MQ=255)
                                              aTAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCt                                   >  1:1503172/1‑90 (MQ=255)
                                              aTAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCt                                   >  2:889206/1‑90 (MQ=255)
                                              aTAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCt                                   >  1:1312397/1‑90 (MQ=255)
                                               tAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCtt                                  <  1:903242/90‑1 (MQ=255)
                                               tAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCtt                                  <  1:612168/90‑1 (MQ=255)
                                               tAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCtt                                  <  2:757857/90‑1 (MQ=255)
                                                  ggCCTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGc                               <  1:301082/90‑1 (MQ=255)
                                                     cTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGt                            >  1:704300/1‑90 (MQ=255)
                                                     cTGGCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGt                            >  2:651703/1‑90 (MQ=255)
                                                       ggCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAg                          <  1:1288494/90‑1 (MQ=255)
                                                       ggCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAg                          >  2:1471460/1‑90 (MQ=255)
                                                       ggCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAg                          <  2:159197/90‑1 (MQ=255)
                                                       ggCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAg                          <  2:79742/90‑1 (MQ=255)
                                                       ggCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAg                          <  1:244003/90‑1 (MQ=255)
                                                        gCAAAATGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGt                         <  2:108308/90‑1 (MQ=255)
                                                            aaTGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGAcc                     >  1:1096202/1‑90 (MQ=255)
                                                             aTGCCCATGCTATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCa                    <  2:658118/90‑1 (MQ=255)
                                                             aTGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCa                    <  2:1579763/90‑1 (MQ=255)
                                                             aTGCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCa                    <  1:766273/90‑1 (MQ=255)
                                                               gCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATg                  <  2:439319/90‑1 (MQ=255)
                                                               gCCCATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATg                  <  2:1312397/90‑1 (MQ=255)
                                                                  cATGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGAcc               <  1:975073/90‑1 (MQ=255)
                                                                   aTGCAATGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCg              >  2:538181/1‑90 (MQ=255)
                                                                        aTGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGTCACt         <  2:498522/90‑1 (MQ=255)
                                                                        aTGTATGAGAACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGTCACt         <  1:68520/90‑1 (MQ=255)
                                                                               agaACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGTCACTGCACaaa  <  1:497041/90‑1 (MQ=255)
                                                                               agaACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGTCACTGCACaaa  <  1:567255/90‑1 (MQ=255)
                                                                                gaACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAAt                                              >  2:770731/1‑45 (MQ=255)
                                                                                gaACGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAAt                                              <  1:770731/45‑1 (MQ=255)
                                                                                  aCGGTGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGTc            <  2:1133814/77‑1 (MQ=255)
                                                                                     |                                                                                   
GTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGTCACTGCACAAA  >  NZ_CP009273/1454689‑1454857

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