Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 161,106 T→G 10.0% 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,1060TG10.0% 1468.5 / 19.4 476intergenic (+85/‑111)hrpB/mrcBputative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base T (176/234);  new base G (46/2);  total (240/236)
Fisher's exact test for biased strand distribution p-value = 1.54e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 479 total aligned reads displayed.

CCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGA  >  CP009273/161040‑161172
                                                                  |                                                                  
ccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTc                                                                   <  1:792881/68‑1 (MQ=255)
 cgtgtgGGAGTGGGCCGGGGGGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCa                                                                  <  1:28293119/66‑1 (MQ=37)
 cctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCa                                                                  <  1:33634880/68‑1 (MQ=255)
 cctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCa                                                                  <  1:27478921/68‑1 (MQ=255)
 cctGTGGGAGAGGGCCGGGGGGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCa                                                                  <  1:7764614/68‑1 (MQ=255)
 cctGGGGGAGCGGGCCGGGGGGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCa                                                                  <  1:6337647/68‑1 (MQ=37)
   tgtgGGAGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                <  1:28076604/68‑1 (MQ=255)
   tggggggGAGGGGGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                <  1:8123813/65‑1 (MQ=25)
     tgGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGCCCCGCCCCCTGACCCTCTCCCGCAAcgg                                                              >  1:33680406/1‑68 (MQ=18)
     tgGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGACCCACACCCGAACCCACGCCCGCAAggg                                                              >  1:5774050/1‑68 (MQ=18)
     tgGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCCCCCTCCCCCTAAACCGCTCCCGCAAggg                                                              >  1:25863346/1‑68 (MQ=21)
     tgGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCGCCCCCGAACCCTCTCCCGCAAggg                                                              >  1:24854600/1‑68 (MQ=21)
      ggtagagGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAggg                                                              <  1:25191798/64‑1 (MQ=255)
      gggCGAGGGCCGGGGGGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAggg                                                              <  1:22175152/67‑1 (MQ=255)
      gggAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAggg                                                              <  1:26495480/67‑1 (MQ=255)
      gggAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAggg                                                              <  1:6792916/67‑1 (MQ=255)
       ggAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGACCCGCCCCCAAAACCACTCCCGCAAggggg                                                            >  1:17316043/1‑66 (MQ=17)
       ggAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGAACCACACACAAAACCACAACCTCAAgggtg                                                            >  1:23757799/1‑66 (MQ=9)
       ggAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCCCCCGCACCCTAACCCAAGCCCGCAAgggg                                                             >  1:4394485/1‑66 (MQ=18)
       ggAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCGCAAGGgag                                                            >  1:6179512/1‑68 (MQ=255)
       ggAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCCCACCCGAACCCACACCCGCAAGGgag                                                            >  1:13753305/1‑68 (MQ=21)
       ggAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCTAACCCTCTCCCGCCAGGgag                                                            >  1:6011363/1‑68 (MQ=25)
       ggAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCCCCCGCCCCCGAACCCCCGCCCCCCAggggg                                                            >  1:13149150/1‑66 (MQ=9)
        gatagGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgag                                                            <  1:12361746/64‑1 (MQ=255)
        gagcGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgag                                                            <  1:28868182/67‑1 (MQ=255)
        gagcGGGCCGGGGGGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgag                                                            <  1:6551491/67‑1 (MQ=38)
        gagaGGGCCGTTGAGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCACCCTCAAGGgag                                                            <  1:4678661/67‑1 (MQ=37)
        gagaGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgag                                                            <  1:14836327/67‑1 (MQ=255)
        gagaGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgag                                                            <  1:19231742/67‑1 (MQ=255)
         agagGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCAACCAAACCCTCTCCCGCCAGGGagag                                                          >  1:9640853/1‑68 (MQ=25)
         agagGGCCGGGGTGAGGGCATCAGCGCGCACGTTAACCCTAACCCAAAACCACTCCCGCAAGGgagcg                                                          >  1:9202360/1‑66 (MQ=18)
           agGGCCGGGGTGAGGGCATCAGCGCGCACGTTGACCCGCACCCAAACCCACGCCCGCAACGGGGAggg                                                        >  1:5324168/1‑68 (MQ=18)
           agGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCGCCCCCGAACCCACGCCCGCAAGGGGGAggg                                                        >  1:30358022/1‑68 (MQ=21)
            gggCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAggg                                                        <  1:747740/67‑1 (MQ=255)
            gggCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAggg                                                        <  1:3632575/67‑1 (MQ=255)
            gggCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAggg                                                        <  1:32245984/67‑1 (MQ=255)
            gggCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAggg                                                        <  1:12910466/67‑1 (MQ=255)
            gggCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAggg                                                        <  1:15886883/67‑1 (MQ=255)
            gggCCGGGGGGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAggg                                                        <  1:12042775/67‑1 (MQ=255)
            gggCCGGGGGGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAggg                                                        <  1:14337027/67‑1 (MQ=255)
              gCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGa                                                      >  1:16767540/1‑67 (MQ=255)
                   ggTGAGGGCATCAGCGCGCACGTTGACCCCCCCCCAAACCCCCTCCCGCCAGGGGGGGGGGGCCGCTc                                                >  1:15508467/1‑68 (MQ=12)
                   ggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGAt                                                 <  1:25595864/67‑1 (MQ=255)
                   ggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAAACCTCTCCCTCAAGGGAGAGGGGCCCGATc                                                >  1:22590170/1‑68 (MQ=255)
                   ggTGAGGGCATCAGCGCGCACGTTCACCCGCCCCCGAACCCACGCCCGCAAGGGCGGGGGGGCCgggg                                                >  1:1810975/1‑65 (MQ=16)
                        gggCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGc                                            <  1:32313277/67‑1 (MQ=255)
                          gCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGca                                           <  1:33911461/66‑1 (MQ=255)
                            aTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACaaa                                       <  1:2961654/68‑1 (MQ=255)
                             tCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAt                                      <  1:32972766/68‑1 (MQ=255)
                                      aCGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGCTCGAGCACAAATTTTgagaga                             <  1:14584826/68‑1 (MQ=255)
                                      aCGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgagaga                             <  1:16008724/68‑1 (MQ=255)
                                      aCGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgagaga                             <  1:24590788/68‑1 (MQ=255)
                                       cGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGat                            <  1:26702804/68‑1 (MQ=255)
                                       cGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGat                            <  1:28625939/68‑1 (MQ=255)
                                          tcaccctcacccTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGatat                          <  1:9550902/67‑1 (MQ=255)
                                          tcaccctcacccTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGatat                          <  1:16408395/67‑1 (MQ=255)
                                          tcaccctcacccTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGatat                          <  1:31125280/67‑1 (MQ=255)
                                          tcaccctcacccTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGatat                          <  1:13448620/67‑1 (MQ=255)
                                          tcaccctcacccTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGatat                          <  1:27903315/67‑1 (MQ=255)
                                          tcaccctcacccTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGatat                          <  1:21425157/67‑1 (MQ=255)
                                          tcaccctcacccTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGatat                          <  1:8949598/67‑1 (MQ=255)
                                             ccctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttc                      <  1:11145606/68‑1 (MQ=255)
                                              cctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttct                     <  1:4953946/68‑1 (MQ=255)
                                              cctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttct                     <  1:24128949/68‑1 (MQ=255)
                                              cctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttct                     <  1:2337036/68‑1 (MQ=255)
                                                   ccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGt                <  1:20973790/68‑1 (MQ=255)
                                                    cctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttct                     >  1:29753685/1‑62 (MQ=255)
                                                    cctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttct                     >  1:34055747/1‑62 (MQ=255)
                                                    cctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTc               <  1:9855852/68‑1 (MQ=255)
                                                    cctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTc               <  1:17123201/68‑1 (MQ=255)
                                                    cctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTc               <  1:22817051/68‑1 (MQ=255)
                                                    cctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTc               <  1:25281051/68‑1 (MQ=255)
                                                       aaCCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTg            <  1:16212727/68‑1 (MQ=255)
                                                        aCCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGt           >  1:3262501/1‑68 (MQ=255)
                                                        aCCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGt           >  1:20477390/1‑68 (MQ=255)
                                                        aCCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGt           >  1:27099694/1‑68 (MQ=255)
                                                         cccTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTa          >  1:10052077/1‑68 (MQ=255)
                                                         cccTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTa          <  1:8591541/68‑1 (MQ=255)
                                                          ccTCTCCCTCAAGGGAGAGGGGCCCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTa          >  1:11465463/1‑67 (MQ=255)
                                                          ccTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTa          >  1:31350763/1‑67 (MQ=255)
                                                          ccTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTa          >  1:26254279/1‑67 (MQ=255)
                                                          ccTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTa          >  1:10420774/1‑67 (MQ=255)
                                                          ccTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTa          >  1:10947818/1‑67 (MQ=255)
                                                          ccTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTa          >  1:943201/1‑67 (MQ=255)
                                                          ccTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTa          >  1:11639603/1‑67 (MQ=255)
                                                          ccTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTa          >  1:34288943/1‑67 (MQ=255)
                                                          ccTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTa          >  1:14015257/1‑67 (MQ=255)
                                                          ccTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTa          >  1:31772016/1‑67 (MQ=255)
                                                          ccTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTa          >  1:2449019/1‑67 (MQ=255)
                                                          ccTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTa          >  1:17791675/1‑67 (MQ=255)
                                                          ccTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTa          >  1:18320301/1‑67 (MQ=255)
                                                          ccTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTa          >  1:1972643/1‑67 (MQ=255)
                                                          ccTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTa          >  1:20675122/1‑67 (MQ=255)
                                                          ccTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTa          >  1:21702787/1‑67 (MQ=255)
                                                          ccTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTa          >  1:29182628/1‑67 (MQ=255)
                                                            tctcCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACa       <  1:33185615/68‑1 (MQ=255)
                                                              tcCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACaga     <  1:23579760/68‑1 (MQ=255)
                                                               cccTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaa    >  1:22291726/1‑68 (MQ=255)
                                                                ccTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaag   >  1:18713670/1‑68 (MQ=255)
                                                                  tCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTc               <  1:12553965/54‑1 (MQ=255)
                                                                  |                                                                  
CCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGA  >  CP009273/161040‑161172

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: