Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,264,241 A→G 12.0% K63K (AAA→AAG ssb → single‑stranded DNA‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,264,2410AG12.0% 1545.3 / 147.9 511K63K (AAA→AAGssbsingle‑stranded DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (209/239);  new base G (25/36);  total (234/277)
Fisher's exact test for biased strand distribution p-value = 4.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01

Only 100 of 542 total aligned reads displayed.

GTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGT  >  CP009273/4264177‑4264307
                                                                |                                                                  
  gATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTgg                                                                <  1:29697912/67‑1 (MQ=255)
  gATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTgg                                                                <  1:12077509/67‑1 (MQ=255)
     aaaGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGa                                                            >  1:4217054/1‑68 (MQ=255)
     aaaGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGa                                                            >  1:13628153/1‑68 (MQ=255)
     aaaGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGa                                                            <  1:24720537/68‑1 (MQ=255)
      aaGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGaa                                                           >  1:23907868/1‑68 (MQ=255)
       aGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAg                                                          >  1:5128122/1‑68 (MQ=255)
       aGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAg                                                          >  1:21081585/1‑68 (MQ=255)
       aGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAg                                                          >  1:1420614/1‑68 (MQ=255)
       aGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAg                                                          >  1:23687491/1‑68 (MQ=255)
        gCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGt                                                         >  1:19717357/1‑68 (MQ=255)
        gCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGt                                                         >  1:27044479/1‑68 (MQ=255)
        gCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGt                                                         >  1:339159/1‑68 (MQ=255)
         cGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTg                                                        >  1:12295616/1‑68 (MQ=255)
                cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc                                                  <  1:28061592/67‑1 (MQ=255)
                cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc                                                  <  1:11908294/67‑1 (MQ=255)
                cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc                                                  <  1:31205351/67‑1 (MQ=255)
                cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc                                                  <  1:20277087/67‑1 (MQ=255)
                cGAGATGAAAGAACACACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc                                                  <  1:1035288/67‑1 (MQ=255)
                 gagaTGATAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagcga                                                <  1:19491950/68‑1 (MQ=255)
                 gagaTGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagcga                                                <  1:4763370/68‑1 (MQ=255)
                         aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc                                        <  1:9767761/68‑1 (MQ=255)
                         aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc                                        <  1:12823893/68‑1 (MQ=255)
                         aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc                                        <  1:8235481/68‑1 (MQ=255)
                         aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc                                        <  1:2881346/68‑1 (MQ=255)
                         aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc                                        <  1:18811762/68‑1 (MQ=255)
                          gagcagacagaATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGc                                        <  1:15731443/64‑1 (MQ=255)
                          gagcagacagaATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGc                                        <  1:27787880/64‑1 (MQ=255)
                          gagcagacagaATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGc                                        <  1:30645784/64‑1 (MQ=255)
                          gAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCg                                       >  1:11308206/1‑68 (MQ=255)
                          gAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCg                                       >  1:29087770/1‑68 (MQ=255)
                          gAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCg                                       >  1:16616924/1‑68 (MQ=255)
                           aaCAGACTGAATGGCCCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGt                                      >  1:15468293/1‑68 (MQ=255)
                           aaCAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGt                                      >  1:14560464/1‑68 (MQ=255)
                            aCAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTTGCAGAAGTGGCGAGAGAATATCTGCGTa                                     >  1:26516218/1‑68 (MQ=255)
                              aGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTaaa                                   >  1:14335600/1‑68 (MQ=255)
                                 cagAATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCaag                                   <  1:22168340/63‑2 (MQ=255)
                                 cagAATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCaag                                   <  1:13047647/63‑2 (MQ=255)
                                 cTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGt                                <  1:15383895/68‑1 (MQ=255)
                                 cTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGt                                >  1:30297757/1‑68 (MQ=255)
                                 cTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGt                                <  1:6505199/68‑1 (MQ=255)
                                 cTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGt                                <  1:7623715/68‑1 (MQ=255)
                                 cTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGt                                <  1:26179686/68‑1 (MQ=255)
                                 cTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGt                                >  1:25204030/1‑68 (MQ=255)
                                     aTGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTtctc                            >  1:26659666/1‑68 (MQ=255)
                                     aTGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTtctc                            >  1:10491883/1‑68 (MQ=255)
                                      tGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCa                           >  1:29428446/1‑68 (MQ=255)
                                      tGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCa                           <  1:12572386/68‑1 (MQ=255)
                                      tGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCa                           >  1:57741/1‑68 (MQ=255)
                                      tGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCa                           >  1:25036822/1‑68 (MQ=255)
                                      tGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGGAAAGGTTCTCa                           >  1:28466679/1‑68 (MQ=255)
                                        gCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                         >  1:23463525/1‑68 (MQ=255)
                                         cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                         <  1:20661061/67‑1 (MQ=255)
                                         cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                         <  1:5339199/67‑1 (MQ=255)
                                         cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                         <  1:7892882/67‑1 (MQ=255)
                                         cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                         <  1:24392046/67‑1 (MQ=255)
                                         cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                         <  1:6814192/67‑1 (MQ=255)
                                         cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                         <  1:8572502/67‑1 (MQ=255)
                                         cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                         <  1:27405863/67‑1 (MQ=255)
                                         cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                         <  1:30006921/67‑1 (MQ=255)
                                         cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                         <  1:6298669/67‑1 (MQ=255)
                                          atcgcgTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGtc                       <  1:739098/66‑2 (MQ=255)
                                          aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                       >  1:21237361/1‑68 (MQ=255)
                                          aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                       >  1:5613531/1‑68 (MQ=255)
                                          aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                       >  1:4964936/1‑68 (MQ=255)
                                          aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                       >  1:10687452/1‑68 (MQ=255)
                                          aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                       >  1:28285190/1‑68 (MQ=255)
                                          aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                       >  1:21695721/1‑68 (MQ=255)
                                          aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                       >  1:2268506/1‑68 (MQ=255)
                                          aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                       >  1:15135169/1‑68 (MQ=255)
                                          aCCGCGATCTGCTGTTCGGCAAACTCGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                       >  1:8395150/1‑68 (MQ=255)
                                               gttgtgttgtTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTtata                   >  1:594416/1‑67 (MQ=255)
                                               gTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTtata                   >  1:14826724/1‑67 (MQ=255)
                                               gTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTtata                   >  1:17448595/1‑67 (MQ=255)
                                               gTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTatat                  <  1:2237273/68‑1 (MQ=255)
                                                    gCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAg             >  1:30971078/1‑68 (MQ=255)
                                                    gCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAg             >  1:31661833/1‑68 (MQ=255)
                                                    gCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAg             >  1:11047399/1‑68 (MQ=255)
                                                    gCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAg             >  1:25611183/1‑68 (MQ=255)
                                                    gCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAg             >  1:18575059/1‑68 (MQ=255)
                                                    gCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAg             >  1:12210184/1‑68 (MQ=255)
                                                    gCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAg             >  1:21500634/1‑68 (MQ=255)
                                                      tGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAgg            >  1:23136197/1‑67 (MQ=255)
                                                         tCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCa         >  1:9417694/1‑67 (MQ=255)
                                                         tCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCa         >  1:16306500/1‑67 (MQ=255)
                                                         tCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCa         >  1:29846497/1‑67 (MQ=255)
                                                         tCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCa         >  1:25984386/1‑67 (MQ=255)
                                                         tCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCa         >  1:16819308/1‑67 (MQ=255)
                                                          cGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGc       <  1:21405567/68‑1 (MQ=255)
                                                           ggCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCt      <  1:9128921/68‑1 (MQ=255)
                                                           ggCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCt      <  1:22826456/68‑1 (MQ=255)
                                                           ggCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCt      <  1:32013124/68‑1 (MQ=255)
                                                           ggCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCt      <  1:19195472/68‑1 (MQ=255)
                                                              aaaCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATTTCGAAGGTCAGCTGCg   >  1:17109572/1‑68 (MQ=255)
                                                              aaaCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCg   >  1:7208726/1‑68 (MQ=255)
                                                              aaaCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCg   >  1:10054196/1‑68 (MQ=255)
                                                              aaaCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCg   >  1:13360538/1‑68 (MQ=255)
                                                              aaaCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCg   >  1:8768603/1‑68 (MQ=255)
                                                              aaaCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCg   >  1:4042809/1‑68 (MQ=255)
                                                               aaCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGt  <  1:1807712/68‑1 (MQ=255)
                                                                |                                                                  
GTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGT  >  CP009273/4264177‑4264307

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