Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,264,202 A→G 6.7% K50K (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,2020AG 6.7% 1911.1 / 73.8 547K50K (AAA→AAG) ‡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 = 1.00e+00

Only 100 of 586 total aligned reads displayed.

GCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGC  >  CP009273/4264137‑4264269
                                                                 |                                                                   
gcAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAg                                                                    <  1:24639751/67‑1 (MQ=255)
gcAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGa                                                                   <  1:14019978/68‑1 (MQ=255)
gcAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGa                                                                   <  1:26524112/68‑1 (MQ=255)
  aGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAAc                                                                 >  1:35161815/1‑68 (MQ=255)
     tGCCAACATTACGCTGGCTACTTCCGAATCCTGGGGTGATAAAGGTACCGGCGAGATGAAAGAACAGa                                                              <  1:28675600/68‑1 (MQ=255)
     tGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAg                                                               <  1:27882697/67‑1 (MQ=255)
     tGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGa                                                              <  1:33251841/68‑1 (MQ=255)
      gCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGAc                                                             <  1:3838335/68‑1 (MQ=255)
        cAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACt                                                            <  1:33516063/67‑1 (MQ=255)
         aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                          <  1:30519216/68‑1 (MQ=255)
         aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                          <  1:22231068/68‑1 (MQ=255)
         aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                          <  1:2889593/68‑1 (MQ=255)
         aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                          <  1:13747569/68‑1 (MQ=255)
         aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                          <  1:14458879/68‑1 (MQ=255)
         aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                          <  1:4196123/68‑1 (MQ=255)
          aCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGCCCGGCGAGATGAAAGAACAGACTGaa                                                         >  1:15959532/1‑68 (MQ=255)
          aCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGaa                                                         >  1:19994384/1‑68 (MQ=255)
          aCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGaa                                                         >  1:12213227/1‑68 (MQ=255)
            aTTACGCTGGCTACTTCCGAATCCTGGCGTGGTAAAGCGACCGGCGAGATGAAAGAACAGACTGAATg                                                       >  1:34281471/1‑68 (MQ=255)
            aTTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGCCCGGCGAGATGAAAGAACAGACTGAATg                                                       >  1:12975867/1‑68 (MQ=255)
            aTTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATg                                                       <  1:31269432/68‑1 (MQ=255)
             ttACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATg                                                       >  1:14766916/1‑67 (MQ=255)
              tACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGc                                                     <  1:7675679/68‑1 (MQ=255)
              tACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGc                                                     <  1:9296576/68‑1 (MQ=255)
                 gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                                  >  1:27717901/1‑68 (MQ=255)
                 gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                                  >  1:7298494/1‑68 (MQ=255)
                 gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                                  >  1:10334437/1‑68 (MQ=255)
                 gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                                  >  1:4951063/1‑68 (MQ=255)
                  cTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                                  >  1:23239558/1‑67 (MQ=255)
                  cTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                                  >  1:33043171/1‑67 (MQ=255)
                  cTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                                  >  1:1568808/1‑67 (MQ=255)
                   tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                                <  1:22873230/68‑1 (MQ=255)
                   tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                                <  1:26725921/68‑1 (MQ=255)
                   tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                                <  1:14917654/68‑1 (MQ=255)
                   tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                                <  1:34001122/68‑1 (MQ=255)
                   tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                                <  1:17996752/68‑1 (MQ=255)
                   tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                                <  1:22599659/68‑1 (MQ=255)
                   tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                                <  1:18178612/68‑1 (MQ=255)
                   tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                                <  1:8372916/68‑1 (MQ=255)
                   tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                                <  1:1865464/68‑1 (MQ=255)
                   tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                                <  1:31820521/68‑1 (MQ=255)
                   tGGCTACTTCCGAATCCTGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                                <  1:17126992/67‑1 (MQ=255)
                   tGGCCACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                                <  1:34696529/68‑1 (MQ=255)
                      cTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGtt                                             >  1:19574497/1‑68 (MQ=255)
                       tACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGtt                                             <  1:10729172/67‑1 (MQ=255)
                        aCTTCCTAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                           <  1:1015993/68‑1 (MQ=255)
                        aCTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                           <  1:5625263/68‑1 (MQ=255)
                        aCTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                           <  1:16476936/68‑1 (MQ=255)
                        aCTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                           <  1:6503395/68‑1 (MQ=255)
                          ttCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGc                                         <  1:26148843/68‑1 (MQ=255)
                           tCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGaa                                                         >  1:36775632/1‑51 (MQ=255)
                           tCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                           >  1:28423686/1‑65 (MQ=255)
                           tCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                           >  1:29989562/1‑65 (MQ=255)
                           tCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                           >  1:17258643/1‑65 (MQ=255)
                           tCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCt                                        <  1:11401322/68‑1 (MQ=255)
                           tCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCt                                        <  1:8174399/68‑1 (MQ=255)
                             cGATTCCTGGCGTGATAAAGCGTCCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                      <  1:3640564/68‑1 (MQ=255)
                             cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                      <  1:29617819/68‑1 (MQ=255)
                             cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                      <  1:32220864/68‑1 (MQ=255)
                             cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                      <  1:31467502/68‑1 (MQ=255)
                             cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                      <  1:4602395/68‑1 (MQ=255)
                                aTCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCg                                   >  1:18984329/1‑68 (MQ=255)
                                aTCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCg                                   >  1:25283920/1‑68 (MQ=255)
                                aTCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCg                                   >  1:25564675/1‑68 (MQ=255)
                                  gcTGGCGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGc                                 >  1:28236892/2‑68 (MQ=255)
                                  ccTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGc                                 >  1:132370/1‑68 (MQ=255)
                                       cGTGATTAAGCGACCGGCCAGAGCAAAGAACAGACTCAAAGGCACCGCGTTGTGCTGTTCGGGaa                               >  1:3563691/1‑65 (MQ=255)
                                       cGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCaa                               >  1:35990880/1‑65 (MQ=255)
                                       cGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCaagc                             <  1:721689/67‑3 (MQ=255)
                                        gTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTg                           <  1:9582837/68‑1 (MQ=255)
                                        gTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTg                           <  1:21742812/68‑1 (MQ=255)
                                        gTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCAAGCtc                           <  1:30312965/68‑2 (MQ=255)
                                          gATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTgg                          <  1:15413461/67‑1 (MQ=255)
                                          gATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTgg                          <  1:17474175/67‑1 (MQ=255)
                                             aaaGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGa                      <  1:27242049/68‑1 (MQ=255)
                                             aaaGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGa                      <  1:19772547/68‑1 (MQ=255)
                                              aaGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGaa                     >  1:16752959/1‑68 (MQ=255)
                                              aaGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGaa                     >  1:22092452/1‑68 (MQ=255)
                                               acagACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAg                    <  1:32842243/65‑1 (MQ=255)
                                               aGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAg                    >  1:32945558/1‑68 (MQ=255)
                                               aGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAg                    <  1:11515882/68‑1 (MQ=255)
                                               aGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAg                    >  1:20899939/1‑68 (MQ=255)
                                                gCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGt                   >  1:1113593/1‑68 (MQ=255)
                                                 agACGGGGGAGATGCGGGAGCAGAAAGAATGGCATCGCGTGGTGCTGTTCGGCAAGCACGCCGAAGTg                  >  1:6080121/2‑68 (MQ=255)
                                                  gACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTgg                 <  1:28968911/68‑1 (MQ=255)
                                                  gACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTg                  >  1:21266614/1‑67 (MQ=255)
                                                     cGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcga              <  1:5160478/68‑1 (MQ=255)
                                                        cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc            <  1:24135675/67‑1 (MQ=255)
                                                        cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc            <  1:10000779/67‑1 (MQ=255)
                                                        cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc            <  1:32611754/67‑1 (MQ=255)
                                                        cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc            <  1:12523483/67‑1 (MQ=255)
                                                        cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc            <  1:2974748/67‑1 (MQ=255)
                                                        cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGAGCTGTTCGGCAAACTGGCAGAAGTGgcgagc            <  1:26887547/67‑1 (MQ=255)
                                                         gagaTGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagcga          >  1:6663306/1‑68 (MQ=255)
                                                         gagaTGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagcga          <  1:24958518/68‑1 (MQ=255)
                                                           gaTGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCTaa         <  1:20699979/67‑1 (MQ=255)
                                                                aaGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCTAGCGAATATCTg   <  1:20563111/68‑1 (MQ=255)
                                                                 aGAACAGACTGATTGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc  >  1:23866499/1‑68 (MQ=255)
                                                                 aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc  <  1:31024564/68‑1 (MQ=255)
                                                                 aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc  <  1:8909076/68‑1 (MQ=255)
                                                                 |                                                                   
GCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGC  >  CP009273/4264137‑4264269

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: