Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,458,375 G→T 25.3% E194* (GAA→TAA)  paaX → transcriptional repressor of phenylacetic acid degradation paa operon, phenylacetyl‑CoA inducer

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,458,3750GT25.3% 31.0 / 3.6 19E194* (GAA→TAA) paaXtranscriptional repressor of phenylacetic acid degradation paa operon, phenylacetyl‑CoA inducer
Reads supporting (aligned to +/- strand):  ref base G (7/7);  new base T (2/3);  total (9/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.48e-01

GAAGGTTTAGATAAATCAACGCTGGCTGATGTCAAAAAGCAGTTGATCTGGCAAGGTTTTGGCGCACTGGCACCCAGCCTGATGGCATCGCCGTCGCAAAAACTGGCCGATGTACAGACACTTTTGCATGAAGCGGGTGTGGCGGATAACGTGATTTGTTTTGAAGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGCCAGAGTAGAAGAGTGCTGGCATTTAACTGAACAAAATGCCATGTACGAAACCTTTATTCAGTCATTCCGCCCGCTGGTGCCGCTTTTAAAAGAGGCGGCAGACGAGTTAACCCCGGAGCGGGCATTTCATATTCAGCTTTTACTGATCCATTTTTATCGCCGTGTCGTCCTTAAAGACCCATTGTTGCCGGAGGAGTTGCTTCCGGCACACTGGGC  >  CP009273/1458141‑1458562
                                                                                                                                                                                                                                          |                                                                                                                                                                                           
gcctgTTTAGAAAAAGCAACGCTTGCTGATCTCAAACAGCAGTTGGTCTGGCAAGGTTTTGGCGCACTGGTACCCAGCCTGATGGTATCGCCGTCGCAAAAACTGGCCGATGTACAGACACTTTTGAATGAAGCGGGTGTGGCGGATAACGTGATTTGTTTTGAAGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGCCAGAGTAGAAGAGGGCTGGCATTTAACTGAACAAAATGCCATGTACGAAACCTTTATTCAGTCATTccg                                                                                                                                                     <  2:82196/271‑1 (MQ=255)
                                                                     gCACCCAGCCTGATAGCATCGCCGTCGCAAAAACTGGCCGATGTACAGACACTTTTGCATGAAGCGGGTGTGGCGGATAACGTGATTTGTTTTGAAGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGCCAGAGTAGAAGAGTGCTGGCATTTAACTGAACAAAATGCCa                                                                                                                                                                                 <  1:373340/178‑1 (MQ=255)
                                                                                     cATCGCCGTCGCAAAAACTGGCCGATGTACAGACACTTTTGCATGAAGCGGGTGTGGCGGATAACGTGATTTGTTTTGAAGCGCAAATACCACTGGCGCTTTCTTGCGCAGCACTGCGTGCCAGAGTAGAAGAGTGCTGGCATTTAACTGAACAAAATGCCATGTACGAAACCTTTATTCAGTCATTCCGCCCGCTGGTGCCGCTTTTAAAAGAGGCGGCAGAc                                                                                                                   >  1:201067/1‑224 (MQ=255)
                                                                                         ggcttCGCAAAAGCTGGTCGATGGACAGACACTTTTGCATGACGCGGGTGTGGCGGATAACGTGATTTGTTTTGAATCGCAAATACCACTGGCGCTTTCTGGCGCAGCACTGCGGGGCAGAGTAGAAGAGTGCTGGCATTTAACTGAACAAAATGCCATGTACGAAACCTTTATTCAGTCATTCCGCCCGCTGGTGCCGCTTTTAAAAGAGGCGGCAGACGAGt                                                                                                               <  2:201061/220‑1 (MQ=255)
                                                                                                                           ttGCATGAAGCGGGTGTGGCGGATAACGTGATTTGTTTTGAAGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGCCAGAGTAGAAGAGTGCTGGCATTTAACTTAACAAAATGCCATGTACGAAACCTTTATTCAGTCATTCCGCCCGCTGGTGCCGCTTTTAAAAGAGGCGGCAGACGAGTTAACCCCGGAGCGGGCATTTCATATTCAGCTTTTACTGATCCATTTTTATCGCCGTGTCGTCCTTAAAGACCCATTGTTGCCgg                           >  1:89803/1‑274 (MQ=255)
                                                                                                                                              cGGTTAACGTGGTTTGTTTTTAAGCGCAAATAACACTGGCGCTTTCTCGCGCAGCACTGCGTGCCAGAGTAGAAGAGTGCTGGCATTTAACTTAACAAAATGCCATGTACGAAACCTTTATTCAGTCAATCCGCCCGCTGGTGCCGCTTTTAAAAGAGGCGGCAGACGAGTTAACCCCGGAGCGGGCATTTCATATTCAGCTTTTACTGATCCATTTTTATCGCCGTGTCGTCCTTAAAGACCCATTGTTGCCGGAGGAGTTGCTTCCGGcaca        <  2:89801/274‑1 (MQ=255)
                                                                                                                                                   aaCGTGATTTGTTTTAAAGCGCAAATAAAACTGGCGCTTTCTCGCGCAGAACTGCGGGCCAGAGTAGAAGAGTGCTGGCATTTGAATGAACAAAATGAAATGTACGAAACATTTATTCAGTCATTAAGCCCGCTGGTGCAGCTTGTAAACGAGGCGACAGACGAGTTAACACCGGAGCGGGCATTTCATATTCAGCAGTTACAACTCCACTTTTATCGCCGTGACGTCCTGAAAGCACCCTTGTTTCCGGAGGCGTTGGTTCCGGCACACtaaga  >  2:70224/1‑271 (MQ=255)
                                                                                                                                                                                aCTGGCGCTTTCTCGCGCAGCACTGCGTGCCAGAGTAGAAGAGTGCTGGCATTTAACTGAACAAAATGCCATGTACGAAACCTTTATTCAGTCATTCCGCCCGCTGGTGCCGCTTTTAAAAGAGGCGGCAGACGAGTTAACCCCGGAGCGGGCATTTCATATTCAGCTTTTACTGATCCATTTTTATCGcc                                                         >  1:216268/1‑191 (MQ=255)
                                                                                                                                                                                    gcgcTTTCTCGCGCAGCACTGCGTGCCAGAGTAGAAGAGTGCTGGCATTTAACTGAACAAAATGCCATGTACGAAACCTTTATTCAGTCATTCCGCCCGCTGGTGCCGCTTTTAAAAGAGGCGGCAGACGAGTTAACCCCGGAGCGGGCATTTCATATTCAGCTTTTACTGATCCATTTTTATCGCCgtgt                                                     <  2:216261/191‑1 (MQ=255)
                                                                                                                                                                                                    caCTGCGTGCCAGAGTAGAAGAGTGCTGGCATTTAACTTAACAAAATGCCATGTACGAAACCTTTATTCAGTCATTCCGCCCGCTgg                                                                                                                                             >  2:401371/1‑87 (MQ=255)
                                                                                                                                                                                                        gCGTGCCAGAGTAGAAGAGTGCTGGCATTTAACTTAACAAAATGCCATGTACGAAACCTTTATTCAGTCATTCCGCCCGCTGGTGcc                                                                                                                                         <  1:401384/87‑1 (MQ=255)
                                                                                                                                                                                                                 agTAGAAGAGTGCTGGCATTTAACTGAACAAAATGTCATGTACGAAACCTTTATTCAGTCATTCCGCCCGCTGGTGCAGCTTTTAAAAGaggagg                                                                                                                        >  1:61738/1‑95 (MQ=255)
                                                                                                                                                                                                                     gaagaGTGCTGGCATTTAACTTGACAAAATGCCATGTACGAAACCTGTATTCAGTCATTCCGCCCGCTGGTGCCGCTTTTAAAAGAGGCGGCAGa                                                                                                                    <  2:61736/95‑1 (MQ=255)
                                                                                                                                                                                                                      aagaGTGCTGGCATTTAACTGAACAAAATGCCATGTACGAAACCTTTa                                                                                                                                                                  >  1:227136/1‑48 (MQ=255)
                                                                                                                                                                                                                          gTGCTGGCATTTAACTGAACAAAATGCCATGTACGAAACCTTTATTCa                                                                                                                                                              <  2:227129/48‑1 (MQ=255)
                                                                                                                                                                                                                           tGCTGGCATTTAACTGAACAAAATGCCATGTACGAAACCTTTATTCAGTCATTCCGCCCGCTGGTGCCGCTTTTAAAAGAGGCGGCAGACGAGTTAACCCCGGAGCGGGCATTTCATATTCAGCTTTTACTGATCCATTTTTATCGCCGTGTCGTCCTTAAAGACCCATTGTTGCCGGAGGAGTTGCTTCCGGcaca        >  1:173257/1‑197 (MQ=255)
                                                                                                                                                                                                                               ggCATTTAACTGAACAAAATGCCATGTACGAAACCTTTATTCAGTCATTCCGCCCGCTGGTGCCGCTTTTAAAAGAGGCGGCAGACGAGTTAACCCCGGAGCGGGCATTTCATATTCAGCTTTTACTGATCCATTTTTATCGCCGTGTCGTCCTTAAAGACCCATTGTTGCCGGAGGAGTTGCTTCCGGCACACTgg    <  2:173251/197‑1 (MQ=255)
                                                                                                                                                                                                                                gCATTTAACTGAACAAAATGCCATGTACGAAACCTTTATTCAGTCATTCCGCCCGCTgg                                                                                                                                             >  2:312952/1‑59 (MQ=255)
                                                                                                                                                                                                                                    ttAACTGAACAAAATGCCATGTACGAAACCTTTATTCAGTCATTCCGCCCGCTGGTGcc                                                                                                                                         <  1:312961/59‑1 (MQ=255)
                                                                                                                                                                                                                                          |                                                                                                                                                                                           
GAAGGTTTAGATAAATCAACGCTGGCTGATGTCAAAAAGCAGTTGATCTGGCAAGGTTTTGGCGCACTGGCACCCAGCCTGATGGCATCGCCGTCGCAAAAACTGGCCGATGTACAGACACTTTTGCATGAAGCGGGTGTGGCGGATAACGTGATTTGTTTTGAAGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGCCAGAGTAGAAGAGTGCTGGCATTTAACTGAACAAAATGCCATGTACGAAACCTTTATTCAGTCATTCCGCCCGCTGGTGCCGCTTTTAAAAGAGGCGGCAGACGAGTTAACCCCGGAGCGGGCATTTCATATTCAGCTTTTACTGATCCATTTTTATCGCCGTGTCGTCCTTAAAGACCCATTGTTGCCGGAGGAGTTGCTTCCGGCACACTGGGC  >  CP009273/1458141‑1458562

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: