Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,851,693 A→G 7.8% intergenic (+85/+34) apt → / ← recR adenine phosphoribosyltransferase/recombination protein RecR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,851,6930AG 7.8% 139.5 / 2.9 51intergenic (+85/+34)apt/recRadenine phosphoribosyltransferase/recombination protein RecR
Reads supporting (aligned to +/- strand):  ref base A (23/24);  new base G (2/2);  total (25/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 = 4.26e-01

TCGGTTGCAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGTGGCACACCATGGGCAATGCGAGACGCCACCAAG  >  NC_002947/4851545‑4851838
                                                                                                                                                    |                                                                                                                                                 
tCGGTTGCAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGGATGTTTTGTAGGCGCGGGTTTa                                                                                                                                                   >  8:220405/1‑149 (MQ=255)
      gCAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCcgc                                                                                                                                              <  4:47497/148‑1 (MQ=255)
        aaGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGCAAGGCGCGTATGGTTTGTAGGCGCGGGTTTACCCGCGaa                                                                                                                                           >  3:7220/1‑149 (MQ=255)
        aaGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGaa                                                                                                                                           <  7:69064/149‑1 (MQ=255)
                  ggTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCa                                                                                                                                 <  7:220405/149‑1 (MQ=255)
                          cctgcagCCTCACCGAGTCTTCCCTCAGCGAGTATTGAGGGGGCTTTGGGCGCTGCTGTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGc                                                                                                                           <  8:232864/141‑1 (MQ=38)
                          ccTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGccc                                                                                                                         <  4:125614/149‑1 (MQ=255)
                            ttCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGcccc                                                                                                                        >  5:4948/1‑148 (MQ=255)
                                  cTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAAt                                                                                                                 <  8:16180/149‑1 (MQ=255)
                                   tGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAAt                                                                                                                 >  1:305922/1‑148 (MQ=255)
                                     aCCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTAGCGGGTAATCCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTAACCGCGAATGGGACGCCAAGCGGCCCCAATGAAt                                                                                                                 >  5:27842/1‑146 (MQ=255)
                                       cGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTGCCCGCGAATgg                                                                                                                                        <  5:88041/121‑1 (MQ=255)
                                       cGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCAGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGGTTGCCCGCGAATgg                                                                                                                                        >  6:88041/1‑121 (MQ=255)
                                                       gAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGgcgc                                                                                            >  3:118717/1‑149 (MQ=255)
                                                        aGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGcc                                                                                           <  1:38104/149‑1 (MQ=255)
                                                            ttGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTg                                                                                        >  8:287429/1‑148 (MQ=255)
                                                                   ggCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTGCCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGc                                                                                                  <  2:139749/131‑1 (MQ=255)
                                                                   ggCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTGCCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGc                                                                                                  >  1:139749/1‑131 (MQ=255)
                                                                      ttttGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATgg                                                                             >  5:213955/1‑149 (MQ=255)
                                                                           gggCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGGGTATGTTTTGTAGGCGCGGGATTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCCGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAg                                                                        >  8:205126/1‑149 (MQ=255)
                                                                             gctgctTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAggg                                                                      >  3:275584/1‑149 (MQ=255)
                                                                                 ctTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTa                                                                    <  1:194779/147‑1 (MQ=255)
                                                                                        gCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTcc                                                           <  4:275584/149‑1 (MQ=255)
                                                                                         cccATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCCGGGTACCGCCGTcca                                                          >  4:314236/1‑149 (MQ=255)
                                                                                              tCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCAtc                                                     >  1:172224/1‑149 (MQ=255)
                                                                                                gcgGGTAAACGCGCTCCTACGAAAGCCGCGTTTGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCAtctc                                                   <  4:218290/149‑1 (MQ=255)
                                                                                                gcgGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCAtctc                                                   <  6:213955/149‑1 (MQ=255)
                                                                                                   ggTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTACAAAGGCATGGGGCAGGGTACCGCCGTCCACCATCTCCAg                                                >  1:310498/1‑149 (MQ=255)
                                                                                                      aaaCCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAg                                                <  6:321842/146‑1 (MQ=255)
                                                                                                      aaaCCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAg                                                >  5:321842/1‑146 (MQ=255)
                                                                                                          gatttCCTACGATATCCGCTTATGTTTTTTAGGCGATGGTTTACCCGCGTATGGGCCGCGAAGGGGCGCCTATTAATCAGAGCGAAATGGGCCGGCGCCCTGCAAATGCATGGGCCAGGGTACCGCCGTCCACCTTCTCCAGCTcgccg                                         <  6:260059/145‑1 (MQ=255)
                                                                                                             ctacTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGccc                                      <  4:7220/146‑1 (MQ=255)
                                                                                                                cTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCAcgccgcagatc                                   >  5:260059/1‑144 (MQ=255)
                                                                                                                  aCGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGTgg                                 >  3:282146/1‑149 (MQ=255)
                                                                                                                   cGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGTgg                                 <  4:282146/148‑1 (MQ=255)
                                                                                                                    gAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGTGGca                               >  4:254684/1‑149 (MQ=255)
                                                                                                                    gAAAGCCGCGTATGTATTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGCGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGCGTACCGCCGTCCACCATCTCCAGCTCGCCCCCCAGTGGca                               >  8:207007/1‑149 (MQ=255)
                                                                                                                       aGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGTGGCacac                            >  3:255887/1‑149 (MQ=255)
                                                                                                                         ccGCGTATGTTTTGTAGGCGCTGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCAtc                                                     <  8:228637/122‑1 (MQ=255)
                                                                                                                         ccGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCAtc                                                     >  7:228637/1‑122 (MQ=255)
                                                                                                                         ccGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGTGGCACACCa                          >  8:24389/1‑149 (MQ=255)
                                                                                                                          cgcgTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGTGGCACACCAt                         <  4:255887/149‑1 (MQ=255)
                                                                                                                           gcgTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGTGGCACACCATg                        <  6:4948/149‑1 (MQ=255)
                                                                                                                               aTGTGTTGTAGGCGCTGGTTTACCCGCGAATGGGCCGACAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCACGCCGCCCAGTGGCACACCATGGGc                     <  6:27842/148‑1 (MQ=255)
                                                                                                                                      gTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGTGGCACACCATGGGCAATGCga              <  7:287429/148‑1 (MQ=255)
                                                                                                                                        aGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGTGGCACACCATGGGCAATGCgaga            <  2:305922/148‑1 (MQ=255)
                                                                                                                                             cgGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTccacc                                                        <  1:327381/99‑1 (MQ=255)
                                                                                                                                             cgGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTccacc                                                        >  2:327381/1‑99 (MQ=255)
                                                                                                                                                gTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGTGGCACACCATGGGCAATGCGAGACGccacc     >  4:264555/1‑147 (MQ=255)
                                                                                                                                                 tttACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGTGGCACACCATGGGCAATGCGAGACGCCACCAAg  <  4:118717/149‑1 (MQ=255)
                                                                                                                                                  ttACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGTGGCACACCATggg                      <  6:296727/128‑1 (MQ=255)
                                                                                                                                                  ttACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGTGGCACACCATggg                      >  5:296727/1‑128 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
TCGGTTGCAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGTGGCACACCATGGGCAATGCGAGACGCCACCAAG  >  NC_002947/4851545‑4851838

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