Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,264,181 T→C 5.3% D43D (GAT→GAC ssb → single‑stranded DNA‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,264,1810TC 5.3% 2146.5 / 65.1 601D43D (GAT→GACssbsingle‑stranded DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (309/258);  new base C (22/10);  total (333/268)
Fisher's exact test for biased strand distribution p-value = 1.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01

Only 100 of 608 total aligned reads displayed.

TACGCTACATGCCAAATGGTGGCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAG  >  CP009273/4264114‑4264248
                                                                   |                                                                   
   gCTACATGCCAAATGGTGGCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATaaa                                                                  >  1:2248395/1‑68 (MQ=255)
   gCTACATGCCAAATGGTGGCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATaaa                                                                  >  1:17181284/1‑68 (MQ=255)
   gCTACATGCCAAATGGTGGCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATaaa                                                                  >  1:10958547/1‑68 (MQ=255)
   gCTACATGCCAAATGGTGGCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGAGaaa                                                                  >  1:18330414/1‑68 (MQ=255)
     tACATGCCAAATGGTGGCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGc                                                                >  1:1559674/1‑68 (MQ=255)
      aCATGCCAAATGGTGGCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCg                                                               >  1:20940919/1‑68 (MQ=255)
           ccAAATGGTGGCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCg                                                           >  1:39509157/1‑67 (MQ=255)
           ccAAATGGTGGCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCg                                                           >  1:24271584/1‑67 (MQ=255)
           ccAAATGGTGGCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCg                                                           >  1:15178732/1‑67 (MQ=255)
           ccAAATGGTGGCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCg                                                           >  1:28568959/1‑67 (MQ=255)
           ccAAATGGTGGCGCAGTTACCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCg                                                           >  1:27860055/1‑67 (MQ=255)
            cAAATGGTGGCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCgg                                                          <  1:1931634/67‑1 (MQ=255)
             aaaTGGTGGCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCg                                                        >  1:22930889/1‑68 (MQ=255)
             aaaTGGTGGCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCg                                                        >  1:23197023/1‑68 (MQ=255)
             aaaTGGTGGCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCg                                                        >  1:23766016/1‑68 (MQ=255)
             aaaTGGTGGCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCg                                                        >  1:25482497/1‑68 (MQ=255)
              aaTGGTGGCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCga                                                       <  1:36458588/68‑1 (MQ=255)
               aTGGTGGCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCgag                                                      <  1:4389610/68‑1 (MQ=255)
               aTGGTGGCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCgag                                                      >  1:31738034/1‑68 (MQ=255)
                tggtggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCgaga                                                     >  1:9591661/1‑68 (MQ=255)
                tggtggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCgaga                                                     >  1:32195816/1‑68 (MQ=255)
                tggtggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCgaga                                                     >  1:32754415/1‑68 (MQ=255)
                tggtggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCgaga                                                     <  1:3336862/68‑1 (MQ=255)
                  gtggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATg                                                   >  1:21980547/1‑68 (MQ=255)
                   tggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATg                                                   <  1:32405502/67‑1 (MQ=255)
                   tggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATg                                                   <  1:26257038/67‑1 (MQ=255)
                   tggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATg                                                   <  1:19499231/67‑1 (MQ=255)
                   tggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATg                                                   <  1:19916321/67‑1 (MQ=255)
                   tggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATg                                                   <  1:13570760/67‑1 (MQ=255)
                       gcAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAg                                               <  1:22137225/67‑1 (MQ=255)
                       gcAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAg                                               <  1:16880758/67‑1 (MQ=255)
                       gcAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAg                                               <  1:8944853/67‑1 (MQ=255)
                       gcAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGa                                              <  1:38860280/68‑1 (MQ=255)
                            tGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAg                                          <  1:34790226/67‑1 (MQ=255)
                            tGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAg                                          <  1:16103367/67‑1 (MQ=255)
                                aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                     <  1:4183329/68‑1 (MQ=255)
                                aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                     <  1:37320614/68‑1 (MQ=255)
                                aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                     <  1:27449720/68‑1 (MQ=255)
                                aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                     <  1:20678180/68‑1 (MQ=255)
                                aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                     <  1:16243843/68‑1 (MQ=255)
                                 aCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGaa                                    >  1:29864249/1‑68 (MQ=255)
                                 aCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGaa                                    >  1:28111572/1‑68 (MQ=255)
                                 aCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGaa                                    >  1:7428330/1‑68 (MQ=255)
                                 aCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGAATGaa                                    >  1:17823267/1‑68 (MQ=255)
                                   aTTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATg                                  >  1:29517767/1‑68 (MQ=255)
                                   aTTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATg                                  >  1:25019740/1‑68 (MQ=255)
                                     tACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGCGATGAAAGAACAGACTGAATGGc                                <  1:5022056/68‑1 (MQ=255)
                                     tACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGc                                <  1:5197201/68‑1 (MQ=255)
                                     tACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGc                                <  1:72072/68‑1 (MQ=255)
                                        gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                             >  1:15811337/1‑68 (MQ=255)
                                        gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                             >  1:22284550/1‑68 (MQ=255)
                                        gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                             >  1:1151942/1‑68 (MQ=255)
                                        gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                             >  1:13783529/1‑68 (MQ=255)
                                        gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                             >  1:34960628/1‑68 (MQ=255)
                                        gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                             >  1:6672988/1‑68 (MQ=255)
                                        gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                             >  1:35269315/1‑68 (MQ=255)
                                        gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                             >  1:13275816/1‑68 (MQ=255)
                                        gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                             >  1:38204684/1‑68 (MQ=255)
                                        gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                             >  1:29321980/1‑68 (MQ=255)
                                        gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                             >  1:25836654/1‑68 (MQ=255)
                                        gCTGGCTACTTCCGAATCATGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                             >  1:4006943/1‑68 (MQ=255)
                                         cTGGCTACTTCCGAATCCTGGCGTGGTAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                             >  1:26600610/1‑67 (MQ=255)
                                         cTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                             >  1:8725799/1‑67 (MQ=255)
                                         cTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                             >  1:7750946/1‑67 (MQ=255)
                                          tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                           <  1:34409498/68‑1 (MQ=255)
                                          tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                           <  1:39111412/68‑1 (MQ=255)
                                          tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                           <  1:12999671/68‑1 (MQ=255)
                                          tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                           <  1:12108898/68‑1 (MQ=255)
                                          tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                           <  1:5395130/68‑1 (MQ=255)
                                          tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                           <  1:30649473/68‑1 (MQ=255)
                                          tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                           <  1:28853430/68‑1 (MQ=255)
                                              tACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGtt                        <  1:620189/67‑1 (MQ=255)
                                               aCTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                      <  1:18640835/68‑1 (MQ=255)
                                               aCTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                      <  1:9377131/68‑1 (MQ=255)
                                               aCTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                      <  1:37649123/68‑1 (MQ=255)
                                                  tCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                      >  1:35568303/1‑65 (MQ=255)
                                                  tCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                      >  1:31318927/1‑65 (MQ=255)
                                                  tCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCt                   <  1:20127233/68‑1 (MQ=255)
                                                  tCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCt                   <  1:10445079/68‑1 (MQ=255)
                                                    cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                 <  1:7980017/68‑1 (MQ=255)
                                                    cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                 <  1:21576722/68‑1 (MQ=255)
                                                    cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                 <  1:20776007/68‑1 (MQ=255)
                                                       aTCCTTGCGTGATAAAGCGACCGGCGATATGAAAGAACAGACTGAATGGAACCGCGTTGTGCTGTTCg              >  1:30307426/1‑68 (MQ=255)
                                                       aTCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCg              >  1:30912126/1‑68 (MQ=255)
                                                       aTCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCg              >  1:10236313/1‑68 (MQ=255)
                                                       aTCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCg              >  1:12548441/1‑68 (MQ=255)
                                                       aTCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCg              >  1:4801264/1‑68 (MQ=255)
                                                       aTCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCg              >  1:33673883/1‑68 (MQ=255)
                                                        tCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCg              >  1:14228075/1‑67 (MQ=255)
                                                         gcTGGCGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGc            >  1:11946437/2‑68 (MQ=255)
                                                         gcTGGCGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGc            >  1:10078067/2‑68 (MQ=255)
                                                         gcTGGCGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGc            >  1:21267713/2‑68 (MQ=255)
                                                         ccTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGc            >  1:26837144/1‑68 (MQ=255)
                                                              cGTTATAAAGCGACCGGCGAGATGAACCAACATACTGAATGGAAAAGCGTTGTGCTGTTCGGCAAACt       <  1:3701552/68‑1 (MQ=255)
                                                              cGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCaa          >  1:30159347/1‑65 (MQ=255)
                                                               gTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTg      <  1:14679753/68‑1 (MQ=255)
                                                               gTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTg      <  1:18832458/68‑1 (MQ=255)
                                                               gTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTg      <  1:35796866/68‑1 (MQ=255)
                                                                 gATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTgg     <  1:17426523/67‑1 (MQ=255)
                                                                 gATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTgg     <  1:24584936/67‑1 (MQ=255)
                                                                   |                                                                   
TACGCTACATGCCAAATGGTGGCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAG  >  CP009273/4264114‑4264248

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: