Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,264,226 T→G 12.8% 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,2260TG12.8% 1582.3 / 182.9 522V58V (GTT→GTGssbsingle‑stranded DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (207/248);  new base G (32/35);  total (239/283)
Fisher's exact test for biased strand distribution p-value = 7.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-01

Only 100 of 573 total aligned reads displayed.

CTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCG  >  CP009273/4264159‑4264293
                                                                   |                                                                   
cTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGtt                                                                     >  1:6666083/1‑68 (MQ=255)
  aCTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                                                   <  1:32507637/68‑1 (MQ=255)
  aCTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                                                   <  1:468061/68‑1 (MQ=255)
  aCTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                                                   <  1:31614834/68‑1 (MQ=255)
     tCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                                                   >  1:20733567/1‑65 (MQ=255)
     tCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                                                   >  1:16666101/1‑65 (MQ=255)
     tCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                                                   >  1:1535441/1‑65 (MQ=255)
     tCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                                                   >  1:25712693/1‑65 (MQ=255)
     tCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCt                                                                <  1:17469864/68‑1 (MQ=255)
     tCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACGGAATGGCACCGCGTtgt                                                                   >  1:30666284/1‑65 (MQ=255)
       cGAATCCTTGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                                              <  1:16218838/68‑1 (MQ=255)
       cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATTGCACCGCGTTGTGCTGt                                                              <  1:10839482/68‑1 (MQ=255)
       cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                                              <  1:3454107/68‑1 (MQ=255)
       cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                                              <  1:4822432/68‑1 (MQ=255)
       cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                                              <  1:969478/68‑1 (MQ=255)
          aTCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCg                                                           >  1:929075/1‑68 (MQ=255)
          aTCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCg                                                           >  1:29566775/1‑68 (MQ=255)
            gcTGGCGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGc                                                         >  1:13479826/2‑68 (MQ=255)
            gcTGGCGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGc                                                         >  1:33863545/2‑68 (MQ=255)
            ccTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGc                                                         >  1:11377089/1‑68 (MQ=255)
            ccTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGc                                                         >  1:21032520/1‑68 (MQ=255)
            ccTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGc                                                         >  1:18213998/1‑68 (MQ=255)
              tGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCaa                                                       <  1:30458430/68‑1 (MQ=255)
                  gTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTg                                                   <  1:23498670/68‑1 (MQ=255)
                    gATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTgg                                                  <  1:27106477/67‑1 (MQ=255)
                     aTAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCa                                                >  1:9041397/1‑68 (MQ=255)
                         acagACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAg                                            >  1:15003546/4‑68 (MQ=255)
                         aGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAg                                            >  1:7131403/1‑68 (MQ=255)
                         aGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAg                                            >  1:24685439/1‑68 (MQ=255)
                         aGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAg                                            >  1:2775/1‑68 (MQ=255)
                          gCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGt                                           >  1:8752722/1‑68 (MQ=255)
                          gCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGt                                           >  1:28056846/1‑68 (MQ=255)
                           cGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTg                                          >  1:32926038/1‑68 (MQ=255)
                           agACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTg                                          >  1:4172381/2‑68 (MQ=255)
                               cGGCGAGATGAAAGATCAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcga                                      <  1:26319181/68‑1 (MQ=255)
                                  cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc                                    <  1:32723880/67‑1 (MQ=255)
                                  cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc                                    <  1:13863843/67‑1 (MQ=255)
                                  cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc                                    <  1:13088647/67‑1 (MQ=255)
                                  cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc                                    <  1:22814366/67‑1 (MQ=255)
                                  cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc                                    <  1:12737091/67‑1 (MQ=255)
                                   gagaTGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagcga                                  <  1:14659905/68‑1 (MQ=255)
                                   gagaTGAAAGAACAGACTGAATGGCACCGCGGTGGGCGGTTCGGCAAACTGGCAGAAGTGgcgagcga                                  <  1:3849406/68‑1 (MQ=255)
                                           aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc                          <  1:31275561/68‑1 (MQ=255)
                                           aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc                          <  1:22973896/68‑1 (MQ=255)
                                           aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc                          <  1:1237528/68‑1 (MQ=255)
                                           aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc                          <  1:10258563/68‑1 (MQ=255)
                                           aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc                          <  1:21538708/68‑1 (MQ=255)
                                            gagcagacagaATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGc                          <  1:1199234/64‑1 (MQ=255)
                                            gagcagacagaATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGc                          <  1:4966601/64‑1 (MQ=255)
                                            gagcagacagaATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGc                          >  1:26505836/4‑67 (MQ=255)
                                            gAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCg                         >  1:30221139/1‑68 (MQ=255)
                                            gAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCg                         >  1:19304692/1‑68 (MQ=255)
                                             aaCAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGt                        >  1:21672850/1‑68 (MQ=255)
                                              gcagacagaATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCgca                       >  1:25392652/2‑66 (MQ=255)
                                              aCAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCg                         <  1:5928976/66‑1 (MQ=255)
                                              aCAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTa                       >  1:19869826/1‑68 (MQ=255)
                                              aCAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTa                       >  1:34033152/1‑68 (MQ=255)
                                                  aCTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAgg                   >  1:6428702/1‑68 (MQ=255)
                                                  aCTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAgg                   >  1:7706872/1‑68 (MQ=255)
                                                  aCTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAgg                   >  1:15495873/1‑68 (MQ=255)
                                                   cTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTaaa                     <  1:14196168/65‑1 (MQ=255)
                                                   cTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGt                  <  1:3086799/68‑1 (MQ=255)
                                                   cTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGt                  <  1:3973668/68‑1 (MQ=255)
                                                   cTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGt                  >  1:5423292/1‑68 (MQ=255)
                                                   cTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGt                  >  1:4436586/1‑68 (MQ=255)
                                                    agAATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAgggc                  >  1:2200578/2‑66 (MQ=255)
                                                        tGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGATTATCTGCGTAAAGGTTCTCa             >  1:26084057/1‑68 (MQ=255)
                                                        tGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCa             >  1:17326165/1‑68 (MQ=255)
                                                        tGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCa             >  1:25028347/1‑68 (MQ=255)
                                                        tGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCa             >  1:2445571/1‑68 (MQ=255)
                                                          gCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg           >  1:28734029/1‑68 (MQ=255)
                                                          gCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg           >  1:19086044/1‑68 (MQ=255)
                                                           cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg           <  1:9824227/67‑1 (MQ=255)
                                                           cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg           <  1:9478163/67‑1 (MQ=255)
                                                           cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg           <  1:1165277/67‑1 (MQ=255)
                                                           cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg           <  1:14428438/67‑1 (MQ=255)
                                                           cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg           <  1:7874789/67‑1 (MQ=255)
                                                           cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg           <  1:16950088/67‑1 (MQ=255)
                                                           cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg           <  1:18506735/67‑1 (MQ=255)
                                                           cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg           <  1:29037630/67‑1 (MQ=255)
                                                           cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg           <  1:10013213/67‑1 (MQ=255)
                                                           cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg           <  1:3357193/67‑1 (MQ=255)
                                                            atcgcgTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGt          <  1:33224163/65‑1 (MQ=255)
                                                            acagcgTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt         >  1:16055092/4‑68 (MQ=255)
                                                            aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt         >  1:17154643/1‑68 (MQ=255)
                                                            aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt         >  1:10541844/1‑68 (MQ=255)
                                                            aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt         >  1:8519465/1‑68 (MQ=255)
                                                            aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt         >  1:8238380/1‑68 (MQ=255)
                                                            aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt         >  1:7365535/1‑68 (MQ=255)
                                                            aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt         >  1:26941511/1‑68 (MQ=255)
                                                            aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt         >  1:17673843/1‑68 (MQ=255)
                                                            aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt         >  1:17959751/1‑68 (MQ=255)
                                                            aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt         >  1:23887033/1‑68 (MQ=255)
                                                            aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt         >  1:19631412/1‑68 (MQ=255)
                                                            aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt         >  1:20431963/1‑68 (MQ=255)
                                                            aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt         >  1:2331885/1‑68 (MQ=255)
                                                             ccGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGttt        <  1:28328194/68‑1 (MQ=255)
                                                               gcgTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTtat      <  1:22348769/68‑1 (MQ=255)
                                                                 gTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTtata     <  1:3195315/67‑1 (MQ=255)
                                                                   tgtgCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATc   <  1:29383648/67‑1 (MQ=255)
                                                                   |                                                                   
CTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCG  >  CP009273/4264159‑4264293

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: