Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,166,231 T→C 9.3% D550G (GAC→GGC)  hrpA ← ATP‑dependent helicase HrpA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,166,2310TC 9.3% 114.6 / 2.8 43D550G (GAC→GGC) hrpAATP‑dependent helicase HrpA
Reads supporting (aligned to +/- strand):  ref base T (22/17);  new base C (2/2);  total (24/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

CGATGAGCATCGCGCCACTCGCGCAGCCGCAAGTAGTTCAGGAACTGCTTGCGACACCAGTTGCGCAGCGGGTTGGCGGTCAGCGCCTGGCGTTGTTCCTCGAAGCCACGCCACAGGTTGACCAGCGCAGCGAAGTCCGAGTCCACGTCCTTCCACTGTGCGTGCGCCTGATCGGCGGCCTGCTGGCGTTCTGGGGGGCGTTCACG‑CGGGTCCTGCACCGACAGCGCGCTGGTGACAATCAGCACTTCCTGCAGGCTGCCCAGGCGCGCGCCTTCAAGCAGCATGCGGCCCAGA  >  NC_002947/5166084‑5166377
                                                                                                                                                   |                                                                                                                                                   
cGATGAGCATCGCGCCACTCGCGCAGCCGCAAGTAGTTCAGGAACTGCTTGCGACACCAGTTGCGCAGCGGGTTGGCGGTCAGCGCCTGGCGTTGTTCCTCGAAGCCACGCCACAGGTTGACCAGCGCAGCGAAGTCCGAGTCCACGt                                                                                                                                                     <  6:309545/148‑1 (MQ=255)
      gCATCGCGCCACTCGCGCAGCCGCAAGTAGTTCAGGAACTGCTTGCGACACCAGTTGCTCAGCGGGTTGGCGGTCAGCGCCTGGCGTTGTTCCTCGAAGCCACGCCACAGGTTGACCAGCGCAGCGAAGTCCGAGTCCACGTCCTTCCa                                                                                                                                              <  6:95489/149‑1 (MQ=255)
      gCATCGCGCCACTCGCGCAGCCGCAAGTAGTTCAGGAACTGCTTGCGACACCAGTTGCGCAGCGGGTTGGCGGTCAGCGCCTGGCGTTGTTCCTCGAAGCCACGCCACAGGTTGACCAGCGCAGCGAAGTCCGAGTCCACGTCCTTCCa                                                                                                                                              >  7:212817/1‑149 (MQ=255)
            cgcCACTCGCGCAGCCGCAAGTAGTTCAGGAACTGCTTGCGACACCAGTTGCGCAGCGGGTTGGCGGTCAGCGCCTGGCGTTGTTCCTCGAAGCCACGCCACAGGTTGACCAGCGCAGCGAAGTCCGAGTCGACGTCCTTCCACTgtgc                                                                                                                                        >  7:295267/1‑149 (MQ=255)
                 cTCGCGCAGCCGCAAGTAGTTCAGGAACTGCTTGCGACACCAGTTGCGCAGCGGGTTGGCGGTCAGCGCCTGGCGTTGTTCCTCGAAGCCACGCCACAGGTTGACCAGCGCAGCGAAGTCCGAGTCCACGCCCTTCCACTGtgcgtgcg                                                                                                                                   >  7:196971/1‑149 (MQ=255)
                 cTCGCGCAGCCGCAAGTAGTTCAGGAACTGCTTGCGACACCAGTTGCGCAGCGGGTTGGCGGTCAGCGCCTGGCGTTGTTCCTCGAAGCCACGCCACAGGTTGACCAGCGCAGCGAAGTCCGAGTCCACGCCCTTCCACTGtgcgtgcg                                                                                                                                   >  5:51900/1‑149 (MQ=255)
                       cAGCCGCAAGTAGTTCAGGAACTGCTTGCGACACCAGTTGCGCAGCGGGTTGGCGGTCAGCGCCTGGCGTTGTTCCTCGAAGCCACGCCACAGGTTGACCAGCGCAGCGAAGTCCGAGTCCACGCCCTTCCACTGTGCGTGCGCCTGAt                                                                                                                             <  6:51900/149‑1 (MQ=255)
                       cAGCCGCAAGTAGTTCAGGAACTGCTTGCGACACCAGTTGCGCAGCGGGTTGGCGGTCAGCGCCTGGCGTTGTTCCTCGAAGCCACGCCACAGGTTGACCAGCGCAGCGAAGTCCGAGTCCACGCCCTTCCACTGTGCGTGCGCCTGAt                                                                                                                             <  8:196971/149‑1 (MQ=255)
                                      cAGGAACTGCTTGCGACACCAGTTGCGCAGCGGGTTGGCGGTCAGCGCCTGGCGTTGTTCCTCGAAGCCACGCCACAGGTTGACCAGCGCAGCGATGTCCGAGTCCACGTCCTTCCACtgtg                                                                                                                                         >  1:296937/1‑122 (MQ=255)
                                      cAGGAACTGCTTGCGACACCAGTTGCGCAGCGGGTTGGCGGTCAGCGCCTGGCGTTGTTCCTCGAAGCCACGCCACAGGTTGACCAGCGCAGCGAAGTCCGAGTCCACGTCCTTCCACtgtg                                                                                                                                         <  2:296937/122‑1 (MQ=255)
                                       aGGAACTGCTTGCGACACCAGTTGCGCAGCGGGTTGGCGGTCAGCGCCTGGCGTTGTTCCTCGAAGCCACGCCACAGGTTGACCAGCGCAGCGAAGTCCGAGTCCACGTCCTTCCACTGTGCGTGCGCCTGATCGGCGGCCTGCTcgcg                                                                                                             >  4:296278/1‑149 (MQ=255)
                                            cTGCTTGCGACACCAGTTGCGCAGCGGGTTGGCGGTCAGCGCCTGGCGTTGTTCCTCGAAGCCACGCCACAGGTTGACCAGCGCAGCGAAGTCCGAGTCCACGTCCTTCCACTGTGCGTGCGCCTGATCGGCGGCCTGCTGGCGTTCTg                                                                                                        <  5:290156/149‑1 (MQ=255)
                                               cTTGCGACACCAGTTGCGCAGCGGGTTGGCGGTCAGCGCCTGGCGTTGTTCCTCGAAGCCACGCCACAGGTTGACCAGCGCAGCGAAGTCCGAGTCCACGTCCTTCCACTGTGCGTGCGCCTGATCGGCGGCCTGCTGGCGTTCTggg                                                                                                      <  2:187890/148‑1 (MQ=255)
                                                      cacCAGTTGCGCAGCGGGTTGGCGGTCAGCGCCTGGCGTTGTTCCTCGAAGCCACGCCACAGGTTGACCAGCGCAGCGAAGTCCGAGTCCACGTCCTTCCACTGTGCGTGCGCCTGATCGGCGGCctgct                                                                                                                 >  8:101006/1‑130 (MQ=255)
                                                         cAGTTGCGCAGCGGGTTGGCGGTCAGCGCCTGGCGTTGTTCCTCGAAGCCACGCCACAGGTTGACCAGCGCAGCGAAGTCCGAGTCCACGTCCTTCCACTGTGCGTGCGCCTGATCGGCGGCCTGCTGGCGTTCTGGGGGGCGTTCAcg                                                                                           >  8:2413/1‑149 (MQ=255)
                                                         cAGTTGCGCAGCGGGTTGGCGGTCAGCGCCTGGCGTTGTTCCTCGAAGCCACGCCACAGGTTGACCAGCGCAGCGAAGTCCGAGTCCACGTCCTCCCACTGTGCGTGCGCCTGATCGGCGGCctgctg                                                                                                                >  2:348164/1‑128 (MQ=255)
                                                         cAGTTGCGCAGCGGGTTGGCGGTCAGCGCCTGGCGTTGTTCCTCGAAGCCACGCCACAGGTTGACCAGCGCAGCGAAGTCCGAGTCCACGTCCTCCCACTGTGCGTGCGCCTGATCGGCGGCctgctg                                                                                                                <  1:348164/128‑1 (MQ=255)
                                                          aGTTGCGCAGCGGGTTGGCGGTCAGCGCCTGGCGTTGTTCCTCGAAGCCACGCCACAGGTTGACCAGCGCAGCGCAGTCCGAGTCCACGTCCTTCCACTGTGCGTGCGCCTGATCGGCGGCCTGCTGGCGTTCTGGGGCGCGTTCAcg‑c                                                                                         >  5:342634/1‑149 (MQ=255)
                                                            ttGCGCAGCGGGTTGGCGGTCAGCGCCTGGCGTTGTTCCTCGAAGCCACGCCACAGGTTGACCAGCGCAGCGAAGTCCGAGTCCACGTCCTTCCACTGTGCGTGCGCCTGATCGGCGGCCTGCTGGCGTTCTGGGGGGCGTTCACG‑Cgg                                                                                       <  8:212817/149‑1 (MQ=255)
                                                              gcgcAGCGGGTTGGCGGTCAGCGCCTGGCGTTGTTCCTCGAAGCCACGCCACAGGTTGACCAGCGCAGCGAAGTCCGAGTCCACGTCCTTCCACTGTGCGTGCGCCTGATCGGCGGCCTGCTGGCGTTCTGGGGGGCGTTCACG‑Cggg                                                                                      <  2:282046/148‑1 (MQ=255)
                                                                                 aGCGCCTGGCGTTGTTCCTCGAAGCCACGCCACAGGTTGACCAGCGCAGCGAAGTCCGAGTCCACGTCCTTCCACTGTGCGTGCGCCTGATCGGCGGCCTGCTGGCGTTCTGGGGGGCGTTCACG‑CGGGTCCTGCACCGACAGCGCGCt                                                                  <  1:341442/149‑1 (MQ=255)
                                                                                  gcgcCTGGCGTTGTTCCTCGAAGCCACGCCACAGGTTGACCAGCGCAGCGAAGTCCGAGTCCACGTCCTTCCACTGTGCGTGCGCCTGATCGGCGGCctgctg                                                                                                                >  2:27543/1‑103 (MQ=255)
                                                                                  gcgcCTGGCGTTGTTCCTCGAAGCCACGCCACAGGTTGACCAGCGCAGCGAAGTCCGAGTCCACGTCCTTCCACTGTGCGTGCGCCTGATCGGCGGCctgctg                                                                                                                <  1:27543/103‑1 (MQ=255)
                                                                                              gttCCTCGAAGCCACGCCACAGGTTGACCAGCGCAGCGAAGTCCGAGTCCACGTCCTTCCACTGTGCGTGCGCCTGATCGGCGGCCTGCTGGCGTTCTGGGGGGCGTTCACG‑CGGGTCCTGCACCGACAGCGCGCTGGTGACAATCAGc                                                     >  6:295312/1‑149 (MQ=255)
                                                                                                   tCGAAGCCACGCCACAGGTTGACCAGCGCAGCGAAGTCCGAGTCCACGTCCTTCCACTGTGCGTGCGCCTGATCGGCGGCCTGCTGGCGTTCTGGGGGGCGTTCACG‑CGGGTCCTGCACCGACAGCGCGCTGGTGACAATCAGCACTTc                                                >  6:181796/1‑149 (MQ=255)
                                                                                                           acgccacAGGTTGACCAGAGCAGCGAAGTCCGAGTCCACGTCCTTCCACTGTGCGTGCGCCTGATCCGCGGCCTGCTGGCCTTCTGGGGGGCGTTCACG‑CGGGTCCTGCACCGACAGCGCGCTGGTGACAATCAGCACTTCCTGCAGGc                                        >  5:125430/1‑149 (MQ=255)
                                                                                                                    gTTGACCAGCGCAGCGAAGTCCGAGTCCACGTCCTTCCACTGTGCGTGCGCCTGATCGGCGGCCTGCTGGCGTTCTGGGGGGCGTTCACG‑CGGGTCCTGCACCGACAGCGCGCTGGTGACAATCAGCACTTCCTGCAGGCTGCCCAGgc                               >  4:171437/1‑149 (MQ=255)
                                                                                                                        aCCAGCGCAGCGAAGTCCGAGTCCACGTCCTTCCACTGTGCGTGCGCCTGATCGGCGGCCTGCTGGCGTTCTGGGGGGCGTTCACGCCGGGTCCTGCACCGACAGCGCGCTGGt                                                               <  1:187907/114‑1 (MQ=255)
                                                                                                                        aCCAGCGCAGCGAAGTCCGAGTCCACGTCCTTCCACTGTGCGTGCGCCTGATCGGCGGCCTGCTGGCGTTCTGGGGGGCGTTCACG‑CGGGTCCTGCACCGACAGCGCGCTGGt                                                               >  2:187907/1‑113 (MQ=255)
                                                                                                                         ccAGCGCAGCGAAGTCCGAGTCCACGTCCTTCCACTGTGCGTGCGCCTGATCGGCGGCCTGCTGGCGTTCTGGGGGGCGTTCACG‑CGGGTCCTGCTCCGACAGCGCGCTGGTGACAATCAGCACTTCCTGCAGGCTGCCCAGgcgcgcg                          >  1:295558/1‑149 (MQ=255)
                                                                                                                           agcgcagcgAAGTCCGAGTCCACGTCCTTCCACTGTGCGTGCGCCTGATCGGCGGCCTGCTGGCGTTCTGGGGGGCGTTCACG‑CGGGTCCTGCACCGACAGCGCGCTGGTGACAATCAGCACTTCCTGCAGGCTGCCCAGGCGCGCGcc                        <  7:2413/149‑1 (MQ=255)
                                                                                                                            gcgcagcgAAGTCCGAGTCCACGTCCTTCCACTGTGCGTGCGCCTGATCGGCGGCCTGCTGGCGTTCTGGGGGGCGTTCACG‑CGGGTCCTGCACCGACAGCGCGCTGGTGACAATCAGCACTTCCTGCAGGCTGCCCAGGCGCGCGCCt                       >  7:3775/1‑149 (MQ=255)
                                                                                                                              gcagcgAAGTCCGAGTCCACGTCCTTCCACTGTGCGTGCGCCTGATCGGCGGCCTGCTGGCGTTCTGGGGGGCGTTCACG‑CGGGTCCTGCACCGACAGCGCGCTGGTGACAATCAGCACTTCCTGCAGGCTGCCCAGGCGCGCGCCTTc                     <  5:181796/149‑1 (MQ=255)
                                                                                                                               cagcgAAGTCCGAGTCCACGTCCTTCCACTGTGCGTGCGCCTGATCGGCGGCCTGCTGGCGTTCTGGGGGGCGTTCACG‑CGGGTCCTGCACCGACAGCGCGCTGGTGACAATCAGCACTTCCTGCAGGCTGCCCAGGCGCGCACCTTCa                    <  6:119716/149‑1 (MQ=255)
                                                                                                                                agcgAAGTCCGAGTCCACGTCCTTCCACTGTGCGTGCACCTGATCGGCGGCCTGCTGGCGTTCTGGGGGGCGTTCACG‑CGGGTCCTGCACCGACAGCGCGCTGGTGACAATCAGCACTTCCTGCAGGCTGCCCAGGCGCGCGCCTTc                     >  6:182548/1‑147 (MQ=255)
                                                                                                                                  cgAAGTCCGAGTCCACGTCCTTCCACTGTGCGTGCGCCTGATCGGCGGCCTGCTGGCGTTCTGGGGGGCGTTCAcg‑c                                                                                         <  1:135265/77‑1 (MQ=255)
                                                                                                                                  cgAAGTCCGAGTCCACGTCCTTCCACTGTGCGTGCGCCTGATCGGCGGCCTGCTGGCGTTCTGGGGGGCGTTCAcg‑c                                                                                         >  2:135265/1‑77 (MQ=255)
                                                                                                                                  cgAAGTCCGAGTCCACGTCCTTCCACTGTGCGTGCGCCTGATCGGCGGCCTGCTGGCGTTCTGGGGGGCGTTCACG‑CGGGTc                                                                                    <  2:170221/82‑1 (MQ=255)
                                                                                                                                  cgAAGTCCGAGTCCACGTCCTTCCACTGTGCGTGCGCCTGATCGGCGGCCTGCTGGCGTTCTGGGGGGCGTTCACG‑CGGGTc                                                                                    >  1:170221/1‑82 (MQ=255)
                                                                                                                                         cGAGTCCACGTCCTTCCACTGTGCGTGCGCCTGATCGGCGGCCTGCTGGCGTTCTGGGGGGCGTTCACG‑CGGGTCCTGCACCGACAGCGCGCTGGTGACAATCAGCACTTCCTGCAGGCTGCCCAGGCGCGCGCCTTCAAGCAGCATGc          >  7:260556/1‑149 (MQ=255)
                                                                                                                                           aGTCCACGTCCTTCCACTCTGCGGGCGCCTGATCGGCGGCCTGCTGGCGTTCTGGGGGGCGTTCACG‑CGGGGCCTGCACCGACAGCGCGCTGGTGACAATCAGCACTTCCTGCAGGCTGCCCAGGCGCGCGCCTTCAAGCAGCATGCgg        <  3:296278/149‑1 (MQ=255)
                                                                                                                                                aCGTCCTTCCACTGTGCGTGCGCCTGATCGGCCGCCTGCTCGCGTTCTGGGGGGCGTTCACG‑CGGGTCCTGCACCGACAGCGCCCTGGTGAAAATCAGCACTTCCTGCAGGCTGCCCAGGCGCGCGCCTTCAAGCAGCATGCGGcc      >  7:65291/1‑146 (MQ=255)
                                                                                                                                                 cGTCCTTCCACTGTGCGTGCGCCTGATCGGCGGCCTGCTGGCGTTCTGGGGGGCGTTCACG‑CGGGTCCTGCACCGACAGCGCGCTGGTGACAATCAGCACTTCCTGCAGGCTGCCCAGGCGCGCGCCTTCAAGCAGCATGCGGCCCAGa  <  5:165821/149‑1 (MQ=255)
                                                                                                                                                 cGTCCTTCCACTGTGCGTGCGCCTGATCGGCGGCCTGCTGGCGTTCTGGGGGGCGTTCACG‑CGGGTCCTGCACCGACAGCGCGCTGGTGACAATCAGCACTTCCTGCAGGCTGCCCAGGCGCGCGCCTTCAAGCAGCATGCGGCCCAGa  >  1:95776/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CGATGAGCATCGCGCCACTCGCGCAGCCGCAAGTAGTTCAGGAACTGCTTGCGACACCAGTTGCGCAGCGGGTTGGCGGTCAGCGCCTGGCGTTGTTCCTCGAAGCCACGCCACAGGTTGACCAGCGCAGCGAAGTCCGAGTCCACGTCCTTCCACTGTGCGTGCGCCTGATCGGCGGCCTGCTGGCGTTCTGGGGGGCGTTCACG‑CGGGTCCTGCACCGACAGCGCGCTGGTGACAATCAGCACTTCCTGCAGGCTGCCCAGGCGCGCGCCTTCAAGCAGCATGCGGCCCAGA  >  NC_002947/5166084‑5166377

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