Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 163,992 G→A A630T (GCG→ACG)  hrpB → putative ATP‑dependent RNA helicase HrpB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913163,9920GA92.3% 68.1 / ‑2.5 26A630T (GCG→ACG) hrpBputative ATP‑dependent RNA helicase HrpB
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (13/11);  total (14/12)
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.75e-01

CTGGCGTCGGCTACAAATCGGTCAGTTGACGGTGAAAGTGCAGCCGCTGGCGAAACCGTCAGAAGACGAGTTGCATCAGGCGATGCTTAATGGCATCCGTGATAAAGGTTTAAGCGTGCTCAACTGGACGGCGGAAGCGGAACAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCA  >  NC_000913/163856‑164126
                                                                                                                                        |                                                                                                                                      
cTGGCGTCGGCTACAAATCGGTCAGTTGACGGTGAAAGTGCAGCCGCTGGCGAAACCGTCAGAAGACGAGTTGCATCAGGCGATGCTTAATGGCATCCGTGATAAAGGTTTAAGCGTGCTCAACTGGACGgcggaaacgg                                                                                                                                     <  1:427631/140‑1 (MQ=255)
 tGGCGTCGGCTACAAATCGGTCAGTTGACGGTGAAAGTGCAGCCGCTGGCGAAACCGTCAGAAGACGAGTTGCATCAGGCGATGCTTAATGGCATCCGTGATAAAGGTTTAAGCGTGCTCAACTGGACGgcggaaacgga                                                                                                                                    <  1:13162/140‑1 (MQ=255)
         gCTACAAATCGGTCAGTTGACGGTGAAAGTGCAGCCGCTGGCGAAACCGTCAGAAGACGAGTTGCATCAGGCGATGCTTAATGGCATCCGTGATAAAGGTTTAAGCGTGCTCAACTGGACGGCGGAAACGGAACAGCTAc                                                                                                                            <  1:5182/140‑1 (MQ=255)
                     tCAGTTGACGGTGAAAGTGCAGCCGCTGGCGAAACCGTCAGAAGACGAGTTGCATCAGGCGATGCTTAATGGCATCCGTGATAAAGGTTTAAGCGTGCTCAACTGGACGGCGGAAACGGAACAGCTACGCTTGCGTTTGt                                                                                                                >  2:467434/1‑140 (MQ=255)
                      cAGTTGACGGTGAAATTGCAGCTGCGGGCTAAACCGTGAGAAGACGAGTTGCATCAGGCGATGCTTAATGGCATCCGTGATAAAGGTTTAAGCGTGCTCAACTGGACGGCGGAAACGGAACAGCTACGCTTGCGTTTGtt                                                                                                               <  1:287030/140‑1 (MQ=255)
                        gTTGACGGTGAAAGTGCAGCCGCTGGCGAAACCGTCAGAAGACGAGTTGCATCAGGCGATGCTTAATGGCATCCGTGATAAAGGTTTAAGCGTGCTCAACTGGACGGCGGAAACGGAACAGCTACGCt                                                                                                                         <  2:67431/128‑1 (MQ=255)
                        gTTGACGGTGAAAGTGCAGCCGCTGGCGAAACCGTCAGAAGACGAGTTGCATCAGGCGATGCTTAATGGCATCCGTGATAAAGGTTTAAGCGTGCTCAACTGGACGGCGGAAACGGAACAGCTACGCt                                                                                                                         >  1:67431/1‑128 (MQ=255)
                                    aGTGCAGCCGCTGGCGAAACCGTCAGAAGACGAGTTGCATCAGGCGATGCTTAATGGCATCCGTGATAAAGGTTTAAGCGTGCTCAACTGGACGGCGGAAACGGAACAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGt                                                                                                 <  1:745565/140‑1 (MQ=255)
                                          gCCGCTGGCGAAACCGTCAGAAGACGAGTTGCATCAGGCGATGCTTAATGGCATCCGTGATAAAGGTTTAAGCGTGCTCAACTGGACGGCGGAAACGGAACAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGc                                                                                           <  1:802213/140‑1 (MQ=255)
                                          gCCGCTGGCGAAACCGTAAGAATACGAGTTGCATCTGGCCATGCTTAATGGCATCCGTGATAAAGGTTTAAGCGTGCTCAACTGGACGGCGGAAACGGAACAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGc                                                                                           <  2:332307/140‑1 (MQ=255)
                                                           cAGAAGACGAGTTGCATCAGGCGATGCTTAATGGCATCCGTGATAAAGGTTTAAGCGTGCTCAACTGGACGGCGGAAACGGAACAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCa                                                                          <  1:467434/140‑1 (MQ=255)
                                                               gacagAGTTGCATCATGCGATGATTCATGGCATCCGTGATAAAGGTTTAAGCGTGCTCAACTGGACGGCGGAAGCGGAACAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCgg                                                                      <  1:387006/136‑1 (MQ=255)
                                                                          aTCAGGCGATGCTTAATGGCATACGTGATAAAGGTTTAAGCGTGCTCAACTGGACGGCGGAAACGGAACAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGATTATGACTGGCCAGCGGTTGATGATGaa                                                           <  1:402474/140‑1 (MQ=255)
                                                                           tCAGGCGATGCTTAATGGCATCCGTGATAAAGGTTTAAGCGTGCTCAACTGGACGGCGGAAACGGAACAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGTATACGATTGGCCAGCGGTTGATGATGaaa                                                          >  1:362662/1‑140 (MQ=255)
                                                                              ggCGATGCTTAATGGCATCCGTGATAAAGGTTTAAGCGTGCTCAACTGGACGGCGGAAACGGAACAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGtt                                                       >  1:703926/1‑140 (MQ=255)
                                                                               gCGATGCTTAATGGCATCCGTGATAAAGGTTTAAGCGTGCTCAACTGGACGGCGGAAGCGGAACAGCTACGCTTGCGTTTGTTATGCGCCGTAAAGTGGTAGCCAGCATATGACTGGACAGCGGTTGATGAGGaactgtc                                                      >  2:387006/1‑135 (MQ=255)
                                                                                    gCTTAATGGCATCCGCGATAAATGTCTAAGCGTGCTCTACTGGACGGCGGAAACGGAACAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTgg                                                 <  1:254519/140‑1 (MQ=255)
                                                                                               tCCGTGATAAAGGTTTAAGCGTGCTCAACTGGACGGCGGAAACGGAACAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGaa                                      >  2:716967/1‑140 (MQ=255)
                                                                                               tCCGTGATAAAGGTTTAAGCGTGCTCAACTGGACGGCGGAAACGGAACAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGaa                                      >  1:567290/1‑140 (MQ=255)
                                                                                                           gTTTAAGCGTGCTCAACTGGACGGCGGAAACGGAACAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAATTTTATTGGCAGCGCTGGAAACGTGGctgctg                          >  2:357600/1‑140 (MQ=255)
                                                                                                                     gCTCAACTGGACGGCGGAAACGGAACAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGa                >  1:775577/1‑140 (MQ=255)
                                                                                                                      cTCAACTGGACGGCGGAAACGGAACAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGAc               >  1:777632/1‑140 (MQ=255)
                                                                                                                      cTCAACTGGACGGCGGAAACGGAACAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGAc               >  1:609933/1‑140 (MQ=255)
                                                                                                                      cTCAACTGGACGGCGGAAACGGAACAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTAGTGGCAGCGCTGTAAACGTGGCGGCTGCCACATATGAc               >  1:303863/1‑140 (MQ=255)
                                                                                                                       tCAACTGGACGGCGGAAACGGAACAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACt              <  1:187469/140‑1 (MQ=255)
                                                                                                                           cTGGACGGCGGAAACGGAACAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCg          >  2:12038/1‑140 (MQ=255)
                                                                                                                           cTGGACGGCGGAAACGGAACAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCg          <  2:501716/140‑1 (MQ=255)
                                                                                                                                cggcggAAACGGAACAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACAt     >  2:103440/1‑140 (MQ=255)
                                                                                                                                   cggaaacggaaCAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCa  >  1:54365/1‑140 (MQ=255)
                                                                                                                                   cggaaacggaaCAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCa  >  1:444270/1‑140 (MQ=255)
                                                                                                                                        |                                                                                                                                      
CTGGCGTCGGCTACAAATCGGTCAGTTGACGGTGAAAGTGCAGCCGCTGGCGAAACCGTCAGAAGACGAGTTGCATCAGGCGATGCTTAATGGCATCCGTGATAAAGGTTTAAGCGTGCTCAACTGGACGGCGGAAGCGGAACAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCA  >  NC_000913/163856‑164126

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