Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,264,271 T→C 8.0% R73R (CGT→CGC ssb → single‑stranded DNA‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,264,2710TC 8.0% 2083.4 / 108.6 623R73R (CGT→CGCssbsingle‑stranded DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (206/365);  new base C (14/36);  total (222/401)
Fisher's exact test for biased strand distribution p-value = 2.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

Only 100 of 630 total aligned reads displayed.

AACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGG  >  CP009273/4264204‑4264338
                                                                   |                                                                   
aaCAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGt                                                                     >  1:8886086/1‑68 (MQ=255)
 aCAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTa                                                                    >  1:18859394/1‑68 (MQ=255)
 aCAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTa                                                                    >  1:12986393/1‑68 (MQ=255)
  cAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTaa                                                                   >  1:14239267/1‑68 (MQ=255)
      cTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGt                                                               >  1:25204030/1‑68 (MQ=255)
      cTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGt                                                               >  1:13785942/1‑68 (MQ=255)
       agAATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTTGCAGGTGAATGTCTGCGCAAgggc                                                               >  1:29985587/2‑66 (MQ=255)
       agAATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAgggc                                                               >  1:23901951/2‑66 (MQ=255)
           tGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCa                                                          >  1:30726652/1‑68 (MQ=255)
           tGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCa                                                          >  1:31467162/1‑68 (MQ=255)
           tGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCa                                                          >  1:19563796/1‑68 (MQ=255)
           tGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCa                                                          >  1:20749756/1‑68 (MQ=255)
           tGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCa                                                          >  1:7905510/1‑68 (MQ=255)
           tGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGGAAAGGTTCTCa                                                          >  1:28466679/1‑68 (MQ=255)
              cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                                        <  1:21665714/67‑1 (MQ=255)
              cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                                        <  1:8227062/67‑1 (MQ=255)
              cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                                        <  1:1484510/67‑1 (MQ=255)
              cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                                        <  1:6814192/67‑1 (MQ=255)
              cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                                        <  1:15595315/67‑1 (MQ=255)
              cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                                        <  1:27405863/67‑1 (MQ=255)
              cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                                        <  1:16753144/67‑1 (MQ=255)
              cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                                        <  1:26155253/67‑1 (MQ=255)
              cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                                        <  1:11963393/67‑1 (MQ=255)
              cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                                        <  1:24652625/67‑1 (MQ=255)
               aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                                                      >  1:15135169/1‑68 (MQ=255)
               aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                                                      >  1:28285190/1‑68 (MQ=255)
               aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                                                      >  1:22797061/1‑68 (MQ=255)
               aCCGCGTTGTGATGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                                                      >  1:4172347/1‑68 (MQ=255)
                  gcgTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTtat                                                   <  1:7674510/68‑1 (MQ=255)
                  gcgTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTtat                                                   <  1:23437762/68‑1 (MQ=255)
                    gttgtgttgtTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTtata                                                  >  1:594416/1‑67 (MQ=255)
                    gTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTtata                                                  >  1:28971180/1‑67 (MQ=255)
                    gTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTtata                                                  >  1:2147523/1‑67 (MQ=255)
                      tgtgCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATc                                                <  1:8648463/67‑1 (MQ=255)
                         gCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAg                                            >  1:30971078/1‑68 (MQ=255)
                         gCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAg                                            >  1:17530948/1‑68 (MQ=255)
                         gCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAg                                            >  1:25611183/1‑68 (MQ=255)
                         gCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAg                                            >  1:12210184/1‑68 (MQ=255)
                         gCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAg                                            >  1:17212995/1‑68 (MQ=255)
                         gCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATACCTGCGTAAAGGTTCTCAGGTTTATATCGAAg                                            >  1:6147347/1‑68 (MQ=255)
                           tGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAgg                                           <  1:14535916/67‑1 (MQ=255)
                           tGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAgg                                           <  1:10657598/67‑1 (MQ=255)
                              tCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCa                                        >  1:21211491/1‑67 (MQ=255)
                              tCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCa                                        >  1:22299497/1‑67 (MQ=255)
                               cGGCAATCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGc                                      <  1:24375733/68‑1 (MQ=255)
                                ggCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCt                                     <  1:9128921/68‑1 (MQ=255)
                                ggCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCt                                     <  1:25960866/68‑1 (MQ=255)
                                ggCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCt                                     <  1:2137404/68‑1 (MQ=255)
                                 gCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCt                                     >  1:10956333/1‑67 (MQ=255)
                                   aaaCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGc                                   <  1:1019996/67‑1 (MQ=255)
                                   aaaCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGc                                   <  1:15404581/67‑1 (MQ=255)
                                   aaaCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCg                                  >  1:5357399/1‑68 (MQ=255)
                                   aaaCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCg                                  >  1:10054196/1‑68 (MQ=255)
                                       ttgCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc                              <  1:11713835/66‑1 (MQ=255)
                                       tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc                              <  1:29387535/68‑1 (MQ=255)
                                       tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc                              <  1:5190392/68‑1 (MQ=255)
                                       tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc                              <  1:4966475/68‑1 (MQ=255)
                                       tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc                              <  1:4474829/68‑1 (MQ=255)
                                       tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc                              <  1:3883491/68‑1 (MQ=255)
                                       tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc                              <  1:32006993/68‑1 (MQ=255)
                                       tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc                              <  1:6376378/68‑1 (MQ=255)
                                       tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc                              <  1:665200/68‑1 (MQ=255)
                                       tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc                              <  1:18606576/68‑1 (MQ=255)
                                       tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc                              <  1:29807700/68‑1 (MQ=255)
                                       tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc                              <  1:7438333/68‑1 (MQ=255)
                                       tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc                              <  1:28093146/68‑1 (MQ=255)
                                       tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc                              <  1:26922702/68‑1 (MQ=255)
                                        ggCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAccc                             <  1:23239725/68‑1 (MQ=255)
                                        ggCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAccc                             <  1:9423352/68‑1 (MQ=255)
                                        ggCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAccc                             <  1:22011428/68‑1 (MQ=255)
                                              aGTGGCGAGCTAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTa                          <  1:18482518/65‑1 (MQ=255)
                                              aGTGGCGAGCGAATATCTTCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTa                          <  1:28791874/65‑1 (MQ=255)
                                              aGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTa                          <  1:1905790/65‑1 (MQ=255)
                                              aGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTa                          <  1:16985651/65‑1 (MQ=255)
                                              aGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTa                          <  1:26494501/65‑1 (MQ=255)
                                              aGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTa                          <  1:1632237/65‑1 (MQ=255)
                                              aGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTa                          <  1:31764019/65‑1 (MQ=255)
                                              aGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTa                          <  1:12608509/65‑1 (MQ=255)
                                              aGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTa                          <  1:30399429/65‑1 (MQ=255)
                                              aGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTa                          <  1:30265592/65‑1 (MQ=255)
                                                      gcgaATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGa                <  1:11332903/67‑1 (MQ=255)
                                                          atatCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCaa            <  1:9702192/67‑1 (MQ=255)
                                                          atatCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCaa            <  1:23564286/67‑1 (MQ=255)
                                                          atatCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCaa            <  1:24967350/67‑1 (MQ=255)
                                                          atatCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAAt           <  1:22918190/68‑1 (MQ=255)
                                                          atatCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAAt           <  1:7161701/68‑1 (MQ=255)
                                                             tCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCg        >  1:3379844/1‑68 (MQ=255)
                                                             tCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCg        >  1:22500898/1‑68 (MQ=255)
                                                              cTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCgg       <  1:19363302/68‑1 (MQ=255)
                                                              cTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCgg       <  1:18197833/68‑1 (MQ=255)
                                                              cTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCgg       <  1:2005980/68‑1 (MQ=255)
                                                              cTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCgg       <  1:14133308/68‑1 (MQ=255)
                                                              cTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCgg       <  1:1321169/68‑1 (MQ=255)
                                                              cTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCgg       <  1:27747196/68‑1 (MQ=255)
                                                              cTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCgg       <  1:20537081/68‑1 (MQ=255)
                                                              cTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCgg       <  1:26637131/68‑1 (MQ=255)
                                                              cTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCgg       <  1:15843785/68‑1 (MQ=255)
                                                              cTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCgg       <  1:5644897/68‑1 (MQ=255)
                                                               tGCATAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCgg       >  1:17869415/1‑67 (MQ=255)
                                                                  gTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAg   <  1:8459697/68‑1 (MQ=255)
                                                                   |                                                                   
AACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGG  >  CP009273/4264204‑4264338

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: