Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,264,193 C→G 5.6% 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 5.6% 1845.7 / 64.9 515G47G (GGC→GGGssbsingle‑stranded DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (217/269);  new base G (20/9);  total (237/278)
Fisher's exact test for biased strand distribution p-value = 1.24e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

Only 100 of 550 total aligned reads displayed.

AAATGGTGGCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGC  >  CP009273/4264127‑4264259
                                                                  |                                                                  
aaaTGGTGGCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCg                                                                   <  1:5825287/68‑1 (MQ=255)
  aTGGTGGCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCgag                                                                 >  1:23798940/1‑68 (MQ=255)
  aTGGTGGCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCgag                                                                 <  1:34673000/68‑1 (MQ=255)
  aTGGTGGCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCgag                                                                 <  1:3054208/68‑1 (MQ=255)
   tggtggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCgaga                                                                >  1:22062223/1‑68 (MQ=255)
   tggtggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCgaga                                                                <  1:7186103/68‑1 (MQ=255)
   tggtggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCgaga                                                                >  1:25388138/1‑68 (MQ=255)
    ggtggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGAt                                                               >  1:435419/1‑68 (MQ=255)
    ggtggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGAt                                                               >  1:3622549/1‑68 (MQ=255)
    ggtggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGAt                                                               >  1:6584324/1‑68 (MQ=255)
      tggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATg                                                              <  1:20183854/67‑1 (MQ=255)
      tggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATg                                                              <  1:13180117/67‑1 (MQ=255)
       ggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGaa                                                            <  1:18924593/68‑1 (MQ=255)
         cgcAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAg                                                          >  1:5011976/1‑68 (MQ=255)
          gcAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGa                                                         <  1:24134109/68‑1 (MQ=255)
          gcAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGa                                                         <  1:1621238/68‑1 (MQ=255)
          gcAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGa                                                         >  1:25664673/1‑68 (MQ=255)
               tGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAg                                                     <  1:30946748/67‑1 (MQ=255)
               tGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAg                                                     >  1:7790893/1‑67 (MQ=255)
               tGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGa                                                    <  1:24576218/68‑1 (MQ=255)
                   aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                <  1:18153485/68‑1 (MQ=255)
                   aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                <  1:17774279/68‑1 (MQ=255)
                   aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                <  1:5500547/68‑1 (MQ=255)
                   aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                <  1:3215744/68‑1 (MQ=255)
                    aCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGaa                                               >  1:3939688/1‑68 (MQ=255)
                    aCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGaa                                               >  1:18813774/1‑68 (MQ=255)
                    aCATTACGCTGGCTACTTCCGAATCCTGGCCTTATAAAGCGAACGGCGAGATGAAAGAAAAGACTGaa                                               >  1:31665371/1‑68 (MQ=255)
                      aTTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATg                                             >  1:16945182/1‑68 (MQ=255)
                      aTTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATg                                             >  1:15811115/1‑68 (MQ=255)
                      aTTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATg                                             >  1:11434968/1‑68 (MQ=255)
                        tACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGc                                           <  1:30302924/68‑1 (MQ=255)
                         aCGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGc                                           >  1:27116974/1‑67 (MQ=255)
                           gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                        >  1:32770600/1‑68 (MQ=255)
                           gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                        >  1:4753360/1‑68 (MQ=255)
                           gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                        >  1:23523982/1‑68 (MQ=255)
                           gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                        >  1:22857492/1‑68 (MQ=255)
                           gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                        >  1:33352520/1‑68 (MQ=255)
                           gCTGGCTACTTCCGAATACTGGCGTGAGAAAGCCACCGGCGAGATCAAAGACCAGACAGAATGGCAcc                                        >  1:28162048/1‑68 (MQ=255)
                            cTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACTGGCGAGATGAAAGAACAGACTGAATGGCAcc                                        >  1:29446417/1‑67 (MQ=255)
                            cTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                        >  1:2491281/1‑67 (MQ=255)
                            cTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                        >  1:15382960/1‑67 (MQ=255)
                            cTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACGGAATGGcccc                                        >  1:11034454/1‑67 (MQ=255)
                             tGGCTACTTCCGAATCCTGGCGTTATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                      <  1:737522/68‑1 (MQ=255)
                             tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                      <  1:28777399/68‑1 (MQ=255)
                             tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                      <  1:26761238/68‑1 (MQ=255)
                             tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                      <  1:31216866/68‑1 (MQ=255)
                             tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                      <  1:31418260/68‑1 (MQ=255)
                             tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                      <  1:19821630/68‑1 (MQ=255)
                             tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                      <  1:17546501/68‑1 (MQ=255)
                             tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                      <  1:5990440/68‑1 (MQ=255)
                             tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                      <  1:10665092/68‑1 (MQ=255)
                             tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                      <  1:1355701/68‑1 (MQ=255)
                             tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                      <  1:12934564/68‑1 (MQ=255)
                             tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                      <  1:12499743/68‑1 (MQ=255)
                                cTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGtt                                   >  1:2093269/1‑68 (MQ=255)
                                cTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGtt                                   >  1:27897230/1‑68 (MQ=255)
                                cTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGt                                    <  1:1830407/67‑1 (MQ=255)
                                 tACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGtt                                   <  1:9216060/67‑1 (MQ=255)
                                  aCTTCCGAATCCTGGCGTGATAAAGCGACCGGCTAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                 <  1:14716436/68‑1 (MQ=255)
                                  aCTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                 <  1:21702984/68‑1 (MQ=255)
                                  aCTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                 <  1:26139559/68‑1 (MQ=255)
                                  aCTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                 <  1:14961657/68‑1 (MQ=255)
                                     tCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCCCCGCGTtgt                                 >  1:9621274/1‑65 (MQ=255)
                                     tCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                 >  1:9074065/1‑65 (MQ=255)
                                       cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                            <  1:10002105/68‑1 (MQ=255)
                                       cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                            <  1:15157255/68‑1 (MQ=255)
                                       cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                            >  1:16070383/1‑68 (MQ=255)
                                       cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                            >  1:17130355/1‑68 (MQ=255)
                                       cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                            <  1:10516775/68‑1 (MQ=255)
                                       cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                            <  1:34455310/68‑1 (MQ=255)
                                       cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                            <  1:19643619/68‑1 (MQ=255)
                                       cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                            <  1:21272150/68‑1 (MQ=255)
                                       cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                            <  1:28566842/68‑1 (MQ=255)
                                       cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                            <  1:24026286/68‑1 (MQ=255)
                                       cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                            <  1:29747044/68‑1 (MQ=255)
                                       cGAATCCTGGCGTGATAAACCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                            <  1:14262217/68‑1 (MQ=255)
                                           tCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                <  1:892233/44‑1 (MQ=255)
                                            gcTGGCGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGc                       >  1:27495003/2‑68 (MQ=255)
                                            gcTGGCGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGc                       >  1:16441812/2‑68 (MQ=255)
                                            gcTGGCGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGc                       >  1:20365222/2‑68 (MQ=255)
                                            gcTGGCGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGc                       >  1:8291167/2‑68 (MQ=255)
                                                 cGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCaa                     >  1:25960545/1‑65 (MQ=255)
                                                 cGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCaa                     >  1:33206964/1‑65 (MQ=255)
                                                 cGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAAc                   <  1:14019295/67‑1 (MQ=255)
                                                 cGTGACAAACAGACGGGGGAGATGCGGGAGCAGACATAATGGCATCGCGTGGTGCTGTTCGGCaagc                   <  1:16703126/67‑3 (MQ=255)
                                                  gagATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTg                 <  1:33565786/66‑1 (MQ=255)
                                                    gATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTgg                <  1:29188117/67‑1 (MQ=255)
                                                    gATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTgg                <  1:10245283/67‑1 (MQ=255)
                                                       aaaGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGa            >  1:11710066/1‑68 (MQ=255)
                                                         aGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAg          >  1:33875702/1‑68 (MQ=255)
                                                          gCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGt         >  1:12046515/1‑68 (MQ=255)
                                                           agACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTg        >  1:18637688/2‑68 (MQ=255)
                                                           agACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTg        >  1:6376794/2‑68 (MQ=255)
                                                               cGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcga    <  1:20672433/68‑1 (MQ=255)
                                                                ggCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgag   <  1:21422650/68‑1 (MQ=255)
                                                                  cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc  <  1:15622644/67‑1 (MQ=255)
                                                                  cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc  <  1:18499032/67‑1 (MQ=255)
                                                                  cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc  <  1:32583949/67‑1 (MQ=255)
                                                                  cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc  <  1:22513487/67‑1 (MQ=255)
                                                                  cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc  <  1:29388604/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: