Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,264,205 A→G 5.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 5.7% 1786.0 / 58.2 507E51E (GAA→GAGssbsingle‑stranded DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (172/304);  new base G (20/9);  total (194/313)
Fisher's exact test for biased strand distribution p-value = 6.22e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.81e-01

Only 100 of 534 total aligned reads displayed.

AGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTA  >  CP009273/4264139‑4264272
                                                                  |                                                                   
aGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAAc                                                                    >  1:3144243/1‑68 (MQ=255)
   tGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAg                                                                  >  1:7790893/1‑67 (MQ=255)
   tGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAg                                                                  <  1:7495242/67‑1 (MQ=255)
   tGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGa                                                                 <  1:5636833/68‑1 (MQ=255)
    gCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGAc                                                                <  1:27439815/68‑1 (MQ=255)
    gCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGAc                                                                <  1:17692994/68‑1 (MQ=255)
      cAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACt                                                               <  1:22350/67‑1 (MQ=255)
      cAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACt                                                               <  1:13847937/67‑1 (MQ=255)
      cAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTg                                                              >  1:18362908/1‑68 (MQ=255)
      cAACATTACGCTGGCTACTTCCGAATCCAGGCGTGATAAAGCGACCGGCGAGATGAAAGCCCAGACt                                                               >  1:18574791/1‑67 (MQ=255)
       aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                             <  1:23733593/68‑1 (MQ=255)
       aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                             <  1:11344631/68‑1 (MQ=255)
       aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                             <  1:2156173/68‑1 (MQ=255)
       aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                             <  1:34431413/68‑1 (MQ=255)
       aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                             <  1:5500547/68‑1 (MQ=255)
       aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                             <  1:20805190/68‑1 (MQ=255)
        aCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGaa                                                            >  1:33003338/1‑68 (MQ=255)
        aCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGaa                                                            >  1:21838426/1‑68 (MQ=255)
        aCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGaa                                                            >  1:27697624/1‑68 (MQ=255)
        aCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGaa                                                            >  1:13311044/1‑68 (MQ=255)
          aTTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATg                                                          >  1:4063573/1‑68 (MQ=255)
          aTTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATg                                                          >  1:15811115/1‑68 (MQ=255)
           ttACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATg                                                          >  1:22649804/1‑67 (MQ=255)
            tACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGc                                                        <  1:33262204/68‑1 (MQ=255)
            tACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGc                                                        <  1:4897855/68‑1 (MQ=255)
            tACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGc                                                        <  1:28064676/68‑1 (MQ=255)
               gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                                     >  1:3164128/1‑68 (MQ=255)
               gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                                     >  1:3265280/1‑68 (MQ=255)
               gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                                     >  1:17978769/1‑68 (MQ=255)
               gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                                     >  1:3689844/1‑68 (MQ=255)
                cTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                                     >  1:2491281/1‑67 (MQ=255)
                cTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                                     >  1:31879918/1‑67 (MQ=255)
                cTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                                     >  1:29272740/1‑67 (MQ=255)
                cTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                                     >  1:15382960/1‑67 (MQ=255)
                 tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                                   <  1:28675548/68‑1 (MQ=255)
                 tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                                   <  1:30213121/68‑1 (MQ=255)
                 tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                                   <  1:19821630/68‑1 (MQ=255)
                 tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                                   <  1:2599277/68‑1 (MQ=255)
                 tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                                   <  1:1883155/68‑1 (MQ=255)
                 tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                                   <  1:11614195/68‑1 (MQ=255)
                    cTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGtt                                                >  1:1916084/1‑68 (MQ=255)
                    cTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGt                                                 <  1:1830407/67‑1 (MQ=255)
                    cTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGt                                                 <  1:24582282/67‑1 (MQ=255)
                     tACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGtt                                                <  1:28422/67‑1 (MQ=255)
                      aCTTCCGAATCCTGGCGTGATAAAGCGACCGGCTAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                              <  1:14716436/68‑1 (MQ=255)
                      aCTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                              <  1:2024318/68‑1 (MQ=255)
                      aCTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                              <  1:16880382/68‑1 (MQ=255)
                         tCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                              >  1:9074065/1‑65 (MQ=255)
                         tCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGcc                                           >  1:6555769/1‑67 (MQ=255)
                         tCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCt                                           <  1:25770994/68‑1 (MQ=255)
                         tCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCt                                           <  1:9205831/68‑1 (MQ=255)
                         tCCGAATACTGGCGTGTTAATGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCt                                           <  1:34615603/68‑1 (MQ=255)
                           cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                         <  1:33482468/68‑1 (MQ=255)
                           cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                         <  1:8168556/68‑1 (MQ=255)
                           cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                         <  1:29747044/68‑1 (MQ=255)
                           cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                         <  1:10848187/68‑1 (MQ=255)
                           cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                         <  1:10516775/68‑1 (MQ=255)
                           cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                         <  1:10002105/68‑1 (MQ=255)
                           cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                         <  1:15157255/68‑1 (MQ=255)
                           cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                         <  1:24365608/68‑1 (MQ=255)
                           cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                         <  1:24026286/68‑1 (MQ=255)
                           cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                         >  1:17130355/1‑68 (MQ=255)
                           cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                         <  1:19643619/68‑1 (MQ=255)
                           cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                         <  1:21272150/68‑1 (MQ=255)
                           cGAATCCTGGCGTGATAAACCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                         <  1:14262217/68‑1 (MQ=255)
                                ccTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGc                                    >  1:7053174/1‑68 (MQ=255)
                                ccTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGc                                    >  1:29679317/1‑68 (MQ=255)
                                     cGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCaa                                  >  1:33721299/1‑65 (MQ=255)
                                     cGTGACAAACAGACGGGGGAGATGCGGGAGCAGACATAATGGCATCGCGTGGTGCTGTTCGGCaagc                                <  1:16703126/67‑3 (MQ=255)
                                     cGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCaagc                                <  1:20519404/67‑3 (MQ=255)
                                     cGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCaagc                                <  1:12058533/67‑3 (MQ=255)
                                      gTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTg                              <  1:13044213/68‑1 (MQ=255)
                                      gTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCAAGCtc                              <  1:32480274/68‑2 (MQ=255)
                                        gATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTgg                             <  1:10245283/67‑1 (MQ=255)
                                        gATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTgg                             <  1:8565948/67‑1 (MQ=255)
                                          tAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCa                           <  1:30468968/67‑1 (MQ=255)
                                          tAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCa                           <  1:20963901/67‑1 (MQ=255)
                                             acagACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAg                       >  1:29403090/4‑68 (MQ=255)
                                             aGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAg                       <  1:31096775/68‑1 (MQ=255)
                                                      cGAGTTGAAAGAACAGACTGAATGGCACCTCGTTGTGGTGTGCGGCATAGTGGCAGAAGTGgcgagc               <  1:25410492/67‑1 (MQ=255)
                                                      cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc               <  1:31225672/67‑1 (MQ=255)
                                                      cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc               <  1:15565318/67‑1 (MQ=255)
                                                      cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc               <  1:6334586/67‑1 (MQ=255)
                                                      cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc               <  1:59718/67‑1 (MQ=255)
                                                       gagaTGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCACACTGGCAGAAGTGgcgagcga             <  1:21348995/68‑1 (MQ=255)
                                                       gagaTGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagcga             <  1:9667264/68‑1 (MQ=255)
                                                       gagaTGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagcga             <  1:27089475/68‑1 (MQ=255)
                                                           tGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAata          <  1:4488751/67‑1 (MQ=255)
                                                               aGAACAGACTGAATGTCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc     <  1:12611698/68‑1 (MQ=255)
                                                               aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc     <  1:23232351/68‑1 (MQ=255)
                                                               aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc     <  1:17341922/68‑1 (MQ=255)
                                                               aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc     <  1:16376104/68‑1 (MQ=255)
                                                               aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc     <  1:1618193/68‑1 (MQ=255)
                                                               aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc     <  1:12445117/68‑1 (MQ=255)
                                                               aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc     <  1:9434727/68‑1 (MQ=255)
                                                               aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGAGAATATCTGc     <  1:29015214/68‑1 (MQ=255)
                                                                 aaCAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAat           >  1:15027333/1‑60 (MQ=255)
                                                                  aCAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTa  >  1:18692156/1‑68 (MQ=255)
                                                                  aCAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTa  <  1:2666744/68‑1 (MQ=255)
                                                                  aCAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTa  >  1:15991380/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: