Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 161,088 T→G 17.9% intergenic (+67/‑129) 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,0880TG17.9% 1410.2 / 63.9 560intergenic (+67/‑129)hrpB/mrcBputative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base T (35/374);  new base G (80/2);  total (183/377)
Fisher's exact test for biased strand distribution p-value = 0.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 807 total aligned reads displayed.

GGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTC  >  CP009273/161022‑161155
                                                                  |                                                                   
  ttggttTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccctca                                                                   <  1:16305182/67‑1 (MQ=255)
  ttggttTTCTCCCTCTCCCTGGGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccctca                                                                   <  1:32456521/67‑1 (MQ=255)
  ttggttTTCCCCCTCCCCCTGGGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccctca                                                                   <  1:19960506/67‑1 (MQ=18)
   tgtttttCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccctcacc                                                                 <  1:27083048/65‑1 (MQ=255)
   tgtttttCTCCCGCGCCGTGGGGTAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccctcacc                                                                 <  1:17533641/65‑1 (MQ=18)
   tggttTTCTCCCTCTCCCTGTGGTAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccctcacc                                                                 <  1:8665346/68‑1 (MQ=255)
   tggttTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccctcacc                                                                 <  1:11210923/68‑1 (MQ=255)
   tggttTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccctcacc                                                                 <  1:3615222/68‑1 (MQ=255)
   tggttTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccctcacc                                                                 <  1:13067398/68‑1 (MQ=255)
   tggttTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccctcacc                                                                 <  1:31642731/68‑1 (MQ=255)
   tggttTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccctcacc                                                                 <  1:44910937/68‑1 (MQ=255)
   tggttTTCTCCCTCTCCCTGGGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccctcacc                                                                 <  1:54025515/68‑1 (MQ=255)
   tggttTTCTCCCTCTCCCGGGGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccctcacc                                                                 <  1:10019478/68‑1 (MQ=255)
   tggttTTCACCCTCTCCCTGTGGGTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccctcacc                                                                 <  1:49577446/68‑1 (MQ=255)
   ggtttttCCCCCTCCCCTGGGGGGAGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTtcaccctcacc                                                                 <  1:33560183/65‑1 (MQ=16)
    ggttTTCGCCCACCCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccctcacc                                                                 <  1:13981906/67‑1 (MQ=18)
    ggttTGCCCCCTGGCCCTGTGGGTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccctcacc                                                                 <  1:26769114/67‑1 (MQ=18)
       tttCTCCCTCTCCCTGTGGGAGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTaa                                                             <  1:43941153/68‑1 (MQ=255)
        ttCTCCCTCCCCGTGGGGTAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAAc                                                            <  1:40989725/68‑1 (MQ=21)
        ttCTCCCTCCCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAAc                                                            <  1:19407217/68‑1 (MQ=255)
         tctccctctcccTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAAcc                                                           <  1:28776606/68‑1 (MQ=255)
          ccccctccccgtGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAAccc                                                          <  1:55479886/66‑1 (MQ=21)
           tccctctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCCCCCGCACCCTAAacc                                                          >  1:57199959/1‑67 (MQ=17)
           accctctccctGTGGGTGTGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCTCCCTCACCCTAACCct                                                         <  1:34131737/67‑1 (MQ=21)
            ccctctccctGTGGGAGAGGGCCGGGGGGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCctc                                                        <  1:41742147/68‑1 (MQ=255)
              ctctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCtctc                                                      <  1:27658041/68‑1 (MQ=255)
              ctctccctGGGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCtctc                                                      <  1:40398047/68‑1 (MQ=255)
              ctcccccgGGGGGAGAGGGCCGGGGGGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCtctc                                                      <  1:13715738/68‑1 (MQ=25)
               tctccctGGGGGTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTcc                                                     <  1:44688655/68‑1 (MQ=25)
               tccccctGCGGGCGCGGGCCGGGCGGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTcc                                                     <  1:50848033/65‑1 (MQ=25)
               ccgccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTcc                                                     <  1:32956349/65‑1 (MQ=255)
                ctccctGTGGGATAGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTccc                                                    <  1:3942063/68‑1 (MQ=25)
                ctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTccc                                                    <  1:25752778/68‑1 (MQ=255)
                ctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTccc                                                    <  1:9188567/68‑1 (MQ=255)
                ctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTccc                                                    <  1:27970052/68‑1 (MQ=255)
                ctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTccc                                                    <  1:22003097/68‑1 (MQ=255)
                ctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTccc                                                    <  1:34972372/68‑1 (MQ=255)
                ctccctGGGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTccc                                                    <  1:57662142/68‑1 (MQ=255)
                 tccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCCCCCGCCCCCGAACCCTCTccc                                                    >  1:51390398/1‑67 (MQ=18)
                 tccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTccc                                                    <  1:14327416/67‑1 (MQ=255)
                 tccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCAACCCCACCCGAACCCTCGccc                                                    >  1:39636327/1‑67 (MQ=17)
                 cccctGTGGGCTTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTccc                                                    <  1:45392980/66‑1 (MQ=37)
                  ccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTc                                                  <  1:23974117/68‑1 (MQ=255)
                  ccctGGGGGTGTGGGCCGGGGGGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTc                                                  <  1:5194452/68‑1 (MQ=37)
                  ccctGGGGGTGCGGGCCGGGGGGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTc                                                  <  1:31066132/68‑1 (MQ=37)
                  ccctGGGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTc                                                  <  1:7527665/68‑1 (MQ=255)
                  ccctGGGGGAGCGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTc                                                  <  1:49255451/68‑1 (MQ=37)
                  ccctGGGGGAGAGGGCCGGGGGGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTc                                                  <  1:52667300/68‑1 (MQ=255)
                   cctGTGGGTTAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCa                                                 <  1:22997381/68‑1 (MQ=255)
                   cctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGACCCGCACCCAAACCCCCGcccgcc                                                 >  1:47006186/1‑65 (MQ=17)
                   cctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTc                                                  <  1:30457075/67‑1 (MQ=255)
                   cctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCa                                                 <  1:54461886/68‑1 (MQ=255)
                     tgtgGGAGCGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                               <  1:18057716/68‑1 (MQ=255)
                     tgtgGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                               <  1:38114346/68‑1 (MQ=255)
                     tgtgGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCGCACCCAAACCCACTCCCGCAAg                                               >  1:42998927/1‑68 (MQ=21)
                     tgggggAGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                               <  1:40577203/65‑1 (MQ=255)
                        gggAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAggg                                             <  1:12067833/67‑1 (MQ=255)
                         ggAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGCCCCCCCCCCCAAACCACACCCCCAAggggg                                           >  1:17029432/1‑66 (MQ=9)
                         ggAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGCCCCCCACCCCAACCCCCACCCCCAAcgggg                                           >  1:50075674/1‑66 (MQ=12)
                         ggAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGACCCACACCCAAAACCACACCCCCAAagggg                                           >  1:48400868/1‑66 (MQ=12)
                          gagaGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgag                                           <  1:10829547/67‑1 (MQ=255)
                          gagaGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgag                                           <  1:36788158/67‑1 (MQ=255)
                          gagaGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCGCACCCGAACCCTCACCCGCAAgggggg                                          >  1:37430021/1‑65 (MQ=21)
                           agagGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCAGACCCAAAACCACTCCACCAAGGgatgg                                         >  1:4777478/1‑65 (MQ=18)
                           agagGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACCCCCGAACCCGCTCCCACAATGgtgac                                         >  1:22465269/1‑67 (MQ=18)
                           agagGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACCCCCAAACCCACACCCGCAACGgagcg                                         >  1:7032333/1‑66 (MQ=18)
                             agGGCCGGGGTGAGGGCATCAGCGCGCACGTTGCCCCCCACCCAGACCCACTCCCGGCAGGGAGCggg                                       >  1:42488505/1‑68 (MQ=12)
                             agGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCGCCCCCAAACCCTGTCCCGGAAGGGGGCggg                                       >  1:30054769/1‑68 (MQ=18)
                             agGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCAAAACCAATACCCAAAGGGAGAggg                                       >  1:52335157/1‑68 (MQ=18)
                              gggCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCCCACCCTAACCCTCTCCCTCAAGGGAGAggg                                       <  1:56365414/67‑1 (MQ=38)
                              gggCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAggg                                       <  1:56788521/67‑1 (MQ=255)
                                gcgggggggAGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAgggg                                      <  1:43297937/63‑1 (MQ=25)
                                gCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGa                                     <  1:41537988/67‑1 (MQ=255)
                                gCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGa                                     <  1:47507478/67‑1 (MQ=255)
                                gCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGa                                     <  1:15589007/67‑1 (MQ=255)
                                gCCGGGGGGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGa                                     <  1:6613756/67‑1 (MQ=255)
                                     ggTGAGGGCATCAGCGCGCACGTTGACCCACACCCAAAACCACACCCCCAAAGGGGAGGGGGCCGCTc                               >  1:29292751/1‑68 (MQ=12)
                                     ggTGAGGGCATCAGCGCGCACGTTCACCCTCCCCCTAACCCTCTCCCGCAAGGGAGCGGGGGCCGGTc                               >  1:26200270/1‑68 (MQ=37)
                                     ggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAgg                                              >  1:25509557/1‑53 (MQ=255)
                                     ggTGAGGGCATCAGCGCGCACGTTCACCCACCCCCGAACCCCCTACCACAAGGGGGGGGGGGCCgctg                               >  1:35148657/1‑65 (MQ=16)
                                          ggccATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGca                          <  1:20539197/65‑1 (MQ=255)
                                           ggCATCAGCGCGCACGTTCACCCGCACCCTAACCCACTCCCTCAAGGGAGAGGGGGCCGCTCGAGccc                         >  1:29785433/1‑66 (MQ=25)
                                            gCATCAGCGCGCACGTTCACCCTCACCCTAACCCACTCCCTCAAGGGAGAGGGGCCCGATCGAGcaca                        >  1:11665968/1‑68 (MQ=255)
                                            gCATCAGCGCGCACGTTCACCCGCCCCCGAACCCTCTACCGCAAGGGAGAGGGGGCCGCTAGCGcaca                        >  1:23340641/1‑68 (MQ=255)
                                              aTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACaaa                      <  1:40027805/68‑1 (MQ=255)
                                              aTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACaaa                      <  1:53402233/68‑1 (MQ=255)
                                                cAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtt                    <  1:14813030/68‑1 (MQ=255)
                                                 aGCGCGCACGTTCACCCTCACCCAAACCCACACCCACAACGGGGAGGGGCCCGCTCGAGCACAAAatt                   >  1:38913739/1‑68 (MQ=255)
                                                        aCGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgagaga            <  1:54858341/68‑1 (MQ=255)
                                                        aCGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgagaga            <  1:18716367/68‑1 (MQ=255)
                                                         cGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGat           <  1:21550493/68‑1 (MQ=255)
                                                            tcaccctcacccTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGatat         <  1:24673378/67‑1 (MQ=255)
                                                            tcaccctcacccTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGatat         <  1:25173173/67‑1 (MQ=255)
                                                            tcaccctcacccTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGatat         <  1:46080425/67‑1 (MQ=255)
                                                            tcaccctcacccTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGatat         <  1:3771598/67‑1 (MQ=255)
                                                             caccctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtct       <  1:15142206/68‑1 (MQ=255)
                                                              accctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtctt      <  1:28419539/68‑1 (MQ=255)
                                                               ccctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttc     <  1:12498250/68‑1 (MQ=255)
                                                                cctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttct    <  1:38479746/68‑1 (MQ=255)
                                                                actcaccatAACCCTCTCCCGCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttct    <  1:20992604/67‑1 (MQ=38)
                                                                  |                                                                   
GGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTC  >  CP009273/161022‑161155

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: