Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,264,205 A→G 6.7% E51E (GAA→GAG ssb → single‑stranded DNA‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,264,2050AG 6.7% 1961.1 / 77.8 568E51E (GAA→GAGssbsingle‑stranded DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (193/334);  new base G (25/13);  total (218/350)
Fisher's exact test for biased strand distribution p-value = 4.92e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 593 total aligned reads displayed.

AGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTA  >  CP009273/4264139‑4264272
                                                                  |                                                                   
   tGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAg                                                                  <  1:27882697/67‑1 (MQ=255)
   tGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGa                                                                 <  1:28820285/68‑1 (MQ=255)
   tGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGa                                                                 <  1:4076132/68‑1 (MQ=255)
   tGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGa                                                                 <  1:33251841/68‑1 (MQ=255)
      cAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACt                                                               <  1:31483462/67‑1 (MQ=255)
      cAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACt                                                               <  1:31279825/67‑1 (MQ=255)
      cAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACt                                                               <  1:35104009/67‑1 (MQ=255)
      cAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACt                                                               <  1:33516063/67‑1 (MQ=255)
       aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAGCAGACTGa                                                             <  1:24479179/68‑1 (MQ=255)
       aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                             <  1:1818302/68‑1 (MQ=255)
       aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                             <  1:30519216/68‑1 (MQ=255)
       aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                             <  1:11737707/68‑1 (MQ=255)
       aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                             <  1:25214471/68‑1 (MQ=255)
        aCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGGGAGATGAAAGAACAGACTGaa                                                            >  1:17503371/1‑68 (MQ=255)
        aCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGaa                                                            >  1:7621000/1‑68 (MQ=255)
          aTTACGCTGGCTACTTCCGAATCCTGGCGTGGTAAAGCGACCGGCGAGATGAAAGAACAGACTGAATg                                                          >  1:34281471/1‑68 (MQ=255)
          aTTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATg                                                          >  1:19068085/1‑68 (MQ=255)
          aTTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACAGAATg                                                          >  1:25472802/1‑68 (MQ=255)
            tACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGc                                                        <  1:13887516/68‑1 (MQ=255)
            tACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGc                                                        <  1:24116944/68‑1 (MQ=255)
            tACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGc                                                        <  1:11446029/68‑1 (MQ=255)
            tACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGc                                                        <  1:9296576/68‑1 (MQ=255)
               gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                                     >  1:24954678/1‑68 (MQ=255)
               gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                                     >  1:19692313/1‑68 (MQ=255)
               gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                                     >  1:26320800/1‑68 (MQ=255)
               gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                                     >  1:27717901/1‑68 (MQ=255)
               gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                                     >  1:29710023/1‑68 (MQ=255)
               gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                                     >  1:14896950/1‑68 (MQ=255)
                cTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                                     >  1:1263191/1‑67 (MQ=255)
                cTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                                     >  1:33043171/1‑67 (MQ=255)
                 tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                                   <  1:34001122/68‑1 (MQ=255)
                 tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                                   <  1:22599659/68‑1 (MQ=255)
                 tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                                   <  1:22233825/68‑1 (MQ=255)
                 tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                                   <  1:21753280/68‑1 (MQ=255)
                 tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                                   <  1:15286694/68‑1 (MQ=255)
                 tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                                   <  1:6404987/68‑1 (MQ=255)
                 tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                                   <  1:20609729/68‑1 (MQ=255)
                 tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                                   <  1:5585478/68‑1 (MQ=255)
                 tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                                   <  1:5156641/68‑1 (MQ=255)
                 tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                                   <  1:25880472/68‑1 (MQ=255)
                 tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                                   <  1:18003373/68‑1 (MQ=255)
                 tGGCTACTTCCGAATCCTGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                                   <  1:17126992/67‑1 (MQ=255)
                 tGGCCACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                                   <  1:34696529/68‑1 (MQ=255)
                    cTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGCACAGACTGAATGGCACCGCGtt                                                >  1:4421484/1‑68 (MQ=255)
                    cTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAAGAGACTGAATGGCACCGCGtt                                                >  1:31838679/1‑68 (MQ=255)
                    cTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGtt                                                >  1:7278486/1‑68 (MQ=255)
                    cTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGtt                                                >  1:2376056/1‑68 (MQ=255)
                    cTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGt                                                 <  1:10397247/67‑1 (MQ=255)
                      aCTTCCTAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                              <  1:1015993/68‑1 (MQ=255)
                      aCTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                              <  1:16476936/68‑1 (MQ=255)
                      aCTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                              <  1:8094970/68‑1 (MQ=255)
                      aCTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                              <  1:28668314/68‑1 (MQ=255)
                         tCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                              >  1:13088736/1‑65 (MQ=255)
                         tCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCt                                           <  1:6817644/68‑1 (MQ=255)
                         tCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCt                                           <  1:10889454/68‑1 (MQ=255)
                           cGATTCCTGGCGTGATAAAGCGTCCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                         <  1:3640564/68‑1 (MQ=255)
                           cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                         <  1:26852143/68‑1 (MQ=255)
                           cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                         <  1:14238947/68‑1 (MQ=255)
                           cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                         <  1:19852452/68‑1 (MQ=255)
                              aTCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTc                                       <  1:13425435/67‑1 (MQ=255)
                                gcTGGCGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGc                                    >  1:12261465/2‑68 (MQ=255)
                                     cGTGATTAAGCGACCGGCCAGAGCAAAGAACAGACTCAAAGGCACCGCGTTGTGCTGTTCGGGaa                                  >  1:3563691/1‑65 (MQ=255)
                                     cGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCaa                                  >  1:11217160/1‑65 (MQ=255)
                                     cGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCaa                                  >  1:35990880/1‑65 (MQ=255)
                                     cGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCaagc                                >  1:23195732/1‑65 (MQ=255)
                                      gTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTg                              <  1:9582837/68‑1 (MQ=255)
                                      gTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTg                              <  1:18658727/68‑1 (MQ=255)
                                      gTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCAAGCtc                              <  1:30312965/68‑2 (MQ=255)
                                        gATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTgg                             <  1:6632202/67‑1 (MQ=255)
                                          tAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCa                           <  1:16070140/67‑1 (MQ=255)
                                           aaaGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGTAGa                         >  1:17286433/1‑68 (MQ=255)
                                           aaaGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGa                         >  1:371390/1‑68 (MQ=255)
                                           aaaGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGa                         <  1:27242049/68‑1 (MQ=255)
                                            aaGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGaa                        >  1:16752959/1‑68 (MQ=255)
                                            aaGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGaa                        >  1:2213796/1‑68 (MQ=255)
                                            aaGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGaa                        >  1:29968255/1‑68 (MQ=255)
                                             acagACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAg                       <  1:32842243/65‑1 (MQ=255)
                                             aGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAg                       >  1:32945558/1‑68 (MQ=255)
                                                gACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTgg                    <  1:2128167/68‑1 (MQ=255)
                                                gACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTgg                    <  1:28968911/68‑1 (MQ=255)
                                                gACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTg                     >  1:4839598/1‑67 (MQ=255)
                                                      cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc               <  1:14700323/67‑1 (MQ=255)
                                                      cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc               <  1:15535597/67‑1 (MQ=255)
                                                      cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc               <  1:592854/67‑1 (MQ=255)
                                                      cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc               <  1:29492483/67‑1 (MQ=255)
                                                      cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc               <  1:32292608/67‑1 (MQ=255)
                                                      cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc               <  1:10000779/67‑1 (MQ=255)
                                                      cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc               <  1:20192782/67‑1 (MQ=255)
                                                      cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc               <  1:20734894/67‑1 (MQ=255)
                                                      cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc               <  1:32611754/67‑1 (MQ=255)
                                                       gagaTGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagcga             <  1:22819566/68‑1 (MQ=255)
                                                               aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc     <  1:31024564/68‑1 (MQ=255)
                                                               aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc     <  1:8335246/68‑1 (MQ=255)
                                                               aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc     <  1:8909076/68‑1 (MQ=255)
                                                               aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc     <  1:26721471/68‑1 (MQ=255)
                                                               aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc     <  1:28423555/68‑1 (MQ=255)
                                                               aGAACAGACTGAATGGAACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc     <  1:20912912/68‑1 (MQ=255)
                                                                 aaCAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGt   >  1:3676653/1‑68 (MQ=255)
                                                                  aCAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTa  >  1:15804879/1‑68 (MQ=255)
                                                                  aCAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTa  >  1:28198128/1‑68 (MQ=255)
                                                                  |                                                                   
AGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTA  >  CP009273/4264139‑4264272

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: