Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,828,849 A→C 12.5% D342A (GAC→GCC)  PP_2480 → xanthine dehydrogenase accessory factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,828,8490AC12.5% 144.0 / 7.4 53D342A (GAC→GCC) PP_2480xanthine dehydrogenase accessory factor
Reads supporting (aligned to +/- strand):  ref base A (22/24);  new base C (3/4);  total (25/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 = 1.00e+00

ACCCTGCTATGTCGGTCAGTTGGGGCCGCGTGAGCGCACCGAGCGCTTGTTGGCCGGCATCCACGAACAATTGGCTAAACCACAGGACGAGCTGCCGGGGCTGGAATGTTTGCATTATCCGATCGGCCTGGACCTGGGTGGCGACACCCCGGAAAGCGTGGCCATGGCGGTGCTGGCGGAAATTCAGGCAGTGCTCAACGGCCGCAATGGCGGCAGCCTGCGCTTCAGATCTGCGAGCATCCATGACAGCGACCCTGTGCTGAAGGCTGGCCGGGATCAAACGGA  >  NC_002947/2828706‑2828990
                                                                                                                                               |                                                                                                                                             
aCCCTGCTATGTCGGTCAGTTGGGGCCGCGTGAGCGCACCGAGCGCTTGTTGGCCGGCATCCACGAACAATTGGCTAAACCACAGGACGAGCTGCCGGGGCTGGAATGTTTGCATTATCCGATCGGCCTGGACCGGGGTGTCGCCAccc                                                                                                                                          <  6:274860/149‑1 (MQ=255)
 cccTGCTATGTCGGTCAGTTGGGGCCGCGTGAGCGCACCGAGCGCTTGTTGGCCGGCATCCACGAACAATTGGCTAAACCACAGGACGAGCTGCCGGGGCTGGAATGTTTGCATTATCCGATCGGCCTGGACCTGGGTGGCGACAcccc                                                                                                                                         >  7:1044/1‑149 (MQ=255)
  ccTGCTATGTCGGTCAGTTGGGGCCGCGTGAGCGCACCGAGCGCTTGTTGGCCGGCATCCACGAACAATTGGCTAAACCACAGGACGAGCTGCCGGGGCTGGAATGTTTGCATTATCCGATCGGCCTGGACCTGGGTGGCGACACCCCg                                                                                                                                        >  1:254512/1‑149 (MQ=255)
   cttcTATGTCGGGCAGTTGGGGCCGCGTGAGCGCACCGAGCGCTTGTTGGCCGGCATCCACGAACAATTGGCTAAACCACAGGACGAGCTGCCGGGGCTTGAATGTTTGCATTATCCGATCGGCCTGGACCTGGGTGTCGCCACCCAgg                                                                                                                                       <  5:215734/146‑1 (MQ=255)
         tGTCGGTCAGTTGGGGCCGCGTGAGCGCACCGAGCGCTTGTTGGCCGGCATCCACGAACAATTGGCTAAACCACAGGACGAGCTGCCGGGGCTGGAATGTTTGCATTATCCGATCGGCCTGGACCTGGGTGGCGACACCCCGGAAAGc                                                                                                                                  <  7:289144/148‑1 (MQ=255)
           tCGGTCAGTTGGGGCCGCGTGAGCGCACCGCGCGCTTGTTGGCCGGCATCCACGAACAATTGGCTAAACCACAGGACGAGCTGCCGGGGCTGGAATGTTTGCATTATCCGATCGGCCTGGACCTGGGTGGCGACACCCCGGAAAGCGTg                                                                                                                               >  8:221011/1‑149 (MQ=255)
           tCGGTCAGTTGGGGCCGCGTGAGCGCACCGAGCGCTTGTTGGCCGGCATCCACGAACAATTGGCTAAACCACAGGACGAGCTGCCGGGGCTGGAATGTTTGCATTATCCGATCGGCCTGGACCTGGGTGGCGACACCCCGGAAAGCGt                                                                                                                                <  2:254512/148‑1 (MQ=255)
           tCGGTCAGTTGGGGCCGCGTGAGCGCACCGAGCGCTTGTTGGCCGGCATCCACGAACAATTGGCTAAACCACAGGACGAGCTGCCGGGGCTGGAATGTTTGCATTATCCGATCGGCCTGGACCTGGGTGGCGACACCCCGGAAAGCGTg                                                                                                                               <  3:34/149‑1 (MQ=255)
              gTCAGTTGGGGCCGCGTGAGCGCACCGAGCGCTTGTTGGCCGGCATCCACGAACAATTGGCTAAACCACAGGACGAGCTGCCGGGGCTGGAATGTTTGCATTATCCGATCGGCCTGGACCTGGGTGGCGCCACCCCGGAAAGCGTGGc                                                                                                                             >  7:49689/1‑148 (MQ=255)
                     ggggCCGCGTGAGCGCACCGAGCGCTTGTTGGCCGGCATCCACGAACAATTGGCTAAACCACAGGACGAGCTGCCGGGGCTGGAATGTTTGCATTATCCGATCGGCCTGGACCTGGGTGGCGACACCCCGGAAAGCGTGGCCATGGCgg                                                                                                                     >  5:13437/1‑149 (MQ=255)
                       ggCCGCGTGAGCGCACCGAGCGCTTGTTGGCCGGCATCCACGAACAATTGGCTAAACCACAGGACGAGCTGCCGGGGCTGGAATGTTTGCATTATCCGATCGGCCTGGACCTGGGTGGCGACACCCCGGAAAGCGTGGCCATGGCGGTg                                                                                                                   <  8:1044/149‑1 (MQ=255)
                        gcagcgTGAGCGCACCGAGCGCTTGTTGGCCGGCATCCACGAACAATTGGCTAAACCACAGGACGAGCTTCCGGGGCTGGAATGTTTGCATTATCCGATCGGCCTGGACCTGGGTGGCGACACCCCGGAAAGCGTGGCCATGGCGGTGc                                                                                                                  <  8:280214/146‑1 (MQ=255)
                           gcgTGAGCGCACCGAGCGCTTGTTGGCCGGCATCCACGAACAATTGGCTAAACCACAGGACGAGCTGCCGGGGCTGGAATGTTTGCATTATCCGATCGGCCTGGACCTGGGTGGCGACACCCCGGAAAGCGTGGCCATg                                                                                                                         <  3:130550/139‑1 (MQ=255)
                           gcgTGAGCGCACCGAGCGCTTGTTGGCCGGCATCCACGAACAATTGGCTAAACCACAGGACGAGCTGCCGGGGCTGGAATGTTTGCATTATCCGATCGGCCTGGACCTGGGTGGCGACACCCCGGAAAGCGTGGCCATg                                                                                                                         >  4:130550/1‑139 (MQ=255)
                             gTGAGCGCACCGAGCGCTTGTTGGCCGTCATCCACGAACAATTGGCTAAACCACAGGACGAGCTGCCGGGGCTGGAATGTTTGCATTATCCGATCTGCCTGGCCCTGGGTGGCTCCACCACGGAAACCGTGGCCATGGCGGTGCTGGCg                                                                                                             >  8:232626/1‑149 (MQ=255)
                                aGCGCACCGAGCGCTTGTTGGCCGGCATCCACGAACAATTGGCTAAACCACAGGACGAGCTGCCGGGGCTGGAATGTTTGCATTATCCGATCGGCCTGGACCTGGGTGGCGACACCCCGGAAAGCGTGGCCATGGCGGTGCTGGCGGa                                                                                                           >  6:183534/1‑148 (MQ=255)
                                     aCCGAGCGCTTGTTGGCCGGCATCCACGAACAATTGGCTAAACCACAGGACGAGCTGCCGGGGCTGGAATGTTTGCATTATCCGATCGGCCTGGACCTGGGTGGCGa                                                                                                                                               >  6:6325/1‑107 (MQ=255)
                                     aCCGAGCGCTTGTTGGCCGGCATCCACGAACAATTGGCTAAACCACAGGACGAGCTGCCGGGGCTGGAATGTTTGCATTATCCGATCGGCCTGGACCTGGGTGGCGa                                                                                                                                               <  5:6325/107‑1 (MQ=255)
                                      ccGAGCGCTTGTTGGCCGGCATACACGAACAATTGGCTAAACCACAGGACGAGCTGCCGGGGCTGGAATGTTTGCATTATCCGATCGGCCTGGACCTGGGTGGCGACACCCCGGAAAGCGTGGCCATGGCGGTGCTGGCGGAAATTCAg                                                                                                    >  6:304578/1‑149 (MQ=255)
                                             cTTGTTGGCCGGCATCCACGAACAATTGGCTAAACCACAGGACGAGCTGCCGGGGCTGGAATGTTTGCATTATCCGATCGGCCTGGACCTGGGTGGCGACACCCCGGAAAGCGTGGCCATGGCGGTGCTGGCGGAAATTCAGGCAGTGc                                                                                             >  1:193791/1‑149 (MQ=255)
                                                   ggCCGGCATCCACGAACAATTGGCTAAACCACAGGACGAGCTGCCGGGGCTGGAATGTTTGCATTATCCGATCGGCCTGGACCTGGGTGGCGACACCCCGGAAAGCGTGGCCATGGCGGTGCTGGCGGAAATTCAGGCAGTGCTCAAc                                                                                        <  8:113588/148‑1 (MQ=255)
                                                   ggCCGGCATCCACGAACAATTGGCTAAACCACAGGACGAGCTGCCGGGGCTGGAATGTTTGCATTATCCGATCGGCCTGGACCTGGGTGGCGACACCCCGGAAAGCGTGGCCATGGCGGTGCTGGCGGAAATTCAGGCAGTGCTCAACg                                                                                       >  3:107282/1‑149 (MQ=255)
                                                    gCCGGCATCCACGAACAATTGGCTAAACCACAGGACGAGCTGCCGGGGCTGGAATGTATGCATTATCCGATCGGCCTGGACCTGGGACGCGACACCCCGGAAAGCGTGGCCATGGCGGTGCTGTCGGAAATTCAGGCAGTGCTCAACgg                                                                                      >  4:257909/1‑149 (MQ=255)
                                                          aTCCACGAACAATTGGCTAAACCACAGGACGAGCTGCCGGGGCTGGAATGTTTGCATTATCCGATCGGCCTGGACCTGGGTGGCGACACCCCGGAAAGCGTGGCCATGGCGGTGCTGGCGGAAATTCAGGCAGTGCTCAACGGCCGCaa                                                                                >  2:287768/1‑149 (MQ=255)
                                                              aCGAACAATTGGCTAAACCACAGGACGAGCTGCCGGGGCTGGAATGTTTGCATTATCCGATCGGCCTGGACCTGGGTGGCGACACCCCGGAAAGCGTGGCCATGGCGGTGCTGGCGGAAATTCAGGCAGTGCTCAACGGCCGCAATggc                                                                            <  6:13437/149‑1 (MQ=255)
                                                                 aaCAATTGGCTAAACCACAGGACGAGCTGCCGGGGCTGGAATGTTTGCATTATCCGATCGGCCTGGACCTGGGTGGCGACACCCCGGAAAGCGTGGCCATGGCGGTGCTGGCGGAAATTCAGGCAGTGCTCAACGGCCGCAATggcggc                                                                         <  7:22051/149‑1 (MQ=255)
                                                                  aCAATTGGCTAAACCACAGGACGAGCTGCCGGGGCTGGAATGTTTGCATTATCCGATCGGCCTGGACCTGGGTGGCGACACCCCGGAAAGCGTGGCCATGGCGGTGCTGGCGGAAATTCAGGCAGTGCTCAACGGCCGCAATggcgg                                                                          >  2:64820/1‑147 (MQ=255)
                                                                             aaCCACAGGACGAGCTGCCGGGGCTTGAATGTTTGCATTATCCGATCGGCCTGGACCTGGGTGGCGACACCCCGGAAAGCGTGGCCATGGCGGTGCTGGCGGAAATTCAGGCAGTGCTCAACGGCCGCAATGGCGGCAGCCTGCGCtt                                                              <  4:107282/148‑1 (MQ=255)
                                                                             aaCCACAGGACGAGCTGCCGGGGCTGGAATGTTTGCATTATCCGATCGGCCTGGACCTGGGTGGCGACACCCCGGAAAGCGTGGCCATGGCGGTGCTGGCGGAAATTCAGGCAGTGCTCAACGGCCGCAATGGCGGCAGCCTGCGCTTc                                                             <  5:183534/149‑1 (MQ=255)
                                                                             aaCCACAGGACGAGCTGCCGGGGCTGGAATGTTTGCATTATCCGATCGGCCTGGACCTGGGTGGCGACACCCCGGAAAGCGTGGCCATGGCGGTGCTGGCGGAAATTCAGGCAGTGCTCAACGGCCGCAATGGCGGCAGCCTGCGCTTc                                                             <  5:236212/149‑1 (MQ=255)
                                                                             aaCCACAGGACGAGCTGCCGGGGCTGGAATGTTTGCATTATCCGATCGGCCTGGACCTGGGTGGCGACACCCCGGAAAGCGTGGCCATGGCGGTGCTGGCGGAAATTCAGGCAGTGCTCAACGGCCGCAATGGCGGCAGCCTGCGCTTc                                                             <  7:232626/149‑1 (MQ=255)
                                                                                 acaGGACGAGCTGCCGGGGCTGGAATGTTTGCATTATCCGATCGGCCTGGACCTGGGTGGCGACACCCCGGAAAGCGTGGCCATGGCGGTGCTGGCGGAAATTCAGGCAGTGCTCAACGGCCGCAATGGCGGCAGCCTGCGCtt                                                              >  1:41140/1‑144 (MQ=255)
                                                                                  caGGACGAGCTGCCGGGGCTGGAATGTTTGCATTATCCGATCGGCCTGGACCTGGGTGGCGACACCCCGGAAAGCGTGGCCATGGCGGTGCTGGCGGAAATTCAGGCAGTGCTCAACGGCCGCAATGGCGGCAGCCTGCGCTTCAGATc                                                        >  2:45043/1‑149 (MQ=255)
                                                                                     gACGAGCTGCCGGGGCTGGAATGTTTGCATTATCCGATCGGCCTGGACCTGGGTGGCGACACCCCGGAAAGCGTGGCCATGGCGGTGCTGGCGGAAATTCAGGCAGTGCTCAACGGCCGCAATGGCGGCAGCCTGCGCTTCAGATCTGc                                                     >  3:299026/1‑149 (MQ=255)
                                                                                       cGAGCTGCCGGGGCTGGAATGTTTGCATTATCCGATCGGCCTGGACCTGGGTGGCGACACCCCGGAAAGCGTGGCCATGGCGGTGCTGGCGGAAATTCAGGCAGTGCTCAACGGCCGCAATGGCGGCAGCCTGCGCTTCAGATCTGCGa                                                   <  1:205471/149‑1 (MQ=255)
                                                                                         aGCTGCCGGGGCTGGAATGTTTGCATTATCCGATCGGCCTGGACCTGGGTGGCGACACCCCGGAAAGCGTGGCCATGGCGGTGCTGGCGGAAATTCAGGCAGTGCTCAACGGCCGCAATGGCGGCAGCCTGCGCTTCAGATCTGCGAGc                                                 <  2:155081/149‑1 (MQ=255)
                                                                                               cGGGGCTGGAATGTTTGCATTATCCGATCGGCCTGGACCTGGGTGGCGACACCCCGGAAAGCGTGGCCATGGCGGTGCTGGCGGAAATTCAGGCAGTGCTCAACGGCCGCAATGGCGGCAGCCTGCGCTTCAGATCTGCGAGCATCCAt                                           <  7:45560/149‑1 (MQ=255)
                                                                                               cGGGGCTGGAATGTTTGCATTATCCGATCGGCCTGGACCTGGGTGGCGACACCCCGGAAAGCGTGGCCATGGCGGTGCTGGCGGAAATTCAGGCAGTGCTCAACGGCCGCAATGGCGGCAGCCTGCGCTTCAGATCTGCGAGCATCCAt                                           <  7:221011/149‑1 (MQ=255)
                                                                                                 gggCTGGAATGTTTGCATTATCCGATCGGCCTGGACCTGGGTGGCGACACCCCGGAAAGCGTGGCCATGGCGGTGCTGGCGGAAATTCAGGCAGTGCTCAACGGCCGCAATGGCGGCAGCCTGCGCTTCAGATCTGCGAGCATCCATGa                                         <  1:287768/149‑1 (MQ=255)
                                                                                                     tGGAATGTTTGCATTATCCGATCGGCCTGGACCTGGGTGGCGACACCCCGGAAAGCGTGGCCATGGCGGTGCTGGCGGAAATTCAGGCAGTGCTCAACGGCCGCAATGGCGGCAGCCTGCGCTTCAGATCTGCGAGCATCCATGACAGc                                     >  7:98958/1‑149 (MQ=255)
                                                                                                            tttGCATTATCCGATCGGCCTGGACCTGGGTGGCGACACCCCGGAAAGCGTGGCCATGGCGGTGCTGGCGGAAATTCAGGCAGTGCTCAACGGCCGCAATGGCGGCAGCCTGCGCTTCAGATCTGCGAGCATCCATGACAGCGACCCtg                              >  5:167322/1‑149 (MQ=255)
                                                                                                                cATTATCCGATCGGCCTGGACCTGGGTGGCGACACCCCGGAAAGCGTGGCCATGGCGGTGCTGGCGGAAATTCAGGCAGTGCTCAACGGCCGCAATGGCGGCAGCCTGCGCTTCAGATCTGCGAGCATCCATGACAGCGACCCTGTGCt                          >  5:30392/1‑149 (MQ=255)
                                                                                                                 aTTATCCGATCGGCCTGGACCTGGGTGGCGACACCCCGGAAAGCGTGGCCATGGCGGTGCTGGCGGAAATTCAGGCAGTGCTCAACGGCCGCAATGGCGGCAGCCTGCGCTTCAGATCTGCGAGCATCCATGACAGCGACCCTGTGCt                          >  4:294170/1‑148 (MQ=255)
                                                                                                                 aTTATCCGATCGGCCTGGACCTGGGTGGCGACACCCCGGAAAGCGTGGCCATGGCGGTGCTGGCGGAAATTCAGGCAGTGCTCAACGGCCGCAATGGCGGCAGCCTGCGCTTCAGATCTGCGAGCATCCATGACAGCGACCCTGTGCt                          >  4:168483/1‑148 (MQ=255)
                                                                                                                    tttcgttCGGGCTGGTCCGGGGTGGGTCCACCCCGGTATCCGGGGCCATGGGGGTGCTTGCGGAAATTCAGGCATTGCTCAACGGCCGCAATGGCGGCAGCCTGCGCTTCTGATCTGCGAGCATCCATGACAGCGACCCTGTGCTGaa                       <  8:49689/142‑1 (MQ=255)
                                                                                                                    aTCCGATCGGCCTGGACCTGGGTGGCGACACCCCGGAAAGCGTGGCCATGGCGGTGCTGGCGGAAATTCAGGCAGTGCTCAACGGCCGCAATGGCGGCAGCCTGCGCTTCAGATCTGCGAGCATCCATGACAGCGACCCTGTGCTGaa                       <  6:30392/148‑1 (MQ=255)
                                                                                                                             gCCTGGACCTGGGTGGCGACACCCCGGAAAGCGTGGCCATGGCGGTGCTGGCGGAAATTCAGGCAGTGCTCAACGGCCGCAATGGCGGCAGCCTGCGCTTCAGATCTGCGAGCATCCATGACAGCGACCCTGTGCTGAAGGCTGGCCgg             <  3:168483/149‑1 (MQ=255)
                                                                                                                                 tttccggggTGTCGCCACCCAGGAAAGCGTGGCCATGGCGGTGCTGGCGGAAATTCAGGCAGTGCTCAACGGTCGCAATGGCGGCAGCCTGCGCTTCAGATCTGCGAGCATCCATGACAGCGACCCtgtg                            <  4:193261/124‑1 (MQ=255)
                                                                                                                                 tttccggggTGTCGCCACCCAGGAAAGCGTGGCCATGGCGGTGCTGGCGGAAATTCAGGCAGTGCTCAACGGTCGCAATGGCGGCAGCCTGCGCTTCAGATCTGCGAGCATCCATGACAGCGACCCtgtg                            >  3:193261/7‑130 (MQ=255)
                                                                                                                                    ccTGGGTGGCGACACCCCGGAAAGCGTGGCCATGGCGGTGCTGGCGGAAATTCAGGCAGTGCTCAACGGCCGCAATGGCGGCAGCCTGCGCTTCAGATCTGCGAGCATCCATGACAGCGACCCTGTGCTGAAGGCTGGCCGGGATCaaa      >  7:90828/1‑149 (MQ=255)
                                                                                                                                     cTGGGTGGCGACACCCCGGAAAGCGTGGCCATGGCGGTGCTGGCGGAAATTCAGGCAGTGCTCAACGGCCGCAATGGCGGCAGCCTGCGCTTCAGATCTGCGAGCATCCATGACAGCGACCCTGTGCTGAAGGCTGGCCGGGATCAAAc     <  2:41140/149‑1 (MQ=255)
                                                                                                                                     cTGGGTGGCGACACCCCGGAAAGCGTGGCCATGGCGGTGCTGGCGGAAATTCAGGCAGTGCTCAACGGCCGCAATGGCGGCAGCCTGCGCTTCAGATCTGCGAGCATCCATGACAGCGACCCTGTGCTGAAGGCTGGCCGGGATCAAAc     <  2:193791/149‑1 (MQ=255)
                                                                                                                                        ggTGGCGACACCCCGGAAAGCGTGGCCATGGCGGTGCTGGCGGAAATTCAGGCAGTGCTCAACGGCCGCAATGGCGGCAGCCTGCGCTTCAGATCTGCGAGCATCCATGACAGCGACCCTGTGCTGAAGGCTGGCCGGGATCAAACGGa  <  8:98958/149‑1 (MQ=255)
                                                                                                                                           ggCGACACCCCGGAAAGCGTGGCCATGGCGGTGCTGGCGGAAATTCAGGCAGTGCTCAACGGCCGCAATGGCGGCAGCCTGCGCTTCAGATCTGCGAGCATCCATGACAGCGACCCTGTGCTGAAGGCTGGCCggg            <  4:105555/136‑1 (MQ=255)
                                                                                                                                           ggCGACACCCCGGAAAGCGTGGCCATGGCGGTGCTGGCGGAAATTCAGGCAGTGCTCAACGGCCGCAATGGCGGCAGCCTGCGCTTCAGATCTGCGAGCATCCATGACAGCGACCCTGTGCTGAAGGCTGGCCggg            >  3:105555/1‑136 (MQ=255)
                                                                                                                                               |                                                                                                                                             
ACCCTGCTATGTCGGTCAGTTGGGGCCGCGTGAGCGCACCGAGCGCTTGTTGGCCGGCATCCACGAACAATTGGCTAAACCACAGGACGAGCTGCCGGGGCTGGAATGTTTGCATTATCCGATCGGCCTGGACCTGGGTGGCGACACCCCGGAAAGCGTGGCCATGGCGGTGCTGGCGGAAATTCAGGCAGTGCTCAACGGCCGCAATGGCGGCAGCCTGCGCTTCAGATCTGCGAGCATCCATGACAGCGACCCTGTGCTGAAGGCTGGCCGGGATCAAACGGA  >  NC_002947/2828706‑2828990

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: