Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,165,415 Δ1 bp 9.1% coding (2513/3906 nt) hrpA ← ATP‑dependent helicase HrpA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,165,4150T. 9.1% 161.5 / 10.4 55coding (2513/3906 nt)hrpAATP‑dependent helicase HrpA
Reads supporting (aligned to +/- strand):  ref base T (23/26);  new base . (3/2);  total (27/28)
Fisher's exact test for biased strand distribution p-value = 6.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

CCGGCAGGCTCGGCAACAGCGGTGCCGGCACCCGCACCGTCACACCATCACGGGGGTGGCCGGGCTCGAAGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTCCGCCGGCAGGCGCGCTTCATAGAAAGCGTACAGCGTTTCTTC  >  NC_002947/5165269‑5165557
                                                                                                                                                   |                                                                                                                                               
ccGGCAGGCTCGGCAACAGCGGTGCCGGCACCCGCACCGTCACACCATCACGGGGGTGGCCGGGCTCGAAGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTcgc                                                                                                                                              >  1:873560/1‑150 (MQ=255)
ccGGCAGGCTCGGCAACAGCGGTGCCGGCACCCGCACCGTCACACCATCACGGGGGTGGCCGGGCTCGAAGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTcgc                                                                                                                                              >  1:663865/1‑150 (MQ=255)
           ggCAACAGCGGTGCCGGCACCCGCACCGTCACACCATCACGGGGGTGGCCGGGCTCGAAGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCCGCGAGGCCAGcgaa                                                                                                                                  <  2:10587/150‑4 (MQ=255)
            gCAACAGCGGTGCCGGCACCCGCACCGTCACACCATCACGGGGGTGGCCGGGCTCGAAGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGt                                                                                                                                 >  2:1400966/1‑150 (MQ=255)
             cAACAGCGGTGCCGGCACCCGCACCGTCACACCATCACGGGGGTGGCCGGGCTCGAAGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGtc                                                                                                                                >  1:282952/1‑150 (MQ=255)
             cAACAGCGGTGCCGGCACCCGCACCGTCACACCATCACGGGGGTGGCCGGGCTCGAAGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCCACCCCCGCGCTCGAGCACGtc                                                                                                                                >  2:748162/1‑150 (MQ=255)
              aaCAGCGGTGCCGGCACCCGCACCGTCACACCATCACGGGGGTGGCCGGGCTCGAAGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCAGTCACTTCGCTGGCCTCGCGGGCCAGCACGtcc                                                                                                                               <  2:217154/150‑1 (MQ=255)
                     gTGCCGGCACCCGCACCGTCACACCATCACGGGGGTGGCCGGGCTCGAAGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTcgc                                                                                                                        <  1:982581/150‑1 (MQ=255)
                         cGGCACCCGCACCGTCACACCATCACGGGGGTGGCCGCGCTCGAAGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTCCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCAt                                                                                                                    <  2:123041/150‑1 (MQ=255)
                               ccGCACCGTCACACCATCACGGGGGTGGCCGGGCTCGAAGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATcag                                                                                                              >  2:1417480/1‑150 (MQ=255)
                                ctcACCGTCACACCATCACGGGGGTGGCCGGGCTCGAAGTTGTCGCTCAGCTGCAGTCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCATCGCGGGCCAGCACGTCCTCCTCGCGCATGATcagc                                                                                                             <  2:1156385/148‑1 (MQ=255)
                                   aCCGTCACACCATCACGGGGGTGGCCGGGCTCGAAGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGAAGAGCACACGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATcagcaag                                                                                                          >  1:1156385/1‑148 (MQ=255)
                                   aCCGTCACACCATCACGGGGGTGGCCGGGCTCGAAGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAgg                                                                                                          <  1:1400966/150‑1 (MQ=255)
                                   aCCGTCACACCATCACGGGGGTGGCCGGGCTCGAAGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCCGCGAGGCCAGcgaa                                                                                                                                  >  2:966529/1‑123 (MQ=255)
                                   aCCGTCACACCATCACGGGGGTGGCCGGGCTCGAAGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCCGCGAGGCCAGcgaa                                                                                                                                  <  1:966529/126‑4 (MQ=255)
                                      gTCACACCATCACGGGGGTGGCCGGGCTCGAAGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTg                                                                                                       <  2:663865/150‑1 (MQ=255)
                                         acacCATCAGGGGGGTGGCCGGGCTCGAAGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGtc                                                                                                                                <  2:284753/122‑1 (MQ=255)
                                         acacCATCACGGGGGTGGCCGGGCTCGAAGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGtc                                                                                                                                >  1:284753/1‑122 (MQ=255)
                                         acacCATCACGGGGGTGGCCGGGCTCGAAGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCg                                                                                                    <  1:979615/150‑1 (MQ=255)
                                           acCATCACGGGGGTGGCCGGGCTCGAAGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTc                                                                                                  >  2:934253/1‑150 (MQ=255)
                                             cATCACGGGGGTGGCCGGGCTCGAAGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCt                                                                                                <  1:1059491/150‑1 (MQ=255)
                                                cACGGGGGTGGCCGGGCTCGAAGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCtttt                                                                                             <  1:376335/150‑1 (MQ=255)
                                                 aCGGGGGTGGCCGGGCTCGAAGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTg                                                                                            >  2:1346482/1‑150 (MQ=255)
                                                                   gAAGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCa                                                                          >  2:968519/1‑150 (MQ=255)
                                                                          tAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCa                                                                   <  1:897258/150‑1 (MQ=255)
                                                                                aGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTc                                                             <  1:1305129/150‑1 (MQ=255)
                                                                                   ggCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGcc                                                          <  1:1346482/150‑1 (MQ=255)
                                                                                    gCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCg                                                         <  1:1417480/150‑1 (MQ=255)
                                                                                      aGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTc                                                       >  1:120296/1‑150 (MQ=255)
                                                                                          gcAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGc                                                                                                                                     <  2:755939/68‑1 (MQ=255)
                                                                                          gcAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGc                                                                                                                                     >  1:755939/1‑68 (MQ=255)
                                                                                          gcAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGACAGCACGTCCTCCTCGCGCCTGATAAGCAGGTTGGCGTCCTTTTGGCTTCCCATGCGGTACCAACTATCAAAGGTCGCCGTCtggt                                                   <  2:873560/150‑1 (MQ=255)
                                                                                          gcAGATCCCCCACCTGTAGACTGTCGGTGTACTGAGCGGCGGTCACTGCGCTGGCCTCGCGGTCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTTGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCtggt                                                   <  2:629663/150‑1 (MQ=255)
                                                                                           cAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCTTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCtggtg                                                  <  1:968519/150‑1 (MQ=255)
                                                                                             gATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGt                                                                                                         >  1:567771/1‑93 (MQ=255)
                                                                                             gATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGt                                                                                                         <  2:567771/93‑1 (MQ=255)
                                                                                                   ccACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGcc                                                          <  1:474756/134‑1 (MQ=255)
                                                                                                   ccACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGcc                                                          >  2:474756/1‑134 (MQ=255)
                                                                                                    cACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTccgc                                         <  1:618590/150‑1 (MQ=255)
                                                                                                     aCCTGCATACTGTCGGGGTACTGTGCGGCGGTCACGTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTccgcc                                        >  1:824948/1‑150 (MQ=255)
                                                                                                           aTACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCCGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTCCGCCGGCAgg                                  >  2:657000/1‑150 (MQ=255)
                                                                                                            tACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTCCGCCGGCAGgc                                 >  1:782166/1‑150 (MQ=255)
                                                                                                               tGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCaa                                                                         <  1:561705/107‑1 (MQ=255)
                                                                                                               tGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCaa                                                                         >  2:561705/1‑107 (MQ=255)
                                                                                                               tGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTCCGCCGGCAGgcgcg                              >  1:1031954/1‑150 (MQ=255)
                                                                                                                  cGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTGCCATCTATCAAAGGTCGCCGTCTGGTGGATTTCCGCCGGCAGGCGCGCtt                           <  1:748162/150‑1 (MQ=255)
                                                                                                                    gggTACTGTGCGGCGGTCACTTCGCTGGCCTTCAGCGTCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTCCGCCGGCAGGCGCGCTTCa                         <  2:824948/150‑1 (MQ=255)
                                                                                                                                        cTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACGCACTATCAAAGGTCGCCGTCTGGTGGATTTCCGCCGGCAGGCGTGCTTCATAGAAAGCGTACAGCGttt      >  1:1022468/1‑150 (MQ=255)
                                                                                                                                         ttCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTCCGCCGGCAGGCGCGCTTCATAGAAAGCGTACa            <  1:1311323/143‑1 (MQ=255)
                                                                                                                                         ttCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTCCGCCGGCAGGCGCGCTTCATAGAAAGCGTACa            >  2:1311323/1‑143 (MQ=255)
                                                                                                                                         ttCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTCCGCCGGCAGGCGCGCTTCATAGAAAGCGTACAGCGTttc     >  1:829699/1‑150 (MQ=255)
                                                                                                                                         ttCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTCCGCCGGCAGGCGCGCTTCATAGAAAGCGTACAGCGTttc     >  2:1382190/1‑150 (MQ=255)
                                                                                                                                         cgtgCTGGCCCGCGAGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTCCGCCGGCAGGCGCGCTTCATAGAAAGCGTACAGCGTttc     >  1:794433/4‑150 (MQ=255)
                                                                                                                                          gtgCTGGCCCGCGAGGCCAGCACGTCGTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCAATGCGGTACCAACTGTCAAAGGTCGCCGTCTGGTGGATTTCCGCCGGCAGGCGCGTTTCATTAAAAGCGTACAGCGTttct    >  2:229757/3‑150 (MQ=255)
                                                                                                                                            gCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTCCGCCGGCAGGCGCGCTTCATAGAAAGCGTACAGCGTttcttc  <  1:934253/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
CCGGCAGGCTCGGCAACAGCGGTGCCGGCACCCGCACCGTCACACCATCACGGGGGTGGCCGGGCTCGAAGTGGTAGCTCAGCGGCAGGCGCAGATCGCCCACCTGCATACTGTCGGGGTACTGTGCGGCGGTCACTTCGCTGGCCTCGCGGGCCAGCACGTCCTCCTCGCGCATGATCAGCAGGTTGGCGTCCTTTTGGCTGCCCATGCGGTACCAACTATCAAAGGTCGCCGTCTGGTGGATTTCCGCCGGCAGGCGCGCTTCATAGAAAGCGTACAGCGTTTCTTC  >  NC_002947/5165269‑5165557

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