Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,264,257 A→G 18.9% 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,2570AG18.9% 1731.8 / inf 694S69G (AGC→GGC) ‡ssbsingle‑stranded DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (221/341);  new base G (33/98);  total (255/439)
Fisher's exact test for biased strand distribution p-value = 2.47e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

Only 100 of 761 total aligned reads displayed.

CGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGA  >  CP009273/4264190‑4264324
                                                                   |                                                                   
cGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcga                                                                     <  1:29761817/68‑1 (MQ=255)
   cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc                                                                   <  1:33517494/67‑1 (MQ=255)
   cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc                                                                   <  1:10000779/67‑1 (MQ=255)
   cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc                                                                   <  1:20734894/67‑1 (MQ=255)
    gagaTGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagcga                                                                 <  1:10523618/68‑1 (MQ=255)
            aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc                                                         >  1:26198178/1‑68 (MQ=255)
            aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc                                                         >  1:5835674/1‑68 (MQ=255)
             gagcagacagaATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGc                                                         <  1:31111888/64‑1 (MQ=255)
             gagcagacagaATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGc                                                         >  1:25656910/4‑67 (MQ=255)
             gAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCg                                                        >  1:22785936/1‑68 (MQ=255)
              aaCAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGt                                                       >  1:15033156/1‑68 (MQ=255)
               aCAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTa                                                      >  1:20576663/1‑68 (MQ=255)
               aCAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTa                                                      >  1:6807260/1‑68 (MQ=255)
               aCAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTa                                                      >  1:13574585/1‑68 (MQ=255)
               aCAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTa                                                      >  1:1698068/1‑68 (MQ=255)
               aCAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTa                                                      >  1:1619048/1‑68 (MQ=255)
                   aCTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAgg                                                  >  1:8026473/1‑68 (MQ=255)
                    cTGAATTGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGt                                                 <  1:15383295/68‑1 (MQ=255)
                    cTGAATGGCACGGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGt                                                 <  1:35120426/68‑1 (MQ=255)
                    cTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTaaa                                                    <  1:33994133/65‑1 (MQ=255)
                    cTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGt                                                 >  1:23473490/1‑68 (MQ=255)
                    cTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGt                                                 <  1:31468493/68‑1 (MQ=255)
                    cTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGt                                                 <  1:21720326/68‑1 (MQ=255)
                    cTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGt                                                 >  1:6907615/1‑68 (MQ=255)
                    cTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGt                                                 <  1:20125744/68‑1 (MQ=255)
                    cTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGt                                                 <  1:18219301/68‑1 (MQ=255)
                    cTGAATGGCACCGCGGTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGt                                                 <  1:34608977/68‑1 (MQ=255)
                         tGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCa                                            >  1:36033114/1‑68 (MQ=255)
                         tGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCa                                            >  1:12965381/1‑68 (MQ=255)
                         tGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCa                                            >  1:36354408/1‑68 (MQ=255)
                            cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                          <  1:19465131/67‑1 (MQ=255)
                            cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                          <  1:30692799/67‑1 (MQ=255)
                            cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                          <  1:9116763/67‑1 (MQ=255)
                            cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                          <  1:6028060/67‑1 (MQ=255)
                            cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                          <  1:17462970/67‑1 (MQ=255)
                            cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                          <  1:4570089/67‑1 (MQ=255)
                            cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                          <  1:37062554/67‑1 (MQ=255)
                            cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                          <  1:25225391/67‑1 (MQ=255)
                             atcgcgTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGt                                         <  1:17866889/65‑1 (MQ=255)
                             atcgcgTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGt                                         <  1:32120458/65‑1 (MQ=255)
                             aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                                        >  1:19827896/1‑68 (MQ=255)
                             aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                                        >  1:32989097/1‑68 (MQ=255)
                             aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                                        >  1:23676761/1‑68 (MQ=255)
                             aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                                        >  1:36684820/1‑68 (MQ=255)
                             aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                                        >  1:2297499/1‑68 (MQ=255)
                             aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                                        >  1:22251127/1‑68 (MQ=255)
                             aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                                        >  1:34895677/1‑68 (MQ=255)
                             aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                                        >  1:18766355/1‑68 (MQ=255)
                             aCCGCGTTGTGCTGTTCAGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                                        >  1:1094593/1‑68 (MQ=255)
                              tcgcgTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGtct                                       <  1:9558709/67‑3 (MQ=255)
                               cgcgTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGtct                                       >  1:25995617/1‑65 (MQ=255)
                                gcgTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTtat                                     <  1:12407490/68‑1 (MQ=255)
                                gcgTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTtat                                     <  1:5066733/68‑1 (MQ=255)
                                  gTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTtata                                    <  1:29077367/67‑1 (MQ=255)
                                    tgtgCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATc                                  <  1:11842288/67‑1 (MQ=255)
                                       gCTGTTCGGCAATCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATAGAAg                              >  1:8749594/1‑68 (MQ=255)
                                       gCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAg                              >  1:21459388/1‑68 (MQ=255)
                                       gCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAg                              >  1:34315332/1‑68 (MQ=255)
                                       gCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAg                              >  1:29309682/1‑68 (MQ=255)
                                         tGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAgg                             <  1:19717475/67‑1 (MQ=255)
                                          gTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATGGAAGGTc                           <  1:26840308/68‑1 (MQ=255)
                                            tCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCa                          >  1:30132108/1‑67 (MQ=255)
                                            tCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCa                          >  1:22010746/1‑67 (MQ=255)
                                            tCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCa                          <  1:3990220/67‑1 (MQ=255)
                                            tCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCa                          >  1:24245326/1‑67 (MQ=255)
                                             cGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGc                        <  1:15834300/68‑1 (MQ=255)
                                             cGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGc                        <  1:27921510/68‑1 (MQ=255)
                                              ggCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCt                       <  1:14826246/68‑1 (MQ=255)
                                              ggCAAACTGGCAGAAGTGGCGAGCGAATATCTGAGTAAAGGTTCTCAGGTTTATATCGAAGGtaa                          >  1:13956095/1‑63 (MQ=255)
                                                 aaaCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCg                    >  1:7165209/1‑68 (MQ=255)
                                                     tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc                <  1:32390668/68‑1 (MQ=255)
                                                     tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc                <  1:22499730/68‑1 (MQ=255)
                                                     tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc                <  1:24991497/68‑1 (MQ=255)
                                                     tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc                <  1:2650082/68‑1 (MQ=255)
                                                     tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc                <  1:27248903/68‑1 (MQ=255)
                                                     tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc                <  1:16709117/68‑1 (MQ=255)
                                                     tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc                <  1:16502386/68‑1 (MQ=255)
                                                      ggCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAccc               <  1:29453100/68‑1 (MQ=255)
                                                      ggCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAccc               <  1:1154618/68‑1 (MQ=255)
                                                         ggAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCTTCGCACCCGt             <  1:12697970/66‑1 (MQ=255)
                                                         ggAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCTTCGCACCCGt             <  1:28319323/66‑1 (MQ=255)
                                                         ggAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCTTCGCACCCGTa            <  1:7558272/67‑1 (MQ=255)
                                                         ggAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCTTCGCACCCGTa            <  1:11297001/67‑1 (MQ=255)
                                                         ggAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCTTCGCACCCGTa            <  1:35594523/67‑1 (MQ=255)
                                                         ggAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCTTCGCACCCGTa            <  1:28640924/67‑1 (MQ=255)
                                                         ggAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCTTCGCACCCGTa            <  1:32596016/67‑1 (MQ=255)
                                                          gAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCTTCGCACCCGTa            >  1:35917869/1‑67 (MQ=255)
                                                            aGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTa            <  1:27630309/65‑1 (MQ=255)
                                                            aGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTa            <  1:5580913/65‑1 (MQ=255)
                                                            aGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTa            <  1:31760154/65‑1 (MQ=255)
                                                            aGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTa            <  1:33380753/65‑1 (MQ=255)
                                                            aGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAAt         <  1:6571567/68‑1 (MQ=255)
                                                            aGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAAt         >  1:21075091/1‑68 (MQ=255)
                                                            aGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAAt         <  1:5250017/68‑1 (MQ=255)
                                                            aGTGCCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTa            <  1:4367469/65‑1 (MQ=255)
                                                             gTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATg        >  1:28794523/1‑68 (MQ=255)
                                                               ggCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATgg       <  1:1213752/67‑1 (MQ=255)
                                                               ggCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATgg       <  1:367890/67‑1 (MQ=255)
                                                                  gagcgaATATCTGCGTAAATGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCg   <  1:24728066/68‑1 (MQ=255)
                                                                  gagcgaATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCg   <  1:33184756/68‑1 (MQ=255)
                                                                   |                                                                   
CGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGA  >  CP009273/4264190‑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: