Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 161,106 T→G 9.6% intergenic (+85/‑111) hrpB → / → mrcB putative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273161,1060TG 9.6% 2746.3 / 38.7 888intergenic (+85/‑111)hrpB/mrcBputative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base T (318/451);  new base G (76/3);  total (434/454)
Fisher's exact test for biased strand distribution p-value = 1.14e-23
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 897 total aligned reads displayed.

TCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAA  >  CP009273/161039‑161173
                                                                    |                                                                   
 ccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTc                                                                    <  1:23974117/68‑1 (MQ=255)
  cctGTGGGTTAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCa                                                                   <  1:22997381/68‑1 (MQ=255)
  cctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTc                                                                    <  1:30457075/67‑1 (MQ=255)
  cctGGGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCa                                                                   <  1:45711963/68‑1 (MQ=255)
   ctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                                                                  <  1:50130280/68‑1 (MQ=255)
   ctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCACTCCCGCaa                                                                  >  1:23383306/1‑68 (MQ=255)
    tgtgGGAGCGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                 <  1:18057716/68‑1 (MQ=255)
    tgtgGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGACCCACCCCCAAACCCACACCCGCAAg                                                                 >  1:28501217/1‑68 (MQ=17)
    tgtgGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                 <  1:56842644/68‑1 (MQ=255)
    tgtgGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                 <  1:55055499/68‑1 (MQ=255)
    tgggggAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                 <  1:54272078/65‑1 (MQ=255)
      tgGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGCCCCGCACCCAAACCCACAACCCCAAggg                                                               >  1:75104/1‑68 (MQ=17)
      tgGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGACCCACACCCGAACCCTCGCCCGCAAggg                                                               >  1:52630681/1‑68 (MQ=21)
      tgGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCCCCCGCCCCCGAACCCTCTCCCTCAAggg                                                               >  1:52400452/1‑68 (MQ=21)
      tgGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCCCCCAAACCCAAACCCACTCCCACAAggg                                                               >  1:27921556/1‑68 (MQ=18)
        ggAGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGga                                                              <  1:8706557/67‑1 (MQ=255)
        ggAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGACCCACACCCAAAACCACACCCCCAAagggg                                                             >  1:48400868/1‑66 (MQ=12)
        ggAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCGCACCCTAAACCGCTCCCGCTAggggg                                                             >  1:34803466/1‑66 (MQ=18)
        ggAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCCCACCCAAACCCACTCCCGCAAtggc                                                              >  1:33012007/1‑66 (MQ=21)
        ggAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCTAACCCCCTACCCCAAGGgag                                                             >  1:12050096/1‑68 (MQ=25)
        ggAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTAACCCACACCCTAAACCACGCCCCCAAggggg                                                             >  1:4285384/1‑66 (MQ=18)
         gagtGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgag                                                             <  1:29465089/67‑1 (MQ=255)
         gagtGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgag                                                             <  1:24720527/67‑1 (MQ=255)
         gagaGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgag                                                             <  1:12592434/67‑1 (MQ=255)
         gagaGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgag                                                             <  1:10894426/67‑1 (MQ=255)
         gagaGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgag                                                             <  1:36550956/67‑1 (MQ=255)
         gagaGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgag                                                             <  1:6077402/67‑1 (MQ=255)
          agagGGGCGGGGTGAGGGCATCTGCGCGCACGTTCACCCTCACCCTAACCCTCTCCGTCAAGGgaga                                                            <  1:32247147/67‑1 (MQ=37)
          agagGGCCGGGGTGAGGGCATCAGCGCGCACGTTGCACCCGCCCCAAAACCCCAGCCCGCAAgggggg                                                            >  1:6748913/1‑65 (MQ=9)
          agagGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCGCACCCGAACCCACTCCCGCAAGGGagag                                                           >  1:29979695/1‑68 (MQ=25)
            agGGCCGGGGTGAGGGCATCAGCGCGCACGTTGACCCGCGCCCGAACCCGCGCCCGCAAGgggggggg                                                         >  1:35129810/1‑68 (MQ=12)
            agGGCCGGGGTGAGGGCATCAGCGCGCACGTTCCCCCGCACCCTAAACCAATCCCGCAAGGGAGAggg                                                         >  1:27600976/1‑68 (MQ=21)
             gggCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAggg                                                         <  1:13227387/67‑1 (MQ=255)
             gggCCGGGGGGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAggg                                                         <  1:37001212/67‑1 (MQ=255)
             gggCCGGGGGGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAgg                                                          <  1:3868330/66‑1 (MQ=255)
              ggCGGGGGTGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAgggg                                                        <  1:33934627/67‑1 (MQ=37)
               gcgggggTGAGGGCATCAGCGCGCACGTTCACGCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGa                                                       <  1:11344035/64‑1 (MQ=255)
               gCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGa                                                       <  1:40438913/67‑1 (MQ=255)
               gCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGa                                                       <  1:47507478/67‑1 (MQ=255)
               gCCGGGGGGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGa                                                       <  1:43158843/67‑1 (MQ=255)
                    ggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGAt                                                  <  1:43643603/67‑1 (MQ=255)
                    ggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGAt                                                  <  1:49402408/67‑1 (MQ=255)
                     gTGAGGGCATCAGCGCGCACGTTCCCCCGCACCCGAACCCACACCCCCCACGGAGAGGGGGCCGCTCg                                                >  1:53103886/1‑68 (MQ=16)
                     gTGAGGGCATCAGCGCGCACGTTCCCCCCCCCCCTAAACCACTCCCGCAAGGGGGCGGGGGCCGATc                                                 >  1:53561441/1‑67 (MQ=19)
                     gTGAGGGCATCAGCGCGCACGTTCCCCCCCCCCCGAACCCGCTCCCACAACGGCGAGGGGGCCGCTc                                                 >  1:37936469/1‑67 (MQ=16)
                       gcgggCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAg                                              <  1:17492242/66‑1 (MQ=255)
                       gAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAg                                              <  1:16495240/68‑1 (MQ=255)
                           gCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGTGAGGGGACCGCTCGAGcaca                                          >  1:54622747/1‑68 (MQ=255)
                             aTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACaaa                                        <  1:37539048/68‑1 (MQ=255)
                             aTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACaaa                                        <  1:55784088/68‑1 (MQ=255)
                                   gcgcACGTTGACCCACACCCCAAACCACACCCGCAAGCGAGAGGGGACCGCGCGCGCACAAAatttta                                  >  1:41325415/1‑66 (MQ=255)
                                   gcgcACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                  <  1:33473859/68‑1 (MQ=255)
                                     gcACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGTCCGATCGAGCACAAATTTTgag                                 <  1:42285499/67‑1 (MQ=255)
                                     gcACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgaga                                <  1:28977552/68‑1 (MQ=255)
                                     gcACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgaga                                <  1:48968108/68‑1 (MQ=255)
                                       aCGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgagaga                              <  1:18716367/68‑1 (MQ=255)
                                       aCGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgagaga                              <  1:14184884/68‑1 (MQ=255)
                                           tcaccctcacccTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGatat                           <  1:4511925/67‑1 (MQ=255)
                                             caccctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtct                         <  1:9370588/68‑1 (MQ=255)
                                              accctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtctt                        <  1:10109390/68‑1 (MQ=255)
                                              accctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtctt                        <  1:15474004/68‑1 (MQ=255)
                                               ccctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttc                       <  1:40159165/68‑1 (MQ=255)
                                                cctcaccctAACCGTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttct                      <  1:52026986/68‑1 (MQ=255)
                                                cctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttct                      <  1:39844982/68‑1 (MQ=255)
                                                cctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttct                      <  1:47976479/68‑1 (MQ=255)
                                                cctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttct                      <  1:46571826/68‑1 (MQ=255)
                                                cctcaacctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttct                      <  1:56385519/68‑1 (MQ=255)
                                                  tcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttcttc                    <  1:21058679/68‑1 (MQ=255)
                                                     ccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGt                 <  1:17210454/68‑1 (MQ=255)
                                                      cctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGatat                           <  1:38981234/57‑1 (MQ=255)
                                                      cctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttct                      >  1:719808/1‑62 (MQ=255)
                                                      cctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttct                      >  1:15158382/1‑62 (MQ=255)
                                                      cctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttct                      >  1:39378177/1‑62 (MQ=255)
                                                      cctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttct                      >  1:20225503/1‑62 (MQ=255)
                                                      cctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTc                <  1:31664028/68‑1 (MQ=255)
                                                      cctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTc                <  1:13823811/68‑1 (MQ=255)
                                                      cctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTc                <  1:25610565/68‑1 (MQ=255)
                                                         aaCCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTg             <  1:22425415/68‑1 (MQ=255)
                                                          aCCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTga            >  1:44613198/1‑67 (MQ=255)
                                                          aCCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGt            >  1:10486859/1‑68 (MQ=255)
                                                          aCCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGt            >  1:2161483/1‑68 (MQ=255)
                                                          aCCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGt            >  1:14647758/1‑68 (MQ=255)
                                                          aCCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGt            >  1:19247220/1‑68 (MQ=255)
                                                          aCCCTCTCCCTCAAGGGAGAGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTa           <  1:27156082/68‑1 (MQ=255)
                                                            ccTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTaa          >  1:57431022/1‑68 (MQ=255)
                                                            ccTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTa           >  1:25203918/1‑67 (MQ=255)
                                                            ccTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTa           >  1:8052173/1‑67 (MQ=255)
                                                            ccTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTa           >  1:36110956/1‑67 (MQ=255)
                                                            ccTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTa           >  1:31235629/1‑67 (MQ=255)
                                                            ccTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTa           >  1:38333656/1‑67 (MQ=255)
                                                            ccTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTa           >  1:46954880/1‑67 (MQ=255)
                                                            ccTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTa           >  1:26535534/1‑67 (MQ=255)
                                                            ccTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTa           >  1:26087871/1‑67 (MQ=255)
                                                            ccTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTa           >  1:41816707/1‑67 (MQ=255)
                                                            ccTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTa           >  1:51633643/1‑67 (MQ=255)
                                                            ccTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTa           >  1:50744448/1‑67 (MQ=255)
                                                              tctcCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACa        <  1:20673081/68‑1 (MQ=255)
                                                                tcCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACaga      >  1:19942280/1‑68 (MQ=255)
                                                                  ccTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaag    >  1:15817505/1‑68 (MQ=255)
                                                                  ccTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaag    >  1:40857483/1‑68 (MQ=255)
                                                                    |                                                                   
TCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAA  >  CP009273/161039‑161173

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: