Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,264,257 A→G 16.2% S69G (AGC→GGC) ‡ ssb → single‑stranded DNA‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,264,2570AG16.2% 1869.3 / 288.8 687S69G (AGC→GGC) ‡ssbsingle‑stranded DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (222/352);  new base G (44/67);  total (268/419)
Fisher's exact test for biased strand distribution p-value = 9.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

Only 100 of 748 total aligned reads displayed.

CGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCG  >  CP009273/4264190‑4264323
                                                                   |                                                                  
cGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcga                                                                    <  1:34786586/68‑1 (MQ=255)
 ggCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgag                                                                   <  1:872502/68‑1 (MQ=255)
   cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc                                                                  <  1:22157582/67‑1 (MQ=255)
   cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc                                                                  <  1:9328039/67‑1 (MQ=255)
   cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc                                                                  <  1:38850468/67‑1 (MQ=255)
   cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc                                                                  <  1:12018109/67‑1 (MQ=255)
    gagaTGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagcga                                                                <  1:27869464/68‑1 (MQ=255)
        tGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAata                                                             <  1:23222216/67‑1 (MQ=255)
            ataaCAGACTGAATGGCACCGCGTTGTGATGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc                                                        <  1:7372627/66‑1 (MQ=255)
            aGAACAGACTGAATGGCTCCGCGGTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc                                                        <  1:28855250/68‑1 (MQ=255)
            aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc                                                        <  1:4223186/68‑1 (MQ=255)
            aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc                                                        <  1:16625589/68‑1 (MQ=255)
             gagcagacagaATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGc                                                        <  1:39125738/64‑1 (MQ=255)
             gagcagacagaATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGc                                                        <  1:22776787/64‑1 (MQ=255)
             gagcagacagaATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGc                                                        >  1:15695116/4‑67 (MQ=255)
              aaCAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGt                                                      >  1:17190829/1‑68 (MQ=255)
               aCAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTa                                                     >  1:14227539/1‑68 (MQ=255)
               aCAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTa                                                     >  1:35668264/1‑68 (MQ=255)
               aCAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTa                                                     >  1:31798951/1‑68 (MQ=255)
               aCAGACTGAATGCCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTa                                                     >  1:373646/1‑68 (MQ=255)
                   aCTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAgg                                                 >  1:28451600/1‑68 (MQ=255)
                    cTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGt                                                <  1:24477115/68‑1 (MQ=255)
                    cTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGt                                                <  1:37640509/68‑1 (MQ=255)
                    cTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGt                                                <  1:11714385/68‑1 (MQ=255)
                    cTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGt                                                <  1:21087274/68‑1 (MQ=255)
                     agAATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAgggc                                                >  1:1900438/2‑66 (MQ=255)
                        aTGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCaagg                                                  <  1:19389923/62‑3 (MQ=255)
                        aTGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAgggcgc                                              <  1:8128127/66‑4 (MQ=255)
                        aTGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAgggcgc                                              <  1:13654546/66‑4 (MQ=255)
                        aTGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAggg                                                 <  1:38142766/63‑1 (MQ=255)
                         tGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCa                                           >  1:5680956/1‑68 (MQ=255)
                         tGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCa                                           >  1:30779608/1‑68 (MQ=255)
                         tGGCACCGCGTTGTGCAGTGCGGCAAACTGGCAGAAGTGCCGAGCGAATATCAGCGTAAAGGTTATCa                                           >  1:32832840/1‑68 (MQ=255)
                          ggCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAg                                          >  1:679535/1‑68 (MQ=255)
                            cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                         <  1:21426311/67‑1 (MQ=255)
                            cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                         <  1:3110424/67‑1 (MQ=255)
                            cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                         <  1:3155609/67‑1 (MQ=255)
                            cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                         <  1:26854498/67‑1 (MQ=255)
                            cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                         <  1:6941494/67‑1 (MQ=255)
                            cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                         <  1:25190609/67‑1 (MQ=255)
                             aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                                       >  1:1767831/1‑68 (MQ=255)
                             aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                                       >  1:27036909/1‑68 (MQ=255)
                             aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                                       >  1:26005444/1‑68 (MQ=255)
                             aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                                       >  1:21836514/1‑68 (MQ=255)
                                gcgTTGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCtac                                    <  1:2557249/68‑2 (MQ=255)
                                gcgTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTtat                                    <  1:1096826/68‑1 (MQ=255)
                                gcgTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTtat                                    >  1:34285048/1‑68 (MQ=255)
                                gcgTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTtat                                    <  1:955573/68‑1 (MQ=255)
                                gcgTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCtac                                    <  1:24610252/68‑2 (MQ=255)
                                  gTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTtata                                   >  1:36107887/1‑67 (MQ=255)
                                    ggtgCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATc                                 <  1:32172478/66‑1 (MQ=255)
                                       gCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAg                             >  1:33790334/1‑68 (MQ=255)
                                       gCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAg                             >  1:34046629/1‑68 (MQ=255)
                                       gCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAg                             >  1:7932004/1‑68 (MQ=255)
                                       gCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAg                             >  1:36915961/1‑68 (MQ=255)
                                       gCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGaa                              <  1:19650440/67‑1 (MQ=255)
                                       gCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGaa                              <  1:16963139/67‑1 (MQ=255)
                                       gCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAg                             >  1:24794705/1‑68 (MQ=255)
                                         tGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAgg                            >  1:5075885/1‑67 (MQ=255)
                                          gTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTc                          <  1:13110347/68‑1 (MQ=255)
                                          gTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTc                          <  1:12811558/68‑1 (MQ=255)
                                            tCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCa                         >  1:18074587/1‑67 (MQ=255)
                                            tCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCa                         >  1:10536137/1‑67 (MQ=255)
                                            tCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAg                        <  1:20762265/68‑1 (MQ=255)
                                              ggCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCt                      <  1:8824306/68‑1 (MQ=255)
                                                 aaaCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCg                   >  1:12476605/1‑68 (MQ=255)
                                                 aaaCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCg                   >  1:22334461/1‑68 (MQ=255)
                                                 aaaCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCg                   >  1:39558249/1‑68 (MQ=255)
                                                 aaaCTGGCAGAAGTGGCGAGCCAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCg                   >  1:2402302/1‑68 (MQ=255)
                                                  aaCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGt                  <  1:1991268/68‑1 (MQ=255)
                                                  aaCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGt                  <  1:35188991/68‑1 (MQ=255)
                                                     tGGCAGATGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc               <  1:30240660/68‑1 (MQ=255)
                                                     tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCCGGTTTATATCGAAGGTCAGCTGCGTAcc               <  1:29968502/68‑1 (MQ=255)
                                                     tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTGTATCGAAGGTCAGCTGCGTAcc               <  1:707441/68‑1 (MQ=255)
                                                     tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc               <  1:5329226/68‑1 (MQ=255)
                                                     tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc               <  1:4607211/68‑1 (MQ=255)
                                                     tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc               <  1:38622338/68‑1 (MQ=255)
                                                     tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc               <  1:1498558/68‑1 (MQ=255)
                                                     tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc               <  1:15998704/68‑1 (MQ=255)
                                                     tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc               <  1:18661757/68‑1 (MQ=255)
                                                     tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc               <  1:33457567/68‑1 (MQ=255)
                                                     tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc               <  1:27779827/68‑1 (MQ=255)
                                                     tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc               <  1:6034364/68‑1 (MQ=255)
                                                     tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc               <  1:295706/68‑1 (MQ=255)
                                                      ggCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAccc              <  1:23610271/68‑1 (MQ=255)
                                                         ggAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCTTCGCACCCGt            <  1:27699385/66‑1 (MQ=255)
                                                         ggAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCTTCGCACCCGTa           <  1:10114356/67‑1 (MQ=255)
                                                         ggAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCTTCGCACCCGTa           >  1:24194420/2‑68 (MQ=255)
                                                         ggAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCTTCGCACCCGTa           <  1:20451093/67‑1 (MQ=255)
                                                         ggAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCTTCGCACCCGTa           <  1:14363507/67‑1 (MQ=255)
                                                          gAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCTTCGCACCCGTa           >  1:34604492/1‑67 (MQ=255)
                                                            aGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTa           <  1:5479940/65‑1 (MQ=255)
                                                            aGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTa           <  1:26283150/65‑1 (MQ=255)
                                                            aGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTa           <  1:2310961/65‑1 (MQ=255)
                                                            aGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTa           <  1:10370770/65‑1 (MQ=255)
                                                             gTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATg       >  1:15248077/1‑68 (MQ=255)
                                                              tGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATgg      >  1:29301712/1‑68 (MQ=255)
                                                                 cgagcgaATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGAcc   >  1:13717772/1‑68 (MQ=255)
                                                                  gagcgaATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCg  <  1:6557964/68‑1 (MQ=255)
                                                                  gagcgaATATCTGCGTAAAGGTTCTCAGGTTTATAACGAAGGTCAGCTGCGTACCCGTAAATGGACCg  <  1:16285581/68‑1 (MQ=255)
                                                                   |                                                                  
CGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCG  >  CP009273/4264190‑4264323

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: