Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,264,226 T→G 15.3% V58V (GTT→GTG ssb → single‑stranded DNA‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,264,2260TG15.3% 1597.3 / 239.7 566V58V (GTT→GTGssbsingle‑stranded DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (211/266);  new base G (37/50);  total (248/318)
Fisher's exact test for biased strand distribution p-value = 8.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 620 total aligned reads displayed.

CTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATC  >  CP009273/4264159‑4264292
                                                                   |                                                                  
cTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGtt                                                                    >  1:27406608/1‑68 (MQ=255)
  aCTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                                                  <  1:17256378/68‑1 (MQ=255)
  aCTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                                                  <  1:11070020/68‑1 (MQ=255)
  aCTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                                                  <  1:6503395/68‑1 (MQ=255)
     tCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                                                  >  1:32059117/1‑65 (MQ=255)
       cGATTCCTGGCGTGATAAAGCGTCCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                                             <  1:3640564/68‑1 (MQ=255)
       cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                                             <  1:21902254/68‑1 (MQ=255)
       cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                                             <  1:27777499/68‑1 (MQ=255)
       cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                                             <  1:11499729/68‑1 (MQ=255)
       cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                                             <  1:5010041/68‑1 (MQ=255)
       cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                                             <  1:32864091/68‑1 (MQ=255)
       cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                                             <  1:34427040/68‑1 (MQ=255)
            gcTGGCGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGc                                                        >  1:13829127/2‑68 (MQ=255)
            gcTGGCGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGc                                                        >  1:26991409/2‑68 (MQ=255)
            gcTGGCGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGc                                                        >  1:1572504/2‑68 (MQ=255)
            ccTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGc                                                        >  1:20878501/1‑68 (MQ=255)
            ccTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGc                                                        >  1:132370/1‑68 (MQ=255)
               ggCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCaaa                                                     <  1:34080872/68‑1 (MQ=255)
               ggCGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCaag                                                     <  1:1612609/68‑2 (MQ=255)
               ggCGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCaag                                                     <  1:14493171/68‑2 (MQ=255)
                 cGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCaa                                                      >  1:33304629/1‑65 (MQ=255)
                 cGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCaagc                                                    <  1:22100162/67‑3 (MQ=255)
                 cGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCaagc                                                    <  1:721689/67‑3 (MQ=255)
                  gTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTg                                                  <  1:15439358/68‑1 (MQ=255)
                  gTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTg                                                  <  1:18658727/68‑1 (MQ=255)
                  gTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTg                                                  <  1:12494034/68‑1 (MQ=255)
                  gTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTg                                                  >  1:24410449/1‑68 (MQ=255)
                     aTAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCa                                               >  1:2928947/1‑68 (MQ=255)
                        aaGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGaa                                            >  1:29968255/1‑68 (MQ=255)
                           cGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTg                                         >  1:3998161/1‑68 (MQ=255)
                            gACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTgg                                        <  1:16679022/68‑1 (MQ=255)
                            gACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTgg                                        <  1:30757722/68‑1 (MQ=255)
                                  cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc                                   <  1:10000779/67‑1 (MQ=255)
                                  cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc                                   <  1:29492483/67‑1 (MQ=255)
                                  cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc                                   <  1:21528928/67‑1 (MQ=255)
                                  cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc                                   <  1:33517494/67‑1 (MQ=255)
                                  cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc                                   <  1:17701029/67‑1 (MQ=255)
                                  cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc                                   <  1:903475/67‑1 (MQ=255)
                                   gagaTGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagcga                                 <  1:22819566/68‑1 (MQ=255)
                                   gagaTGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagcga                                 >  1:6663306/1‑68 (MQ=255)
                                          aaGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCTAGCGAATATCTg                          <  1:20563111/68‑1 (MQ=255)
                                          aaGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTg                          <  1:20684362/68‑1 (MQ=255)
                                           aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc                         <  1:28732222/68‑1 (MQ=255)
                                           aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc                         <  1:28423555/68‑1 (MQ=255)
                                           aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc                         <  1:8335246/68‑1 (MQ=255)
                                           aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc                         <  1:10284450/68‑1 (MQ=255)
                                           aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc                         >  1:26198178/1‑68 (MQ=255)
                                           aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc                         <  1:12184096/68‑1 (MQ=255)
                                           aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc                         <  1:12903423/68‑1 (MQ=255)
                                           aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc                         <  1:23940253/68‑1 (MQ=255)
                                           aGAACAGACTGAATGGCACCGCGTAGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc                         <  1:17392943/68‑1 (MQ=255)
                                            gagcagacagaATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGc                         <  1:31111888/64‑1 (MQ=255)
                                            gagcagacagaATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGc                         <  1:23204716/64‑1 (MQ=255)
                                            gagcagacagaATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGc                         <  1:4356757/64‑1 (MQ=255)
                                            gagcagacagaATGGCATCGCGTGGAGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGc                         <  1:19288403/64‑1 (MQ=255)
                                            gagcagacagaATGGCATCGCGGGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGc                         <  1:14970024/64‑1 (MQ=255)
                                              gcagacagaATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCgca                      >  1:5237860/2‑66 (MQ=255)
                                              aCAGACTGGATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTa                      >  1:5931752/1‑68 (MQ=255)
                                              aCAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCAGCGTa                      >  1:36768721/1‑68 (MQ=255)
                                               cAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTaa                     >  1:33008400/1‑68 (MQ=255)
                                                  aCTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAgg                  >  1:28180886/1‑68 (MQ=255)
                                                  aCTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAgg                  >  1:19563211/1‑68 (MQ=255)
                                                   cTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGt                 >  1:29687065/1‑68 (MQ=255)
                                                   cTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGt                 <  1:29006837/68‑1 (MQ=255)
                                                   cTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGt                 <  1:20125744/68‑1 (MQ=255)
                                                     gAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTtc               <  1:35914176/68‑1 (MQ=255)
                                                       aTGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAgggcgc               <  1:3022927/66‑4 (MQ=255)
                                                        tGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCa            >  1:36033114/1‑68 (MQ=255)
                                                        tGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCa            >  1:11264969/1‑68 (MQ=255)
                                                          gCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg          >  1:32671806/1‑68 (MQ=255)
                                                           cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg          <  1:37134970/67‑1 (MQ=255)
                                                           cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg          <  1:7681113/67‑1 (MQ=255)
                                                           cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg          <  1:7961996/67‑1 (MQ=255)
                                                           cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg          <  1:4570089/67‑1 (MQ=255)
                                                           cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg          <  1:13499712/67‑1 (MQ=255)
                                                           cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg          <  1:26718204/67‑1 (MQ=255)
                                                           cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg          <  1:25225391/67‑1 (MQ=255)
                                                           cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg          <  1:18070681/67‑1 (MQ=255)
                                                           cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg          <  1:33138852/67‑1 (MQ=255)
                                                           cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg          <  1:32294664/67‑1 (MQ=255)
                                                           cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg          <  1:3152534/67‑1 (MQ=255)
                                                           cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg          <  1:22269026/67‑1 (MQ=255)
                                                           cACCGCGTGGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg          <  1:21219046/67‑1 (MQ=255)
                                                           cACCGCGTAGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg          <  1:16972927/67‑1 (MQ=255)
                                                            atcgcgTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGt         <  1:2493176/65‑1 (MQ=255)
                                                            aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt        >  1:35120251/1‑68 (MQ=255)
                                                            aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt        >  1:5510550/1‑68 (MQ=255)
                                                            aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt        >  1:23676761/1‑68 (MQ=255)
                                                            aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt        >  1:22658941/1‑68 (MQ=255)
                                                            aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt        >  1:19827896/1‑68 (MQ=255)
                                                            aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt        >  1:36684820/1‑68 (MQ=255)
                                                            aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt        >  1:35585969/1‑68 (MQ=255)
                                                            aCCGCGTTGTGCGGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt        >  1:34785390/1‑68 (MQ=255)
                                                             tcgcgTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGtct       <  1:9558709/67‑3 (MQ=255)
                                                              cgcgTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGtct       >  1:25995617/1‑65 (MQ=255)
                                                               gcgTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCtac     <  1:25511374/68‑2 (MQ=255)
                                                                 gTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTtata    >  1:6834324/1‑67 (MQ=255)
                                                                 gTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTatat   <  1:16367608/68‑1 (MQ=255)
                                                                 gTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTatat   <  1:10846159/68‑1 (MQ=255)
                                                                   tgtgCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATc  <  1:11842288/67‑1 (MQ=255)
                                                                   |                                                                  
CTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATC  >  CP009273/4264159‑4264292

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: