Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,876,733 G→A 7.2% R74C (CGT→TGT)  pucM ← 5‑hydroxyisourate hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,876,7330GA 7.2% 156.6 / 2.0 55R74C (CGT→TGT) pucM5‑hydroxyisourate hydrolase
Reads supporting (aligned to +/- strand):  ref base G (27/24);  new base A (2/2);  total (29/26)
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.68e-01

GACCAACTAGCTTCCACGATAGGTCGAATAGCTGTACGGCGAAATCAGCAGCGGCACGTGGTAGTGCTCCTGCGACTCGTCGATGCCAAAGCGCAGCACGACAACATCCAGAAACGCCTGGGCCGGCAATTGCACGCCACGGGCGCGGTAGTAATCCCCAGCGCTGAACTGCAACTGGTAAACGCCAGTGCGGTAATCGTCCCCCTGCAGCAGCGGCGCGTCGACACGGCCATCGCTGTTGGTCAGGGCGGTGTTCACCAGTTCCAGCTGCTGGCCTTCGACGC  >  NC_002947/4876593‑4876876
                                                                                                                                            |                                                                                                                                               
gACCAACTAGCTTCCACGATAGGTCGAATAGCTGTACGGCGAAATCAGCAGCGGCACGTGGTAGTGCTCCTGCGACTCGTCGATGCCAAAGCGCAGCACGACAACATCCAGAAACGCCTGGGCCGGCAATTGCACGCCACGGGCGCgg                                                                                                                                          <  7:170063/148‑1 (MQ=255)
 aCCAACTAGCTTCCACGATAGGTCGAATAGCTGTACGGCGAAATCAGCAGCGGCACGTGGTAGTGCTCCTGCGACTCGTCGATGCCAAAGCGCAGCACGACAACATCCAGAAACGCCTGGGCCGGCAATTGCACGCCACGGGCGCGgta                                                                                                                                        <  7:43117/149‑1 (MQ=255)
 aCCAACTAGCTTCCACGATAGGTCGAATAGCTGTACGGCGAAATCAGCAGCGGCACGTGGTAGTGCTCCTGCGACTCGTCGATGCCAAAGCGCAGCACGACAACATCCAGAAACGCCTGGGCCGGCAATTGCACGCCACGGGCGCGgta                                                                                                                                        <  3:138649/149‑1 (MQ=255)
       tAGCTTCCACGATAGGTCGAATAGCTGTACGGCGAAATCAGCAGCGGCACGTGGTAGTGCTCCTGCGACTCGTCGATGCCAAAGCGCAGCACGACAACATCCAGAAACGCCTGGGCCGGCAATTGCACGCCACGGGCGCGGTAGTAATc                                                                                                                                  >  4:239342/1‑149 (MQ=255)
         gCTTCCACGATAGGTCGAATAGCTGTACGGCGAAATCAGCAGCGGCACGTGGTAGTGCTCCTGCGACTCGTCGATGCCAAAGCGCAGCACGACAACATCCAGAAACGCCTGGGCCGGCAATTGCACGCCACGGGCGCGGTAGTAATccc                                                                                                                                >  7:199694/1‑149 (MQ=255)
               aCGATAGGTCGAATAGCTGTACGGCGAAATCAGCAGCGGCACGTGGTAGTGCTCCTGCGACTCGTCGATGCCAAAGCGCAGCACGACAACATCCAGAAACGCCTGGGCCGGCAATTGCACGCCACGGGCGCGGTAGTAATCCCCAgcgc                                                                                                                          >  4:77615/1‑149 (MQ=255)
                    aGGTCGAATAGCTGTACGGCGAAATCAGCAGCGGCACGTGGTAGTGCTCCTGCGACTCGTCGATGCCAAAGCGCAGCACGACAACATCCAGAAACGCCTGGGCCGGCAATTACACGCCACGGGCGCGGTAGTAATCCCCAGCGCTGAAc                                                                                                                     >  7:103394/1‑149 (MQ=255)
                          aaTAGCTGTACGGCGAAATCAGCAGCGGCACGTGGTAGTGCTCCTGCGACTCGTCGATGCCAAAGCGCAGCACGACAACATCCAGAAACGCCTGGGCCGGCAATTACACGCCACGGGCGCGGTAGTAATCCCCAGCGCTGAACTGCAAc                                                                                                               <  8:103394/149‑1 (MQ=255)
                           aTAGCTGTACGGCGAAATCAGCAGCGGCACGTGGTAGTGCTCCTGCGACTCGTCGATGCCAAAGCGCAGCACGACAACATCCAGAAACGCCTGGGCCGGCAATTGCACGCCACGGGCGCGGTAGTAATCCCCAGCGCTGAACTGCAACt                                                                                                              >  5:214527/1‑149 (MQ=255)
                            tAGCTGTACGGCGAAATCAGCAGCGGCACGTGGTAGTGCTCCTGCGACTCGTCGATGCCAAAGCGCAGCACGACAACATCCAGAAACGCCTGGGCCGGCAATTGCACGCCACGGGCGCGGTAGTAATCCCCAGCGCTGAACTGCAACTg                                                                                                             <  8:199694/149‑1 (MQ=255)
                            tAGCTGTACGGCGAAATCAGCAGCGGCACGTGGTAGTGCTCCTGCGACTCGTCGATGCCAAAGCGCAGCACGACAACATCCAGAAACGCCTGGGCCGGCAATTGCACGCCACGGGCGCGGTAGTAATCCCCAGCGCTGAACTGCAACTg                                                                                                             <  7:117587/149‑1 (MQ=255)
                                tGTACGGCGAAATCAGCAGCGGCACGTGGTAGTGCTCCTGCGACTCGTCGATGCCAAAGCGCAGCACGACAACATCCAGAAACGCCTGGGCCGGCAATTGCACGCCACGGGCGCGGTAGTAAt                                                                                                                                   >  8:318463/1‑123 (MQ=255)
                                tGTACGGCGAAATCAGCAGCGGCACGTGGTAGTGCTCCTGCGACTCGTCGATGCCAAAGCGCAGCACGACAACATCCAGAAACGCCTGGGCCGGCAATTGCACGCCACGGGCGCGGTAGTAAt                                                                                                                                   <  7:318463/123‑1 (MQ=255)
                                      gCGAAATCAGCAGCGGCACGTGGTAGTGCTCCTGCGACTCGTCGATGCCAAAGCGCAGCACGACAACATCCAGAAACGCCTGGGCCGGCAATTGCACGCCACGGGCGCGGTAGTAATCCCCAGCGCTGAACTGCAACTGGTAAACGc                                                                                                     >  5:69344/1‑147 (MQ=255)
                                       cGAAATCAGCAGCGGCACGTGGTAGTGCTCCTGCGACTCGTCGATGCCAAAGCGCAGCACGACAACATCCAGAAACGCCTGGGCCGGCAATTGCACGCCACGGGCGCGGTAGTAATCCCCAGCGCTGAACTGCAACTGGTAAACGCCAg                                                                                                  >  3:327137/1‑149 (MQ=255)
                                       cGAAATCAGCAGCGGCACGTGGTAGTGCTCCTGCGACTCGTCGATGCCAAAGCGCAGCACGACAACATCCAGAAACGCCTGGGCCGGCAATTGCACGCCACGGGCGCGGTAGAAATCCCCAGCGCTGAACTGc                                                                                                                  <  4:38871/133‑1 (MQ=255)
                                       cGAAATCAGCAGCGGCACGTGGTAGTGCTCCTGCGACTCGTCGATGCCAAAGCGCAGCACGACAACATCCAGAAACGCCTGGGCCGGCAATTGCACGCCACGGGCGCGGTAGAAATCCCCAGCGCTGAACTGc                                                                                                                  >  3:38871/1‑133 (MQ=255)
                                         aaaTCAGCAGCGGCACGTGGTAGTGCTCCTGCGACTCGTCGATGCCAAAGCGCAGCACGACAACATCCAGAAACGCCTGGGCCGGCAATTGCACGCCACGGGCGCGGTAGTAATCCCCAGCGCTGAACTGCAACTGGTAAACGCCAGTg                                                                                                <  3:239342/149‑1 (MQ=255)
                                         aaaTCAGCAGCGGCACGTGGTAGTGCTCCTGCGACTCGTCGATGCCAAAGCGCAGCACGACAACATCCAGAAACGCCTGGGCCGGCAATTGCACGCCACGGGCGCGGTAGTAATCCCCAGCGCTGAACTGCAACTGGTAAACGCCAGTg                                                                                                <  2:98933/149‑1 (MQ=255)
                                          aaTCAGCAGCGGCACGTGGTAGTGCTCCTGCGACTCGTCGATGCCAAAGCGCAGCACGACAACATCCAGAAACGCCTGGGCCGGCAATTGCACGCCACGGGCGCGGTAGTAATCCCCAGCGCTGAACTg                                                                                                                   >  2:210298/1‑129 (MQ=255)
                                          aaTCAGCAGCGGCACGTGGTAGTGCTCCTGCGACTCGTCGATGCCAAAGCGCAGCACGACAACATCCAGAAACGCCTGGGCCGGCAATTGCACGCCACGGGCGCGGTAGTAATCCCCAGCGCTGAACTg                                                                                                                   <  1:210298/129‑1 (MQ=255)
                                             cagcagCGGCACGTGGTAGTGCTCCTGCGACTCGTCGATGCCAAAGCGCAGCACGACAACATCCAGAAACGCCTGGGCCGGCAATTGCACGCCACGGGCGCgg                                                                                                                                          <  6:224437/103‑1 (MQ=255)
                                             cagcagCGGCACGTGGTAGTGCTCCTGCGACTCGTCGATGCCAAAGCGCAGCACGACAACATCCAGAAACGCCTGGGCCGGCAATTGCACGCCACGGGCGCgg                                                                                                                                          >  5:224437/1‑103 (MQ=255)
                                                 agcGGCACGTGGTAGTGCTCCTGCGACTCGTCGATGCCAAAGCGCAGCACGACAACATCCAGAAACGCCTGGGCCGGCAATTGCACGCCACGGGCGCGGTAGTAATCCCCAGCGCTGTACTGCAACTGGTAAACGCCAGTGCGGTAATc                                                                                        <  4:327137/149‑1 (MQ=255)
                                                    ggCACGTGGTAGTGCTCCTGCGACTCGTCGATGCCAAAGCGCAGCACGACAACATCCAGAAACGCCTGGGCCGGCAATTGCACGCCACGGGCGCGGTAGTAATCCCCAGCGCTGAACTGCAACTGGTAAACGCCAGTGCGGTAATCGTc                                                                                     <  6:214527/149‑1 (MQ=255)
                                                            gTAGTGCTCCTGCGACTCGTCGATGCCAAAGCGCAGCACGACAACATCCAGAAACGCCTGGGCCGGCAATTGCACGCCACGGGCGCGGTAGTAATCCCCAGCGCTGAACTGCAACTGGTAAACGCCAGTGCGGTAATCGTCCCCCTgca                                                                             <  8:89551/149‑1 (MQ=255)
                                                                        cGACTCGTCGATGCCAAAGCGCAGCACGACAACATCCAGAAACGCCTGGGCCGGCAATTGCACGCCACGGGCGCgg                                                                                                                                          >  3:2247/1‑76 (MQ=255)
                                                                        cGACTCGTCGATGCCAAAGCGCAGCACGACAACATCCAGAAACGCCTGGGCCGGCAATTGCACGCCACGGGCGCgg                                                                                                                                          <  4:2247/76‑1 (MQ=255)
                                                                           cTCGGCGATGCCAAAGCGCAGCACGACAACATCCAGAAACGCCTGGGCCGGCAATTGCACGCCACGGGCGCGGTAGTAATCCCCAGCGCTGAACTGCAACTGGTAAACGCCAGTGCGGTAATCGTCCCCCTGCAGCAGCGGCGCGTCGa                                                              <  6:69344/149‑1 (MQ=255)
                                                                               tcgATGCCAAAGCGCAGCACGACAACATCCAGAAACGCCTGGGCCGGCAATTGCACGCCACGGGCGCGGTAGTAATCCCCAGCGCTGAACTGCAACTGGTAAACGCCAGTGCGGTAATCGTCCCCCTGCAGCAGCGGCGCGTCGACACg                                                          >  5:230399/1‑149 (MQ=255)
                                                                               tcgATGCCAAAGCGCAGCACGACAACATCCAGAAACGCCTGGGCCGGCAATTGCACGCCACGGGCGCGGTAGTAATCCCCAGCGCTGAACTGCAACTGGTAAACGCCAGTGCGGTAATCGTCCCCCTGCAGCAGCGGCGCGTCGACACg                                                          >  7:149943/1‑149 (MQ=255)
                                                                                    gCCAAAGCGCAGCACGACAACATCCAGAAACGCCTGGGCCGGCAATTGCACGCCACGGGCGCGGTAGTAATCCCCAGCGCTGAACTGCAACTGGTAAACGCCAGTGCGGTAATCGTCCCCCTGCAGCAGCGGCGCGTCGACACGGCCAt                                                     >  8:149143/1‑149 (MQ=255)
                                                                                    gCCAAAGCGCAGCACGACAACATCCAGAAACGCCTGGGCCGGCAATTGCACGCCACAGGCGCGGTAGTAATCCCCAGCGCTGAACTGCAACTGGTAAACGCCAGTGCGGTAATCGTCCCCCTGCAGCAGCGGCGCGTCGACACGGCCAt                                                     >  4:254858/1‑149 (MQ=255)
                                                                                    gCCAAAGCGCAGCACGACAACATCCAGAAACGCCTGGGCCGGCAATTGCACGCCACAGGCGCGGTAGTAATCCCCAGCGCTGAACTGCAACTGGTAAACGCCAGTGCGGTAATCGTCCCCCTGCAGCAGCGGCGCGTCGACACGGCCAt                                                     >  4:250435/1‑149 (MQ=255)
                                                                                     ccAAAGCGCAGCACGACAACATCCAGAAACGCCTGGGCCGGCAATTGCACGCCACGGGCGCGGTAGTAATCCCCAGCGCTGAACTGCAACTGGTAAACGCCAGTGCGGTAATCGTCCCCCTGCAGCAGCGGCGCGTCGa                                                              <  3:33224/139‑1 (MQ=255)
                                                                                     ccAAAGCGCAGCACGACAACATCCAGAAACGCCTGGGCCGGCAATTGCACGCCACGGGCGCGGTAGTAATCCCCAGCGCTGAACTGCAACTGGTAAACGCCAGTGCGGTAATCGTCCCCCTGCAGCAGCGGCGCGTCGa                                                              >  4:33224/1‑139 (MQ=255)
                                                                                      cAAAGCGCAGCACGACAACATCCAGAAACGCCTGGGCCGGCAATTGCACGCCACGGGCGCGGTAGTAATCCCCAGCGCTGAACTGCAACTGGTAAACGCCAGTGCGGTAATCGTCCCCCTGCAGCAGCGGCGCGTCGACACGGCCATc                                                    >  4:157745/1‑148 (MQ=255)
                                                                                       aaaGCGCAGCACGACAACATCCAGAAACGCCTGGGCCGGCAATTGCACGCCACGGGCGCGGTAGTAATCCCCAGCGCTGAACTGCAACTGGTAAACGCCAGTGCGGTAATCGTCCCCCTGCAGCAGCGGCGCGTCGACACGGCCATCGc                                                  >  5:256341/1‑149 (MQ=255)
                                                                                       aaaGCGCAGCACGACAACATCCAGAAACGCCTGGGCCGGCAATTGCACGCCACGGGCGCGGTAGTAATCCCCAGCGCTGAACTGCAACTGGTAAACGCCAGTGCGGTAATCGTCCCCCTGCAGCAGCGGCGCGTCGACACGGCCATCGc                                                  >  7:250186/1‑149 (MQ=255)
                                                                                       aaaGCGCAGCACGACAACATCCAGAAACGCCTGGGCCGGCAATTGCACGCCACAGGCGCGGTAGTAATCCCCAGCGCTGAACTGCAACTGGTAAACGCCAGTGCGGTAATCGTCCCCCTGCAGCAGCGGCGCGTCGACACGGCCATCGc                                                  <  3:254858/149‑1 (MQ=255)
                                                                                       aaaGCGCAGCACGACAACATCCAGAAACGCCTGGGCCGGCAATTGCACGCCACAGGCGCGGTAGTAATCCCCAGCGCTGAACTGCAACTGGTAAACGCCAGTGCGGTAATCGTCCCCCTGCAGCAGCGGCGCGTCGACACGGCCATCGc                                                  <  3:250435/149‑1 (MQ=255)
                                                                                              agcaCGACAACATCCAGAAACGCCTGGGCCGGCAATTGCACGCCACGGGCGCGGTAGTAATCCCCAGCGCTGAACTGCAACTGGTAAACGCCAGTGCGGTAATCGTCCCCCTGCAGCAGCGGCGCGTCGACACGGCCATCGCTGTTGGt                                           <  8:250186/149‑1 (MQ=255)
                                                                                              agcaCGACAACATCCAGAAACGCCTGGGCCGGCAATTGCACGCCACGGGCGCGGTAGTAATCCCCAGCGCTGAACTGCAACTGGTAAACGCCAGTGCGGTAATCGTCCCCCTGCAGCAGCGGCGCGTCGACACGGCCATCGCTGTTGGt                                           >  7:293462/1‑149 (MQ=255)
                                                                                                  cGACAACATCCAGAAACGCCTGGGCCGGCAATTGCACGCCACGGGCGCGGTAGTAATCCCCAGCGCTGAACTGCAACTGGTAAACGCCAGTGCGGTAATCGTCCCCCTGCAGCAGCGGCGCGTCGACACGGCCATCGCTGTTGGTCAgg                                       <  8:149943/149‑1 (MQ=255)
                                                                                                       acaTCCAGAAACGCCTGGGCCGGCAATTGCACGCCACGGGCGCGGTAGTAATCCCCAGCGCTGAACTGCAACTGGTAAACGCCAGTGCGGTAATCGTCCCCCTGCAGCAGCGGCGCGTCGACACGGCCATCGCTGTTGGTCAGGGCGgt                                  <  7:149143/149‑1 (MQ=255)
                                                                                                          tCCAGAAACGCCTGGGCCGGCAATTGCACGCCACGGGCGCGGTAGTAATCCCCAGCGCTGAACTGCAACTGGTAAACGCCAGTGCGGTAATCGTCCCCCTGCAGCAGCGGCGCGTCGACACGGCCATCGCTGTTGGTCAGGGCGGTGtt                               >  5:158208/1‑149 (MQ=255)
                                                                                                            cAGAAACGCCTGGGCCGGCAATTGCACGCCACGGGCGCGGTAGTAATCCCCAGCGCTGAACTGCaa                                                                                                                >  8:206822/1‑66 (MQ=255)
                                                                                                            cAGAAACGCCTGGGCCGGCAATTGCACGCCACGGGCGCGGTAGTAATCCCCAGCGCTGAACTGCAACTGGTAAACGCCAGTGCGGTAATCGTCCCCCTGCAGCAGCGGCGCGTCGACACGGCCATCGCTGTTGGTCAGGGCGGTGTTc                              <  3:77615/148‑1 (MQ=255)
                                                                                                                     cTGGGCCGGCAATTGCACGCCACGGGCGCGGTAGTAATCCCCAGCGCTGAACTGCAACTGGTAAACGCCAGTGCGGTAATCGTCCCCCTGCAGCAGCGGCGCGTCGACACGGCCATCGCTGTTGGTCAGGGCGGTGTTCACCAGTTc                      >  2:36145/1‑147 (MQ=255)
                                                                                                                                  tGCACGCCACGGGCGCGGTAGTAATCCCCAGCGCTGAACTGCAACTGGTAAACGCCAGTGCGGTAATCGTCCCCCTGCAGCAGCGGCGCGTCGACACGGCCATCGCTGTTGGTCAGGGCGGTGTTCACCAGTTCCAGCTGCTGGCCTTc       >  3:68499/1‑149 (MQ=255)
                                                                                                                                   gCACGCCACGGGCGCGGTAGTAATCCCCAGCGCTGAACTGCAACTGGTAAACGCCAGTGCGGTAATCGTCCCCCTGCAGCAGCGGCGCGTCGACACGGCCATCGCTGTTGGTCAGGGCGGTGTTCACCAGTTCCAgctgct              <  2:63261/141‑1 (MQ=255)
                                                                                                                                   gCACGCCACGGGCGCGGTAGTAATCCCCAGCGCTGAACTGCAACTGGTAAACGCCAGTGCGGTAATCGTCCCCCTGCAGCAGCGGCGCGTCGACACGGCCATCGCTGTTGGTCAGGGCGGTGTTCACCAGTTCCAgctgct              >  1:63261/1‑141 (MQ=255)
                                                                                                                                    cACGCCACGGGCGCGGTAGTAATCCCCAGCGCTGAACTGCAACTGGTAAACGCCAGTGCGGTAATCGTCCCCCTGCAGCAGCGGCGCGTCGACACGGCCATCGCTGTTGGTCAGGGCGGTGTTCACCAGTTCCAGCTGCTGGCCTTCGa     >  2:161887/1‑149 (MQ=255)
                                                                                                                                       gCCACGGGCGCGGTAGTAATCCCCAGCGCTGAACTGCAACTGGTAAACGCCAGTGCGGTAATCGTCCCCCTGCAGCAGCGGCGCGTCGACACGGCCATCGCTGTTGGTCAGGGCGGTGTTCACCAGTTCCAGCTGCTGGCCTTCGAcgc  >  6:150192/1‑149 (MQ=255)
                                                                                                                                       gCCACGGGCGCGGTAGTAATCCCCAGCGCTGAACTGCAACTGGTAAACGCCAGTGCGGTAATCGTCCCCCTGCAGCAGCGGCGCGTCGACACGGCCATCGCTGTTGGTCAGGGCGGTGTTCACCAGTTCCAGCTGCTGGCCTTCGAcgc  <  1:161887/149‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                               
GACCAACTAGCTTCCACGATAGGTCGAATAGCTGTACGGCGAAATCAGCAGCGGCACGTGGTAGTGCTCCTGCGACTCGTCGATGCCAAAGCGCAGCACGACAACATCCAGAAACGCCTGGGCCGGCAATTGCACGCCACGGGCGCGGTAGTAATCCCCAGCGCTGAACTGCAACTGGTAAACGCCAGTGCGGTAATCGTCCCCCTGCAGCAGCGGCGCGTCGACACGGCCATCGCTGTTGGTCAGGGCGGTGTTCACCAGTTCCAGCTGCTGGCCTTCGACGC  >  NC_002947/4876593‑4876876

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: ‑