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

Only 100 of 544 total aligned reads displayed.

GGCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTG  >  CP009273/4264134‑4264268
                                                                   |                                                                   
ggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGaa                                                                     >  1:4927365/1‑68 (MQ=255)
   gcAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAg                                                                   <  1:16880758/67‑1 (MQ=255)
   gcAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGa                                                                  <  1:27672594/68‑1 (MQ=255)
   gcAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGa                                                                  >  1:8286448/1‑68 (MQ=255)
   gcAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGa                                                                  >  1:38620603/1‑68 (MQ=255)
   gcAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGa                                                                  <  1:3701664/68‑1 (MQ=255)
    cAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGa                                                                  >  1:12034315/1‑67 (MQ=255)
    cAGTTGCCAACATTACGCTGGCTACTTACGAATCCTGGCGTGATATAGCGACCGGCGAGATGAAAGaa                                                                 <  1:19925193/68‑1 (MQ=255)
     aGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGaa                                                                 <  1:32613348/67‑1 (MQ=255)
     aGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAAc                                                                >  1:38770743/1‑68 (MQ=255)
        tGCCAGCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGa                                                             <  1:37776094/68‑1 (MQ=255)
        tGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAg                                                              <  1:16103367/67‑1 (MQ=255)
        tGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAg                                                              <  1:13321311/67‑1 (MQ=255)
        tGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGa                                                             <  1:2358398/68‑1 (MQ=255)
           cAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACt                                                           <  1:15657461/67‑1 (MQ=255)
           cAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACt                                                           <  1:39631023/67‑1 (MQ=255)
            aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                         <  1:26688199/68‑1 (MQ=255)
            aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                         <  1:4183329/68‑1 (MQ=255)
            aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                         <  1:4441417/68‑1 (MQ=255)
            aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                         <  1:20770260/68‑1 (MQ=255)
            aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                         <  1:22916722/68‑1 (MQ=255)
            aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                         <  1:5125372/68‑1 (MQ=255)
            aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                         <  1:26046017/68‑1 (MQ=255)
            aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                         <  1:735361/68‑1 (MQ=255)
            aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                         <  1:36413122/68‑1 (MQ=255)
            aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                         <  1:34644862/68‑1 (MQ=255)
             aCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGaa                                                        >  1:29864249/1‑68 (MQ=255)
             aCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGaa                                                        >  1:16227126/1‑68 (MQ=255)
             aCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGaa                                                        >  1:15144807/1‑68 (MQ=255)
             aCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGaa                                                        >  1:11589520/1‑68 (MQ=255)
             aCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGaa                                                        >  1:7564668/1‑68 (MQ=255)
             aCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGAATGaa                                                        >  1:17823267/1‑68 (MQ=255)
               aTTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATg                                                      >  1:30748957/1‑68 (MQ=255)
               aTTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATg                                                      >  1:29517767/1‑68 (MQ=255)
               aTTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATg                                                      >  1:30421097/1‑68 (MQ=255)
                 tACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGTAAGAACAGACTGAATGGc                                                    <  1:30147606/68‑1 (MQ=255)
                 tACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGc                                                    <  1:33551752/68‑1 (MQ=255)
                 tACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGc                                                    <  1:5341481/68‑1 (MQ=255)
                 tACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGc                                                    <  1:19426916/68‑1 (MQ=255)
                 tACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCCAGATGAAAGAACAGACTGAATGGc                                                    <  1:28324822/68‑1 (MQ=255)
                    gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                                 >  1:29321980/1‑68 (MQ=255)
                    gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                                 >  1:22088606/1‑68 (MQ=255)
                    gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                                 >  1:7857073/1‑68 (MQ=255)
                    gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                                 >  1:12984497/1‑68 (MQ=255)
                    gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                                 >  1:31186031/1‑68 (MQ=255)
                     cTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                                 >  1:25594848/1‑67 (MQ=255)
                      tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                               <  1:32271781/68‑1 (MQ=255)
                      tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                               <  1:25173421/68‑1 (MQ=255)
                      tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                               <  1:24577900/68‑1 (MQ=255)
                      tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                               <  1:39111412/68‑1 (MQ=255)
                      tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                               <  1:21693267/68‑1 (MQ=255)
                      tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                               <  1:4772044/68‑1 (MQ=255)
                      tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                               <  1:20346903/68‑1 (MQ=255)
                      tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                               <  1:11001690/68‑1 (MQ=255)
                          tACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGtt                                            <  1:620189/67‑1 (MQ=255)
                           aCTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                          <  1:18640835/68‑1 (MQ=255)
                           aCTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                          <  1:34922212/68‑1 (MQ=255)
                              tCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                          >  1:22357720/1‑65 (MQ=255)
                              tCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                          >  1:13755039/1‑65 (MQ=255)
                              tCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCt                                       <  1:10445079/68‑1 (MQ=255)
                                cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                     <  1:33845979/68‑1 (MQ=255)
                                cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                     <  1:35188106/68‑1 (MQ=255)
                                cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                     >  1:3682611/1‑68 (MQ=255)
                                cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                     <  1:37330179/68‑1 (MQ=255)
                                cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                     <  1:26324195/68‑1 (MQ=255)
                                cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                     <  1:21360490/68‑1 (MQ=255)
                                cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                     <  1:1222150/68‑1 (MQ=255)
                                  aaTCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTc                                   >  1:8921264/1‑68 (MQ=255)
                                   aTCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCg                                  >  1:6660916/1‑68 (MQ=255)
                                   aTCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCg                                  >  1:10236313/1‑68 (MQ=255)
                                   aTCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCg                                  >  1:15825465/1‑68 (MQ=255)
                                    tCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCg                                  >  1:14228075/1‑67 (MQ=255)
                                    agcTGGCGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCgg                                 <  1:9607880/66‑1 (MQ=255)
                                     gcTGGCGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGc                                >  1:10078067/2‑68 (MQ=255)
                                     ccTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGc                                >  1:27169964/1‑68 (MQ=255)
                                        ggCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCaaa                             <  1:28754357/68‑1 (MQ=255)
                                        ggCGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCaa                              >  1:24191398/1‑67 (MQ=255)
                                          cGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCaa                              >  1:22263136/1‑65 (MQ=255)
                                          cGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCaagc                            <  1:12701015/67‑3 (MQ=255)
                                          cGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCaagc                            <  1:27994499/67‑3 (MQ=255)
                                          cGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCaagc                            >  1:38902939/1‑65 (MQ=255)
                                           gTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTg                          <  1:35796866/68‑1 (MQ=255)
                                           gTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTg                          <  1:14679753/68‑1 (MQ=255)
                                           gTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTg                          <  1:18832458/68‑1 (MQ=255)
                                           gTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTg                          >  1:23062370/1‑68 (MQ=255)
                                             gATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTgg                         <  1:18310271/67‑1 (MQ=255)
                                                  aGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAg                   >  1:17400595/1‑68 (MQ=255)
                                                   gCGACCGGCGAGATGAAAGAACAGACTGATTGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGt                  >  1:17316317/1‑68 (MQ=255)
                                                   gCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGt                  >  1:35312853/1‑68 (MQ=255)
                                                     gACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGt                  <  1:2481810/66‑1 (MQ=255)
                                                     gACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTg                 >  1:38031744/1‑67 (MQ=255)
                                                       cgggCGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTg                 <  1:1964432/63‑1 (MQ=255)
                                                        cGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcga             <  1:20996704/68‑1 (MQ=255)
                                                        cGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcga             <  1:20591527/68‑1 (MQ=255)
                                                         ggCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgag            <  1:872502/68‑1 (MQ=255)
                                                           cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc           <  1:5794612/67‑1 (MQ=255)
                                                           cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc           <  1:25756272/67‑1 (MQ=255)
                                                           cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc           <  1:31586514/67‑1 (MQ=255)
                                                           cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc           <  1:9328039/67‑1 (MQ=255)
                                                                   aaGAACAGACTGAATGGCACCGCGTTGTGATGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTg  <  1:33159092/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: