Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,264,200 A→C 6.7% K50Q (AAA→CAA) ‡ ssb → single‑stranded DNA‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,264,2000AC 6.7% 1978.5 / 90.7 547K50Q (AAA→CAA) ‡ssbsingle‑stranded DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (185/325);  new base C (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 = 1.00e+00

Only 100 of 558 total aligned reads displayed.

GGCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATC  >  CP009273/4264134‑4264266
                                                                  |                                                                  
ggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGaa                                                                   <  1:21618851/68‑1 (MQ=255)
ggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGaa                                                                   <  1:585870/68‑1 (MQ=255)
ggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGaa                                                                   <  1:28506933/68‑1 (MQ=255)
   gcAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAg                                                                 <  1:8382636/67‑1 (MQ=255)
   gcAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAg                                                                 <  1:20433307/67‑1 (MQ=255)
   gcAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAg                                                                 <  1:33273667/67‑1 (MQ=255)
   gcAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAg                                                                 <  1:16781108/67‑1 (MQ=255)
   gcAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGa                                                                <  1:3653513/68‑1 (MQ=255)
   gcAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGa                                                                <  1:24123178/68‑1 (MQ=255)
   gcAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGa                                                                <  1:14019978/68‑1 (MQ=255)
   gcAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGa                                                                >  1:5043938/1‑68 (MQ=255)
     aGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAAc                                                              >  1:22217641/1‑68 (MQ=255)
        tGCCAACATTACGCTGGCTACTTCCGAATCCTGGGGTGATAAAGGTACCGGCGAGATGAAAGAACAGa                                                           <  1:28675600/68‑1 (MQ=255)
        tGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGa                                                           <  1:4076132/68‑1 (MQ=255)
           cAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACt                                                         <  1:35104009/67‑1 (MQ=255)
           cAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACt                                                         <  1:31483462/67‑1 (MQ=255)
           cAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACt                                                         <  1:31279825/67‑1 (MQ=255)
            aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                       <  1:26341505/68‑1 (MQ=255)
            aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                       <  1:11508117/68‑1 (MQ=255)
            aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                       <  1:23719703/68‑1 (MQ=255)
            aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                       <  1:14458879/68‑1 (MQ=255)
            aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                       <  1:25542497/68‑1 (MQ=255)
             aCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGaa                                                      >  1:28920261/1‑68 (MQ=255)
               aTTACGCTGGCTACTTCCGAATCCTGGCGTGGTAAAGCGACCGGCGAGATGAAAGAACAGACTGAATg                                                    >  1:34281471/1‑68 (MQ=255)
               aTTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATg                                                    >  1:6145859/1‑68 (MQ=255)
                ttACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATg                                                    >  1:14766916/1‑67 (MQ=255)
                 tACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGc                                                  <  1:16662874/68‑1 (MQ=255)
                 tACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGc                                                  <  1:2667374/68‑1 (MQ=255)
                 tACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGc                                                  <  1:20685991/68‑1 (MQ=255)
                    gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                               >  1:27717901/1‑68 (MQ=255)
                    gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                               >  1:7298494/1‑68 (MQ=255)
                    gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                               >  1:30415521/1‑68 (MQ=255)
                    gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAc                                                >  1:17713685/1‑67 (MQ=255)
                    gCTGGCTACTACCGAATCCTGGCGTGATAAAGCGCCCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                               >  1:29619026/1‑68 (MQ=255)
                     cTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                               >  1:5398170/1‑67 (MQ=255)
                      tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                             <  1:18470441/68‑1 (MQ=255)
                      tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                             <  1:25880472/68‑1 (MQ=255)
                      tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                             <  1:34001122/68‑1 (MQ=255)
                      tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                             <  1:17398987/68‑1 (MQ=255)
                      tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                             <  1:1584104/68‑1 (MQ=255)
                      tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                             <  1:15286694/68‑1 (MQ=255)
                      tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                             <  1:12585987/68‑1 (MQ=255)
                      tGGCTACTTCCGAATCCTGGCGGGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                             <  1:13801628/68‑1 (MQ=255)
                      tGGCCACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                             <  1:34696529/68‑1 (MQ=255)
                         cTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAAGAGACTGAATGGCACCGCGtt                                          >  1:31838679/1‑68 (MQ=255)
                         cTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGtt                                          >  1:27406608/1‑68 (MQ=255)
                         cTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGt                                           <  1:2465055/67‑1 (MQ=255)
                         cTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGt                                           <  1:10397247/67‑1 (MQ=255)
                          tACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGtt                                          <  1:18909475/67‑1 (MQ=255)
                          tACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGtt                                          <  1:10729172/67‑1 (MQ=255)
                           aCTTCCTAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                        <  1:1015993/68‑1 (MQ=255)
                           aCTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                        <  1:11354846/68‑1 (MQ=255)
                              tCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGCACAGACTGAATGGCACCGCGTtgt                                        >  1:32983594/1‑65 (MQ=255)
                              tCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                        >  1:13088736/1‑65 (MQ=255)
                              tCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                        >  1:29989562/1‑65 (MQ=255)
                                cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                   <  1:18074909/68‑1 (MQ=255)
                                cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                   <  1:26852143/68‑1 (MQ=255)
                                cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                   <  1:32864091/68‑1 (MQ=255)
                                   aTCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCg                                >  1:28261817/1‑68 (MQ=255)
                                   aTCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCg                                >  1:25283920/1‑68 (MQ=255)
                                   aTCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCg                                >  1:33083092/1‑68 (MQ=255)
                                     gcTGGCGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGc                              >  1:12261465/2‑68 (MQ=255)
                                     ccTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGc                              >  1:9454697/1‑68 (MQ=255)
                                     ccTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGc                              >  1:13279884/1‑68 (MQ=255)
                                      cTGGCGTGACAAACAGACGGGGGAGATGCGGAAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCa                             >  1:36384751/1‑68 (MQ=255)
                                        ggCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCaa                            >  1:35638760/1‑67 (MQ=255)
                                        ggCGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCaag                           <  1:1612609/68‑2 (MQ=255)
                                          cGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCaa                            >  1:35990880/1‑65 (MQ=255)
                                          cGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCaagc                          <  1:721689/67‑3 (MQ=255)
                                          cGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCaagc                          <  1:17130693/67‑3 (MQ=255)
                                          cGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCaagc                          <  1:7856076/67‑3 (MQ=255)
                                           gTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTg                        <  1:8224133/68‑1 (MQ=255)
                                           gTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCAAGCtc                        <  1:30312965/68‑2 (MQ=255)
                                           gTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCAAGCt                         <  1:4384841/67‑1 (MQ=255)
                                             gATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTgg                       <  1:2303447/67‑1 (MQ=255)
                                              aTAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCa                     >  1:2928947/1‑68 (MQ=255)
                                               tAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAg                    >  1:22552346/1‑68 (MQ=255)
                                                aaaGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGa                   <  1:19772547/68‑1 (MQ=255)
                                                aaaGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGa                   >  1:371390/1‑68 (MQ=255)
                                                 aaGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGaa                  >  1:22092452/1‑68 (MQ=255)
                                                  acagACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAg                 >  1:6525216/4‑68 (MQ=255)
                                                  acagACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAg                 >  1:929101/4‑68 (MQ=255)
                                                  aGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAg                 >  1:32704909/1‑68 (MQ=255)
                                                  aGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAg                 >  1:20899939/1‑68 (MQ=255)
                                                   gCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGt                >  1:1113593/1‑68 (MQ=255)
                                                     gACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTgg              <  1:2128167/68‑1 (MQ=255)
                                                     gACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTg               >  1:4839598/1‑67 (MQ=255)
                                                        cGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcga           <  1:29761817/68‑1 (MQ=255)
                                                          gctagaTGAAAGAACAGACTGAATGGCACGGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc         <  1:19966158/65‑1 (MQ=255)
                                                           cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc         <  1:10000779/67‑1 (MQ=255)
                                                           cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc         <  1:903475/67‑1 (MQ=255)
                                                           cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc         <  1:19289704/67‑1 (MQ=255)
                                                           cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc         <  1:15535597/67‑1 (MQ=255)
                                                           cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc         <  1:33517494/67‑1 (MQ=255)
                                                           cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc         <  1:32292608/67‑1 (MQ=255)
                                                           cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc         <  1:3099029/67‑1 (MQ=255)
                                                            gagaTGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagcga       >  1:6663306/1‑68 (MQ=255)
                                                            gagaTGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagcga       <  1:22819566/68‑1 (MQ=255)
                                                            gagaTGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagcga       <  1:24958518/68‑1 (MQ=255)
                                                            gagaTGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagcga       <  1:27984693/68‑1 (MQ=255)
                                                                  |                                                                  
GGCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATC  >  CP009273/4264134‑4264266

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: