Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,434,888 T→G 10.3% N142H (AAC→CAC)  hrpB ← ATP‑dependent helicase HrpB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,434,8880TG10.3% 90.6 / 3.8 39N142H (AAC→CAC) hrpBATP‑dependent helicase HrpB
Reads supporting (aligned to +/- strand):  ref base T (15/20);  new base G (2/2);  total (17/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.87e-01

CCCCAGCGGATATCCACCGGATGCATGCGGCCTTCGCTGCTGACCACCGGGGCATCGTCCAGCAACCGCGACAGGCGTTCGCCTTCCAGGGTCGCGGACATCAGCAGAATCTTCAGTGGCGGGTCGTCACGCAGCAACTCGCGGCCGTTAAGGCTCAGCGCAAGTGCGAGATCGGCATCCAGGCTGCGCTCGTGAAATTCGTCGAAGATCAGCAGGCCCACACCTTCCAGCGCCGGGTCGGCTTGCAGGCGGCGTGTGAGGATGCCTTCGGTCACCACTTCGA  >  NC_002947/5434740‑5435022
                                                                                                                                                    |                                                                                                                                      
ccccAGCGGATATCCACCGGATGCATGCGGCCTTCGCTGCTGACCACCGGGGCATCGTCCAGCAACCGCGACAGGCGTTCGCCTTCCAGGGTCGCGGACATCAGCAGAATCTTCAGTGGCGGGTCGTCACGCAGCAACTCGCGGCCGTTa                                                                                                                                       >  1:201680/1‑150 (MQ=255)
    aGCGGATATCCACCGGATGCATGCGGCCTTCGCTGCTGACCACCGGGGCATCGTCCAGCAACCGCGACAGGCGTTCGCCTTCCAGGGTCGCGGACATCAGCAGAATCTTCAGTGGCGGGTCGTCACGCAGCAACTCGCGGCCGTTAAGGc                                                                                                                                   >  2:375121/1‑150 (MQ=255)
     gCGGATATCCACCGGATGCATGCGGCCTTCGCTGCTGACCACCGGGGCATCGTCCAGCAACCGCGACAGGCGTTCGCCTTCCAGGGTCGCGGACATCAGCAGAATCTTCAGTGGCGGGTCGTCACGCAGAAACTCGCGGCCGTTAAGGCt                                                                                                                                  >  2:514805/1‑150 (MQ=255)
               aCCGGATGCATGCGGCCTTCGCTGCTGACCACCGGGGCATCGTCCAGCAACCGCGACAGGCGTTCGCCTTCCAGGGTCGCGGACATCAGCAGAATCTTCAGTGGCGGGTCGTCACGCAGCAACTCGCGGCCGTTAAGGCTCAGCGCAAGt                                                                                                                        <  1:759449/150‑1 (MQ=255)
                ccGGATGCATGCGGCCTTCGCTGCTGACCACCGGGGCATCGTCCAGCAACCGCGACAGGCGTTCGCCTTCCAGGGTCGCGGACATCAGCAGAATCTTCAGTGGCGGGTCGTCACGCAGCAACTCGCGGCCGTTAAGGCTCAGCGCAAGTg                                                                                                                       >  2:269382/1‑150 (MQ=255)
                        atgcGGCCTTCGCTGCTGACCACCGGGGCATCGTCCAGCAACCGCGACAGGCGTTCGCCTTCCAGGGTCGCGGACATCAGCAGAATCTTCAGTGGCGGGTCGTCACGCAGCAACTCGCGGCCGTTAAGGCTCAGCGCAAGTGCGAGATCg                                                                                                               <  1:204822/150‑1 (MQ=255)
                                ttCGCTGCTGACCACCGGGGCATCGTCCAGCAACCGCGACAGGCGTTCGCCTTCCAGGGTCGCGGACATCAGCAGAATCTTCAGTGGCGGGTCGTCACGCAGCAACTCGCGGCCGTTAAGGCTCAGCGCAAGTGc                                                                                                                      >  2:641206/1‑135 (MQ=255)
                                ttCGCTGCTGACCACCGGGGCATCGTCCAGCAACCGCGACAGGCGTTCGCCTTCCAGGGTCGCGGACATCAGCAGAATCTTCAGTGGCGGGTCGTCACGCAGCAACTCGCGGCCGTTAAGGCTCAGCGCAAGTGc                                                                                                                      <  1:641206/135‑1 (MQ=255)
                                  cGCTGCTGACCACCGGGGCATCGTCCAGCAACCGCGACAGGCGTTCGCCTTCCAGGGTCGCGGACATCAGCAGAATCTTCAGTGGCGGGTCGTCACGCAGCAACTCGCGGCCGTTAAGGCTCAGCGCAAGTGCGAGATCGGCATCCAGGc                                                                                                     >  1:535025/1‑150 (MQ=255)
                                         gACCACCGGGGCATCGTCCAGCAACCGCACAAAGAGTTCTCCTGTCAGGTTCTCGGACATCAGCAGAACTTTAGGTCGCGGGTCGTCACGCAGCAACTCGCGGCCGTTAAGGCTCAGCGCAAGTGCGAGATCGGCATCCAGGCTGCGCTc                                                                                              <  1:405337/150‑1 (MQ=255)
                                                     aTCGTCCAGCAACCGCGACAGGCGTTCGCCTTCCAGGGTCGCGGACATCAGCAGAATCTTCAGTGGCGGGTCGTCACGCAGCAACTCGCGGCCGTTAAGGCTCAGCGCAAGTGCGAGATCGGCATCCAGGCTGCGCTCGTGAAATtcgtc                                                                                  <  1:243623/150‑1 (MQ=255)
                                                     aTCGTCCAGCAACCGCGACAGGCGTTCGCCTTCCAGGGTCGCGGACATCAGCAGAATCTTCAGTGGCGGGTCGTCACGCAGCAACTCGCGGCCGTTAAGGCTCAGCGCAAGTGCGAGATCGGCATCCAGGCTGCGCTCGTGAAATtcgtc                                                                                  <  1:269382/150‑1 (MQ=255)
                                                     aTCGTCCAGCAACCGCGACAGGCGTTCGCCTTCCAGGGTCGCGGACATCAGCAGAATCTTCAGTGGCGGGTCGTCACGCAGCAACTCGCGGCCGTTAAGGCTCAGCGCAAGTGCGAGATCGGCATCCAGGCCGCGCTCGTGAAATtcgtc                                                                                  <  2:535025/150‑1 (MQ=255)
                                                               aaCCGCGACAGGCGTTCGCCTTCCAGGGTCGCGGACATCAGCAGAATCTTCAGTGGCGGGGCGTCGCGCAGCAACTCGCGGCCGTTAAGGCTCAGCGCAAGTGCGAGATCGGCATCCAGGCTGCGCTCGTGAAATTCGTCGAAGATcagc                                                                        >  1:383151/1‑150 (MQ=255)
                                                                   gcgACAGGCGTTCGCCTTCCAGGGTCGCGGACATCAGCAGAATCTTCAGTGGCGGGTCGTCACGCAGCAACTCGCGGCCGTTAAGGCTCAGCGCAAGTGCGAGATCGGCATCCAGGCTGCGCTCGTGAAATTCGTCGAAGATCAGCAGGc                                                                    <  2:771500/150‑1 (MQ=255)
                                                                   gcgACAGGCGTTCGCCTTCCAGGGTCGCGGACATCAGCAGAATCTTCAGTGGCGGGTCGTCACGCAGCAACTCGCGGCCGTTAAGGCTCAGCGCAAGTGCGAGATCGGCATCCAGGCTGCGCTCGTGAAATTCGTCGAAGATCAGCAGGc                                                                    >  2:643798/1‑150 (MQ=255)
                                                                             ttCGCCTTCCAGGGTCGCGGACATCAGCAGAATCTTCAGTGGCGGGTCGTCACGCAGCAACTCGCGGCCGTTAAGGCTCAGCGCAAGTGCGAGATCGGCATCCAGGCTGCGCTCGTGAAATTCGTCGAAGATCAGCAGGCCCACACCTTc                                                          <  1:779830/150‑1 (MQ=255)
                                                                                    tCCAGGGTCGCGGACATCAGCAGAATCTTCAGTGGCGGGTCGTCACGCAGCAACTCGCGGCCGTTAAGGCTCAGCGCAAGTGCGAGATCGGCATCCAGGCTGCGCTCGTGAAATTCGTCGAAGATCAGCAGGCCCACACCTTCCAGCGcc                                                   <  1:375121/150‑1 (MQ=255)
                                                                                              cgGACATCAGCAGAATCTTCAGTGGCGGGTCGTCACGCAGCAACTCGCGGCCGTTAAGGCTCAGCGCAAGTGCGAGATCGGCATCCAGGCTGCGCTCGTGAAATTCGTCGAAGATCAGCAGGCCCACACCTTCCAGCGCCGGGTCGGCtt                                         <  2:746439/150‑1 (MQ=255)
                                                                                                gACATCAGCAGAATCTTCAGTGGCGGGTCGTCACGCAGCAACTCGCGGCCGTTAAGGCTCAGCGCAAGTGCGAGATCGGCATCCAGGCTGCGCTCGTGAAATTCGTCGAAGATCAGCAGGCCCACACCTTCCAGCGCCGGGTCGGCTTGc                                       <  2:563761/150‑1 (MQ=255)
                                                                                                  cATCAGCAGAATCTTCAGTGGCGGGTCGTCACGCAGCAACTCGCGGCCGTTAAGGCTCAGCGCAAGTGCGAGATCGGCATCCAGGCTGCGCTCGTGAAATTCGTCGAAGATCAGCAGGCCCACACCTTCCAGCGCCGGGTCGGCTTCCAg                                     <  2:461818/150‑1 (MQ=255)
                                                                                                   aTCAGCAGAATCTTCAGTGGCGGGTCGTCACGCAGCAACTCGCGGCCGTTAAGGCTCAGCGCAAGTGCGAGATCGGCATCCAGGCTGCGCTCGTGAAATTCGTCGAAGATCAGCAGGCCCACACCTTCCAGCGCCTGGTCGGCt                                          <  2:423455/144‑1 (MQ=255)
                                                                                                   aTCAGCAGAATCTTCAGTGGCGGGTCGTCACGCAGCAACTCGCGGCCGTTAAGGCTCAGCGCAAGTGCGAGATCGGCATCCAGGCTGCGCTCGTGAAATTCGTCGAAGATCAGCAGGCCCACACCTTCCAGCGCCGGGTCGGCt                                          >  1:423455/1‑144 (MQ=255)
                                                                                                     cagcagAATCTTCAGTGGCGGGTCGTCACGCAGCAACTCGCGGCCGTTAAGGCTCAGCGCAAGTGCGAGATCGGCATCCAGGCTGCGCTCGTGAAATTCGTCGAAGATCAGCAGGCCCACACCTTCCAGCGCCGGGTCGGCTTGCAggcg                                  >  2:454246/1‑150 (MQ=255)
                                                                                                          gAATCTTCAGTGGCGGGTCGTCACGCAGCAACTCGCGGCCGTTAAGGCTCAGCGCAAGTGCGAGATCGGCATCCAGGCTGCGCTCGTGAAATTCGTTGAAGATCAGCAGGCCCACACCTTCCAGCGCCGGGTCGGCTTGCAGGCGGCgtg                             >  1:415559/1‑150 (MQ=255)
                                                                                                                tCAGTGGCGGGTCGTCACTCAGCAACTCGCGGCCGTTAAGGCTCAGCGCAAGTGCGAGATCGGCATCCAGGCTGCGCTGGTGAAATTCGTCGAAAATCAGCAGGCCCACACCTTCCAGCGCCGGGTCGCCTTCCAGGCGGCGAGtgatga                       >  2:348145/1‑150 (MQ=255)
                                                                                                                 cAGTGGCGGGTCGTCACGCAGCAACTCGCGGCCGTTAAGGCTCAGCGCAAGTGCGAGATCGGCATCCAGGCTGCGCTCGTGAAATTCGTCGAAGATCAGCAGGCCCACACCTTCCAGCGCCGGGTCGGCTTGCAGGCGGCGTGTGAGGAt                      <  1:734339/150‑1 (MQ=255)
                                                                                                                  aGTGGCGGGTCGTCACGCAGCAACTCGCGGCCGTTAAGGCTCAGCGCAAGTGCGAGATCGGCATCCAGGCTGCGCTCGTGAAATTCGTCGAAGATCAGCAGGCCCACACCTTCCAGCGCCGGGTCGGCTTGCAGGCGGCGTGTGAGGATg                     >  2:46234/1‑150 (MQ=255)
                                                                                                                      gCGGGTCGTCACGCAGCAACTCGCGGCCGTGAAGGCTCAGCGCAAGTGCGAGATCGGCATCCAGGCTGCGCTCGTGAAATTCGTCGAAGATCAGCAGGCCcaca                                                               <  1:716973/104‑1 (MQ=255)
                                                                                                                      gCGGGTCGTCACGCAGCAACTCGCGGCCGTGAAGGCTCAGCGCAAGTGCGAGATCGGCATCCAGGCTGCGCTCGTGAAATTCGTCGAAGATCAGCAGGCCcaca                                                               >  2:294539/1‑104 (MQ=255)
                                                                                                                      gCGGGTCGTCACGCAGCAACTCGCGGCCGTGAAGGCTCAGCGCAAGTGCGAGATCGGCATCCAGGCTGCGCTCGTGAAATTCGTCGAAGATCAGCAGGCCcaca                                                               >  2:716973/1‑104 (MQ=255)
                                                                                                                      gCGGGTCGTCACGCAGCAACTCGCGGCCGTGAAGGCTCAGCGCAAGTGCGAGATCGGCATCCAGGCTGCGCTCGTGAAATTCGTCGAAGATCAGCAGGCCcaca                                                               <  1:294539/104‑1 (MQ=255)
                                                                                                                        gggTCGTCACGCAGCAACTCGCCGCCGTTAAGGCTCAGCGCAAGTGCGAGATCGGCATCCAGGCTGCGCTCGTGAAATTCGTCGAAGATCAGCAGGCCCACACCTTCCAGCGCCGGGTCGGCTTGCAGGCGGCGTGTGAGGATGCCTTCg               <  1:680955/150‑1 (MQ=255)
                                                                                                                           tcgtcACGCAGCAACTCGCGGCCGTTAAGGCTCAGCGCAAGTGCGAGATCGGCATCCAGGCTGCGCTCGTGAAATTCGTCGAAGATCAGCAGGCCCACACCTTCCAGCGCCgg                                                 <  2:682259/113‑1 (MQ=255)
                                                                                                                           tcgtcACGCAGCAACTCGCGGCCGTTAAGGCTCAGCGCAAGTGCGAGATCGGCATCCAGGCTGCGCTCGTGAAATTCGTCGAAGATCAGCAGGCCCACACCTTCCAGCGCCgg                                                 >  1:682259/1‑113 (MQ=255)
                                                                                                                           tcgtcACGCAGCAACTCGCGGCCGTTAAGGCTCAGCGCAAGTGCGAGATCGGCATCCAGGCTGCGCTCGTGAAATTCGTCGAAGATCAGCAGGCCCACACCTTCCAGCGCCGGGTCGGCTTGCAGGCGGCGTGTGAGGAt                      <  1:570115/140‑1 (MQ=255)
                                                                                                                           tcgtcACGCAGCAACTCGCGGCCGTTAAGGCTCAGCGCAAGTGCGAGATCGGCATCCAGGCTGCGCTCGTGAAATTCGTCGAAGATCAGCAGGCCCACACCTTCCAGCGCCGGGTCGGCTTGCAGGCGGCGTGTGAGGAt                      >  2:570115/1‑140 (MQ=255)
                                                                                                                                  gcagcaACTCGCGGCCGTTAAGGCTCAGCGCAAGTGCGAGATCGGCATCCAGGCTGCGCTCGTGAAATTCGTCGAAGATCAGCAGGCCCACACCTTCCAGCGCCGGGTCGGCTTGCAGGCGGCGTGTGAGGATGCCTTCGGTCACCACtt     <  1:643798/150‑1 (MQ=255)
                                                                                                                                     gcaACTCGCGGCCGTTAAGGCTCAGCGCAAGTGCGAGATCGGCATCCAGGCTGCGCTCGTGAAATTCGTCGAAGATCAGCAGGCCCACACCTTCCAGCGCCGGGTCGGCTTGCAGGCGGCGTGTGAGGATGCCTTCGGTCACCACTTCGa  <  2:644583/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                      
CCCCAGCGGATATCCACCGGATGCATGCGGCCTTCGCTGCTGACCACCGGGGCATCGTCCAGCAACCGCGACAGGCGTTCGCCTTCCAGGGTCGCGGACATCAGCAGAATCTTCAGTGGCGGGTCGTCACGCAGCAACTCGCGGCCGTTAAGGCTCAGCGCAAGTGCGAGATCGGCATCCAGGCTGCGCTCGTGAAATTCGTCGAAGATCAGCAGGCCCACACCTTCCAGCGCCGGGTCGGCTTGCAGGCGGCGTGTGAGGATGCCTTCGGTCACCACTTCGA  >  NC_002947/5434740‑5435022

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