Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,264,256 G→A 16.3% A68A (GCG→GCA ssb → single‑stranded DNA‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,264,2560GA16.3% 1773.1 / inf 693A68A (GCG→GCAssbsingle‑stranded DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (204/375);  new base A (34/79);  total (238/455)
Fisher's exact test for biased strand distribution p-value = 3.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.90e-01

Only 100 of 752 total aligned reads displayed.

CCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCG  >  CP009273/4264189‑4264323
                                                                   |                                                                   
  ggCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgag                                                                   <  1:9733140/68‑1 (MQ=255)
    cGAGATTAAAGAACAGCCTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc                                                                  <  1:34459172/67‑1 (MQ=255)
    cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc                                                                  <  1:1195031/67‑1 (MQ=255)
    cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc                                                                  <  1:4668274/67‑1 (MQ=255)
    cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc                                                                  <  1:17164642/67‑1 (MQ=255)
    cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc                                                                  <  1:15622644/67‑1 (MQ=255)
    cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc                                                                  <  1:15472078/67‑1 (MQ=255)
     gagaTGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagcga                                                                <  1:31232049/68‑1 (MQ=255)
         tGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAata                                                             >  1:9478479/1‑67 (MQ=255)
         tGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAat                                                              <  1:18590341/66‑1 (MQ=255)
          gAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATc                                                           <  1:32206589/68‑1 (MQ=255)
             aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc                                                        <  1:24373280/68‑1 (MQ=255)
             aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc                                                        <  1:23232351/68‑1 (MQ=255)
             aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc                                                        <  1:25511904/68‑1 (MQ=255)
             aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc                                                        <  1:2724244/68‑1 (MQ=255)
              gagcagacagaATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGc                                                        <  1:18960734/64‑1 (MQ=255)
              gagcagacagaATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGc                                                        <  1:20477032/64‑1 (MQ=255)
               agcagacagaATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCgc                                                      >  1:25303926/3‑67 (MQ=255)
                aCAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTa                                                     >  1:19266863/1‑68 (MQ=255)
                aCAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTa                                                     >  1:7795029/1‑68 (MQ=255)
                aCAGACTGAATGGCACCGCGTTGTGCTGTTAGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTa                                                     >  1:29737470/1‑68 (MQ=255)
                 cAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCCGCGTaa                                                    >  1:19546560/1‑68 (MQ=255)
                   gACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAg                                                  >  1:17540030/1‑68 (MQ=255)
                     cTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGt                                                <  1:16056376/68‑1 (MQ=255)
                     cTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGt                                                <  1:24478295/68‑1 (MQ=255)
                     cTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGt                                                <  1:10915913/68‑1 (MQ=255)
                     cTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGt                                                <  1:17828243/68‑1 (MQ=255)
                     cTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGt                                                >  1:11356461/1‑68 (MQ=255)
                     cTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGt                                                >  1:14887391/1‑68 (MQ=255)
                         aTGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAgggcgc                                              <  1:16447086/66‑4 (MQ=255)
                          tGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCa                                           >  1:22633436/1‑68 (MQ=255)
                          tGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCa                                           <  1:28872662/68‑1 (MQ=255)
                           ggCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAg                                          >  1:8063841/1‑68 (MQ=255)
                            gCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                         >  1:2915751/1‑68 (MQ=255)
                             cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                         <  1:7084722/67‑1 (MQ=255)
                             cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                         <  1:18403157/67‑1 (MQ=255)
                             cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                         <  1:6417025/67‑1 (MQ=255)
                             cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                         <  1:25742795/67‑1 (MQ=255)
                             cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                         <  1:31031688/67‑1 (MQ=255)
                             cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                         <  1:21576115/67‑1 (MQ=255)
                             cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                         <  1:30360569/67‑1 (MQ=255)
                             cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                         <  1:542871/67‑1 (MQ=255)
                             cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                         <  1:16340694/67‑1 (MQ=255)
                             cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                         <  1:13625426/67‑1 (MQ=255)
                              aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                                       >  1:31862646/1‑68 (MQ=255)
                              aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                                       >  1:32972912/1‑68 (MQ=255)
                              aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                                       >  1:34242453/1‑68 (MQ=255)
                              aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                                       >  1:874307/1‑68 (MQ=255)
                              aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                                       >  1:20948181/1‑68 (MQ=255)
                              aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGt                                        <  1:14025103/67‑1 (MQ=255)
                                   gTTGTGCTGTTCGGCAAACTGGCAGAAGTAGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTtata                                   <  1:402808/67‑1 (MQ=255)
                                        gCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAg                             >  1:22233537/1‑68 (MQ=255)
                                        gCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAg                             >  1:23528778/1‑68 (MQ=255)
                                        gCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAg                             >  1:10534157/1‑68 (MQ=255)
                                        gCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAg                             >  1:32683342/1‑68 (MQ=255)
                                          tGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAgg                            <  1:15186577/67‑1 (MQ=255)
                                          tGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAgg                            <  1:13297127/67‑1 (MQ=255)
                                          tGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAgg                            <  1:5100659/67‑1 (MQ=255)
                                             tCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCa                         >  1:27792398/1‑67 (MQ=255)
                                             tCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCa                         >  1:21873797/1‑67 (MQ=255)
                                             tCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCa                         >  1:26620560/1‑67 (MQ=255)
                                             tCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCa                         >  1:13101747/1‑67 (MQ=255)
                                              cGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGc                       <  1:2625867/68‑1 (MQ=255)
                                              cGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGc                       <  1:21340836/68‑1 (MQ=255)
                                              cGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGc                       <  1:4368102/68‑1 (MQ=255)
                                               ggCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGc                       <  1:1448402/67‑1 (MQ=255)
                                               ggCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCt                      <  1:17303736/68‑1 (MQ=255)
                                               ggCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCt                      <  1:11543432/68‑1 (MQ=255)
                                                  aaaCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGc                    <  1:6223418/67‑1 (MQ=255)
                                                  aaaCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCg                   >  1:3541100/1‑68 (MQ=255)
                                                  aaaCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCg                   >  1:8945767/1‑68 (MQ=255)
                                                  aaaCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCg                   >  1:841514/1‑68 (MQ=255)
                                                   aaCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGt                  <  1:26942923/68‑1 (MQ=255)
                                                      tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc               <  1:7537118/68‑1 (MQ=255)
                                                      tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc               <  1:5801645/68‑1 (MQ=255)
                                                      tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc               <  1:24954363/68‑1 (MQ=255)
                                                      tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc               <  1:3350850/68‑1 (MQ=255)
                                                      tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc               <  1:10241556/68‑1 (MQ=255)
                                                      tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc               <  1:11683141/68‑1 (MQ=255)
                                                      tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc               <  1:12672491/68‑1 (MQ=255)
                                                      tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc               <  1:20114220/68‑1 (MQ=255)
                                                      tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc               <  1:2030731/68‑1 (MQ=255)
                                                      tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc               <  1:2463271/68‑1 (MQ=255)
                                                      tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc               <  1:27984203/68‑1 (MQ=255)
                                                        gCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGt                                        <  1:29450887/41‑1 (MQ=255)
                                                          ggAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCTTCGCACCCGt            <  1:16779013/66‑1 (MQ=255)
                                                          ggAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCTTCGCACCCGt            <  1:33780386/66‑1 (MQ=255)
                                                          ggAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCTTCGCACCCGTa           <  1:30058415/67‑1 (MQ=255)
                                                          ggAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCTTCGCACCCGTa           <  1:12504382/67‑1 (MQ=255)
                                                          aGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTa           <  1:12300311/68‑1 (MQ=255)
                                                             aGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTa           <  1:2392034/65‑1 (MQ=255)
                                                             aGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTa           <  1:5984903/65‑1 (MQ=255)
                                                             aGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTa           <  1:14300917/65‑1 (MQ=255)
                                                             aGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTa           <  1:8536375/65‑1 (MQ=255)
                                                             aGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTa           <  1:10014983/65‑1 (MQ=255)
                                                             aGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCTTCGCACCCGTagct        <  1:30641158/68‑4 (MQ=255)
                                                             aGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCTTCGCACCCGTagc         <  1:31550759/67‑3 (MQ=255)
                                                              gTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATg       >  1:9280186/1‑68 (MQ=255)
                                                                ggCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATgg      >  1:2744006/1‑67 (MQ=255)
                                                                  cgagcgaATATCTGCGTAAAGGTTCTCAGGTTTATATAGAAGGTCAGCTGCGTGCCCGTAAATGGAcc   >  1:20813124/1‑68 (MQ=255)
                                                                   |                                                                   
CCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCG  >  CP009273/4264189‑4264323

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: