Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,587,965 C→A 7.0% A253D (GCC→GAC)  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,9650CA 7.0% 175.0 / 2.9 57A253D (GCC→GAC) catA‑IIcatechol 1,2‑dioxygenase
Reads supporting (aligned to +/- strand):  ref base C (27/26);  new base A (2/2);  total (29/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.26e-01

GGGGCCAACCCAGGAATGCCTGGACCTGCTGGGCCGTCATGGCCAGCGCCCGGCGCACGTGCACTTCTTTATCTCGGCCCCAGGGTACCGGCACCTGACCACGCAGATAAACCTGTCGGGGGACAAGTACCTGTGGGATGACTTTGCCTATGCCACACGGGATGGGCTGGTCGGGGAGGTGGTGTTCGTCGAAGGGCCGGATGGTCGGCATGCCGAGCTGAAGTTCGACTTCCAGTTGCAGCAGGCCCAGGGCGGTGCCGATGAGCAGCGCAGCGGGCGGCCGCGAGCTTTG  >  NC_002947/3587819‑3588110
                                                                                                                                                  |                                                                                                                                                 
ggggCCAACCCAGGAATGCCTGGACCTGCTGGGCCGTCATGGCCAGCGCCCGGCGCACGTGCACTTCTTTATCTCGGCCCCAGGGTACCGGCACCTGACCACGCAGATAAACCTGTCGGGGGACAAGTACCTGTGGGATGACTTTGCCTa                                                                                                                                                >  2:594820/1‑150 (MQ=255)
 gggCCAACCCAGGAATGCCTGGACCTGCTGGGCCGTCATGGCCAGCGCCCGGCGCACGTGCACTTGTTGATCTCGGCCCCAGGGTACCGGCACCTGACCACGCAGATAAACCTGTCGGGGGACAAGTACCTGTGGGATGACTTTGCCTAt                                                                                                                                               <  1:347750/150‑1 (MQ=255)
  ggCGAACCCAGGAATGCCTGGACCTGCTGGGCCGTCATGGCCAGCGCCCGGCGTACGTGCACTTCTTTATCTCGGCCCCCGGGTACCGGCACCTGACCCCGCAGATAAACCTGTCGGGGGACAAGTACCTGTGGGATGACTTTGCCTATg                                                                                                                                              <  2:44745/150‑1 (MQ=255)
     cAACCCAGGAATGCCTGGACCTGCTGGGCCGTCATGGCCAGCGCCCGGCGCACGTGCACTTCTTTATCTCGGCCCCAGGGTACCGGCACCTGACCACGCAGATAAACCTGTCGGGGGACAAGTACCTGTGGGATGACTTTGCCTATGCca                                                                                                                                           >  2:689760/1‑150 (MQ=255)
         ccAGGAATGCCTGGACCTGCTGGGCCGTCATGGCCAGCGCCCGGCGCACGTGCACTTCTTTATCTCGGCCCCAGGGTACCGGCACCTGACCACGCAGATAAACCTGTCGGGGGACAAGTACCTGTGGGATGACTTTGCCTATGCCACACg                                                                                                                                       >  1:224135/1‑150 (MQ=255)
         ccAGGAATGCCTGGACCTGCTGGGCCGTCATGGCCAGCGCCCGGCGCACGTGCACTTCTTTATCTCGGCCCCAGGGTACCGGCACCTGACCACGCAGATAAACCTGTCGGGGGACAAGTACCTGTGGGATGACTTTGCCTATGCCACACg                                                                                                                                       >  2:321122/1‑150 (MQ=255)
           aGGAATGCCTGGACCCGCTGGGCCGTCATGGCCAGCGCCCGGCGAACGTGCTGTTCTTTATCTCGGCCCCAGGGTACCGGCACCTGACCACGCAGATAAACCTGTCGGGGGACAAGTACCTGTGGGATGACTTTGCCTATGCCACACggg                                                                                                                                     >  1:581248/1‑150 (MQ=255)
                             tgGGCCGTCATGGCCAGCGCCCGGCGCACGTGCACTTCTTTATCTCGGCCCCAGGGTACCGGCACCTGACCACGCAGATAAACCTGTCGGGGGACAAGTACCTGTGGGATGACTTTGCCTATGCCACACGGGATGGGCTGGTCGGGGAgg                                                                                                                   <  1:321122/150‑1 (MQ=255)
                             tgGGCCGTCATGGCCAGCGCCCGGCGCACGTGCACTTCTTTATCTCGGCCCCAGGGTACCGGCACCTGACCACGCAGATAAACCTGTCGGGGGACAAGTACCTGTGGGATGACTTTGCCTATGCCACACGGGATGGGCTGGTCGGGGAgg                                                                                                                   <  2:796788/150‑1 (MQ=255)
                                  tgTCATGGCCAGCGCCCGGCGCACGTGCACTTCTTTATCTCGGCCCCAGGGTACCGGCACCTGACCACGCAGATAAACCTGTCGGGGGACAAGTACCTGTGGGATGACTTTGCCTATGCCACACGGGATGGGCTGGTCGGGGAGgtggtg                                                                                                              <  2:471890/149‑1 (MQ=255)
                                           cAGCGCCCGGCGCACGTGCACTTCTTTATCTCGGCCCCAGGGTACCGGCACCTGACCACGCAGATAAACCTGTCGGGGGACAAGTACCTGTGGGATGACTTTGCCTATGCCACACGGGATGGGCTGGTCGGGGAGGTGGTGTTCGTCGaa                                                                                                     >  1:846236/1‑150 (MQ=255)
                                           cAGCGCCCGGCGCACGTGCACTTCTTTATCTCGGCCCCAGGGTACCGGCACCTGACCACGCAGATAAACCTGTCGGGGGACAAGTACCTGTGGGATGACTTTGCCTATGCCACACGGGATGGGCTGGTCGGGGAGGTGGTGTTCGTCGaa                                                                                                     <  1:755634/150‑1 (MQ=255)
                                                    gcgcACGTGCACTTCTTTATCTCGGCCCCAGGGTACCGGCACCTGACCACGCAGATAAACCTGTCGGGGGACAAGTACCTGTGGGATGACTTTGCCTCTCCCTCACTAGATGGGCTGGTCGTGGAGGTGGTGTTCGTCGAAGGGCCGGAt                                                                                            >  2:22805/1‑150 (MQ=255)
                                                    gcgcACGTGCACTTCTTTATCTCGGCCCCAGGGTACCGGCACCTGACCACGCAGATAAACCTGTCGGGGGACAAGTACCTGTGGGATGACTTTGCCTATGCCACACGGGATGGGCTGGTCGGGGAGGTGGTGTTCGTCGAAGGGCCGGAt                                                                                            <  2:665843/150‑1 (MQ=255)
                                                        aCGTGCACTTCTTTATCTCGGCCCCAGGGTACCGGCACCTGACCACGCAGATAAACCTGTCGGGGGACAAGTACCTGTGGGATGACTTTGCCTATGCCACACGGGATGGGCTGGTCGGGGAGGTGGTGTTCGTCGAAGGGCCGGATGGTc                                                                                        <  1:689760/150‑1 (MQ=255)
                                                             cACTTCTTTATCTCGGCCCCAGGGTACCGGCACCTGACCACGCAGATAAACCTGTCGGGGGACAAGTACCTGTGGGATGACTTTGCCTATGCCACACGGGATGGGCTGGTCGGGGAGGTGGTGTTCGTCGAAGGGCCGGATGGTCGGCAt                                                                                   >  2:659755/1‑150 (MQ=255)
                                                              aCTTCTTTATCTCGGCCCCAGGGTACCGGCACCTGACCACGCAGATAAACCTGTCGGGGGACAAGTACCTGTGGGATGACTTTGCCTATGCCACACGGGATGGGCTGGTCGGGGAGGTGGTGTTCGTCGAAGGGCCGGa                                                                                             <  2:224164/139‑1 (MQ=255)
                                                              aCTTCTTTATCTCGGCCCCAGGGTACCGGCACCTGACCACGCAGATAAACCTGTCGGGGGACAAGTACCTGTGGGATGACTTTGCCTATGCCACACGGGATGGGCTGGTCGGGGAGGTGGTGTTCGTCGAAGGGCCGGa                                                                                             >  1:224164/1‑139 (MQ=255)
                                                              aCTTCTTTATCTCGGCCCCAGGGTACCGGCACCTGACCACGCAGATAAACCTGTCGGGGGACAAGTACCTGTGGGATGACTTTGCCTATGCCACACGGGATGGGCTGGTCGGGGAGGTGGTGTTCGTCGAAGGGCCGGATGGTCGGCATg                                                                                  >  2:553417/1‑150 (MQ=255)
                                                               cttcttTATCTCGGCCCCAGGGTACCGGCACCTGACCACGCAGATAAACCTGTCGGGGGACAAGTACCTGTGGGATGACTTTGCCTATGCCACACGGGATGGGCTGGTCGGGGAGGTGGTGTTCGTCGAAGGGCCGGATGGTCGGCATGc                                                                                 >  1:766177/1‑150 (MQ=255)
                                                                ttcttTATCTCGGCCCCAGGGTACCGGCACCTGACCACGCAGATAAACCTGTCGGGGGACAAGTACCTGTGGGATGACTTTGCCTATGCCACACGGGATGGGCTGGTCGGGGAGGTGgtgt                                                                                                             <  1:810462/121‑1 (MQ=255)
                                                                ttcttTATCTCGGCCCCAGGGTACCGGCACCTGACCACGCAGATAAACCTGTCGGGGGACAAGTACCTGTGGGATGACTTTGCCTATGCCACACGGGATGGGCTGGTCGGGGAGGTGgtgt                                                                                                             >  2:810462/1‑121 (MQ=255)
                                                                ttcttTATCTCGGCCCCAGGGTACCGGCACCTGACCACGCAGATAAACCTGTCGGGGGACAAGTACCTGTGGGATGACTTTGCCTATGCCACACGGGATGGGCTGGTCGGGGAGGTGGTGTTCGTCGAAGGGCCGGATGGTCGGCATGcc                                                                                >  1:359185/1‑150 (MQ=255)
                                                                ttcttTATCTCGGCCCCAGGGTACCGGCACCTGACCACGCAGATAAACCTGTCGGGGGACAAGTACCTGTGGGATGACTTTGCCTATGCCACACGGGATGGGCTGGTCGGGGAGGTGGTGTTCGTCGAAGGGCCGGATGGTCGGCATGcc                                                                                >  1:802416/1‑150 (MQ=255)
                                                                    ttATCTCGGCCCCAGGGTACCGGCACCTGACCACGCAGATAAACCTGTCGGGGGACAAGTACCTGTGGGATGACTTTGCCTATGCCACACGGGATGGGCTGGTCGGGGAGGTGGTGTTCGTCGAAGGGCCGGATGGTCGGCATGCCGAGc                                                                            >  1:457426/1‑150 (MQ=255)
                                                                     tATCTCGGCCCCAGGGTACCGGCACCTGACCACGCAGATAAACCTGTCGGGGGACAAGTACCTGTGGGATGACTTTGCCTATGCCACACGAGATGGGCTGGTCGGGGAGGTGGTGTTAGTCTAAGGGCCGGATGGTCGGCATGCCGAGCt                                                                           <  1:22805/150‑1 (MQ=255)
                                                                      aTCTCGGCCCCAGGGTACCGGCACCTGACCACGCAGATAACCCTGTCGGGGGACAAGTACCTGTGGGATGACTTTGCCTATGCCACACGGGATGGGCTGGTCGGGGAGGGGGTGTTCGTCGAAGGGCCGGATGGTCGGCATGCCGAGCTg                                                                          >  1:812028/1‑150 (MQ=255)
                                                                            gCCCCAGGGTACCGGCACCTGACCACGCAGATAAACCTGTCGGGGGACAAGTACCTGTGGGATGACTTTGCCTATGCCACACGGGATGGGCTGGTCGGGGAGGTGGTGTTCGTCGAAGGGCCGGATGGTCGGCATGCCGAGCTGAAGTTc                                                                    <  1:326513/150‑1 (MQ=255)
                                                                                  gggTACCGGCACCTGACCACGCAGATAAACCTGTCGGGGGACAAGTACCTGTGGGATGACTTTGCCTATGCCACACGGGATGGGCTGGTCGGGGAGGTGGTGTTCGTCGAAGGGCCGGATGGTCGGCATGCCGAGCTGAAGTTCGACTTc                                                              <  1:659755/150‑1 (MQ=255)
                                                                                     tACCGGCACCTGACCACGCAGATAAACCTGTCGGGGGACAAGTACCTGTGGGATGACTTTGCCTATGCCACACGGGATGGGCTGGTCGGGGAGGTGGTGTTCGTCGAAGGGCCGGATGGTCGGCATGCCGAGCTGAAGTTCGACTTCCAg                                                           <  1:317574/150‑1 (MQ=255)
                                                                                               tGACCACGCAGATAAACCTGTCGGGGGACAAGTACCTGTGGGATGACTTTGCCTATGCCACACGGGATGGGCTGGTCGGGGAggtg                                                                                                                 >  1:70359/1‑86 (MQ=255)
                                                                                               tGACCACGCAGATAAACCTGTCGGGGGACAAGTACCTGTGGGATGACTTTGCCTATGCCACACGGGATGGGCTGGTCGGGGAggtg                                                                                                                 <  2:70359/86‑1 (MQ=255)
                                                                                                    aCGCAGATAAACCTGTCGGGGGACAAGTACCTGTGGGATGACTTTGCCTATGCCACACGGGATGGGCTGGTCGGGGAGGTGGTGTTCGTCGAAGGGCCGGATGGTCGGCATGCCGAGCTGAAGTTCGACTTCCAGTTGCAGCAGGCCCAg                                            >  2:93022/1‑150 (MQ=255)
                                                                                                     cGCAGATAAACCTGTCGGGGGACAAGTACCTGTGGGATGACTTTGCCTATGCCACACGGGATGGGCTGGTCGGGGAGGTGGTGTTCGTCGAAGGGCCGGATGGTCGGCATGCCGAGCTGAAGTTCGACTTCCAGTTGCAGCAGGCCCAgg                                           <  1:719878/150‑1 (MQ=255)
                                                                                                        aGATAAACCTGTCGGGGGACAAGTACCTGTGGGATGACTTTGCCTATGCCACACGGGATGGGCTGGTCGGGGAGGTGGTGTTCGTCGAAGGGCCGGATGGTCGGCATGCCGAGCTGAAg                                                                       >  1:137265/1‑119 (MQ=255)
                                                                                                        aGATAAACCTGTCGGGGGACAAGTACCTGTGGGATGACTTTGCCTATGCCACACGGGATGGGCTGGTCGGGGAGGTGGTGTTCGTCGAAGGGCCGGATGGTCGGCATGCCGAGCTGAAg                                                                       <  2:137265/119‑1 (MQ=255)
                                                                                                             aaCCTGTCGGGGGACAAGTACCTGTGGGATGACTTTGCCTATGCCACACGGGATGGGCTGGTCGGGGAGGTGGTGTTCGTCGAAGGGCCGGATGGTCGGCATGCCGAGCTGAAGTTCGACTTCCAGTTGCAGCAGGCCCAGGGCGGTGcc                                   <  1:93022/150‑1 (MQ=255)
                                                                                                              aCCTGTCGGGGGACAAGTACCTGTGGGATGACTTTGCCTATGCcac                                                                                                                                          <  2:691914/46‑1 (MQ=255)
                                                                                                              aCCTGTCGGGGGACAAGTACCTGTGGGATGACTTTGCCTATGCcac                                                                                                                                          >  1:691914/1‑46 (MQ=255)
                                                                                                              aCCTGTCGGGGGACAAGTACCTGTGGGATGACTTTGCCTATGCCACACGGGATGGGCTGGTCGGGGAGGTGGTGTTCGTCGAAGGGCCGGATGGTCGGCATGCCGAGCTGAAGTTCGCCTTCCAGTGGCAGCAGGCCCAGGGCGGTGCCg                                  >  2:534975/1‑150 (MQ=255)
                                                                                                                cTGTCGGGGGACAAGTACCTGTGGGATGACTTTGCCTATGCCACACGGGATGGGCTGGTCGGGGAGGTGGTGTTCGTCGAAGGGCCGGATGGTCGGCATGCCGAGCTGAAGTTCGACTTCCAGTTGCAGCAGGCCCAGGGCGGTGCCGAt                                <  1:553417/150‑1 (MQ=255)
                                                                                                                  gTCGGGGGACAAGTACCTGTGGGATGACTTTGCCTATGCCACACGGGATGGGCTGGTCGGGGAGGTGGTGTTCGTCGAAGGGCCGGATGGTCGGCATGCCGAGCTGAAGTTCGACTTCCAGTTGCAGCAGGCCCAGGGCGGTGCCGATGa                              <  2:224135/150‑1 (MQ=255)
                                                                                                                      ggggACAAGTACCTGTGGGATGACTTTGCCTATGCCACACGGGATGGGCTGGTCGGGGAGGTGGTGTTCGTCGAAGGGCCGGATGGTCGGCATGCCGAGCTGAAGTTCGACTTCCAGTTGCAGCAGGCCCAGGGCGGTGCCGATGagcag                          <  1:594820/150‑1 (MQ=255)
                                                                                                                      ggggACAAGTACCTGTGGGATGACTTTGACTATGCCACACGGGATGGGCTGGTCGGGGAGGTGGTGTTCGTCGAAGGGCCGGATGGTc                                                                                        <  1:679901/88‑1 (MQ=255)
                                                                                                                      ggggACAAGTACCTGTGGGATGACTTTGACTATGCCACACGGGATGGGCTGGTCGGGGAGGTGGTGTTCGTCGAAGGGCCGGATGGTc                                                                                        <  1:679843/88‑1 (MQ=255)
                                                                                                                      ggggACAAGTACCTGTGGGATGACTTTGACTATGCCACACGGGATGGGCTGGTCGGGGAGGTGGTGTTCGTCGAAGGGCCGGATGGTc                                                                                        >  2:679901/1‑88 (MQ=255)
                                                                                                                      ggggACAAGTACCTGTGGGATGACTTTGACTATGCCACACGGGATGGGCTGGTCGGGGAGGTGGTGTTCGTCGAAGGGCCGGATGGTc                                                                                        >  2:679843/1‑88 (MQ=255)
                                                                                                                        ggACAAGTACCTGTGGGATGACTTTGCCTATGCCACACGGGATGGGCTGGTCGGGGAGGTGGTGTTCGTCGAAGGGCCGGATGCTCGGCATGCCGAGCTGAAGTTCGACTTCCAGTTGCAGCAGGCCCAGGGCGGTGCCGATGAgcagcg                        >  2:248979/1‑150 (MQ=255)
                                                                                                                              gTACCTGTGGGATGACTTTGCCTATGCCACACGGGATGGGCTGGTCGGGGAGGTGGTGTTCGTCGAAGGGCCGGATGGTCGGCATGCCGAGCTGAAGTTCGACTTCCAGTTGCAGCAGGCCCAGGGCGGTGcc                                   <  1:822709/133‑1 (MQ=255)
                                                                                                                              gTACCTGTGGGATGACTTTGCCTATGCCACACGGGATGGGCTGGTCGGGGAGGTGGTGTTCGTCGAAGGGCCGGATGGTCGGCATGCCGAGCTGAAGTTCGACTTCCAGTTGCAGCAGGCCCAGGGCGGTGcc                                   >  2:822709/1‑133 (MQ=255)
                                                                                                                              gTACCTGTGGGATGACTTTGCCTATGCCACACGGGATGGGCTGGTCGGGGAGGTGGTGTTCGTCGAAGGGCCGGATGGTCGGCATGCCGAGCTGAAGTTCGACTTCCAGTTGCAGCAGGCCCAGGGCGGTGCCGATGAGCAGCGCAgcgg                  <  2:846236/150‑1 (MQ=255)
                                                                                                                                     tgGGATGACTTTGCCTATTCCACACGGGATGGGCTGGTCGGGGAGGTGGTGTTCGTCGAAGGGCCGGATGGTCGGCATGCCGAGCTGAAGTTCGACTTCCAGTTGCAGCAGGCCCATGGCGGTGCCGATGAGCAGCGCAGCGGGCGGCcg           >  2:393700/1‑150 (MQ=255)
                                                                                                                                       ggATGACTTTGCCTATGCCACACGGGATGGGCTGGTCGGGGGGGTGGTGTTCGTCGAAGGGCCGGATGGTCGGCATGCCGAGCTGAAGTTCGACTTCCAGTTGCAGCAGGCCCAGGGCGGTGCCGATGAGCAGCGCAGCGGGCGGCcgcg         >  2:273991/1‑150 (MQ=255)
                                                                                                                                       ggATGACTTTGCCTATGCCACACGGGATGGGCTGGTCGGGGAGGTGGTGTTCGTCGAAGGGCCGGATGGTCGGCATGCCGAGCTGAAGTTCGACTTCCAGTTGCAGCAGGCCCAGGGCGGTGCCGATGAGCAGCGCAGCGGGCGGCcgcg         >  2:340173/1‑150 (MQ=255)
                                                                                                                                          tGACTTTGCCTATGCCACACGGGATGGGCTGGTCGGGGAGGTGGTGTTCGTCGaa                                                                                                     >  2:298214/1‑55 (MQ=255)
                                                                                                                                          tGACTTTGCCTATGCCACACGGGATGGGCTGGTCGGGGAGGTGGTGTTCGTCGaa                                                                                                     <  1:298214/55‑1 (MQ=255)
                                                                                                                                              tttGCCTATGCCACACGGGATGGGCTGGTCGGGGAGGTGGTGTTCGTCGAAGGGCCGGATGGTCGGCATGCCGAGCTGAAGTTCGACTTCCAGTTGCAGCAGGCCCAGGGCGGTGCCGATGAGCAGCGCAGCGGGCGGCCGCGAGCTTTg  <  2:812028/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GGGGCCAACCCAGGAATGCCTGGACCTGCTGGGCCGTCATGGCCAGCGCCCGGCGCACGTGCACTTCTTTATCTCGGCCCCAGGGTACCGGCACCTGACCACGCAGATAAACCTGTCGGGGGACAAGTACCTGTGGGATGACTTTGCCTATGCCACACGGGATGGGCTGGTCGGGGAGGTGGTGTTCGTCGAAGGGCCGGATGGTCGGCATGCCGAGCTGAAGTTCGACTTCCAGTTGCAGCAGGCCCAGGGCGGTGCCGATGAGCAGCGCAGCGGGCGGCCGCGAGCTTTG  >  NC_002947/3587819‑3588110

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