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

Only 100 of 541 total aligned reads displayed.

GGCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATC  >  CP009273/4264134‑4264266
                                                                  |                                                                  
  cgcAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAg                                                                 >  1:32116219/1‑68 (MQ=255)
  cgcAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAg                                                                 >  1:23385144/1‑68 (MQ=255)
   gcAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAg                                                                 <  1:15833883/67‑1 (MQ=255)
   gcAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGa                                                                >  1:38620603/1‑68 (MQ=255)
   gcAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGa                                                                <  1:27612655/68‑1 (MQ=255)
   gcAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGa                                                                >  1:8286448/1‑68 (MQ=255)
    cAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGa                                                                >  1:12034315/1‑67 (MQ=255)
        tGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAg                                                            <  1:16916219/67‑1 (MQ=255)
        tGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAg                                                            <  1:16103367/67‑1 (MQ=255)
        tGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAg                                                            <  1:13321311/67‑1 (MQ=255)
        tGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGa                                                           <  1:39316797/68‑1 (MQ=255)
           cAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACt                                                         <  1:39631023/67‑1 (MQ=255)
           cAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACt                                                         <  1:4497522/67‑1 (MQ=255)
           cAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACt                                                         <  1:15657461/67‑1 (MQ=255)
            aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                       <  1:35554196/68‑1 (MQ=255)
            aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                       <  1:38128488/68‑1 (MQ=255)
            aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                       <  1:20770260/68‑1 (MQ=255)
            aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                       <  1:35200708/68‑1 (MQ=255)
            aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                       <  1:5125372/68‑1 (MQ=255)
            aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                       <  1:735361/68‑1 (MQ=255)
            aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                       <  1:3013834/68‑1 (MQ=255)
            aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                       <  1:22916722/68‑1 (MQ=255)
            aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                       <  1:32866075/68‑1 (MQ=255)
             aCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGaa                                                      >  1:16227126/1‑68 (MQ=255)
             aCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGaa                                                      >  1:15144807/1‑68 (MQ=255)
             aCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGaa                                                      >  1:11589520/1‑68 (MQ=255)
             aCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGAATGaa                                                      >  1:17823267/1‑68 (MQ=255)
               aTTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATg                                                    >  1:30748957/1‑68 (MQ=255)
               aTTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATg                                                    >  1:32566536/1‑68 (MQ=255)
                 tACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGc                                                  <  1:19959670/68‑1 (MQ=255)
                 tACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGc                                                  <  1:19426916/68‑1 (MQ=255)
                 tACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGc                                                  <  1:33551752/68‑1 (MQ=255)
                 tACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCCAGATGAAAGAACAGACTGAATGGc                                                  <  1:28324822/68‑1 (MQ=255)
                    gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                               >  1:29321980/1‑68 (MQ=255)
                    gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                               >  1:37843650/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)
                     cTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACTGACTGAATGGCACcg                                              >  1:18325875/1‑68 (MQ=255)
                     cTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                               >  1:25594848/1‑67 (MQ=255)
                      tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                             <  1:28705492/68‑1 (MQ=255)
                      tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                             <  1:26663670/68‑1 (MQ=255)
                      tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                             <  1:21693267/68‑1 (MQ=255)
                      tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                             <  1:25173421/68‑1 (MQ=255)
                      tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                             <  1:11001690/68‑1 (MQ=255)
                      tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                             <  1:5395130/68‑1 (MQ=255)
                      tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                             <  1:20346903/68‑1 (MQ=255)
                      tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                             <  1:4772044/68‑1 (MQ=255)
                          tACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGtt                                          <  1:620189/67‑1 (MQ=255)
                           aCTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                        <  1:18640835/68‑1 (MQ=255)
                           aCTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                        <  1:34940712/68‑1 (MQ=255)
                              tCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                        >  1:22357720/1‑65 (MQ=255)
                              tCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                        >  1:14661682/1‑65 (MQ=255)
                              tCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                        >  1:13755039/1‑65 (MQ=255)
                              tCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCt                                     <  1:10445079/68‑1 (MQ=255)
                                cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                   <  1:21360490/68‑1 (MQ=255)
                                cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                   <  1:37330179/68‑1 (MQ=255)
                                cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                   <  1:26324195/68‑1 (MQ=255)
                                cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                   >  1:3682611/1‑68 (MQ=255)
                                cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                   <  1:1222150/68‑1 (MQ=255)
                                cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                   <  1:33845979/68‑1 (MQ=255)
                                  aaTCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTc                                 >  1:8921264/1‑68 (MQ=255)
                                   aTCCTTGCGTGATAAAGCGACCGGCGATATGAAAGAACAGACTGAATGGAACCGCGTTGTGCTGTTCg                                >  1:30307426/1‑68 (MQ=255)
                                   aTCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCg                                >  1:6660916/1‑68 (MQ=255)
                                   aTCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCg                                >  1:10236313/1‑68 (MQ=255)
                                   aTCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCg                                >  1:4199158/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)
                                     gcTGGCGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGc                              >  1:29848451/2‑68 (MQ=255)
                                     ccTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGc                              >  1:27169964/1‑68 (MQ=255)
                                        ggCGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCaa                            >  1:24191398/1‑67 (MQ=255)
                                          cGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCaa                            >  1:22263136/1‑65 (MQ=255)
                                          cGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCaagc                          >  1:34774578/1‑65 (MQ=255)
                                          cGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCaagc                          >  1:38902939/1‑65 (MQ=255)
                                           gTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTg                        >  1:23062370/1‑68 (MQ=255)
                                           gTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTg                        <  1:18832458/68‑1 (MQ=255)
                                           gTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTg                        >  1:5001793/1‑68 (MQ=255)
                                             gATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTgg                       <  1:17426523/67‑1 (MQ=255)
                                               tAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAg                    >  1:37019754/1‑68 (MQ=255)
                                                  acagACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAg                 >  1:26038428/4‑68 (MQ=255)
                                                  acagACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAg                 <  1:35945244/65‑1 (MQ=255)
                                                  aGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAg                 >  1:2452551/1‑68 (MQ=255)
                                                   gCGACCGGCGAGATGAAAGAACAGACTGATTGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGt                >  1:17316317/1‑68 (MQ=255)
                                                   gCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGGAAACTGGCAGAAGt                >  1:36425123/1‑68 (MQ=255)
                                                    agACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTg               >  1:12584852/2‑68 (MQ=255)
                                                     gACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGt                <  1:2481810/66‑1 (MQ=255)
                                                       cgggCGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTg               <  1:1964432/63‑1 (MQ=255)
                                                        cGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcga           <  1:20591527/68‑1 (MQ=255)
                                                        cGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcga           <  1:20996704/68‑1 (MQ=255)
                                                         ggCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgag          <  1:872502/68‑1 (MQ=255)
                                                           cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc         <  1:5794612/67‑1 (MQ=255)
                                                           cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc         <  1:38850468/67‑1 (MQ=255)
                                                           cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc         <  1:7573596/67‑1 (MQ=255)
                                                           cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc         <  1:25756272/67‑1 (MQ=255)
                                                           cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc         <  1:31134358/67‑1 (MQ=255)
                                                           cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc         <  1:9328039/67‑1 (MQ=255)
                                                           cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc         <  1:31586514/67‑1 (MQ=255)
                                                            gagaTGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagcga       <  1:27869464/68‑1 (MQ=255)
                                                            gagaTGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagcga       <  1:29459729/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: