Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,264,259 C→T 18.7% S69S (AGC→AGT) ‡ ssb → single‑stranded DNA‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,264,2590CT18.7% 1621.0 / inf 694S69S (AGC→AGT) ‡ssbsingle‑stranded DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (222/339);  new base T (33/96);  total (255/439)
Fisher's exact test for biased strand distribution p-value = 3.25e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 774 total aligned reads displayed.

CGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGA  >  CP009273/4264193‑4264324
                                                                  |                                                                 
cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc                                                                   <  1:19289704/67‑1 (MQ=255)
cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc                                                                   <  1:21528928/67‑1 (MQ=255)
cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc                                                                   <  1:11315037/67‑1 (MQ=255)
cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc                                                                   <  1:23581191/67‑1 (MQ=255)
cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc                                                                   <  1:10000779/67‑1 (MQ=255)
cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc                                                                   <  1:31339380/67‑1 (MQ=255)
cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc                                                                   <  1:30318348/67‑1 (MQ=255)
 gagaTGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagcga                                                                 <  1:10523618/68‑1 (MQ=255)
 gagaTGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagcga                                                                 <  1:22819566/68‑1 (MQ=255)
   gaTGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGaa                                                                <  1:655092/67‑1 (MQ=255)
      gAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATc                                                            <  1:33047327/68‑1 (MQ=255)
      gAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATc                                                            <  1:6711584/68‑1 (MQ=255)
         aGAACAGACTGATTGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc                                                         >  1:23866499/1‑68 (MQ=255)
         aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc                                                         <  1:31024564/68‑1 (MQ=255)
         aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc                                                         <  1:28732222/68‑1 (MQ=255)
         aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc                                                         <  1:27801428/68‑1 (MQ=255)
          gagcagacagaATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGc                                                         <  1:4791593/64‑1 (MQ=255)
           aaCAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGt                                                       >  1:15033156/1‑68 (MQ=255)
           aaCAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGt                                                       >  1:16177542/1‑68 (MQ=255)
            aCAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTa                                                      <  1:2917891/68‑1 (MQ=255)
            aCAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTa                                                      >  1:15804879/1‑68 (MQ=255)
            aCAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTa                                                      >  1:1698068/1‑68 (MQ=255)
             cAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTaa                                                     >  1:27444165/1‑68 (MQ=255)
                aCTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAgg                                                  >  1:19563211/1‑68 (MQ=255)
                 cTGAATTGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGt                                                 <  1:15383295/68‑1 (MQ=255)
                 cTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGt                                                 <  1:36976390/68‑1 (MQ=255)
                   gAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTtc                                               <  1:27052193/68‑1 (MQ=255)
                     aTGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTtct                                              <  1:29908566/67‑1 (MQ=255)
                      tgccACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCa                                            >  1:17321428/4‑68 (MQ=255)
                      tGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCa                                            <  1:14779232/68‑1 (MQ=255)
                      tGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCa                                            >  1:12965381/1‑68 (MQ=255)
                      tGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCa                                            >  1:21875337/1‑68 (MQ=255)
                        gCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                          >  1:35647660/1‑68 (MQ=255)
                        gCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                          >  1:32671806/1‑68 (MQ=255)
                         cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                          <  1:4292655/67‑1 (MQ=255)
                         cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                          <  1:16445592/67‑1 (MQ=255)
                         cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                          <  1:34123563/67‑1 (MQ=255)
                         cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                          <  1:20774591/67‑1 (MQ=255)
                         cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                          <  1:9116763/67‑1 (MQ=255)
                         cACCGCGTGGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                          <  1:21219046/67‑1 (MQ=255)
                          aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                                        >  1:5510550/1‑68 (MQ=255)
                          aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                                        >  1:24403239/1‑68 (MQ=255)
                          aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                                        >  1:35060207/1‑68 (MQ=255)
                          aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                                        >  1:29400522/1‑68 (MQ=255)
                          aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                                        >  1:35974285/1‑68 (MQ=255)
                          aCCGCGTTGTGCTGTTCAGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                                        >  1:1094593/1‑68 (MQ=255)
                          aCCGCGTTGTGCGGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                                        >  1:34785390/1‑68 (MQ=255)
                           tcgcgTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGtct                                       <  1:9558709/67‑3 (MQ=255)
                           tcgcgTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGtct                                       <  1:35485204/67‑3 (MQ=255)
                             gcgTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTtat                                     <  1:12407490/68‑1 (MQ=255)
                             gcgTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCtac                                     <  1:5805826/68‑2 (MQ=255)
                               gTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTtata                                    >  1:22275552/1‑67 (MQ=255)
                                 tgtgCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATc                                  <  1:11842288/67‑1 (MQ=255)
                                    gCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAg                              >  1:19793008/1‑68 (MQ=255)
                                    gCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAg                              >  1:33335856/1‑68 (MQ=255)
                                    gCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAg                              >  1:31977116/1‑68 (MQ=255)
                                    gCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAg                              >  1:24882906/1‑68 (MQ=255)
                                      tGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAgg                             <  1:19919152/67‑1 (MQ=255)
                                         tCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCa                          >  1:36733363/1‑67 (MQ=255)
                                         tCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCa                          >  1:4444574/1‑67 (MQ=255)
                                         tCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCa                          >  1:28169703/1‑67 (MQ=255)
                                          cGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGc                        <  1:26036910/68‑1 (MQ=255)
                                          cGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGc                        <  1:25050527/68‑1 (MQ=255)
                                           ggCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCt                       <  1:23085503/68‑1 (MQ=255)
                                           ggCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCt                       <  1:17825837/68‑1 (MQ=255)
                                           ggCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCt                       <  1:8744019/68‑1 (MQ=255)
                                           ggCAAACTGGCAGAAGTGGCGAGCGAATATCTGAGTAAAGGTTCTCAGGTTTATATCGAAGGtaa                          >  1:13956095/1‑63 (MQ=255)
                                              aaaCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAGGGTCAGCTGCg                    >  1:31651821/1‑68 (MQ=255)
                                              aaaCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGc                     <  1:5959597/67‑1 (MQ=255)
                                              aaaCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCg                    >  1:18474809/1‑68 (MQ=255)
                                                 cTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTa                  >  1:25306564/1‑67 (MQ=255)
                                                  tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc                <  1:26689224/68‑1 (MQ=255)
                                                  tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc                <  1:6891240/68‑1 (MQ=255)
                                                  tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc                <  1:20637896/68‑1 (MQ=255)
                                                  tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc                <  1:22091793/68‑1 (MQ=255)
                                                  tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc                <  1:32456590/68‑1 (MQ=255)
                                                  tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc                <  1:16709117/68‑1 (MQ=255)
                                                   ggCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAccc               <  1:1154618/68‑1 (MQ=255)
                                                      ggAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCTTCGCACCCGt             <  1:7160517/66‑1 (MQ=255)
                                                      ggAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCTTCGCACCCGt             <  1:12697970/66‑1 (MQ=255)
                                                      ggAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCTTCGCACCCGt             <  1:33414822/66‑1 (MQ=255)
                                                      ggAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCTTCGCACCCGTa            <  1:7558272/67‑1 (MQ=255)
                                                      ggAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCTTCGCACCCGTa            <  1:29011038/67‑1 (MQ=255)
                                                      ggAAGTGGCAGGGGAAAGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCTTCGCACCCGt             <  1:25838749/66‑1 (MQ=255)
                                                      aGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTATATCGAAGGTCAGCTGCGTACCCGTa            <  1:3662105/67‑1 (MQ=255)
                                                         aGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTa            <  1:3610987/65‑1 (MQ=255)
                                                         aGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTa            <  1:18203885/65‑1 (MQ=255)
                                                         aGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTa            <  1:22514816/65‑1 (MQ=255)
                                                         aGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTa            <  1:5179227/65‑1 (MQ=255)
                                                         aGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAAt         <  1:26264668/68‑1 (MQ=255)
                                                         aGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCTTCGCACCCGTagc          <  1:33682347/67‑3 (MQ=255)
                                                         aGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCTTCGCACCCGTagc          <  1:3796741/67‑3 (MQ=255)
                                                         aGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCTTCGCACCCGTagc          <  1:20344507/67‑3 (MQ=255)
                                                          gTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATg        >  1:32253179/1‑68 (MQ=255)
                                                          gTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATg        >  1:2848892/1‑68 (MQ=255)
                                                            ggCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATgg       <  1:1213752/67‑1 (MQ=255)
                                                            ggCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATgg       <  1:34439516/67‑1 (MQ=255)
                                                                 gcgaATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGa  <  1:29531795/67‑1 (MQ=255)
                                                                 gcgaATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGa  <  1:13573559/67‑1 (MQ=255)
                                                                 gcgaATATCTGCGTAAAGGTTCTCAGGGTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGa  <  1:8241349/67‑1 (MQ=255)
                                                                  |                                                                 
CGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGA  >  CP009273/4264193‑4264324

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: