Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,264,201 A→G 6.7% K50R (AAA→AGA) ‡ ssb → single‑stranded DNA‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,264,2010AG 6.7% 1902.8 / 74.9 547K50R (AAA→AGA) ‡ssbsingle‑stranded DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (185/325);  new base G (25/12);  total (210/337)
Fisher's exact test for biased strand distribution p-value = 3.23e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

Only 100 of 562 total aligned reads displayed.

GGCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTG  >  CP009273/4264134‑4264268
                                                                   |                                                                   
   gcAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAg                                                                   <  1:8382636/67‑1 (MQ=255)
   gcAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAg                                                                   <  1:20433307/67‑1 (MQ=255)
   gcAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAg                                                                   <  1:24639751/67‑1 (MQ=255)
   gcAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGa                                                                  >  1:5043938/1‑68 (MQ=255)
   gcAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGa                                                                  <  1:3653513/68‑1 (MQ=255)
        tGCCAACATTACGCTGGCTACTTCCGAATCCTGGGGTGATAAAGGTACCGGCGAGATGAAAGAACAGa                                                             <  1:28675600/68‑1 (MQ=255)
        tGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAActg                                                              <  1:36340122/67‑3 (MQ=255)
        tGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAg                                                              <  1:27882697/67‑1 (MQ=255)
        tGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGa                                                             <  1:33251841/68‑1 (MQ=255)
        tGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGa                                                             <  1:34185940/68‑1 (MQ=255)
           cAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACt                                                           <  1:27651239/67‑1 (MQ=255)
           cAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACt                                                           <  1:33516063/67‑1 (MQ=255)
           cAACATTACGCTGGCTACTTCCGAATCCTGGCGGGATAAAGCGACCGGCGAGATGAAAGAACAGACt                                                           <  1:21088820/67‑1 (MQ=255)
            aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                         <  1:30286986/68‑1 (MQ=255)
            aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                         <  1:20825579/68‑1 (MQ=255)
            aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                         <  1:11508117/68‑1 (MQ=255)
            aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                         <  1:28442823/68‑1 (MQ=255)
            aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                         <  1:2889593/68‑1 (MQ=255)
             aCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGaa                                                        >  1:25750767/1‑68 (MQ=255)
             aCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGaa                                                        >  1:5750486/1‑68 (MQ=255)
               aTTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATg                                                      <  1:31269432/68‑1 (MQ=255)
               aTTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATg                                                      >  1:5340296/1‑68 (MQ=255)
               aTTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATg                                                      >  1:6145859/1‑68 (MQ=255)
                ttACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATg                                                      >  1:14766916/1‑67 (MQ=255)
                 tACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGc                                                    <  1:24116944/68‑1 (MQ=255)
                 tACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGc                                                    <  1:16662874/68‑1 (MQ=255)
                    gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                                 >  1:7298494/1‑68 (MQ=255)
                    gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                                 >  1:30415521/1‑68 (MQ=255)
                    gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                                 >  1:26665441/1‑68 (MQ=255)
                    gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                                 >  1:26320800/1‑68 (MQ=255)
                    gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                                 >  1:24954678/1‑68 (MQ=255)
                    gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAc                                                  >  1:17713685/1‑67 (MQ=255)
                      tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                               <  1:22523112/68‑1 (MQ=255)
                      tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                               <  1:4003283/68‑1 (MQ=255)
                      tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                               <  1:12585987/68‑1 (MQ=255)
                      tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                               <  1:18470441/68‑1 (MQ=255)
                      tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                               <  1:17398987/68‑1 (MQ=255)
                      tGGCTACTTCCGAATCCTGGCGGGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                               <  1:13801628/68‑1 (MQ=255)
                         cTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGtt                                            >  1:31793944/1‑68 (MQ=255)
                         cTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGt                                             <  1:10397247/67‑1 (MQ=255)
                         cTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGt                                             <  1:13290763/67‑1 (MQ=255)
                          tACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGtt                                            <  1:22944316/67‑1 (MQ=255)
                          tACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGtt                                            <  1:10729172/67‑1 (MQ=255)
                          tACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGtt                                            <  1:18909475/67‑1 (MQ=255)
                           aCTTCCTAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                          <  1:1015993/68‑1 (MQ=255)
                           aCTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                          <  1:6503395/68‑1 (MQ=255)
                           aCTTCCGAATCCTGGCGGGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgg                                          <  1:25460899/68‑2 (MQ=255)
                              tCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                          >  1:13088736/1‑65 (MQ=255)
                              tCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                          >  1:15879110/1‑65 (MQ=255)
                              tCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                          >  1:26843569/1‑65 (MQ=255)
                              tCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                          >  1:29989562/1‑65 (MQ=255)
                              tCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCt                                       <  1:11401322/68‑1 (MQ=255)
                              tCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCt                                       <  1:8174399/68‑1 (MQ=255)
                                cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                     <  1:21902254/68‑1 (MQ=255)
                                cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                     <  1:29617819/68‑1 (MQ=255)
                                cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                     <  1:31467502/68‑1 (MQ=255)
                                cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                     <  1:32220864/68‑1 (MQ=255)
                                cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                     <  1:18074909/68‑1 (MQ=255)
                                   aTCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCg                                  >  1:22798564/1‑68 (MQ=255)
                                   aTCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCg                                  >  1:16169810/1‑68 (MQ=255)
                                     gcTGGCGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGc                                >  1:28236892/2‑68 (MQ=255)
                                     gcTGGCGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGc                                >  1:12261465/2‑68 (MQ=255)
                                     gcTGGCGTGACAAACAGACGGGGGAGATGCGGCAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGc                                >  1:34615631/2‑68 (MQ=255)
                                     ccTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGc                                >  1:9454697/1‑68 (MQ=255)
                                        ggCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCaaa                             <  1:32835500/68‑1 (MQ=255)
                                        ggCGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCaag                             <  1:14493171/68‑2 (MQ=255)
                                          cGTGATTAAGCGACCGGCCAGAGCAAAGAACAGACTCAAAGGCACCGCGTTGTGCTGTTCGGGaa                              >  1:3563691/1‑65 (MQ=255)
                                          cGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCaa                              >  1:35002699/1‑65 (MQ=255)
                                          cGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCaagc                            <  1:7856076/67‑3 (MQ=255)
                                          cGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCaagc                            <  1:721689/67‑3 (MQ=255)
                                           gTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCAAGCt                           <  1:4384841/67‑1 (MQ=255)
                                                aaaGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGa                     <  1:27242049/68‑1 (MQ=255)
                                                aaaGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGa                     <  1:19772547/68‑1 (MQ=255)
                                                  acagACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAg                   >  1:929101/4‑68 (MQ=255)
                                                  aGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAg                   >  1:4770634/1‑68 (MQ=255)
                                                  aGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAg                   >  1:37132072/1‑68 (MQ=255)
                                                  aGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAg                   >  1:32945558/1‑68 (MQ=255)
                                                   gCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGt                  >  1:1113593/1‑68 (MQ=255)
                                                    cGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTg                 >  1:14118499/1‑68 (MQ=255)
                                                     gACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTgg                <  1:33810418/68‑1 (MQ=255)
                                                          gctagaTGAAAGAACAGACTGAATGGCACGGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc           <  1:19966158/65‑1 (MQ=255)
                                                           cGAGATGAAAGATCAGACTGATTGGGACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc           <  1:16825893/67‑1 (MQ=255)
                                                           cGAGATGAAAGATCAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc           <  1:15322926/67‑1 (MQ=255)
                                                           cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc           <  1:15535597/67‑1 (MQ=255)
                                                           cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc           <  1:903475/67‑1 (MQ=255)
                                                           cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc           <  1:3099029/67‑1 (MQ=255)
                                                           cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc           <  1:10000779/67‑1 (MQ=255)
                                                           cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc           <  1:19289704/67‑1 (MQ=255)
                                                           cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc           <  1:35949505/67‑1 (MQ=255)
                                                           cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc           <  1:21528928/67‑1 (MQ=255)
                                                           cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc           <  1:2218600/67‑1 (MQ=255)
                                                           cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc           <  1:23581191/67‑1 (MQ=255)
                                                           cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc           <  1:25224734/67‑1 (MQ=255)
                                                           cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc           <  1:32611754/67‑1 (MQ=255)
                                                           cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc           <  1:29004817/67‑1 (MQ=255)
                                                           cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc           <  1:2974748/67‑1 (MQ=255)
                                                            gagaTGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagcga         >  1:6663306/1‑68 (MQ=255)
                                                                 gAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATc    <  1:33047327/68‑1 (MQ=255)
                                                                   aaGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCt   <  1:17177363/67‑1 (MQ=255)
                                                                   aaGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTg  <  1:20684362/68‑1 (MQ=255)
                                                                   |                                                                   
GGCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTG  >  CP009273/4264134‑4264268

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: