Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,447,085 C→T 30.6% A68T (GCT→ACT)  gspB ← part of gsp divergon involved in type II protein secretion

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,447,0850CT30.6% 529.9 / inf 392A68T (GCT→ACT) gspBpart of gsp divergon involved in type II protein secretion
Reads supporting (aligned to +/- strand):  ref base C (74/196);  new base T (35/84);  total (110/282)
Fisher's exact test for biased strand distribution p-value = 7.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

Only 100 of 425 total aligned reads displayed.

CGACGGGGAGTGGCTGTTTTTTGAAGAAATAGTGTTGTGTTGATTGTTTAAAAGCGGGCGTTTCAGCCTCTACAGTGACTTCGGCTTTTACTATCCAAAGCACCCATAGCTGCCTGGCATACCCACCAACCAGC  >  CP009273/3447019‑3447152
                                                                  |                                                                   
cgacgGGGAGTGGCTGTTTTTTGAAGAAATAGTGTTGTGTTGATTGTTTAAAAGCGGGCGTTTCAGcc                                                                    <  1:33741788/68‑1 (MQ=255)
cgacgGGGAGTGGCTGTTTTTTGAAGAAATAGTGTTGTGTTGATTGTTTAAAAGCGGGCGTTTCAGcc                                                                    <  1:18781740/68‑1 (MQ=255)
cgacgGGGAGTGGCTGTTTTTTGAAGAAATAGTGTTGTGTTGATTGTTTAAAAGCGGGCGTTTCAGcc                                                                    <  1:21593358/68‑1 (MQ=255)
   cgGGGAGTGGCTGTTTTTTGAAGAAATAGTGTTGTGTTGATTGTTTAAAAGCGGGCGTTTCAGtctc                                                                  <  1:31291871/67‑1 (MQ=255)
   cgGGGAGTGGCTGTTTTTTGAAGAAATAGTGTTGTGTTGATTGTTTAAAAGCGGGCGTTTCAGTctct                                                                 <  1:7863808/68‑1 (MQ=255)
   cgGGGAGTGGCTGTTTTTTGAAGAAATAGTGTTGTGTTGATTGTTTAAAAGCGGGCGTTTCAGCctct                                                                 <  1:16893963/68‑1 (MQ=255)
   cgGGGAGTGGCTGTTTTTTGAAGAAATAGTGTTGTGTTGATTGTTTAAAAGCGGGCGTTTCAGCctc                                                                  <  1:35341374/67‑1 (MQ=255)
    ggggAGTGGCTGTTTTTTGAAGAAATAGTGTTGTGTTGATTGTTTAAAAGCGGGCGTTTCAGCCTCTa                                                                <  1:25257593/68‑1 (MQ=255)
      ggAGTGGCTGTTTTTTGAAGAAATAGTGTTGTGTTGATTGTTTAAAAGCGGGCGTTTCAGTCTCTAc                                                               <  1:4454912/67‑1 (MQ=255)
      ggAGTGGCTGTTTTTTGAAGAAATAGTGTTGTGTTGATTGTTTAAAAGCGGGCGTTTCAGCCTCTAc                                                               <  1:6286005/67‑1 (MQ=255)
         gTGGCTGTTTTTTGAAGAAATAGTGTTGTGTTGATTGTTTAAAAGCGGGCGTTTCAGCCTCTACAGTg                                                           >  1:3690158/1‑68 (MQ=255)
         gTGGCTGTTTTTTGAAGAAATAGTGTTGTGTTGATTGTTTAAAAGCGGGCGTTTCAGCCTCTACAGTg                                                           >  1:22966184/1‑68 (MQ=255)
           ggCTGTTTTTTGTAGAAATTGTGTTGTGTTGATTGTTTCAAAGCGGGCGTTTCAGCCTCTACAGTGa                                                          <  1:35882801/67‑1 (MQ=255)
           ggCTGTTTTTTGAAGAAATAGTGTTGTGTTGATTGTTTAAAAGCGGGCGTTTCAGTCTCTACAGTGa                                                          <  1:2645752/67‑1 (MQ=255)
           ggCTGTTTTTTGAAGAAATAGTGTTGTGTTGATTGTTTAAAAGCGGGCGTTTCAGCCTCTACAGTGa                                                          <  1:4899193/67‑1 (MQ=255)
              tGTTTTTTGAAGAAATAGTGTTGTGTTGATTGTTTAAAAGCGGGCGTTTCAGTCTCTACAGTGACTTc                                                      >  1:29604207/1‑68 (MQ=255)
              tGTTTTTTGAAGAAATAGTGTTGTGTTGATTGTTTAAAAGCGGGCGTTTCAGCCTCTACAGTGACTTc                                                      >  1:30339852/1‑68 (MQ=255)
                 tttttGAAGAAATAGTGTTGTGTTGATTGTTTAAAAGCGGGCGTTTCAGTCTCTACAGTGACTTCgg                                                    <  1:23235103/67‑1 (MQ=255)
                 tttttGAAGAAATAGTGTTGTGTTGATTGTTTAAAAGCGGGCGTTTCAGTCTCTACAGTGACTTCgg                                                    <  1:31588913/67‑1 (MQ=255)
                 tttttGAAGAAATAGTGTTGTGTTGATTGTTTAAAAGCGGGCGTTTCAGTCTCTACAGTGACTTCgg                                                    <  1:17647871/67‑1 (MQ=255)
                 tttttGAAGAAATAGTGTTGTGTTGATTGTTTAAAAGCGGGCGTTTCAGTCTCTACAGTGACTTCgg                                                    <  1:11753913/67‑1 (MQ=255)
                 tttttGAAGAAATAGTGTTGTGTTGATTGTTTAAAAGCGGGCGTTTCAGCCTCTACAGTGACTTCgg                                                    <  1:18928098/67‑1 (MQ=255)
                 tttttGAAGAAATAGTGTTGTGTTGATTGTTTAAAAGCGGGCGTTTCAGCCTCTACAGTGACTTCgg                                                    <  1:7770720/67‑1 (MQ=255)
                 tttttGAAGAAATAGTGTTGTGTTGATTGTTTAAAAGCGGGCGTTTCAGCCTCTACAGTGACTTCgg                                                    <  1:5986838/67‑1 (MQ=255)
                   tttGAAGAAATAGTGTTGTGTTGATTGTTTAAAAGCGGGCGTTTCAGTCTCTACAGTGACTTCGGCtt                                                 >  1:30612432/1‑68 (MQ=255)
                     tGAAGAAATAGTGTTGTGTTGATTGTTTAAACGCGGGCGTTTCAGGCTCTACAGTGACTTCGGCtttt                                               <  1:8961702/68‑1 (MQ=255)
                     tGAAGAAATAGTGTTGTGTTGATTGTTTAAAAGCGGGCGTTTCAGCCTCTACAGTGACTTCGGCtttt                                               <  1:26032317/68‑1 (MQ=255)
                     tGAAGAAATAGTGTTGTGTTGATTGTTTAAAAGCGGGCGTTTCAGCCTCTACAGTGACTTCGGCtttt                                               <  1:34972399/68‑1 (MQ=255)
                     tGAAGAAATAGTGTTGTGTTGATTGTTTAAAAGCGGGCGTTTCAGCCTCTACAGTGACTTCGGCtttt                                               <  1:34210406/68‑1 (MQ=255)
                     tGAAGAAATAGTGTTGTGTTGATTGTTTAAAAGCGGGCGTTTCAGCCTCTACAGTGACTTCGGCtttt                                               <  1:17488669/68‑1 (MQ=255)
                     tGAAGAAATAGTGTTGTGTTGATTGTTTAAAAGCGGGCGTTTCAGCCTCTACAGTGACTTCGGCtttt                                               <  1:5131674/68‑1 (MQ=255)
                     tGAAGAAATAGTGTGGTGTTGATTGTTTAAAAGCGGGCGTTTCAGTCTCTACAGTGACTTCGGCtttt                                               <  1:20781553/68‑1 (MQ=255)
                        agaaATAGTGTTGTGTTGATTGTTTAAAAGCGGGCGTTTCAGCCTCTACAGTGACTTCGGCTTTTa                                              >  1:25923775/1‑66 (MQ=255)
                        agaaATAGTGTTGTGTTGATTGTTTAAAAGCGGGCGTTTCAGACTCTACATTGACTTCGGCTTTTACt                                            >  1:4147564/1‑68 (MQ=255)
                           aaTAGTGTTGTGTTGATTGTTTAAAAGCGGGCGTTTCAGCCTCTACAGTGACTTCGGCTTTTACTAt                                          <  1:27109330/67‑1 (MQ=255)
                           aaTAGTGTTGTGTTGATTGTTTAAAAGCGGGCGTTTCAGCCTCTACAGTGACTTCGGCTTTTACTAt                                          <  1:13434069/67‑1 (MQ=255)
                             tAGTGTTGTGTTGATTGTTTAAAAGCGGGCGTTTCAGCCTCTACAGTGACTTCGGCTTTTACTATCCa                                       >  1:29276184/1‑68 (MQ=255)
                              aGTGTTGTGTTGATTGTTTAAAAGCGGGCGTTTCAGCCTCTACAGTGACTTCGGCTTTTACTATCCaa                                      <  1:24276507/68‑1 (MQ=255)
                                tgttgtgttgATTGTTTAAAAGCGGGCGTTTCAGCCTCTACAGTGACTTCGGCTTTTACTATCCaaa                                     >  1:2388062/1‑67 (MQ=255)
                                 gttgtgttgATTGTTTAAAAGCGGGCGTTTCAGCCTCTACAGTGACTTCGGCTTTTACTATCCAAAGc                                   >  1:1006385/1‑68 (MQ=255)
                                 gttgtgttgATTGTTTAAAAGCGGGCGTTTCAGCCTCTACAGTGACTTCGGCTTTTACTATCCAAAGc                                   >  1:4360111/1‑68 (MQ=255)
                                   tgtgtttATTGTTTAAAAGCGTTGGTTTCAGTCTCTACAGTGTGTTCGTGTTTTACTTTCCAATGCAc                                 <  1:26715876/68‑1 (MQ=255)
                                   tgtgttgATTGTTTAAAAGCGGGCGTTTTAGTCTCTACAGTGACTTCGGCTTTTACTATCCAAAGCAc                                 <  1:8280850/68‑1 (MQ=255)
                                   tgtgttgATTGTTTAAAAGCGGGCGTTTCAGTCTCTACAGTGACTTCGGCTTTTACTATCCAAAGCa                                  <  1:29143794/67‑1 (MQ=255)
                                   tgtgttgATTGTTTAAAAGCGGGCGTTTCAGCCTCTACAGTGACTTCGGCTTTTACTATCCAAAgaa                                  >  1:14550317/1‑65 (MQ=255)
                                   tgtgttgATTGTTTAAAAGCGGGCGTTTCAGCCTCTACAGTGACTTCGGCTTTTACTATCCAAAGCa                                  >  1:18306849/1‑67 (MQ=255)
                                   tgtgttgATTGTTTAAAAGCGGGCGTTTCAGCCTCTACAGTGACTTCGGCTTTTACTATCCAAAGCa                                  >  1:20433667/1‑67 (MQ=255)
                                   tgtgttgATTGTTTAAAAGCGGGCGTTTCAGCCTCTACAGTGACTTCGGCTTTTACTATCCAAAGCAc                                 <  1:15644814/68‑1 (MQ=255)
                                   tgtgttgATTGTTTAAAAGCGGGCGTTTCAGCCTCTACAGTGACTTCGGCTTTTACTATCCAAAGCAc                                 <  1:22668255/68‑1 (MQ=255)
                                   tgtgttgATTGTTTAAAAGCGGGCGTTTCAGCCTCTACAGTGACTTCGGCTTTTACTATCCAAAGCAc                                 <  1:14300651/68‑1 (MQ=255)
                                   tgtgttgATTGTTTAAAAGCGGGCGTTTCAGCCTCTACAGTGACTTCGGCTTTTACTATCCAAAGCAc                                 <  1:21053361/68‑1 (MQ=255)
                                   tgtgttgATTGTTTAAAAGCGGGCGTTTCAGCCTCTACAGTGACTTCGGCTTTTACTATCCAAAGCAc                                 <  1:14011211/68‑1 (MQ=255)
                                      gttgATTGTTTAAAAGCGGGCGTTTCAGTCTCTACAGTGACTTCGGCTTTTACTATCCAAAGCAccaa                              >  1:24008324/1‑66 (MQ=255)
                                      gttgATTGTTTAAAAGCGGGCGTTTCAGTCTCTACAGTGACTTCGGCTTTTACTATCCAAAGCACCCa                              >  1:35388940/1‑68 (MQ=255)
                                      gttgATTGTTTAAAAGCGGGCGTTTCAGCCTCTACAGTGACTTCGGCTTTTACTATCCAAAGCACCCa                              >  1:8642134/1‑68 (MQ=255)
                                        tgATTGTTTAAAAGCGGGCGTTTCAGCCTCTACAGTGACTTCGGCTTTTACTATCCAAAGCACCCATa                            >  1:32355949/1‑68 (MQ=255)
                                        tgATTGTGTAAAAGCGGGCGTTTCAGCCTCTACAGTGACTTCGGCTTTTACTATCCAAAGCACCCATa                            <  1:16520059/68‑1 (MQ=255)
                                          aTTGTTTAAAAGCGGGCGTTTCAGTCTCTACAGTGACTTCGGCTTTTACTATCCAAAGCACCCATa                            >  1:27436819/1‑66 (MQ=255)
                                           ttGTTTAAAAGCGGGCGTTTCAGTCTCTACAGTGACTTCGGCTTTTACTATCCAAAGCACCCATa                            >  1:11213456/1‑65 (MQ=255)
                                           ttGTTTAAAAGCGGGCGTTTCAGTCTCTACAGTGACTTCGGCTTTTACTATCCAAAGCACCCATa                            >  1:31740160/1‑65 (MQ=255)
                                           ttGTTTAAAAGCGGGCGTTTCAGCCTCTACAGTGACTTCGGCTTTTACTATCCAAAGCACCCATa                            >  1:26891514/1‑65 (MQ=255)
                                           ttGTTTAAAAGCGGGCGTTTCAGCCTCTACAGTGACTTCGGCTTTTACTATCCAAAGCACCCATAGc                          <  1:16734190/67‑1 (MQ=255)
                                            tGTTTAAAAGCGGGCGTTTCAGTCTCTACAGTGACTTCGGCTTTTACTATCCAAAGCACCCATAGCTg                        <  1:1934968/68‑1 (MQ=255)
                                            tGTTTAAAAGCGGGCGTTTCAGCCTCTACAGTGACTTCGGCTTTTACTATCCAAAGCACCCATAGCTg                        <  1:33216876/68‑1 (MQ=255)
                                             gTTTAAAAGCGGGCGTTTCAGCCTCTACAGTGACTTCGGCTTTTACTATCCAAAGCACCCATAGCTGc                       <  1:5638350/68‑1 (MQ=255)
                                             gTTTAAAAGCGGGCGTTTCAGCCTCTACAGTGACTTCGGCTTTTACTATCCAAAGCACCCATAGCTGc                       <  1:18251666/68‑1 (MQ=255)
                                               ttAAAAGCGGGCGTTTCAGCCTCTACAGTGACTTCGGCTTTTACTATCCAAAGCACCCATAGCTGCCt                     <  1:28839135/68‑1 (MQ=255)
                                               ttAAAAGCGGGCGTTTCAGCCTCTACAGTGACTTCGGCTTTTACTATCCAAAGCACCCATAGCTGCCt                     <  1:32005390/68‑1 (MQ=255)
                                               ttAAAAGCGGGCGTTTCAGCCTCTACAGTGACTTCGGCTTTTACTATCCAAAGCACCCATAGCTGCCt                     <  1:15358484/68‑1 (MQ=255)
                                                 aaaaGCGGGCGTTTCAGTCTCTACAGTGACTTCGGCTTTTACTATCCAAAGCACCCATAGCTGCCTgg                   <  1:34511929/68‑1 (MQ=255)
                                                 aaaaGCGGGCGTTTCAGTCTCTACAGTGACTTCGGCTTTTACTATCCAAAGCACCCATAGCTGCCTgg                   <  1:17239690/68‑1 (MQ=255)
                                                 aaaaGCGGGCGTTTCAGTCTCTACAGTGACTTCGGCTTTTACTATCCAAAGCACCCATAGCTGCCTgg                   <  1:1296240/68‑1 (MQ=255)
                                                 aaaaGCGGGCGTTTCAGTCTCTACAGTGACTTCGGCTTTTACTATCCAAAGCACCCATAGCTGCCTgg                   <  1:19466415/68‑1 (MQ=255)
                                                  aaaGCGGGCGTTTCAGCCTCTACAGTGACTTCGGCTTTTACTATCCAAAGCACCCATAGCTGCCTGGc                  <  1:25408082/68‑1 (MQ=255)
                                                   aaGCGGGCGTTTCAGTCTCTACAGTGACTTCGGCTTTTACTATCCAAAGCACCCATAGCTGCCTGGCa                 <  1:3652235/68‑1 (MQ=255)
                                                    aGCGGGCGTTTCAGCCTCTACAGTGACTTCGGCTTTTACTATCCAAAGCACCCATAGCTGCCTGGCAt                >  1:7061780/1‑68 (MQ=255)
                                                    aGCGGGCGTTTCAGCCTCTACAGTGACTTCGGCTTTTACTATCCAAAGCACCCATAGCTGCCTGGCAt                >  1:7562820/1‑68 (MQ=255)
                                                    aGCGGGCGTTTCAGCCTCTACAGTGACTTCGGCTTTTACTATCCAAAGCACCCATAGCTGCCTGGCAt                >  1:10296278/1‑68 (MQ=255)
                                                    aGCGGGCGTTTCAGCCTCTACAGTGACTTCGGCTTTTACTATCCAAAGCACCCATAGCTGCCTGGCAt                >  1:30006581/1‑68 (MQ=255)
                                                    aGCGGGCGTTTCAGCCTCTACAGTGACTTCGGCTTTTACTATCCAAAGCACCCATAGCTGCCTGGCAt                >  1:14924772/1‑68 (MQ=255)
                                                    aGCGGGCGTTTCAGCCTCTACAGTGACTTCGGCTTTTACTATCCAAAGCACCCATAGCTGCCTGGCAt                >  1:23469775/1‑68 (MQ=255)
                                                    aGCGGGCGTTTCAGCCTCTACAGTGACTTCGGCTTTTACTATCCAAAGCACCCATAGCTGCCTGGCAt                >  1:158778/1‑68 (MQ=255)
                                                    aGCGGGCGTTTCAGCCTCTACAGTGACTTCGGCTTTTACTATCCAAAGCACCCATAGCTGCCTGGCAt                >  1:1849167/1‑68 (MQ=255)
                                                       gggCGTTTCAGCCTCTACAGTGACTTCGGCTTTTACTATCCAAAGCACCCATAGCTGCCTGGCATAcc             <  1:12243945/68‑1 (MQ=255)
                                                         gCGTTTCAGTCTCTACAGTGACTTCGGCTTTTACTATCCAAAGCACCCATAGCTGCCTGGCATAccc            <  1:32860596/67‑1 (MQ=255)
                                                         gCGTTTCAGTCTCTACAGTGACTTCGGCTTTTACTATCCAAAGCACCCATAGCTGCCTGGCATAccc            <  1:24612025/67‑1 (MQ=255)
                                                         gCGTTTCAGCCTCTACAGTGACTTCGGCTTTTACTATCCAAAGCACCCATAGCTGCCTGGCATAccc            <  1:35544212/67‑1 (MQ=255)
                                                         gCGTTTCAGCCTCTACAGTGACTTCGGCTTTTACTATCCAAAGCACCCATAGCTGCCTGGCATAccc            <  1:14763460/67‑1 (MQ=255)
                                                         gCGTTTCAGCCTCTACAGTGACTTCGGCTTTTACTATCCAAAGCACCCATAGCTGCCTGGCATAccc            <  1:22019965/67‑1 (MQ=255)
                                                             ttCAGTCTCTACAGTGACTTCGGCTTTTACTATCCAAAGCACCCATAGCTGCCTGGCATACCCaccaa       <  1:29692260/68‑1 (MQ=255)
                                                             ttCAGTCTCTACAGTGACTTCGGCTTTTACTATCCAAAGCACCCATAGCTGCCTGGCATACCCaccaa       <  1:1980854/68‑1 (MQ=255)
                                                             ttCAGCCTCTACAGTGACTTCGGCTTTTACTATCCAAAGCACCCATAGCTGCCTGGCATACCCaccaa       <  1:9375607/68‑1 (MQ=255)
                                                             ttCAGCCTCTACAGTGACTTCGGCTTTTACTATCCAAAGCACCCATAGCTGCCTGGCATACCCaccaa       <  1:10698598/68‑1 (MQ=255)
                                                             ttCAGCCTCTACAGTGACTTCGGCTTTTACTATCCAAAGCACCCATAGCTGCCTGGCATACCCaccaa       <  1:28589842/68‑1 (MQ=255)
                                                             ttCAGCCTCTACAGTGACTTCGGCTTTTACTATCCAAAGCACCCATAGCTGCCTGGCATACCCaccaa       <  1:29040191/68‑1 (MQ=255)
                                                             ttCAGCCTCTACAGTGACTTCGGCTTTTACTATCCAAAGCACCCATAGCTGCCTGGCATACCCaccaa       <  1:33505096/68‑1 (MQ=255)
                                                              tCAGCCTCTACAGTGACTTCGGCTTTTACTATCCAAAGCACCCATAGCTGCCTGGCATACCCaccaac      <  1:28118350/68‑1 (MQ=255)
                                                                  ccTCTACAGTGACTTCGGCTTTTACTATCCAAAGCACCCATAGCTGCCTGGCATACCCACCAACCAGc  >  1:19119459/1‑68 (MQ=255)
                                                                  ccTCTACAGTGACTTCGGCTTTTACTATCCAAAGCACCCATAGCTGCCTGGCATACCCACCAACCAGc  >  1:19613465/1‑68 (MQ=255)
                                                                  ccTCTACAGTGACTTAGGCTTTTACTATCCAAAGCACCCATACCTGCCTGGCATACCCaccaacctg   >  1:35087331/1‑65 (MQ=255)
                                                                  |                                                                   
CGACGGGGAGTGGCTGTTTTTTGAAGAAATAGTGTTGTGTTGATTGTTTAAAAGCGGGCGTTTCAGCCTCTACAGTGACTTCGGCTTTTACTATCCAAAGCACCCATAGCTGCCTGGCATACCCACCAACCAGC  >  CP009273/3447019‑3447152

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: