Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,264,202 A→G 8.8% 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 8.8% 1777.1 / 103.1 535K50K (AAA→AAG) ‡ssbsingle‑stranded DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (211/275);  new base G (32/15);  total (243/292)
Fisher's exact test for biased strand distribution p-value = 1.88e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 564 total aligned reads displayed.

CGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGC  >  CP009273/4264136‑4264269
                                                                  |                                                                   
cgcAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAg                                                                    >  1:18595065/1‑68 (MQ=255)
 gcAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAg                                                                    <  1:22137225/67‑1 (MQ=255)
 gcAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAg                                                                    <  1:15833883/67‑1 (MQ=255)
 gcAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAg                                                                    <  1:16880758/67‑1 (MQ=255)
 gcAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGa                                                                   >  1:38620603/1‑68 (MQ=255)
   aGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAAc                                                                 >  1:20697663/1‑68 (MQ=255)
   aGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAAc                                                                 >  1:37429647/1‑68 (MQ=255)
      tGCCAGCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGa                                                              <  1:37776094/68‑1 (MQ=255)
      tGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAg                                                               <  1:11866544/67‑1 (MQ=255)
       gCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGAc                                                             <  1:2951978/68‑1 (MQ=255)
         cAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAATGAACAGACt                                                            <  1:22298790/67‑1 (MQ=255)
         cAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACt                                                            <  1:39631023/67‑1 (MQ=255)
         cAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACt                                                            <  1:17772068/67‑1 (MQ=255)
         cAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACt                                                            <  1:8107808/67‑1 (MQ=255)
         cAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACt                                                            <  1:15657461/67‑1 (MQ=255)
          aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                          <  1:4183329/68‑1 (MQ=255)
          aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                          <  1:5125372/68‑1 (MQ=255)
          aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                          <  1:16243843/68‑1 (MQ=255)
          aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                          <  1:38128488/68‑1 (MQ=255)
          aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                          <  1:35558852/68‑1 (MQ=255)
          aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                          <  1:35200708/68‑1 (MQ=255)
          aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                          <  1:34644862/68‑1 (MQ=255)
           aCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGaa                                                         >  1:2876267/1‑68 (MQ=255)
           aCATTACGCTGGCTACTTCCGAATCCTCGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGaa                                                         >  1:12068442/1‑68 (MQ=255)
               tACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGc                                                     <  1:30775358/68‑1 (MQ=255)
               tACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGc                                                     <  1:19959670/68‑1 (MQ=255)
               tACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGc                                                     <  1:33551752/68‑1 (MQ=255)
                  gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                                  >  1:2934243/1‑68 (MQ=255)
                  gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                                  >  1:21030707/1‑68 (MQ=255)
                  gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                                  >  1:29940034/1‑68 (MQ=255)
                  gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                                  >  1:7121466/1‑68 (MQ=255)
                  gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                                  >  1:4405805/1‑68 (MQ=255)
                  gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                                  >  1:23665074/1‑68 (MQ=255)
                  gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                                  >  1:31186031/1‑68 (MQ=255)
                   cTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                                  >  1:10545012/1‑67 (MQ=255)
                   cTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                                  >  1:7750946/1‑67 (MQ=255)
                   cTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                                  >  1:25594848/1‑67 (MQ=255)
                    tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                                <  1:30625252/68‑1 (MQ=255)
                    tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                                <  1:11001690/68‑1 (MQ=255)
                    tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                                <  1:12392292/68‑1 (MQ=255)
                    tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                                <  1:26663670/68‑1 (MQ=255)
                    tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                                <  1:25173421/68‑1 (MQ=255)
                    tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                                <  1:12999671/68‑1 (MQ=255)
                    tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                                <  1:22465396/68‑1 (MQ=255)
                    tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                                <  1:32271781/68‑1 (MQ=255)
                    tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                                <  1:39111412/68‑1 (MQ=255)
                    tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                                <  1:5395130/68‑1 (MQ=255)
                     ggCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgcg                                               <  1:2833609/68‑1 (MQ=255)
                       cTACTTCCGTATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACATACTTAATGGCACCGCGt                                              <  1:28058027/67‑1 (MQ=255)
                       cTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGt                                              <  1:31677443/67‑1 (MQ=255)
                        tACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGtt                                             <  1:620189/67‑1 (MQ=255)
                         aCTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                           <  1:34940712/68‑1 (MQ=255)
                         aCTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                           <  1:21987861/68‑1 (MQ=255)
                         aCTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                           <  1:17353130/68‑1 (MQ=255)
                         aCTTCCGAATCCTGGCGGGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                           <  1:8600622/68‑1 (MQ=255)
                            tCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                           >  1:13755039/1‑65 (MQ=255)
                            tCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                           >  1:16109189/1‑65 (MQ=255)
                              cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                      <  1:27772391/68‑1 (MQ=255)
                              cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                      <  1:26324195/68‑1 (MQ=255)
                              cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                      <  1:33845979/68‑1 (MQ=255)
                              cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                      >  1:3682611/1‑68 (MQ=255)
                              cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                      <  1:14312547/68‑1 (MQ=255)
                                aaTCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTc                                    >  1:8921264/1‑68 (MQ=255)
                                 aTCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCg                                   >  1:10236313/1‑68 (MQ=255)
                                 aTCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCg                                   >  1:25726945/1‑68 (MQ=255)
                                 aTCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCg                                   >  1:23024125/1‑68 (MQ=255)
                                 aTCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGGTGTGCTGTTc                                    <  1:15475318/67‑1 (MQ=255)
                                  agcTGGCGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCgg                                  <  1:9607880/66‑1 (MQ=255)
                                   gcTGGCGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGc                                 >  1:10078067/2‑68 (MQ=255)
                                   gcTGGCGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGc                                 >  1:12772543/2‑68 (MQ=255)
                                   gcTGGCGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGc                                 >  1:36431113/2‑68 (MQ=255)
                                   gcTGGCGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGc                                 >  1:21387644/2‑68 (MQ=255)
                                   ccTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGc                                 >  1:27169964/1‑68 (MQ=255)
                                        cGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCaa                               >  1:30159347/1‑65 (MQ=255)
                                        cGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAAc                             <  1:13358568/67‑1 (MQ=255)
                                        cGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCaagc                             >  1:38902939/1‑65 (MQ=255)
                                         gTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTg                           >  1:5001793/1‑68 (MQ=255)
                                         gTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTg                           <  1:18832458/68‑1 (MQ=255)
                                         gTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTg                           <  1:14679753/68‑1 (MQ=255)
                                           gATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTgg                          <  1:24584936/67‑1 (MQ=255)
                                           gATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTgg                          <  1:18310271/67‑1 (MQ=255)
                                             tAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAg                       >  1:37019754/1‑68 (MQ=255)
                                                acagACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAg                    >  1:26038428/4‑68 (MQ=255)
                                                 gCGACCGGCGAGATGAAAGAACAGACTGATTGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGt                   >  1:17316317/1‑68 (MQ=255)
                                                 gCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGt                   >  1:20432469/1‑68 (MQ=255)
                                                 gCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGt                   >  1:4743179/1‑68 (MQ=255)
                                                  cGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAACCTGGCAGAAGTg                  >  1:7517820/1‑68 (MQ=255)
                                                  cGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTg                  >  1:32682846/1‑68 (MQ=255)
                                                  agACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTg                  >  1:6514919/2‑68 (MQ=255)
                                                  agACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTg                  >  1:9327093/2‑68 (MQ=255)
                                                     cgggCGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTg                  <  1:1964432/63‑1 (MQ=255)
                                                         cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc            <  1:5794612/67‑1 (MQ=255)
                                                         cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc            <  1:38850468/67‑1 (MQ=255)
                                                              tGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAata       <  1:23222216/67‑1 (MQ=255)
                                                                 aaGAACAGACTGAATGGCACCGCGTTGTGATGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTg   <  1:33159092/68‑1 (MQ=255)
                                                                  aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc  <  1:19391251/68‑1 (MQ=255)
                                                                  aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc  <  1:20801491/68‑1 (MQ=255)
                                                                  aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc  <  1:24451890/68‑1 (MQ=255)
                                                                  aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc  <  1:35947118/68‑1 (MQ=255)
                                                                  aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc  <  1:24813003/68‑1 (MQ=255)
                                                                  |                                                                   
CGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGC  >  CP009273/4264136‑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: