Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,264,247 A→G 14.2% A65A (GCA→GCG ssb → single‑stranded DNA‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,264,2470AG14.2% 1915.1 / 240.0 663A65A (GCA→GCGssbsingle‑stranded DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (204/365);  new base G (37/57);  total (241/422)
Fisher's exact test for biased strand distribution p-value = 5.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01

Only 100 of 668 total aligned reads displayed.

ATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTA  >  CP009273/4264180‑4264314
                                                                   |                                                                   
 tAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCa                                                                     <  1:20963901/67‑1 (MQ=255)
  aaaGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGa                                                                   >  1:11710066/1‑68 (MQ=255)
    aGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAg                                                                 >  1:18888650/1‑68 (MQ=255)
    aGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAg                                                                 >  1:21326380/1‑68 (MQ=255)
           ggCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgag                                                          <  1:9733140/68‑1 (MQ=255)
           ggCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgag                                                          <  1:29705758/68‑1 (MQ=255)
             cGAGATGATAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc                                                         <  1:12238258/67‑1 (MQ=255)
             cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc                                                         <  1:15565318/67‑1 (MQ=255)
             cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc                                                         <  1:18499032/67‑1 (MQ=255)
             cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc                                                         <  1:32583949/67‑1 (MQ=255)
             cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc                                                         <  1:4917045/67‑1 (MQ=255)
             cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc                                                         <  1:5835644/67‑1 (MQ=255)
              gagaTGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagcga                                                       <  1:29976805/68‑1 (MQ=255)
              gagaTGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagcga                                                       <  1:27089475/68‑1 (MQ=255)
                      aGAACATACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc                                               <  1:34025802/68‑1 (MQ=255)
                      aGAACAGACTGAATGTCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc                                               <  1:12611698/68‑1 (MQ=255)
                      aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc                                               <  1:29381215/68‑1 (MQ=255)
                      aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc                                               <  1:27602749/68‑1 (MQ=255)
                      aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc                                               <  1:24530803/68‑1 (MQ=255)
                      aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc                                               <  1:16376104/68‑1 (MQ=255)
                      aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc                                               <  1:24763888/68‑1 (MQ=255)
                      aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc                                               <  1:15113208/68‑1 (MQ=255)
                      aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGAGAATATCTGc                                               <  1:29015214/68‑1 (MQ=255)
                       gagcagacagaATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGc                                               <  1:2084107/64‑1 (MQ=255)
                       gAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCg                                              >  1:10423969/1‑68 (MQ=255)
                       gAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCg                                              >  1:9414345/1‑68 (MQ=255)
                         aCAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTa                                            >  1:15991380/1‑68 (MQ=255)
                         aCAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTa                                            >  1:4207152/1‑68 (MQ=255)
                         aCAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTa                                            >  1:16737129/1‑68 (MQ=255)
                         aCAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTa                                            >  1:32165101/1‑68 (MQ=255)
                         aCAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTa                                            <  1:20361756/68‑1 (MQ=255)
                           aGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTaaa                                          >  1:22881378/1‑68 (MQ=255)
                              cagAATGGCATCTCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCaag                                          <  1:23461949/63‑2 (MQ=255)
                              cTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGt                                       >  1:21995679/1‑68 (MQ=255)
                              cTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGt                                       >  1:19718289/1‑68 (MQ=255)
                              cTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGt                                       <  1:12900785/68‑1 (MQ=255)
                               agAATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAgggc                                       >  1:10822061/2‑66 (MQ=255)
                                gAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTtc                                     <  1:23849260/68‑1 (MQ=255)
                                  aTGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAgggcgc                                     <  1:18096176/66‑4 (MQ=255)
                                  aTGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAggg                                        <  1:25895092/63‑1 (MQ=255)
                                   tGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCa                                  >  1:7137413/1‑68 (MQ=255)
                                      cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                <  1:7854677/67‑1 (MQ=255)
                                      cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                <  1:30566257/67‑1 (MQ=255)
                                      cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                <  1:27873629/67‑1 (MQ=255)
                                      cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                <  1:31223090/67‑1 (MQ=255)
                                      cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                <  1:822604/67‑1 (MQ=255)
                                      cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                <  1:10984139/67‑1 (MQ=255)
                                      cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                <  1:20132047/67‑1 (MQ=255)
                                      cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                <  1:8865775/67‑1 (MQ=255)
                                      cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                <  1:15356096/67‑1 (MQ=255)
                                      cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                <  1:21446254/67‑1 (MQ=255)
                                      cACCGCGTTGTGCCGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                <  1:10028724/67‑1 (MQ=255)
                                       aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATGTCTGCGTAAAGGTTCTCAGGtt                              >  1:21611743/1‑68 (MQ=255)
                                       aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                              >  1:13542757/1‑68 (MQ=255)
                                       aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                              >  1:60862/1‑68 (MQ=255)
                                       aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                              >  1:31862646/1‑68 (MQ=255)
                                       aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                              >  1:3338580/1‑68 (MQ=255)
                                       aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                              >  1:11368053/1‑68 (MQ=255)
                                       aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                              >  1:8432168/1‑68 (MQ=255)
                                          gcgTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTtat                           <  1:30325913/68‑1 (MQ=255)
                                            gTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGTGCGAATATCTGCGTAAAGGTTCTCAGGTTTatat                         <  1:22551871/68‑1 (MQ=255)
                                            gTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTtata                          <  1:14197361/67‑1 (MQ=255)
                                            gTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTtata                          <  1:31027888/67‑1 (MQ=255)
                                            gTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTtata                          <  1:13873916/67‑1 (MQ=255)
                                              ggtgCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATc                        <  1:3790119/66‑1 (MQ=255)
                                                 gCTTTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGaa                     <  1:869849/67‑1 (MQ=255)
                                                 gCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAg                    >  1:6556156/1‑68 (MQ=255)
                                                 gCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAg                    >  1:33594999/1‑68 (MQ=255)
                                                 gCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGaa                     <  1:18654898/67‑1 (MQ=255)
                                                 gCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGaa                     <  1:25546321/67‑1 (MQ=255)
                                                 gCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAg                    >  1:26659158/1‑68 (MQ=255)
                                                 gCTATTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATGGAAg                    >  1:1731951/1‑68 (MQ=255)
                                                   tGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAgg                   <  1:19190973/67‑1 (MQ=255)
                                                   tGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAgg                   <  1:1481122/67‑1 (MQ=255)
                                                   tGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAgg                   <  1:13158059/67‑1 (MQ=255)
                                                   tGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAgg                   <  1:34705128/67‑1 (MQ=255)
                                                      tCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCa                >  1:12452683/1‑67 (MQ=255)
                                                      tCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCa                >  1:31484894/1‑67 (MQ=255)
                                                       cGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGc              <  1:6288502/68‑1 (MQ=255)
                                                       cGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGc              <  1:4450707/68‑1 (MQ=255)
                                                           aaaCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGc           <  1:16252320/67‑1 (MQ=255)
                                                           aaaCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGc           <  1:24414012/67‑1 (MQ=255)
                                                           aaaCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCg          >  1:32851493/1‑68 (MQ=255)
                                                           aaaCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCg          >  1:11915551/1‑68 (MQ=255)
                                                           aaaCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCg          >  1:28360385/1‑68 (MQ=255)
                                                               tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc      <  1:17169972/68‑1 (MQ=255)
                                                               tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc      <  1:34412415/68‑1 (MQ=255)
                                                               tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc      <  1:17677291/68‑1 (MQ=255)
                                                               tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc      <  1:914772/68‑1 (MQ=255)
                                                               tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc      <  1:20606109/68‑1 (MQ=255)
                                                               tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc      <  1:14408892/68‑1 (MQ=255)
                                                               tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc      <  1:24261703/68‑1 (MQ=255)
                                                               tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc      <  1:25074787/68‑1 (MQ=255)
                                                               tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc      <  1:7537118/68‑1 (MQ=255)
                                                               tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc      <  1:25380624/68‑1 (MQ=255)
                                                               tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc      <  1:27261497/68‑1 (MQ=255)
                                                               tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc      <  1:5585948/68‑1 (MQ=255)
                                                               tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc      <  1:5382521/68‑1 (MQ=255)
                                                               tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc      <  1:11558867/68‑1 (MQ=255)
                                                                 gCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGt                               <  1:29450887/41‑1 (MQ=255)
                                                                   |                                                                   
ATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTA  >  CP009273/4264180‑4264314

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: