Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,237,547 Δ1 bp 9.8% coding (703/1122 nt) catB ← muconate cycloisomerase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,237,5470G. 9.8% 125.4 / 7.5 41coding (703/1122 nt)catBmuconate cycloisomerase 1
Reads supporting (aligned to +/- strand):  ref base G (19/18);  new base . (2/2);  total (21/20)
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.56e-01

GCGCGCGGGCCACCATTCTTGGCGATTTTCAGGGCGAAGATGCTGGCGGCGCCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTCGACCCGCACGCTGGCGCTGT  >  NC_002947/4237405‑4237678
                                                                                                                                              |                                                                                                                                   
gcgcgcgGGCCACCATTCTTGGCGATTTTCAGGGCGAAGATGCTGGCGGCGCCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAgc                                                                                                                               >  4:72411/1‑149 (MQ=255)
 cgcgcgGGCCACCATTCTTGGCGATTTTCAGGGCGAAGATGCTGGCGGCGCCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGTCTCATCGGCCATAATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAgcg                                                                                                                              <  7:298511/149‑1 (MQ=255)
  gcgcgGGCCACCATTCTTGGCGATTTTCAGGGCGAAGATGCTGGCGGCGCCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAgcgc                                                                                                                             >  4:142144/1‑149 (MQ=255)
      gggCCACCATTCTTGGCGATTTTCAGGGCGAAGATGCTGGCGGCGCCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGtt                                                                                                                         <  3:221232/149‑1 (MQ=255)
         ccaccaTTCTTGGCGATTTTCAGGGCGAAGATGCTGGCGGCGCCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGAt                                                                                                                      <  4:5950/149‑1 (MQ=255)
         ccaccaTTCTTGGCGATTTTCAGGGCGAAGATGCTGGCGGCGCCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGAt                                                                                                                      <  6:202708/149‑1 (MQ=255)
                   tGGCGATTTTCAGGGCGAAGATGCTGGCGGCGCCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATc                                                                                                            >  8:1610/1‑149 (MQ=255)
                    ggCGATTTTCAGGGCGAAGATGCTGGCGGCGCCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGaaa                                                                                                              <  5:265856/146‑1 (MQ=255)
                       gATTTTCAGGGCGAAGATGCTGGCGGCGCCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGgct                                                                                                        >  8:103202/1‑149 (MQ=255)
                       gATTTTCAGGGCGAAGATGCTGGCGGCGCCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGgct                                                                                                        >  1:41882/1‑149 (MQ=255)
                           ttCAGGGCGAAGATGCTGGCGGCGCCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTc                                                                                                    >  7:140811/1‑149 (MQ=255)
                                    aaGATGCTGGCGGCGCCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTc                                                                                           <  2:41882/149‑1 (MQ=255)
                                     aGATGCTGGCGGCGCCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTc                                                                                           >  7:290750/1‑148 (MQ=255)
                                         gCTGGCGGCGCCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAAt                                                                                                             <  6:293191/126‑1 (MQ=255)
                                         gCTGGCGGCGCCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAAt                                                                                                             >  5:293191/1‑126 (MQ=255)
                                             gcggcgCCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCg                                                                                    <  8:290750/147‑1 (MQ=255)
                                                gcgcCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTc                                                                               <  3:72411/149‑1 (MQ=255)
                                                     gTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGccc                                                                           >  8:7106/1‑148 (MQ=255)
                                                             ccAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTgg                                                                  <  8:140811/149‑1 (MQ=255)
                                                                     aaGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGACAGGcgcg                                                          >  3:83451/1‑149 (MQ=255)
                                                                                gACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTgg                                                <  7:103202/148‑1 (MQ=255)
                                                                                     tttCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCg                                          >  4:188600/1‑149 (MQ=255)
                                                                                          atcgaCTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTcc                                       >  8:73773/1‑147 (MQ=255)
                                                                                            cgaCTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAg                                    >  5:68815/1‑148 (MQ=255)
                                                                                             gaCTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTa                                  <  3:142144/149‑1 (MQ=255)
                                                                                             gaCTCATCGGCCATGATCGGAGCCGCACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGAGCGAAATAGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGt                                   >  8:146215/1‑148 (MQ=255)
                                                                                              aCTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTAc                                 <  3:99174/149‑1 (MQ=255)
                                                                                                      gCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGAGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGAc                         <  7:7106/149‑1 (MQ=255)
                                                                                                        cATGATCGGAGCCGTACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAg                                    <  1:124113/136‑1 (MQ=255)
                                                                                                          tGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTCg                     <  7:1610/149‑1 (MQ=255)
                                                                                                                gAGCCGGACTGCGCTGGTTCAGGCGCACCT‑GCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTcc                                       >  2:41147/1‑124 (MQ=255)
                                                                                                                gAGCCGGACTGCGCTGGTTCAGGCGCACCT‑GCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTcc                                       >  2:36386/1‑124 (MQ=255)
                                                                                                                gAGCCGGACTGCGCTGGTTCAGGCGCACCT‑GCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTcc                                       <  1:41147/124‑1 (MQ=255)
                                                                                                                gAGCCGGACTGCGCTGGTTCAGGCGCACCT‑GCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTcc                                       <  1:36386/124‑1 (MQ=255)
                                                                                                                 aGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTCGACCCGCa              <  6:68815/149‑1 (MQ=255)
                                                                                                                      gACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTCGACCCGCACGCTg         >  6:275687/1‑149 (MQ=255)
                                                                                                                           cgcTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTCGAc                   <  4:203726/134‑1 (MQ=255)
                                                                                                                           cgcTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTCGAc                   >  3:203726/1‑134 (MQ=255)
                                                                                                                           cgcTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAAAACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTCGCCCCGCACGCTGGCGCt    >  8:82223/1‑149 (MQ=255)
                                                                                                                            gcTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTCGACCCGCACGCTGGCGCTg   >  7:226922/1‑149 (MQ=255)
                                                                                                                             cTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTCGACCCGCACGCTGGCGCTGt  >  1:140282/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                   
GCGCGCGGGCCACCATTCTTGGCGATTTTCAGGGCGAAGATGCTGGCGGCGCCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTCGACCCGCACGCTGGCGCTGT  >  NC_002947/4237405‑4237678

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