Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,264,193 C→G 6.5% G47G (GGC→GGG ssb → single‑stranded DNA‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,264,1930CG 6.5% 2064.3 / 84.6 588G47G (GGC→GGGssbsingle‑stranded DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (240/308);  new base G (27/12);  total (267/321)
Fisher's exact test for biased strand distribution p-value = 2.48e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 630 total aligned reads displayed.

AAATGGTGGCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGC  >  CP009273/4264127‑4264259
                                                                  |                                                                  
aaaTGGTGGCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCg                                                                   >  1:1934547/1‑68 (MQ=255)
  aTGGTGGCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCgag                                                                 <  1:7008798/68‑1 (MQ=255)
   tggtggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCgaga                                                                >  1:13576838/1‑68 (MQ=255)
   tggtggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCgaga                                                                >  1:26080257/1‑68 (MQ=255)
    ggtggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGAt                                                               >  1:15727605/1‑68 (MQ=255)
    ggtggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGAt                                                               >  1:22711231/1‑68 (MQ=255)
     gtggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATg                                                              >  1:25366693/1‑68 (MQ=255)
     gtggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATg                                                              >  1:31770094/1‑68 (MQ=255)
      tggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATg                                                              <  1:18299174/67‑1 (MQ=255)
      tggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATg                                                              <  1:37084681/67‑1 (MQ=255)
      tggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATg                                                              <  1:10785982/67‑1 (MQ=255)
      tggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATg                                                              <  1:31241881/67‑1 (MQ=255)
      tggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATg                                                              <  1:21354442/67‑1 (MQ=255)
       ggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGaa                                                            <  1:25595027/68‑1 (MQ=255)
       ggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGaa                                                            >  1:9265993/1‑68 (MQ=255)
          gcAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAg                                                          <  1:16781108/67‑1 (MQ=255)
          gcAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGa                                                         <  1:3653513/68‑1 (MQ=255)
               tGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAActg                                                     <  1:36340122/67‑3 (MQ=255)
               tGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAg                                                     <  1:32474696/67‑1 (MQ=255)
               tGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGa                                                    <  1:29588343/68‑1 (MQ=255)
               tGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGa                                                    <  1:33251841/68‑1 (MQ=255)
                  cAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACt                                                  <  1:14734194/67‑1 (MQ=255)
                   aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                <  1:35603362/68‑1 (MQ=255)
                   aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                <  1:23463372/68‑1 (MQ=255)
                   aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                <  1:25214471/68‑1 (MQ=255)
                   aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                <  1:26341505/68‑1 (MQ=255)
                   aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                <  1:30519216/68‑1 (MQ=255)
                   aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                <  1:11508117/68‑1 (MQ=255)
                   aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                <  1:15278901/68‑1 (MQ=255)
                   aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                <  1:9426570/68‑1 (MQ=255)
                   aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                <  1:16507469/68‑1 (MQ=255)
                    aCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGaa                                               >  1:5750486/1‑68 (MQ=255)
                    aCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGaa                                               >  1:12213227/1‑68 (MQ=255)
                    aCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGaa                                               >  1:22383902/1‑68 (MQ=255)
                      aTTACGCTGGCTACTTCCGAATCCTGGCGTGGTAAAGCGACCGGCGAGATGAAAGAACAGACTGAATg                                             >  1:34281471/1‑68 (MQ=255)
                      aTTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATg                                             >  1:5340296/1‑68 (MQ=255)
                        tACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCTACCGGCGAGATGAAAGAACAGACTGAATGGc                                           <  1:28034248/68‑1 (MQ=255)
                        tACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGc                                           <  1:24723293/68‑1 (MQ=255)
                        tACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGc                                           <  1:14003121/68‑1 (MQ=255)
                        tACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGc                                           <  1:7675679/68‑1 (MQ=255)
                        tACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGc                                           <  1:6681060/68‑1 (MQ=255)
                          cGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAc                                         <  1:1433625/68‑1 (MQ=255)
                           gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGCCCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                        >  1:27512973/1‑68 (MQ=255)
                           gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                        >  1:29710023/1‑68 (MQ=255)
                           gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                        >  1:28554266/1‑68 (MQ=255)
                           gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                        >  1:1046036/1‑68 (MQ=255)
                           gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                        >  1:8234755/1‑68 (MQ=255)
                           gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                        >  1:8984102/1‑68 (MQ=255)
                            cTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                        >  1:32998456/1‑67 (MQ=255)
                            cTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                        >  1:28388481/1‑67 (MQ=255)
                             tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                      <  1:28813218/68‑1 (MQ=255)
                             tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                      <  1:1865464/68‑1 (MQ=255)
                             tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                      <  1:26725921/68‑1 (MQ=255)
                             tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                      <  1:32899657/68‑1 (MQ=255)
                             tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                      <  1:34001122/68‑1 (MQ=255)
                             tGGCTACTTCCGAATCCTGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                      <  1:17126992/67‑1 (MQ=255)
                                cTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGtt                                   >  1:7278486/1‑68 (MQ=255)
                                cTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGt                                    <  1:13290763/67‑1 (MQ=255)
                                  aCTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                 <  1:8094970/68‑1 (MQ=255)
                                  aCTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                 <  1:6503395/68‑1 (MQ=255)
                                  aCTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                 <  1:26915659/68‑1 (MQ=255)
                                  aCTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                 <  1:10230131/68‑1 (MQ=255)
                                     tCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                 >  1:29989562/1‑65 (MQ=255)
                                     tCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                 >  1:15508981/1‑65 (MQ=255)
                                     tCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                 >  1:21032103/1‑65 (MQ=255)
                                     tCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                 >  1:32059117/1‑65 (MQ=255)
                                     tCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCt                              <  1:11401322/68‑1 (MQ=255)
                                       cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                            <  1:31467502/68‑1 (MQ=255)
                                       cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                            <  1:18074909/68‑1 (MQ=255)
                                       cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                            <  1:5020670/68‑1 (MQ=255)
                                          aTCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCg                         >  1:23967926/1‑68 (MQ=255)
                                          aTCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCg                         >  1:18984329/1‑68 (MQ=255)
                                          aTCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCg                         >  1:19734732/1‑68 (MQ=255)
                                          aTCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCg                         >  1:20435600/1‑68 (MQ=255)
                                            gcTGGCGTGACAAACAGACGGGGGAGATGCGGCAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGc                       >  1:34615631/2‑68 (MQ=255)
                                               ggCGTGATAAATCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCaa                     >  1:22907186/1‑67 (MQ=255)
                                                 cGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCaa                     >  1:35002699/1‑65 (MQ=255)
                                                 cGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCaagc                   <  1:22100162/67‑3 (MQ=255)
                                                  gTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTg                 >  1:24410449/1‑68 (MQ=255)
                                                  gTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTg                 <  1:21742812/68‑1 (MQ=255)
                                                  gTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCAAGCtc                 <  1:30312965/68‑2 (MQ=255)
                                                  gTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCAAGCt                  <  1:4384841/67‑1 (MQ=255)
                                                    gATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTgg                <  1:19978430/67‑1 (MQ=255)
                                                      tAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCa              <  1:16070140/67‑1 (MQ=255)
                                                         aGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAg          >  1:32704909/1‑68 (MQ=255)
                                                         aGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAg          >  1:4770634/1‑68 (MQ=255)
                                                          gCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGt         >  1:13052103/1‑68 (MQ=255)
                                                          gCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGt         >  1:1113593/1‑68 (MQ=255)
                                                           cGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTg        >  1:3998161/1‑68 (MQ=255)
                                                           agACGGGGGAGATGCGGGAGCAGAAAGAATGGCATCGCGTGGTGCTGTTCGGCAAGCACGCCGAAGTg        >  1:6080121/2‑68 (MQ=255)
                                                            gACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTg        >  1:17311417/1‑67 (MQ=255)
                                                            gACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTgg       <  1:28968911/68‑1 (MQ=255)
                                                                  cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc  <  1:35949505/67‑1 (MQ=255)
                                                                  cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc  <  1:33517494/67‑1 (MQ=255)
                                                                  cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc  <  1:20734894/67‑1 (MQ=255)
                                                                  cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc  <  1:17701029/67‑1 (MQ=255)
                                                                  cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc  <  1:3099029/67‑1 (MQ=255)
                                                                  cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc  <  1:12523483/67‑1 (MQ=255)
                                                                  cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc  <  1:10000779/67‑1 (MQ=255)
                                                                  cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc  <  1:27827135/67‑1 (MQ=255)
                                                                  |                                                                  
AAATGGTGGCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGC  >  CP009273/4264127‑4264259

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