Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 568,254 G→C 8.2% R414R (CGC→CGG)  uvrA ← uvrABC excinuclease protein A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947568,2540GC 8.2% 139.6 / 4.7 49R414R (CGC→CGG) uvrAuvrABC excinuclease protein A
Reads supporting (aligned to +/- strand):  ref base G (24/21);  new base C (2/2);  total (26/23)
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.45e-01

GCAGATTTCCTTGAGGATTTTTGCCGCGATTTCGCCACGCCGGCCAGTAAGGGTCAGTTGACCAAAGTAGTTGCTGGCTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGC‑AGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTCGCGGTAGCGGCGCTCAAGGTTCGGCACAATGCCTTCGAACGGGTGCGATCGCTTGACGATG  >  NC_002947/568107‑568401
                                                                                                                                                   |                                                                                                                                                    
gCAGATTTCCTTGAGGATTTTTGCCGCGATTTCGCCACGCCGGCCAGTAAGGGTCAGTTGACCAAAGTAGTTGCTGGCTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCgcggc                                                                                                                                                     >  2:206540/1‑149 (MQ=255)
  aGATTTCCTTGAGGATTTTTGCCGCGATTTCGCCACGCCCGCCAGTAACGGTCAGATGACCAAAGTAGTTGATGGCTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCgcggcg                                                                                                                                                    >  3:134569/1‑148 (MQ=255)
       tCCTTGAGGATTTTTGCCGCGATTTCGCCACGCCGGCCAGTAAGGGTCAGTTGACCAAAGTAGTTGCTGGCTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGC‑AGgc                                                                                                                                              <  4:102592/148‑1 (MQ=255)
         cTTGAGGATTTTTGCCGCGATTTCGCCACGCCGGCCAGTAAGGGTCAGTTGACCAAAGTAGTTGCTGGCTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGC‑AGgcgc                                                                                                                                            <  2:76349/148‑1 (MQ=255)
            gAGGATTTTTGCCGCGATTTCGCCACGCCGGCCAGTAAGGGTCAGTTGACCAAAGTAGTTGCTGGCTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGC‑AGGCGCGt                                                                                                                                          <  2:251680/147‑1 (MQ=255)
             aGGATTTTTGCCGCGATTTCGCCACGCCGGCCAGTAAGGGTCAGTTGACCAAAGTAGTTGCTGGCTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCT‑GCGCCGCGAGGcgcgcca                                                                                                                                        <  4:244322/148‑4 (MQ=255)
                 tttttGCCGCGATTTCGCCACGCCGGCCAGTAAGGGTCAGTTGACCAAAGTAGTTGCTGGCTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGC‑AGGCGCGTACCCCgg                                                                                                                                   >  1:125865/1‑149 (MQ=255)
                  ttttGCCGCGATTTCGCCACGCCTGCCAGTAAGGTTCAGTTGACCAATGTAGTTGCTGGCTTCACGGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCTCCTCGCGGCGC‑AGGCGCGTACCCAgg                                                                                                                                   <  4:134569/148‑1 (MQ=255)
                             tttCGCCACGCCGGCCAGTAAGGGTCAGTTGACCAAAGTAGTTGCTGGCTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGC‑AGGCGCGTACCCCGGCAATCCGGGCAc                                                                                                                       >  5:163150/1‑149 (MQ=255)
                                cGCCACGCCGGCCAGTAAGGGTCAGTTGACCAAAGTAGTTGCTGGCTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGC‑AGgcg                                                                                                                                             <  6:157548/124‑1 (MQ=255)
                                cGCCACGCCGGCCAGTAAGGGTCAGTTGACCAAAGTAGTTGCTGGCTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGC‑AGgcg                                                                                                                                             >  5:157548/1‑124 (MQ=255)
                                 gCCACGCCGGCCAGTAAGGGTCAGTTGACCAAAGTAGTTGCTGGCTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGC‑AGGCGCGTACCCCGGc                                                                                                                                  <  1:97400/134‑1 (MQ=255)
                                 gCCACGCCGGCCAGTAAGGGTCAGTTGACCAAAGTAGTTGCTGGCTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGC‑AGGCGCGTACCCCGGc                                                                                                                                  >  2:97400/1‑134 (MQ=255)
                                  ccACGCCGGCCAGTAAGGGTCAGTTGACCAAAGTAGTTGCTGGCTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGTGC‑AGGCGCGTACCCCGGCAATCCGGGCACGGTTg                                                                                                                  >  7:120355/1‑149 (MQ=255)
                                          gCCAGTAAGGGTCAGTTGACCAAAGTAGTTGCTGGCTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGTGC‑AGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCaa                                                                                                          <  8:120355/149‑1 (MQ=255)
                                             aGTAAGGGTCAGTTGACCAAAGTAGTTGCTGGCTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGC‑AGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAga                                                                                                        >  1:95990/1‑147 (MQ=255)
                                                       aGTTGACCAAAGTAGTTGCTGGCTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGT‑AGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAg                                                                                             <  1:206540/149‑1 (MQ=255)
                                                       aGTTGACCAAAGTAGTTGCTGGCTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGC‑AGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGTCCAg                                                                                             >  6:256106/1‑149 (MQ=255)
                                                            aCCAAAGTAGTTGCTGGCTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGC‑AGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGttct                                                                                         >  4:99332/1‑148 (MQ=255)
                                                             ccAAAGTAGTTGCTGGCTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGC‑AGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGttctt                                                                                        >  3:194427/1‑148 (MQ=255)
                                                              cAAAGTAGTTGCTGGCTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGC‑AGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTca                                                                                      >  8:37595/1‑149 (MQ=255)
                                                              cAAAGTAGTTGCTGGCTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGC‑AGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTca                                                                                      <  2:28351/149‑1 (MQ=255)
                                                               aaaGTAGTTGCTGGCTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGC‑AGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTcac                                                                                     >  1:36470/1‑149 (MQ=255)
                                                                  gtagtTGCTGGCTTCACCGATCGGCAAATTGGGTACCGCCGGCAGGGTTTTCTCGCCCACCCACAGATGGCGCGCCTCGCGGCGC‑AGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTcacgca                                                                                  <  1:177786/149‑1 (MQ=255)
                                                                                 aCCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGC‑AGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGc                                                                                                             >  1:33579/1‑107 (MQ=255)
                                                                                 aCCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGC‑AGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGc                                                                                                             <  2:33579/107‑1 (MQ=255)
                                                                                 aCCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGC‑AGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCgg                                                                   >  3:121712/1‑149 (MQ=255)
                                                                                 aCCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGC‑AGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGGGGCCGATTCgg                                                                   >  8:208831/1‑149 (MQ=255)
                                                                                         gCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGC‑AGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTCGCgg                                                           >  6:183523/1‑149 (MQ=255)
                                                                                         gCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGC‑AGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTCGCgg                                                           >  4:144269/1‑149 (MQ=255)
                                                                                          cAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGC‑AGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTCGCgg                                                           >  6:12141/1‑148 (MQ=255)
                                                                                                 gTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGC‑AGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTCGCGGTAgcggcg                                                   <  5:12141/149‑1 (MQ=255)
                                                                                                  ttACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGC‑AGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTCGCGGTAGCGgcgc                                                  <  1:72680/149‑1 (MQ=255)
                                                                                                      cgccgGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGC‑AGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTCGCGGTAGCGGCGCTCaa                                              <  4:194427/149‑1 (MQ=255)
                                                                                                         cgGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGC‑AGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTCGCGGTAGCGGCGCTCAAGGt                                           <  7:37595/149‑1 (MQ=255)
                                                                                                               ggTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGC‑AGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTCGCGGTAGCGGCGCTCAAGGTTCGGca                                     >  5:252406/1‑149 (MQ=255)
                                                                                                                         cccacccaCACATGGCGCGCCTCGCGGCGC‑AGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTCGCGGTAGCGGCGCTCAAGGTTCGGCACAATGCCTTc                           >  4:24436/1‑149 (MQ=255)
                                                                                                                         cccacccaCACATGGCGCGCCTCGCGGCGC‑AGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTCGCGGTAGCGGCGCTCAAGGTTCGGCACAATGCCTTc                           <  1:72658/149‑1 (MQ=255)
                                                                                                                         cccacccaCACATGGCGCGCCTCGCGGCGC‑AGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCAAGCACGTTGGCCGATTCGGTTTCGCTGTCTCGGCGCTCAAGCTACGGCACACTGCCTTc                           >  8:114303/1‑149 (MQ=255)
                                                                                                                              ccacACATGGCGCGCCTCGCGGCGC‑AGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTCGCGGTAGCGGCGCTCAAGGTTCGGCACAATGCCTTCGAACg                      >  2:115421/1‑149 (MQ=255)
                                                                                                                               cacacaTGGCGCGCCTCGCGGCGC‑AGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTCGCGGTAGCGGCGCTCAAGGTTCGGCACAATGCCTTCGAACgg                     <  3:144269/149‑1 (MQ=255)
                                                                                                                                acacaTGGCGCGCCTCGCGGCGC‑AGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCTCGGTGGCCGATTCGGTTTCGCGGTAGCGGAGCTCAAGGTTCGGCACAATGCCTTCGAACggg                    >  6:188890/1‑149 (MQ=255)
                                                                                                                                acacaTGGCGCGCCTCGCGGCGC‑AGGAGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTCGCGGTAGCGGCGCTCAAGGTTCGGCACAATGCCTTCGAACggg                    >  1:11399/1‑149 (MQ=255)
                                                                                                                                 cacaTGGCGCGCCTCGCGGCGC‑AGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTCGCGGTAGCGGCGCTCAAGGTTCGGCACAATGCCTTCGAACGGGt                   <  3:99332/149‑1 (MQ=255)
                                                                                                                                  acaTGGCGCGCCTCGCGGCGC‑AGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTCGCGGTAGCGGCGCTCAAGGTTCGGCACAATGCCTTCGAACGGGTg                  <  2:36470/149‑1 (MQ=255)
                                                                                                                                   gggtacgcgcCT‑GCGCCGCGAGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGAGTCGGTTTCGCGGTAGCGGCGCTCAAGGTTCGGCACAATGCCTTCGAACGGGTGc                 >  7:147705/6‑149 (MQ=255)
                                                                                                                                    ggtacgcgcCT‑GCGCCGCGAGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTCGCGGTAGCGGCGCTCAAGGTTCGGCACAATGCCTTCGAAc                       <  7:125732/138‑1 (MQ=255)
                                                                                                                                    ggtacgcgcCT‑GCGCCGCGAGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTCGCGGTAGCGGCGCTCAAGGTTCGGCACAATGCCTTCGAAc                       >  8:125732/5‑142 (MQ=255)
                                                                                                                                    aTGGCGCGCCTCGCGGCGC‑AGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTCGCGGTAGCGGCGCTCAAGGTTCGGCACAATGCCTTCGAACGGGTGc                 >  1:169521/1‑148 (MQ=255)
                                                                                                                                            ccTCGCGGCGC‑AGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTCGCGGTAGCGGCGCTCAAGGTTCGGCACAATGCCTTCGAACGGGTGCGATCGCTTg        <  5:256106/149‑1 (MQ=255)
                                                                                                                                             cTCGCGGCGC‑AGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTCGCGGTAGCGGCGCTCAAGGTTCGGCACAATGCCTTCGAACGGGTGCGATCGCTTGa       <  6:163150/149‑1 (MQ=255)
                                                                                                                                                  ggcgC‑AGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTCGCGGTAGCGGCGCTCAAGGTTCGGCACAATGCCTTCGAACGGGTGCGATCGCTTGACGATg  >  2:192383/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GCAGATTTCCTTGAGGATTTTTGCCGCGATTTCGCCACGCCGGCCAGTAAGGGTCAGTTGACCAAAGTAGTTGCTGGCTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGC‑AGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTCGCGGTAGCGGCGCTCAAGGTTCGGCACAATGCCTTCGAACGGGTGCGATCGCTTGACGATG  >  NC_002947/568107‑568401

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