Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,587,631 Δ1 bp 8.5% coding (467/915 nt) catA‑II → catechol 1,2‑dioxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,587,6310T. 8.5% 139.0 / 7.6 47coding (467/915 nt)catA‑IIcatechol 1,2‑dioxygenase
Reads supporting (aligned to +/- strand):  ref base T (23/20);  new base . (2/2);  total (25/22)
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.24e-01

CCGTTGTACGTGGCTGGCGCACCGATTGCCCAAGGTGAAGTGCGCATGGACGACGGCAGCGAGGAGGGCGTGGCCACGGTGATGTTCCTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTG‑CGTC‑GGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGC  >  NC_002947/3587486‑3587767
                                                                                                                                                  |                                                                                                                                          
ccGTTGTACGTGGCTGGCGCACCGATTGCCCAAGGTGAAGTGCGCATGGACGACGGCAGCGAGGAGGGCGTGGCCACGGTGATGTTCCTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCtgtg                                                                                                                                         >  1:16766/1‑149 (MQ=255)
      tACGTGGCTGGCGCACCGATTGCCCAAGGTGAAGTGCGCATGGACGACGGCAGCGAGGAGGGCGTGGCCACGGTGATGTTCCTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATGc                                                                                                                                   <  2:16766/149‑1 (MQ=255)
        cGTGGCTGGCGCACCGATTGCCCAAGGTGAAGTGCGCATGGACGACGGCAGCGAGGAGGGCGTGGCCACGGTGATGTTCCTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCACAGGTCGACC‑GTGGCacccgg                                                                                                                                 <  2:3231/149‑6 (MQ=255)
                gcgcACCGATTGCCCAAGGTGAAGTGCGCATGGACGACGGCAGCGAGGAGGGCGTGGCCACGGTGATGTTCCTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATGCCAATACCCGt                                                                                                                         <  7:246643/149‑1 (MQ=255)
                      cGATTGCCCAAGGTGAAGTGCGCATGGACGACGGCAGCGAGGAGGGCGTGGCCACGGTGATGTTCCTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATGCCAATACCCGTGGTAcc                                                                                                                   <  1:123460/149‑1 (MQ=255)
                            cccAAGGTGAAGTGCGCATGGACGACGGCAGCGAGGAGGGCGTGGCCACGGTGATGTTCCTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATGc                                                                                                                                   <  1:271478/127‑1 (MQ=255)
                            cccAAGGTGAAGTGCGCATGGACGACGGCAGCGAGGAGGGCGTGGCCACGGTGATGTTCCTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATGc                                                                                                                                   >  2:271478/1‑127 (MQ=255)
                              cAAGGTGAAGTGCGCATGGACGACGGCAGCGAGGAGGGCGTGGCCACGGTGATGTTCCTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGtt                                                                                                           >  5:252767/1‑149 (MQ=255)
                               aaGGTGAAGTGCGCATGGACGACGGCAGCGAGGAGGGCGTGGCCACGGTGATGTTCCTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATg                                                                                                                                    <  4:108126/123‑1 (MQ=255)
                               aaGGTGAAGTGCGCATGGACGACGGCAGCGAGGAGGGCGTGGCCACGGTGATGTTCCTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATg                                                                                                                                    >  3:108126/1‑123 (MQ=255)
                                aGGAGAAGTGCGCATGGACGACGGCAGCGAGGAGGGCGTGGCCACGCTGATGATCCTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGCGTGCCAC‑GGTCGACCTCTGGCATGcc                                                                                                                                  >  5:110903/1‑124 (MQ=255)
                                 ggTGAATTGCGCATGGACGACGGCAGCGAGGAGGGCGTGGCCACGGTGACGTTCCGGGAAGGCCAGGTGCTGGACCGGCACGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGttctt                                                                                                        <  4:82131/149‑1 (MQ=255)
                                  gTGAAGTGCGCATGGACGACGGCAGCGAGGAGGGCGTGGCCACGGGGATGTTCCTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGttcttc                                                                                                       <  5:261040/149‑1 (MQ=255)
                                     aaGTGCGCATGGACGACGGCAGCGAGGAGGGCGTGGCCACGGTGATGTTCCTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGa                                                                                                     <  6:252767/148‑1 (MQ=255)
                                             aTGGACGACGGCAGCGAGGAGGGCGTGGCCACGGTGATGTTCCTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCtg                                                                                                                                           >  7:255742/1‑102 (MQ=255)
                                             aTGGACGACGGCAGCGAGGAGGGCGTGGCCACGGTGATGTTCCTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCtg                                                                                                                                           <  8:255742/102‑1 (MQ=255)
                                               ggACGACGGCAGCGAGGAGGGCGTGGCCACGGTGATGTTCCTGGAAGGCCAGGTGCTGGACCCGCACTGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGt                                                                                          >  7:188438/1‑149 (MQ=255)
                                               ggACGACGGCAGCGAGGAGGGCGTGGCCACGGTGATGTTCCTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGa                                                                                                     >  7:220973/1‑138 (MQ=255)
                                               ggACGACCGCAGCGAGCAGGGCGTGGCCACGGTGATGCTCCTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCGGCCGGGTGCCAC‑GGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGa                                                                                                     <  8:220973/138‑1 (MQ=255)
                                                  cgacgGCAGCGAGGAGGGCGTGGCCACGGTGATGTTCCTGGAAGGCCAGGTGCTCGACCCGCACGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATGCAAATACCCGTGGTACCTACTCGTTCTTCGACCAACGCCAGTCgg                                                                                       >  5:190610/1‑149 (MQ=255)
                                                   gacgGCAGCGAGGAGGGCGTGGCCACGGTGATGTTCCTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGc                                                                                      >  8:145049/1‑149 (MQ=255)
                                                        cAGCGAGGAGGGCGTGGCCACGGTGATGTTCCTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTAc                                                                                  >  8:108525/1‑148 (MQ=255)
                                                                      tGGCCACGGTGATGTTCCTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTG‑CGTC‑Ggcg                                                                  >  6:252464/1‑148 (MQ=255)
                                                                       ggCCACGGTGATGTACCTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCAA‑GGTCGACCTGTGGCATGCCAATACCCGAGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTG‑CGTC‑GGCGCa                                                                >  5:113578/1‑149 (MQ=255)
                                                                        gcgactgTGATGTTGCTGGAAGGCCAGGTGCTGTACCCGCACGGACGCCCGCTGCCGGGTGCCAC‑GGTCTACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTG‑CGTC‑TGGGCAt                                                               <  6:113578/143‑1 (MQ=255)
                                                                        gCCACGGTGATGTTCCTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCg                                                                                                             >  8:256339/1‑105 (MQ=255)
                                                                        gCCACGGTGATGTTCCTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCg                                                                                                             <  7:256339/105‑1 (MQ=255)
                                                                        gCCACGGTGATGTTCCTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTG‑CGTC‑GGCGCAt                                                               >  1:186848/1‑149 (MQ=255)
                                                                         ccACGGTGATGTTCCTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTG‑CGTC‑GGCGCATc                                                              <  3:187863/149‑1 (MQ=255)
                                                                           aCGGTGATGTTCCTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTG‑CGTC‑GGCGCATCGt                                                            >  2:245680/1‑149 (MQ=255)
                                                                                    ttCCTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTG‑CGTC‑GGCGCATCGTTACCGATGc                                                   >  6:108806/1‑149 (MQ=255)
                                                                                        tGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTG‑CGTC‑GGCGCATCGTTACCGATGCCCAg                                               >  2:66003/1‑149 (MQ=255)
                                                                                          gAAGGCCAGGTGCTGGACCCGCACTGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACATCCTG‑CGTC‑GGCGCATCGTTACCGATGCCCAgg                                              <  8:188438/148‑1 (MQ=255)
                                                                                          gAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATGTCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGCCGTACAACCTG‑CGTC‑GGCGCATCGTTACCGATGccaacgg                                             >  8:304135/1‑144 (MQ=255)
                                                                                              gCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTG‑CGTC‑GGCGCATCGTTACCGATGCCCAGGGgcgc                                         <  6:65848/149‑1 (MQ=255)
                                                                                                cAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTG‑CGTC‑GGCGCATCGTTACCGATGCCCAGGGGCGCTa                                       >  2:195688/1‑149 (MQ=255)
                                                                                                      tgtgACCCGCACGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATGCCAATAGCCGTGGTACCTACTAGTTCTTCGACCATAGCGAGTCGGCGTACAACCTG‑CGTC‑GGCGCATCGTTACCGATGCCCAGGGGCGCTAgcgcgc                                 <  6:190610/146‑1 (MQ=255)
                                                                                                                aCGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTG‑CGTC‑GGCGCATCGt                                                            >  8:271566/1‑112 (MQ=255)
                                                                                                                aCGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTG‑CGTC‑GGCGCATCGt                                                            <  7:271566/112‑1 (MQ=255)
                                                                                                                aCGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTG‑CGTC‑GGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATc                       >  3:241540/1‑149 (MQ=255)
                                                                                                                   gACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTG‑CGTC‑GGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTg                    <  5:108806/149‑1 (MQ=255)
                                                                                                                   gACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGGGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGTCGTC‑GGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGt                     <  1:245680/149‑1 (MQ=255)
                                                                                                                      gCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCCGCAAACAACCTG‑AATCGGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGc                   >  1:103676/1‑148 (MQ=255)
                                                                                                                                 catGCCACAGGTCGACC‑GTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTG‑CGTC‑GGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCt                            >  3:306385/3‑127 (MQ=255)
                                                                                                                                 catGCCACAGGTCGACC‑GTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTG‑CGTC‑GGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCt                            <  4:306385/125‑1 (MQ=255)
                                                                                                                                 catGCCACAGGTCGACC‑GTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTG‑CGTC‑GGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATgg      >  1:285259/3‑149 (MQ=255)
                                                                                                                                     ccAC‑GGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTG‑CGTC‑GGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGc  <  5:252464/149‑1 (MQ=255)
                                                                                                                                     ccAC‑GGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTG‑CGTC‑GGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGc  >  4:88699/1‑149 (MQ=255)
                                                                                                                                       aC‑GGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTG‑CGTC‑GGCGCATCGTTACCGATGCCCAggg                                             >  1:313672/1‑104 (MQ=255)
                                                                                                                                       aC‑GGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTG‑CGTC‑GGCGCATCGTTACCGATGCCCAggg                                             <  2:313672/104‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                          
CCGTTGTACGTGGCTGGCGCACCGATTGCCCAAGGTGAAGTGCGCATGGACGACGGCAGCGAGGAGGGCGTGGCCACGGTGATGTTCCTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTG‑CGTC‑GGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGC  >  NC_002947/3587486‑3587767

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