Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,264,259 C→T 16.3% 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,2590CT16.3% 1764.3 / 279.3 693S69S (AGC→AGT) ‡ssbsingle‑stranded DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (203/375);  new base T (34/79);  total (238/455)
Fisher's exact test for biased strand distribution p-value = 3.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 764 total aligned reads displayed.

CGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGA  >  CP009273/4264193‑4264324
                                                                  |                                                                 
cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc                                                                   <  1:17164642/67‑1 (MQ=255)
cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc                                                                   <  1:10523864/67‑1 (MQ=255)
cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc                                                                   <  1:31225672/67‑1 (MQ=255)
cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc                                                                   <  1:32583949/67‑1 (MQ=255)
cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc                                                                   <  1:15472078/67‑1 (MQ=255)
cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc                                                                   <  1:34777169/67‑1 (MQ=255)
cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc                                                                   <  1:1195031/67‑1 (MQ=255)
cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc                                                                   <  1:6334586/67‑1 (MQ=255)
cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc                                                                   <  1:20563328/67‑1 (MQ=255)
 gagaTGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagcga                                                                 <  1:5771008/68‑1 (MQ=255)
         aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc                                                         <  1:31793489/68‑1 (MQ=255)
         aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc                                                         <  1:491760/68‑1 (MQ=255)
          gagcagacagaATGGTATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGc                                                         <  1:4258229/64‑1 (MQ=255)
          gagcagacagaATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGc                                                         <  1:18960734/64‑1 (MQ=255)
          gagcagacagaATGGCATCGCGGGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGc                                                         <  1:30114492/64‑1 (MQ=255)
          gAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCg                                                        >  1:21343676/1‑68 (MQ=255)
           agcagacagaATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCgc                                                       >  1:25303926/3‑67 (MQ=255)
           aaCAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGt                                                       >  1:2108529/1‑68 (MQ=255)
            gcagacagaATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCgca                                                      >  1:27362984/2‑66 (MQ=255)
            aCAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCg                                                        <  1:7050252/66‑1 (MQ=255)
            aCAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTa                                                      >  1:32165101/1‑68 (MQ=255)
            aCAGACTGAATGGCACCGCGTTGTGCTGTTAGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTa                                                      >  1:29737470/1‑68 (MQ=255)
               gACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAg                                                   >  1:17540030/1‑68 (MQ=255)
                aCTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAgg                                                  >  1:26487677/1‑68 (MQ=255)
                 cTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGt                                                 >  1:9465019/1‑68 (MQ=255)
                 cTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGt                                                 <  1:17828243/68‑1 (MQ=255)
                 cTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGt                                                 <  1:24478295/68‑1 (MQ=255)
                 cTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGt                                                 <  1:16056376/68‑1 (MQ=255)
                 cTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGt                                                 >  1:24060596/1‑68 (MQ=255)
                   gAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTtc                                               <  1:8585503/68‑1 (MQ=255)
                     aTGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAgggcgc                                               <  1:16447086/66‑4 (MQ=255)
                       ggCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAg                                           >  1:8063841/1‑68 (MQ=255)
                        gCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                          >  1:15607531/1‑68 (MQ=255)
                         cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGTATATCTGCGTAAAGGTTCTCAgg                                          <  1:4620472/67‑1 (MQ=255)
                         cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                          <  1:16340694/67‑1 (MQ=255)
                         cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                          <  1:9274725/67‑1 (MQ=255)
                         cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                          <  1:18403157/67‑1 (MQ=255)
                          aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                                        >  1:20832512/1‑68 (MQ=255)
                          aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                                        >  1:3338580/1‑68 (MQ=255)
                          aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                                        >  1:13648938/1‑68 (MQ=255)
                          aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                                        >  1:11328819/1‑68 (MQ=255)
                          aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                                        >  1:874307/1‑68 (MQ=255)
                          aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGt                                         <  1:10902783/67‑1 (MQ=255)
                          aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGt                                         <  1:14025103/67‑1 (MQ=255)
                          aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCCGCGTAAAGGTTCTCAGGtt                                        >  1:15174577/1‑68 (MQ=255)
                               gTTGTGCTGTTCGGCACACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTatat                                   <  1:14307417/68‑1 (MQ=255)
                               gTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTtata                                    >  1:22295033/1‑67 (MQ=255)
                               gTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTtata                                    <  1:21601532/67‑1 (MQ=255)
                               gTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTtata                                    >  1:26801243/1‑67 (MQ=255)
                                 ggtgCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATc                                  <  1:2284043/66‑1 (MQ=255)
                                    gCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAg                              >  1:3829422/1‑68 (MQ=255)
                                    gCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAg                              >  1:33594999/1‑68 (MQ=255)
                                    gCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGaa                               <  1:18654898/67‑1 (MQ=255)
                                    gCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAg                              >  1:6152050/1‑68 (MQ=255)
                                    gCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAg                              >  1:34145957/1‑68 (MQ=255)
                                      tGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAgg                             <  1:13297127/67‑1 (MQ=255)
                                      tGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAgg                             <  1:1481122/67‑1 (MQ=255)
                                         tCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCa                          >  1:31484894/1‑67 (MQ=255)
                                         tCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCa                          >  1:10156750/1‑67 (MQ=255)
                                         tCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCa                          <  1:25927457/67‑1 (MQ=255)
                                         tCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGAGTAAAGGTTCTCAGGTTTATATCGAAGGTCa                          >  1:27138207/1‑67 (MQ=255)
                                          cGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGc                        <  1:30414308/68‑1 (MQ=255)
                                           ggCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGc                        <  1:1448402/67‑1 (MQ=255)
                                           ggCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCt                       <  1:28782158/68‑1 (MQ=255)
                                           ggCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCt                       <  1:17303736/68‑1 (MQ=255)
                                           ggCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCt                       <  1:11655264/68‑1 (MQ=255)
                                           ggCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCt                       <  1:11543432/68‑1 (MQ=255)
                                             cataCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTATCAGGTTTATATCGAAGGTCAGCTg                      <  1:21990407/64‑1 (MQ=255)
                                              aaaCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCg                    >  1:841514/1‑68 (MQ=255)
                                              aaaCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCg                    >  1:32851493/1‑68 (MQ=255)
                                              aaaCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCg                    >  1:5973378/1‑68 (MQ=255)
                                                  tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc                <  1:24954363/68‑1 (MQ=255)
                                                  tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc                <  1:5382521/68‑1 (MQ=255)
                                                  tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc                <  1:20114220/68‑1 (MQ=255)
                                                  tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc                <  1:2030731/68‑1 (MQ=255)
                                                  tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc                <  1:12672491/68‑1 (MQ=255)
                                                  tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc                <  1:23585551/68‑1 (MQ=255)
                                                  tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc                <  1:27648305/68‑1 (MQ=255)
                                                  tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc                <  1:2439412/68‑1 (MQ=255)
                                                  tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc                <  1:2463271/68‑1 (MQ=255)
                                                    gCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGt                                         <  1:29450887/41‑1 (MQ=255)
                                                      ggAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCTTCGCACCCGt             <  1:16779013/66‑1 (MQ=255)
                                                      ggAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCTTCGCACCCGTa            <  1:12504382/67‑1 (MQ=255)
                                                      ggAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCTTCGCACCCGTa            <  1:7788409/67‑1 (MQ=255)
                                                      aGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTa            <  1:12300311/68‑1 (MQ=255)
                                                         aGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTa            <  1:25676468/65‑1 (MQ=255)
                                                         aGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTa            <  1:34420563/65‑1 (MQ=255)
                                                         aGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTa            <  1:27935702/65‑1 (MQ=255)
                                                         aGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTa            <  1:29154238/65‑1 (MQ=255)
                                                         aGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTa            <  1:9686011/65‑1 (MQ=255)
                                                         aGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTa            <  1:13118307/65‑1 (MQ=255)
                                                         aGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTa            <  1:7485621/65‑1 (MQ=255)
                                                         aGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTa            <  1:19270560/65‑1 (MQ=255)
                                                         aGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCTTCGCACCCGTagct         <  1:30641158/68‑4 (MQ=255)
                                                            ggCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATgg       <  1:23385408/67‑1 (MQ=255)
                                                                 gcgaATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGa  >  1:19592581/1‑67 (MQ=255)
                                                                 gcgaATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGa  <  1:8912447/67‑1 (MQ=255)
                                                                 gcgaATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGa  <  1:25505974/67‑1 (MQ=255)
                                                                 gcgaATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGa  <  1:30992830/67‑1 (MQ=255)
                                                                 gcgaATATCTGCGTAAAGGCTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGa  <  1:10006388/67‑1 (MQ=255)
                                                                  |                                                                 
CGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGA  >  CP009273/4264193‑4264324

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