Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,234,948 G→C 8.9% A39A (GCC→GCG)  PP_3711 ← diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,234,9480GC 8.9% 172.2 / 5.6 56A39A (GCC→GCG) PP_3711diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base G (24/27);  new base C (2/3);  total (26/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

CGCCCGTTTCTTCAAAGCGTCCGGTGCTGCGCGCCTGCTCCAGGTTGAATGCCGGCAACCCGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGATCAACAGCCGGTACATGATGTCGGGGTTTGCTGCGAGCGACATGAACTTTACCTGAAC  >  NC_002947/4234806‑4235079
                                                                                                                                              |                                                                                                                                   
cgcCCGTTTCTTCAAAGCGTCCGGTGCTGCGCGCCTGCTCCAGGTTGAATGCCGGCAACCCGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCa                                                                                                                               <  3:342522/149‑1 (MQ=255)
cgcCCGTTTCTTCAAAGCGTCCGGTGCTGCGCGCCTGCTCCAGGTTGAATGCCGGCAACCCGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCa                                                                                                                               <  2:144655/149‑1 (MQ=255)
  cccGTTTCTTCAAAGCGTCCGGTGCTGCGCGCCTGCTCCAGGTTGAATGCCGGCAACCCGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCCCAGCGCGCCaag                                                                                                                             <  7:37474/149‑3 (MQ=255)
    cGTTTCTTCAAAGCGTCCGGTGCTGCGCGCCTGCTCCAGGTTGAATGCCGGCAACCCGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCCTGGCGCCCAGCGCGCCAaggg                                                                                                                           <  1:117372/149‑2 (MQ=255)
          ttcAAAGCGTCCGGTGCTGCGCGCCTGCTCCAGGTTGAATGCCGGCAACCCGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGt                                                                                                                     >  8:209740/1‑149 (MQ=255)
           tcAAAGCGTCCGGTGCTGCGCGCCTGCTACAGGTTGAATGCCTGCAACCAGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGAAGATTTCATCGGAACGGTACCCCTTGGCGCGCTGGgcgc                                                                                                                                 >  5:194315/1‑136 (MQ=255)
           tcAAAGCGTCAGGTTCGGCGCGCCTGCTCCAGGTTTAATGCCGGCAAACCTGCTGAGCTCTCGGAATCGGTGTAGAACAAGGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGgcgc                                                                                                                                 <  6:194315/136‑1 (MQ=255)
            cAAAGCGTCCGGTGCTGCGCGCCTGCTCCAGGTTGAATGCCGGCAACCCGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCACCAGGGTTCCAGTTg                                                                                                                   >  1:4402/1‑149 (MQ=255)
              aaGCGTCCGGTGCTGCGCGCCTGCTCCAGGTTGAATGCCGGCAACCCTGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGc                                                                                                                 <  5:18997/149‑1 (MQ=255)
                   tCCGGTGCTGCGCGCCTGCTCCAGGTTGAATGCCGGCAACCCGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCACCAGGGTTCCAGTTGGcgacga                                                                                                            <  2:4402/149‑1 (MQ=255)
                           tgcGCGCCTGCTCCAGGTTGAATGCCGGCAACCCGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCtt                                                                                                      >  4:77323/1‑147 (MQ=255)
                                  cTGCTCCAGGTTGAATGCCGGCAACCCGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGa                                                                                             <  3:77323/149‑1 (MQ=255)
                                  cTGCTCCAGGTTGAATGCCGGCAACCCGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGa                                                                                             <  7:209740/149‑1 (MQ=255)
                                     cTCCAGGTTGAATGCCGGCAACCCGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAg                                                                                            <  1:284968/147‑1 (MQ=255)
                                     cTCCAGGTTGAATGCCGGCAACCCGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCCCAGCGCGCCAaggg                                                                                                                           >  7:233029/1‑115 (MQ=255)
                                     cTCCAGGTTGAATGCCGGCAACCCGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCCCAGCGCGCCAaggg                                                                                                                           <  8:233029/116‑2 (MQ=255)
                                        cAGGTTGAATGCCGGCAACCCGGGTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAg                                                                                         <  8:310174/147‑1 (MQ=255)
                                        cAGGTTGAATGCCGGCAACCCGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGt                                                                                       <  4:117662/149‑1 (MQ=255)
                                              gAATGCCGGCAACCCGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATg                                                                                  <  4:66774/148‑1 (MQ=255)
                                              gAATGCCGGCAACCCGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATg                                                                                  >  3:66774/1‑148 (MQ=255)
                                                aTGCCGGCAACCCGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCa                                                                               >  6:22144/1‑149 (MQ=255)
                                                       cAACCCGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATcca                                                                        >  1:21384/1‑149 (MQ=255)
                                                          cccGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATccacca                                                                     <  5:22144/149‑1 (MQ=255)
                                                           ccGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCcaccac                                                                    >  1:81134/1‑149 (MQ=255)
                                                                   gcgcTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGa                                                            <  5:178833/149‑1 (MQ=255)
                                                                     gcTCGGCATCGGTGTAGAACAACGATTATTGCTTGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTTAGCAGGTAGATgg                                                                                 <  6:351397/126‑1 (MQ=255)
                                                                     gcTCGGCATCGGTGTAGAACAACGAGTAGTGCTGTCCGACGATTTACTAGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGAt                                                           >  6:316333/1‑148 (MQ=255)
                                                                     gcTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAg                                                                                                   <  4:187475/108‑1 (MQ=255)
                                                                     gcTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGACGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATgg                                                                                 >  5:351397/1‑126 (MQ=255)
                                                                           cATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGa                                                                                    <  3:178328/117‑1 (MQ=255)
                                                                           cATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGa                                                                                    >  4:178328/1‑117 (MQ=255)
                                                                              cGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGATCAACAGc                                                    <  7:109628/146‑1 (MQ=255)
                                                                              cGGTGTAGAACAACGAGTAGTGCTGGCCCACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGATCAACAGc                                                    >  8:109628/1‑146 (MQ=255)
                                                                                     gAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGATCAACAGCCGGTACatga                                          >  2:242785/1‑149 (MQ=255)
                                                                                         aacGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGa                                                            >  5:186151/1‑127 (MQ=255)
                                                                                         aacGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGa                                                            <  6:186151/127‑1 (MQ=255)
                                                                                              gtagtGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTgg                                                             >  1:6969/1‑121 (MQ=255)
                                                                                              gtagtGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTgg                                                             <  2:6969/121‑1 (MQ=255)
                                                                                                 gtGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGATCAACAGCCGGTACATGATGTCGGGGttt                               <  1:242785/148‑1 (MQ=255)
                                                                                                     gtgCCGAGGAAATCGTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTCGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATACACCACGCTCTGGATCAACAGCCGGTACATGATGTCGGGGTTtgctgc                          >  4:44117/3‑149 (MQ=255)
                                                                                                      tgCCGAGGAAATCGTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCGGCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGATCACCAGCCGGTACATGATGTCGGGGTTTGCTgcg                         >  8:246026/2‑149 (MQ=255)
                                                                                                      ggCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGATCAACAGCCGGTACATGATGTCGGGGTTtgctgc                          >  5:355113/1‑148 (MQ=255)
                                                                                                       gCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGATCAACAGCCGGTACATGATGTCGGGGTTTGCTgcga                        <  5:316333/149‑1 (MQ=255)
                                                                                                        ccGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGATCAACAGCCGGTACATGATGTCGGGGTTtgctgc                          >  4:92944/1‑146 (MQ=255)
                                                                                                        ccGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGATCAACAGCCGGTACATGATGTCGGGGTTTGCTgcga                        <  3:92944/148‑1 (MQ=255)
                                                                                                         cgacgaTTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGATCAACAGCCGGTACATGATGTCGGGGTTTGCTgcgagc                      >  2:97397/1‑149 (MQ=255)
                                                                                                            cgaTTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGATCAACAGCCGGTACATGATGTCGGGGTTTGCTGCGAGCGAc                   <  2:21384/149‑1 (MQ=255)
                                                                                                             gaTTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGATCAACAGCCGGTACATGATGTCGGGGTTTGCTGCGAGCGACa                  >  5:341507/1‑149 (MQ=255)
                                                                                                                  ccTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGATCAACAGCCGGTACATGATGTCGGGGTTTGCTGCGAGCGACATGaa              <  6:355113/148‑1 (MQ=255)
                                                                                                                             tACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGAt                                                                                   >  7:197040/1‑68 (MQ=255)
                                                                                                                             tACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGAt                                                                                   <  8:197040/68‑1 (MQ=255)
                                                                                                                             ggaaccctGGCGCCCAGCGCGCCAAGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGATCAACAGCCGGTACATGATGTCGGGGTTTGCTGCGAGCGACATGAACTTTACCTGAAc  >  8:28230/8‑149 (MQ=255)
                                                                                                                                cccTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGATCAACAGCCGGTACATGATGTc                                      >  6:181402/1‑110 (MQ=255)
                                                                                                                                cccTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGATCAACAGCCGGTACATGATGTc                                      <  5:181402/110‑1 (MQ=255)
                                                                                                                                  cTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGATCAACAGCCGGTACATGATGTCGGGGTTTGCTgcgagcga                    <  2:239399/126‑1 (MQ=255)
                                                                                                                                  cTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGATCAACAGCCGGTACATGATGTCGGGGTTTGCTgcgagcga                    >  1:239399/1‑126 (MQ=255)
                                                                                                                                              |                                                                                                                                   
CGCCCGTTTCTTCAAAGCGTCCGGTGCTGCGCGCCTGCTCCAGGTTGAATGCCGGCAACCCGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGATCAACAGCCGGTACATGATGTCGGGGTTTGCTGCGAGCGACATGAACTTTACCTGAAC  >  NC_002947/4234806‑4235079

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