Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,246,437 G→C 9.2% G231A (GGC→GCC)  PP_2841 → site‑specific recombinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,246,4370GC 9.2% 154.4 / 3.1 51G231A (GGC→GCC) PP_2841site‑specific recombinase
Reads supporting (aligned to +/- strand):  ref base G (20/26);  new base C (3/2);  total (23/28)
Fisher's exact test for biased strand distribution p-value = 6.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGCAGGCCTGGCTTGAACTGCGCCGCCAGGACCTGCCTGCCGGGGCGGAAGATGACCCGTTCCTGTTCAACCGCATCCGTCGCGGCAACC‑ACATTACCCGGGCACGCATTACCAAGCATGCGATCTACTACATTGCCCGCCAGCGCGGCGCGCAGGTGGGGGTGAAGATCATGCCGCACGATTTTCGCCGGGCATTCATTACCCGGGTGATCGAAGAGCACGACCTGTCGATTGCGCAGAAGCTGGCGCACCATGCCAATATCCAGACCACGGCCATCTATGACCGGCGTGAT  >  NC_002947/3246289‑3246582
                                                                                                                                                     |                                                                                                                                                 
cTGCAGGCCTGGCTTGAACTGCGCCGCCAGGACCTGCCTGCCGGGGCGGAAGATGACCCGTTCCTGTTCAACCGCATCCGTCGCGGCAACC‑ACATTACCCGGGCACGCATTACCAAGCCTGCGATCTACTACATTGCCCGCCAgcgcgg                                                                                                                                                   >  7:328393/1‑149 (MQ=255)
  gCAGGCCTGGCTTGAACTGCGCCGCCAGGACCTGCCTGCCGGGGCGGAAGATGACCCGTTCCTGTTCAACCGCATCCGTCGCGGCAACC‑ACATTACCCGGGCACGCATTACCAAGCATGCGATCTACTACATTGCCCGCCAgcgcggcg                                                                                                                                                 >  3:41161/1‑149 (MQ=255)
  gCAGGCCTGGCTTGAACTGCGCCGCCAGGACCTGCCTGCCGGGGCGGAAGATGACCCGTTCCTGTTCAACCGCATCCGTCGCGGCAACC‑ACATTACCCGGGCACGCATTACCAAGCATGCGATCTACTACATTGCCCGCCAgcgcggcg                                                                                                                                                 >  3:142889/1‑149 (MQ=255)
     ggCCTGGCTTGAACTGCGCCGCCAGGACCTGCCTGCCGGGGCGGAAGATGACCCGTTCCTGTTCAACCGCATCCGTCGCGGCAACC‑ACATTACCCGGGCACGCATTACCAAGCATGCGATCTACTACATTGCCCGCCAgcgcggcgcg                                                                                                                                               >  4:133234/1‑148 (MQ=255)
        cTGGCTTGAACTGCGCCGCCAGGACCTGCCTGCCGGGGCGGGAGATGACCCGTTCCTGTTCAACCACATCCGTCGCGGCAACCTCCATTACCCGGGCACGCATTACCAAGCATGCGATCTACTACATTGCCCGCCAGCGCGGCGCGCAg                                                                                                                                            <  1:111454/149‑1 (MQ=255)
           gCTTGAACTGCGCCGCCAGGACCTGCCTGCCGGGGCGGAAGATGACCCGTTCCTGTTCAACCGCATCCGTCGCGGCAACC‑ACATTACCCGGGCACGCATTACCAAGCATGCGATCTACTACATTGCCCGCCAGCGCGGCGCGCAGGt                                                                                                                                          <  8:85936/147‑1 (MQ=255)
           gCTTGAACTGCGCCGCCAGGACCTGCCTGCCGGGGCGGAAGATGACCCGTTCCTGTTCAACCGCATCCGTCGCGGCAACC‑ACATTACCCGGGCACGCATTACCAAGCATGCGATCTACTACATTGCCCACCTGCGCGCCGCGCTGGCgg                                                                                                                                        <  6:13114/149‑1 (MQ=255)
           gCTTGAACTGCGCCGCCAGGACCTGCCTGCCGGGGCGGAAGATGACCCGTTCCTGTGCAACCGCATCCGTCGCGGCAACC‑ACATTACCCGGGCACGCATTACCAAGCATGCGATCTACTACATTGCCCGCCAGCGCGGCGCGCAGGt                                                                                                                                          <  4:53943/147‑1 (MQ=255)
                  cTGCGCCGCCAGGACCTGCCTGCCGGGGCGGAAGATGACCCGTTCCTGTTCAACCGCATCCGTCGCGGCAACC‑ACATTACCCGGGCACGCATTACCAAGCATGCGATCTACTACATTGCCCGCCAGCGCGGCGCGCAGGTGGGGGTGaa                                                                                                                                 <  4:174925/149‑1 (MQ=255)
                       gcgccAGGACCTGCCTGCCGGGGCGGAAGATGACCCTCTCCTCTTCACCCGCCTCCCTCGCCGCAACC‑ACATTCCCCGGCCACGCACTACCAAGCAAGAGATCTCCTACAATGCCCGCCTCCCCGCGGCGCAGGTGGTGGTGAACATCa                                                                                                                            <  6:54998/148‑1 (MQ=255)
                       ccgccAGGACCTGCCTGCCGGGGCGGAAGATGACCCGTTCCTGTTCAACCGCATCCGTCGCGGCAACC‑ACATTACCCGGGCACGCATTACCAAGCATGCGATCTACTACATTGCCCGCCAGCGCGGCGCGCAGGTGGGGGTGAAGATCa                                                                                                                            >  4:307409/1‑149 (MQ=255)
                          ccAGGACCTGCCTGCCGGGGCGGAAGATGACCCGTTCCTGTTTAACCGCATCCGTCGCGGCAACC‑ACATTACCCGGGCACGCATTACCAAGCATGCGATCTACTACATTGCCCGCCAGCGCGGCGCGCAGGTGGGGGTGAAGATCATGc                                                                                                                         <  1:124829/149‑1 (MQ=255)
                           cAGGACCTGCCTGCCGGGGCGGAAGATGACCCGTTCCTGTTCAACCGCATCCGTCGCGGCAACC‑ACATTACCCGGGCACGCATTACCAAGCATGCGATCTACTACATTGCCCGCCAGCGCGGCGCGCAGGTGGGGGTGAAGATCATGcc                                                                                                                        <  4:41161/149‑1 (MQ=255)
                                      tgccGGGGCGGAAGATGACCCGTTCCTGTTCAACCGCATCCGTCGCGGCAACC‑ACATTACCCGGGCACGCATTACCAAGCATGCGATCTACTACATTGCCCGCCAGCGCGGCGCGCAGGTGGGGGTGAAGATCATGCCGCACGATTTTc                                                                                                             >  4:232658/1‑149 (MQ=255)
                                        ccGGGGCGGAAGATGACCCGTTCCTGTTCAACCGCATCCGTCGCGGCAACC‑ACATTACCCGGGCACGCATTACCAAGCATGCGATCTACTACATTGCCCGCCAGCTCGGCGCGCAGGTGGGGGTGAAGATCATGCCTCACGATTTTc                                                                                                             <  6:214414/147‑1 (MQ=255)
                                         cGGGGCGGAAGATGACCCGTTCCTGTTCAACCGCATCCGTCGCGGCAACC‑ACATTACCCGGGCACGCATTACCAAGCATGCGATCTACTACATTGCCCGCCAGCGCGGCGCGCAGGTGGGGGTGAAGATCATGCCGCACGATTTTCGcc                                                                                                          >  2:325581/1‑149 (MQ=255)
                                          ggggCGGAAGATGACCCGTTCCTGTTCAACCGCATCCGTCGCGGCAACC‑ACATTACCCGGGCACGCATTACCAAGCATGCGATCTACTACATTGCCCGCCAGCGCGGCGCGCAGGTGGGGGTGAAGATCATGCCGCACGATTTTCGCCg                                                                                                         <  3:232658/149‑1 (MQ=255)
                                           gggCGGAAGATGACCCGTTCCTGTTCAACCGCATCCGTCGCGGCAACC‑ACATTACCCGGGCACGCATTACCAAGCATGCGATCTACTACATTGCCCGCCAGCGCGGCGCGCAGGTGGGGGTGAAGATCATGCCGCACGATTTTCGCCgg                                                                                                        <  4:288241/149‑1 (MQ=255)
                                                     tGACCCGTTCCTGTTCAACCGCATCCGTCGCGGCAACC‑ACATTACCCGGGCACGCATTACCAAGCATGCGATCTACTACATTGCCCGCCAGCGCGGCGCGCAGGTGCGGGTGAAGATCATGCCGCACGATTTTCGCCGGGCATTCATTa                                                                                              >  1:155725/1‑149 (MQ=255)
                                                         ccGTTCCTGTTCAACCGCATCCGTCGCGGCAACC‑ACATTACCCGGGCACGCATTACCAAGCATGCGATCTACTACATTGCCCGCCAGCGCGGCGCGCAGGTGGGGGTGAAGATCATGCCGCACGATTTTCGCCGGGCATTCATTACCCg                                                                                          <  2:155725/149‑1 (MQ=255)
                                                          cGTTCCTGTTCAACCGCATCCGTCGCGGCAACC‑ACATTACCCGGGCACGCATTACCAAGCATGCGATCTACTACATTGCCCGCCAGCGCGGCGCGCAGGTGGGGGTGAAGATCATGCCGCACGATTTTCGCCGGGCATTCATTACCCgg                                                                                         <  6:106714/149‑1 (MQ=255)
                                                            ttCCTGTTCCACCGCATCCGTCGCGTCAACC‑ACATTCCCCGGGCACGCATTACCAAGCATGCGATCTACTACATTGCCCGCCAGCGCGGCGCGCAGGTGGGGGTGAAGATCATGCCGCACGATTTTCGCCGGGCATTCATTACCCGGGt                                                                                       <  8:328393/149‑1 (MQ=255)
                                                                     aaCCGCATCCGTCGCGGCAACC‑ACATTACCCGGGCACGCATTACCAAGCATGCGATCTACTACATTGCCCGCCAGCGCGGCGCGCAGGTGGGGGTGAAGATCATGCCGCACGATTTTCGCCGGGCATTCATTAccc                                                                                           <  5:207225/136‑1 (MQ=255)
                                                                     aaCCGCATCCGTCGCGGCAACC‑ACATTACCCGGGCACGCATTACCAAGCATGCGATCTACTACATTGCCCGCCAGCGCGGCGCGCAGGTGGGGGTGAAGATCATGCCGCACGATTTTCGCCGGGCATTCATTAccc                                                                                           >  6:207225/1‑136 (MQ=255)
                                                                     aaCCGCATCCGTCGCGGCAACC‑ACATTACCCGGGCACGCATTACCAAGCATGCGATCTACTACATTGCCCGCCAGCGCGGCGCGCAGGTGGGGGTGAAGATCATGCCGCACGATTTTCGCCGGGCATTCATTACCCGGGTGATCGAaga                                                                              <  4:260157/149‑1 (MQ=255)
                                                                        cccATCCGTCGCGGCAACC‑ACATTACCCGGGCACGCATTACCAAGCATGCGATCTACTACATTGCCCGCCAGCGCGGCGCGCAGGTGGGGGTGAAGATCATGCCGCACGATTTTCGCCGGGCATTCATTACCCGGGTGATCGAAGAGCa                                                                           >  4:136935/3‑149 (MQ=255)
                                                                         gCATCCGTCGCGGCAACC‑ACATTACCCGGGCACGCATTACCAAGCATGCGATCTACTACATTGCCCGCCAGCGCGGCGCGCAGGTGGGGGTGAAGATCATGCCGCACGATTTTCGCCGGGCATTCATTACCCGGGTGATCGAAGAGCAc                                                                          <  3:307409/149‑1 (MQ=255)
                                                                         gCATCCGTCGCGGCAACC‑ACATTACCCGGGCACGCATTACCAAGCATGCGATCTACTACATTGCCCGCCAGCGCGGCGCGCAGGTGGGGGTGAAGATCATGCCGCACGATTTTCGCCGGGCATTCATTACCCGGGTGATCGAAGAGCAc                                                                          >  6:210554/1‑149 (MQ=255)
                                                                           aTCCGTCGCGGCAACC‑ACATTACCCGGGCACGCATTACCAAGCATGCGATCTACTACATTGCCCGCCAGCGCGGCGCGCAGGTGGGGGTGAAGATCATGCCGCACGATTTTCGCCGGGCATTCATTACCCGGGTGATCGAAGAGCACGa                                                                        <  4:342975/149‑1 (MQ=255)
                                                                                          c‑acaTTACCCGGGCACGCATTACCAAGCATGCGATCTACTACATTGCCCGCCAGCGCGGCGCGCAGGTGGGGGTGAAGATCATGCCGCACGATTTTCGCCGGGCATTCATTACCCGGGTGATCGAaga                                                                              <  3:13078/128‑1 (MQ=255)
                                                                                          c‑acaTTACCCGGGCACGCATTACCAAGCATGCGATCTACTACATTGCCCGCCAGCGCGGCGCGCAGGTGGGGGTGAAGATCATGCCGCACGATTTTCGCCGGGCATTCATTACCCGGGTGATCGAaga                                                                              >  4:13078/1‑128 (MQ=255)
                                                                                             caTTACCCGGGCACGCATTACCAAGCATGCGATCTACTACATTGCCCGCCAGCGCGGCGCGCAGGTGGGGGTGAAGATCATGCCGCACGATTTTCGCCGGGCATTCATTACCCGGGTGATCGAAGAGCACGACCTGTCGATTGCGCAg                                                        >  4:314585/1‑148 (MQ=255)
                                                                                                 aCCCGGGCACGCATTACCAAGCATGCGATCTACTACATTGCCCGCCAGCGCGGCGCGCAGGTGGGGGTGAAGATCATGCCGCACGATTTTCGCCGGGCATTCATTACCCGGGTGATCGAAGAGCACGACCTGTCGATTGCGCAGAAGc                                                    <  1:325581/148‑1 (MQ=255)
                                                                                                    cGGGCACGCATTACCAAGCATGCGATCTACTACATTGCCCGCCAGCGCGGCGCGCAGGTGGGGGTGAAGATCATGCCGCACGATTTTCGCCGGGCATTCATTACCCGGGTGATCGAAGAGCACGACCTGTCGATTGCGCAGAAGCTGgc                                                <  3:133234/149‑1 (MQ=255)
                                                                                                        cACGCATTACCAAGCATGCGATCTACTACATTGCCCGCCAGCGCGGCGCGCAGGTGGGGGTGAAGATCATGCCGCACGATTTTCGCCGGGCATTCATTACCCGGGTGATCGAAGAGCACGACCTGTCGATTGCGCAGAAGCTGGCGCAc                                            >  6:269910/1‑149 (MQ=255)
                                                                                                                   aaGCATGCGATCTACTACATTGCCCGCCAGCGCGGCGCGCAGGTGGGGGTGAAGATCATGCCGCACGATTTTCGCCGGGCATTCATTACCCGGGTGATCGAAGAGCACGACCTGTCGATTGCGCAGAAGCTGGCGCACCATGc                                       <  8:263681/143‑1 (MQ=255)
                                                                                                                   aaGCATGCGATCTACTACATTGCCAGCCAGCGCGCCGCGCACGTGGGGGTGAAGATCATGCCGCACGATTTTCGCCGGGCATTCATTACCCGGGTGATCGAAGAGCACGACCTGTCGATTGCGCAGAAGCTGGCGCACCATGc                                       >  7:263681/1‑143 (MQ=255)
                                                                                                                               tactacATTGCCCGCCAGCGCGGCGCGCAGGTGGGGGTGAAGATCATGCCGCACGATTTTCGCCGGGCATTCATTACCCGGGTGATCGAAGAGCACGACCTGTCGATTGCGCAGAAGCTGGCGCACCATGCCAATATCCAGACCACGGc                     >  5:128546/1‑149 (MQ=255)
                                                                                                                               tactacATTGCCCGCCAGCGCGGCGCGCAGGTGGGGGTGAAGATCATGCCGCACGATTTTCGCCGGGCATTCATTACACGGGTGATCGAAGAGCACGACCTGTCGATTGCGCAGAAGCTGGCGCACCATGCCAATATCCAGACCCCGGc                     >  8:101850/1‑149 (MQ=255)
                                                                                                                                 ctacATTGCCCGCCAGCGCGGCGCGCAGGTGGGGGTGAAGATCATGCCGCACGATTTTCGCCGGGCATTCATTACCCGGGTGATCGAAGAGCACGACCTGTCGATTGCGCAGAAGCTGGCGCACCATGCCAATATCCAGACCACGGcc                    >  8:282095/1‑148 (MQ=255)
                                                                                                                                 ctacATTGCCCGCCAGCGCGGCGCGCAGGTGGGGGTGAAGATCATGCCGCACGATTTTCGCCGGGCATTCATTACCCGGGTGATCGAAGAGCACGACCTGTCGATTGCGCAGAAGCTGGCGCACCATGCCAATATCCAGACCAAGGCCa                   >  3:218581/1‑149 (MQ=255)
                                                                                                                                  tacATTGCCCGCCAGCGCGGCGCGCAGGTGGGGGTGAAGATCATGCCGCACGATTTTCGCCGGGCATTCATTACCCGGGTGATCGAAGAGCACGACCTGTCGATTGCGCAGAAGCTGGCGCACCATGCCAATATCCAGACCACGGCCAt                  >  2:215757/1‑149 (MQ=255)
                                                                                                                                    cATTGCCCGCCAGCGCGGCGCGCAGGTGGGGGTGAAGATCATGCCGCACGATTTTCGCCGGGCATTCATTACCCGGGTGATCGAAGAGCACGACCTGTCGATTGCGCAGAAGCTGGCGCACCATGCCAATATCCAGACCACGGCCATc                 <  5:269910/148‑1 (MQ=255)
                                                                                                                                     aTTGCCCGCCAGCGCGGCGCGCAGGTGGGGGTGAAGATCATGCCGCACGATTTTCGCCGGGCATTCATTACCCGGGTGATCGAAGAGCACGACCTGTCGATTGCGCAGAAGCTGGCGCACCATGCCAATATCCAGACCACGGCCATCTa               >  8:169809/1‑149 (MQ=255)
                                                                                                                                      ttGCCCGCCAGCGCGGCGCGCAGGTGGGGGTGAAGATCATGCCGCACGATTTTCGCCgg                                                                                                        >  2:324465/1‑59 (MQ=255)
                                                                                                                                      ttGCCCGCCAGCGCGGCGCGCAGGTGGGGGTGAAGATCATGCCGCACGATTTTCGCCgg                                                                                                        <  1:324465/59‑1 (MQ=255)
                                                                                                                                       tGCCCGCCAGCGCGGCGCGCAGGTGGGGGTGAAGATCATGCCGCACGATTTTCGCCGGGCATTCATTACCCGGGTGATCGAAGAGCACGACCTGTCGATTGCGCAGAAGCTGGCGCACCATGCCAATATCCAGACCACGGCCATCTAt              >  7:187103/1‑148 (MQ=255)
                                                                                                                                       cccccACCTGCGCGCCGCGCTGGCGGGGGTGAAGATCATGCCGCACGATTTTCGCCGGGCATTCATTACCCGGGTGATCGAAGAGCACGACCTGTCGATTGCGCAGAAGCTGGCGCACCATGCCAATATCCAGACCACGGCCATCTATg             >  5:83798/3‑149 (MQ=255)
                                                                                                                                        ccccACCTGCGCGCCGCGCTGGCGGGGGTGAAGATCATGCCGCACGATTTTCGCCGGGCATTCATTACCCGGGTGATCGAAGAGCACGACCTGTCGATTGCGCAGAAGCTGGCGCACCATGCCAATATCCAGACCACGGCCATCTATGa            >  5:173384/2‑149 (MQ=255)
                                                                                                                                         cccGCCAGCGCGGCGCGCAGGTGGGGGTGAAGATCATGCCGCACGATTTTCGCCGGGCATTCATTACCCGGGTGATCGAAGAGCACGACCTGTCGATTGCGCAGAAGCTGGCGCACCATGCCAATATCCAGACCACGGCCATCTATGAc           <  7:37382/149‑1 (MQ=255)
                                                                                                                                          ccGCCAGCGCGGCGCGCAGGTGGGGGTGAAGATCATGCCGCACGATTTTCGCCGGGCATTCATTACCCGGGTGATCGAAGAGCACGACCTGTCGATTGCGCAGAAGCTGGCGCACCATGCCAATATCCAGACCACGGCCATCTATGAcc          >  5:308579/1‑149 (MQ=255)
                                                                                                                                           cGCCAGCGCGGCGCGCAGGTGGGGGTGAAGATCATGCCGCACGATTTTCGCCGGGCATTCATTACCCGGGTGATCGAAGAGCACGACCTGTCGATTGCGCAGAAGCTGGCGCACCATGCCAATATCCAGACCACGGCCATCTATGACCg         >  2:117263/1‑149 (MQ=255)
                                                                                                                                             ccAGCGCGGCGCGCAGGTGGGGGTGAAGATCATGCCGCACGATTTTCGCCGGGCATTCATTACCCGGGTGATCGAAGAGCACGACCTGTCGATTGCGCAGAAGCTGGCGCACCATGCCAATATCCAGACCACGGCCATCTATGACCGGc       <  4:142889/149‑1 (MQ=255)
                                                                                                                                               aGCGCGGCGCGCAGGTGGGGGTGAAGATCATGCCGCACGATTTTCGCCGGGCATTCATTACCCGGGTGATCGAAGAGCACGACCTGTCGATTGCGCAGAAGCTGGCGCACCATGCCAATATCCAGACCACGGCCATCTATGACCGGCGt     <  8:187103/149‑1 (MQ=255)
                                                                                                                                               aGCGCGGCGCGCAGGTGGGGGTGAAGATCATGCCGCACGATTTTCGCCGGGCATTCATTACCCGGGTGATCGAAGAGCACGACCTGTCGATTGCGCAGAAGCTGGCGCACCATGCCAATATCCAGACCACGGCCATCTATGACCGGCGt     <  3:314585/149‑1 (MQ=255)
                                                                                                                                                  gcggcgcgCAGGTGGGGGTGAAGATCATGCCGCACGATTTTCGCCGGGCATTCATTACCCGGGTGATCGAAGAGCACGACCTGTCGATTGCGCAGAAGCTGGCGCACCATGCCAATATCCAGACCACGGCCATCTATGACCGGCGtgat  <  1:215757/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
CTGCAGGCCTGGCTTGAACTGCGCCGCCAGGACCTGCCTGCCGGGGCGGAAGATGACCCGTTCCTGTTCAACCGCATCCGTCGCGGCAACC‑ACATTACCCGGGCACGCATTACCAAGCATGCGATCTACTACATTGCCCGCCAGCGCGGCGCGCAGGTGGGGGTGAAGATCATGCCGCACGATTTTCGCCGGGCATTCATTACCCGGGTGATCGAAGAGCACGACCTGTCGATTGCGCAGAAGCTGGCGCACCATGCCAATATCCAGACCACGGCCATCTATGACCGGCGTGAT  >  NC_002947/3246289‑3246582

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