Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,264,200 A→C 6.0% 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 6.0% 1778.3 / 74.7 483K50Q (AAA→CAA) ‡ssbsingle‑stranded DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (177/277);  new base C (20/9);  total (197/286)
Fisher's exact test for biased strand distribution p-value = 2.80e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

Only 100 of 487 total aligned reads displayed.

GGCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATC  >  CP009273/4264134‑4264266
                                                                  |                                                                  
ggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGaa                                                                   >  1:5070397/1‑68 (MQ=255)
ggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGaa                                                                   <  1:30817474/68‑1 (MQ=255)
     aGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGaa                                                               <  1:14102237/67‑1 (MQ=255)
     aGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAAc                                                              >  1:30147866/1‑68 (MQ=255)
     aGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAAc                                                              >  1:9054624/1‑68 (MQ=255)
        tGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAg                                                            <  1:8557460/67‑1 (MQ=255)
        tGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGa                                                           <  1:15188599/68‑1 (MQ=255)
         gCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGAc                                                          <  1:17692994/68‑1 (MQ=255)
            aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                       <  1:18772828/68‑1 (MQ=255)
            aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                       <  1:34431413/68‑1 (MQ=255)
            aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                       <  1:7506456/68‑1 (MQ=255)
            aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                       <  1:23733593/68‑1 (MQ=255)
            aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                       <  1:23885756/68‑1 (MQ=255)
             aCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGaa                                                      >  1:24441308/1‑68 (MQ=255)
             aCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGaa                                                      >  1:14744699/1‑68 (MQ=255)
             aCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGaa                                                      >  1:12102787/1‑68 (MQ=255)
             aCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGaa                                                      >  1:21838426/1‑68 (MQ=255)
               aTTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATg                                                    >  1:1833101/1‑68 (MQ=255)
               aTTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATg                                                    >  1:31137758/1‑68 (MQ=255)
                ttACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATg                                                    >  1:22649804/1‑67 (MQ=255)
                 tACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGc                                                  <  1:21558965/68‑1 (MQ=255)
                 tACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGc                                                  <  1:4894833/68‑1 (MQ=255)
                 tACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGc                                                  <  1:20654705/68‑1 (MQ=255)
                 tACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGc                                                  <  1:19802788/68‑1 (MQ=255)
                 tACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGc                                                  <  1:33262204/68‑1 (MQ=255)
                 tACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGc                                                  <  1:17990796/68‑1 (MQ=255)
                  aCGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGc                                                  >  1:29679329/1‑67 (MQ=255)
                    gCTGGCTACTTCCGAATCCTGGCTTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                               >  1:11057662/1‑68 (MQ=255)
                    gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGCCCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                               >  1:11836937/1‑68 (MQ=255)
                    gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                               >  1:3164128/1‑68 (MQ=255)
                    gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                               >  1:22362609/1‑68 (MQ=255)
                    gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                               >  1:10178410/1‑68 (MQ=255)
                    gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                               >  1:3689844/1‑68 (MQ=255)
                    gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                               >  1:16966560/1‑68 (MQ=255)
                    gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                               >  1:13671440/1‑68 (MQ=255)
                    gCTGGCTACTTCCGAATACTGGCGTGAGAAAGCCACCGGCGAGATCAAAGACCAGACAGAATGGCAcc                                               >  1:28162048/1‑68 (MQ=255)
                     cTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                               >  1:32727606/1‑67 (MQ=255)
                     cTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                               >  1:915357/1‑67 (MQ=255)
                      tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGACTGGCACcgc                                             <  1:24036417/68‑1 (MQ=255)
                      tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                             <  1:17261594/68‑1 (MQ=255)
                      tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                             <  1:26761238/68‑1 (MQ=255)
                      tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                             <  1:2599277/68‑1 (MQ=255)
                      tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                             <  1:25610623/68‑1 (MQ=255)
                      tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                             <  1:12499743/68‑1 (MQ=255)
                      tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                             <  1:12934564/68‑1 (MQ=255)
                      tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                             <  1:13162582/68‑1 (MQ=255)
                      tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                             <  1:23216985/68‑1 (MQ=255)
                      tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                             <  1:15727981/68‑1 (MQ=255)
                      tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                             <  1:15847285/68‑1 (MQ=255)
                      tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                             <  1:31418260/68‑1 (MQ=255)
                      tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                             <  1:20047783/68‑1 (MQ=255)
                       ggCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgcg                                            >  1:15391712/1‑68 (MQ=255)
                         cTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGtt                                          >  1:2093269/1‑68 (MQ=255)
                         cTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGtt                                          >  1:27897230/1‑68 (MQ=255)
                           aCTTCCTAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                        <  1:25812892/68‑1 (MQ=255)
                           aCTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                        <  1:7872464/68‑1 (MQ=255)
                           aCTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                        <  1:16212988/68‑1 (MQ=255)
                           aCTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                        <  1:16880382/68‑1 (MQ=255)
                              tCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                        >  1:5568466/1‑65 (MQ=255)
                              tCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                        >  1:24659146/1‑65 (MQ=255)
                              tCCGAATACTGGCGTGTTAATGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCt                                     <  1:34615603/68‑1 (MQ=255)
                                cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                   >  1:16070383/1‑68 (MQ=255)
                                cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                   <  1:10516775/68‑1 (MQ=255)
                                cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                   <  1:11441363/68‑1 (MQ=255)
                                cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                   <  1:28566842/68‑1 (MQ=255)
                                cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                   <  1:10002105/68‑1 (MQ=255)
                                cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                   <  1:33482468/68‑1 (MQ=255)
                                cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                   <  1:21135399/68‑1 (MQ=255)
                                cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                   <  1:31943442/68‑1 (MQ=255)
                                cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                   <  1:33164669/68‑1 (MQ=255)
                                   aTCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCg                                >  1:29037958/1‑68 (MQ=255)
                                   aTCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCg                                >  1:29332193/1‑68 (MQ=255)
                                   aTCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCg                                >  1:15069121/1‑68 (MQ=255)
                                   aTCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCg                                >  1:19531464/1‑68 (MQ=255)
                                     gcTGGCGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGc                              >  1:16441812/2‑68 (MQ=255)
                                     gcTGGCGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGc                              >  1:6736041/2‑68 (MQ=255)
                                     gcTGGCGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGc                              >  1:20365222/2‑68 (MQ=255)
                                     gcTGGCGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGc                              >  1:27495003/2‑68 (MQ=255)
                                          cGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCaa                            >  1:32544425/1‑65 (MQ=255)
                                          cGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCaa                            >  1:33721299/1‑65 (MQ=255)
                                          cGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCaagc                          >  1:30361170/1‑65 (MQ=255)
                                             gATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTgg                       <  1:25105863/67‑1 (MQ=255)
                                               tAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAg                    >  1:6015530/1‑68 (MQ=255)
                                                aaaGCGACCGGCGATATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGa                   <  1:28766220/68‑1 (MQ=255)
                                                  aGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAATCTGGCAGAAg                 >  1:7198177/1‑68 (MQ=255)
                                                  aGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAg                 >  1:18888650/1‑68 (MQ=255)
                                                  aGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAg                 >  1:21326380/1‑68 (MQ=255)
                                                    agACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTg               >  1:4028232/2‑68 (MQ=255)
                                                    agACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTg               >  1:6376794/2‑68 (MQ=255)
                                                     gACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTgg              <  1:29426409/68‑1 (MQ=255)
                                                         ggCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgag          <  1:9733140/68‑1 (MQ=255)
                                                           cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc         <  1:18499032/67‑1 (MQ=255)
                                                           cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc         <  1:5835644/67‑1 (MQ=255)
                                                           cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc         <  1:31225672/67‑1 (MQ=255)
                                                           cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc         <  1:14316076/67‑1 (MQ=255)
                                                           cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc         <  1:10841853/67‑1 (MQ=255)
                                                            gagaTGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagcga       <  1:27089475/68‑1 (MQ=255)
                                                                tGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAata    <  1:4488751/67‑1 (MQ=255)
                                                                tGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAata    >  1:9478479/1‑67 (MQ=255)
                                                                tGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAat     <  1:18590341/66‑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: