Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,264,274 A→G 10.3% K74K (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,2740AG10.3% 2265.5 / 179.4 712K74K (AAA→AAGssbsingle‑stranded DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (205/431);  new base G (24/49);  total (229/483)
Fisher's exact test for biased strand distribution p-value = 8.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.95e-01

Only 100 of 828 total aligned reads displayed.

AGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATC  >  CP009273/4264207‑4264341
                                                                   |                                                                   
   cTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGt                                                                  <  1:24478295/68‑1 (MQ=255)
   cTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGt                                                                  >  1:19718289/1‑68 (MQ=255)
   cTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGt                                                                  >  1:19296234/1‑68 (MQ=255)
     gAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTtc                                                                <  1:23849260/68‑1 (MQ=255)
     gAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTtc                                                                <  1:8585503/68‑1 (MQ=255)
       aTGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAggg                                                                   <  1:25895092/63‑1 (MQ=255)
        tGGCACCGTGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCa                                                             >  1:11450043/1‑68 (MQ=255)
        tGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCa                                                             >  1:28570333/1‑68 (MQ=255)
           cACCGCGTTGTTCTGTTCGGCAAACTTTCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                                           <  1:14101625/67‑1 (MQ=255)
           cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                                           <  1:17021203/67‑1 (MQ=255)
           cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                                           <  1:7401722/67‑1 (MQ=255)
           cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                                           <  1:8455844/67‑1 (MQ=255)
           cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                                           <  1:542871/67‑1 (MQ=255)
           cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                                           <  1:27851307/67‑1 (MQ=255)
           cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                                           <  1:34204148/67‑1 (MQ=255)
            atcgcgTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGt                                                          <  1:24276615/65‑1 (MQ=255)
            aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                                                         >  1:3338580/1‑68 (MQ=255)
            aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                                                         >  1:31260688/1‑68 (MQ=255)
            aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                                                         >  1:16163512/1‑68 (MQ=255)
            aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                                                         >  1:10881199/1‑68 (MQ=255)
            aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGt                                                          <  1:11874929/67‑1 (MQ=255)
             tcgcgTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGtct                                                        >  1:2345396/2‑66 (MQ=255)
               gcgTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTtat                                                      <  1:14921673/68‑1 (MQ=255)
                 gTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTtata                                                     <  1:21601532/67‑1 (MQ=255)
                      gCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAg                                               >  1:601136/1‑68 (MQ=255)
                      gCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAg                                               >  1:3788071/1‑68 (MQ=255)
                      gCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAg                                               >  1:15215350/1‑68 (MQ=255)
                      gCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAg                                               >  1:1279531/1‑68 (MQ=255)
                      gCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGaa                                                <  1:16676754/67‑1 (MQ=255)
                      gCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGaa                                                <  1:18654898/67‑1 (MQ=255)
                      gCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGaa                                                <  1:9156473/67‑1 (MQ=255)
                      gCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAg                                               >  1:13571387/1‑68 (MQ=255)
                      gCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAg                                               >  1:24644386/1‑68 (MQ=255)
                      gCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAg                                               >  1:12573419/1‑68 (MQ=255)
                      gCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAg                                               >  1:10534157/1‑68 (MQ=255)
                        tGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAgg                                              <  1:8876724/67‑1 (MQ=255)
                           tCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGG‑‑CGCAGGTCTACATCGAAGGTCAGCt                                        >  1:20331265/1‑68 (MQ=255)
                           tCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCa                                           >  1:10156750/1‑67 (MQ=255)
                           tCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCa                                           >  1:13101747/1‑67 (MQ=255)
                            cGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGc                                         <  1:30414308/68‑1 (MQ=255)
                            cGGAAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGc                                         <  1:34501026/68‑1 (MQ=255)
                             ggCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCt                                        <  1:17303736/68‑1 (MQ=255)
                             ggCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCt                                        <  1:8746213/68‑1 (MQ=255)
                             ggCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCt                                        <  1:12342087/68‑1 (MQ=255)
                             ggCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCt                                        <  1:29928393/68‑1 (MQ=255)
                                aaaCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGc                                      <  1:32458263/67‑1 (MQ=255)
                                aaaCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGc                                      <  1:8024950/67‑1 (MQ=255)
                                aaaCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCg                                     >  1:29124981/1‑68 (MQ=255)
                                    tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc                                 <  1:3187780/68‑1 (MQ=255)
                                    tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc                                 <  1:11618794/68‑1 (MQ=255)
                                    tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc                                 <  1:30983352/68‑1 (MQ=255)
                                    tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc                                 <  1:31600735/68‑1 (MQ=255)
                                    tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc                                 <  1:21336786/68‑1 (MQ=255)
                                    tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc                                 <  1:25074787/68‑1 (MQ=255)
                                    tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc                                 <  1:13188900/68‑1 (MQ=255)
                                    tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc                                 <  1:649307/68‑1 (MQ=255)
                                    tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc                                 <  1:5801645/68‑1 (MQ=255)
                                    tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc                                 <  1:20114220/68‑1 (MQ=255)
                                    tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc                                 <  1:4189101/68‑1 (MQ=255)
                                     ggCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAccc                                <  1:25397556/68‑1 (MQ=255)
                                     ggCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAccc                                <  1:29588697/68‑1 (MQ=255)
                                      gCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGt                                                          <  1:29450887/41‑1 (MQ=255)
                                        ggAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCTTCGCACCCGTa                             <  1:9663687/67‑1 (MQ=255)
                                           aGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTa                             <  1:4676627/65‑1 (MQ=255)
                                           aGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTa                             <  1:25676468/65‑1 (MQ=255)
                                           aGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTa                             <  1:16367210/65‑1 (MQ=255)
                                           aGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTa                             <  1:7757200/65‑1 (MQ=255)
                                           aGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTa                             <  1:14545232/65‑1 (MQ=255)
                                           aGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCTTCGCACCCGTagc                           <  1:27052569/67‑3 (MQ=255)
                                           aGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCTTCGCACCCGTagc                           <  1:30606212/67‑3 (MQ=255)
                                            gTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATg                         >  1:15733794/1‑68 (MQ=255)
                                                cgagcgaATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGAcc                     >  1:26730169/1‑68 (MQ=255)
                                                   gcgaATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGa                   >  1:14340862/1‑67 (MQ=255)
                                                   gcgaATATCTGCGTAAAGGCTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGa                   <  1:10006388/67‑1 (MQ=255)
                                                       atatCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTa                             <  1:21114889/53‑1 (MQ=255)
                                                       atatCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCaa               <  1:23109324/67‑1 (MQ=255)
                                                         atCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATcc            >  1:18923561/1‑68 (MQ=255)
                                                           cTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCgg          <  1:17869295/68‑1 (MQ=255)
                                                           cTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCgg          <  1:15479995/68‑1 (MQ=255)
                                                           cTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCgg          <  1:32970523/68‑1 (MQ=255)
                                                           cTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCgg          <  1:32694174/68‑1 (MQ=255)
                                                           cTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCgg          <  1:25250685/68‑1 (MQ=255)
                                                           cTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCgg          <  1:27293433/68‑1 (MQ=255)
                                                           cTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGCACCGATCAATCCgg          <  1:6979548/68‑1 (MQ=255)
                                                               gTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAg      >  1:9422297/1‑68 (MQ=255)
                                                                 aaaGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAgg     <  1:4380639/67‑1 (MQ=255)
                                                                 aaaGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAgg     <  1:22269405/67‑1 (MQ=255)
                                                                 aaaGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAgg     <  1:27529090/67‑1 (MQ=255)
                                                                 aaaGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGa    <  1:13808670/68‑1 (MQ=255)
                                                                  aaGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGa    >  1:11134593/1‑67 (MQ=255)
                                                                   aGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATc  <  1:1826299/68‑1 (MQ=255)
                                                                   aGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATc  <  1:6278250/68‑1 (MQ=255)
                                                                   aGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATc  <  1:17520262/68‑1 (MQ=255)
                                                                   aGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATc  <  1:5146403/68‑1 (MQ=255)
                                                                   aGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATc  <  1:28128195/68‑1 (MQ=255)
                                                                   aGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATc  <  1:2078786/68‑1 (MQ=255)
                                                                   aGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATc  <  1:21967323/68‑1 (MQ=255)
                                                                   aGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATc  <  1:33621255/68‑1 (MQ=255)
                                                                   aGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATc  <  1:22674084/68‑1 (MQ=255)
                                                                   aGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATc  <  1:26362865/68‑1 (MQ=255)
                                                                   |                                                                   
AGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATC  >  CP009273/4264207‑4264341

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: