Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,385,477 Δ1 bp 13.5% coding (401/576 nt) PP_2093 ← two‑component system response regulator NasT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,385,4770T.13.5% 93.4 / 12.0 37coding (401/576 nt)PP_2093two‑component system response regulator NasT
Reads supporting (aligned to +/- strand):  ref base T (16/15);  new base . (3/2);  total (19/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.84e-01

TTCGGCAACCTGGATCAGCTTCTGCTGCCGGCTCATGGCCTGGCGACGCATCAGGGTGTAGGCCTGCTCTTCGTTGCAGTCCTTCATTTTCATCAGCAAGCCCTTGGCCTGCTCGATGCGCTTGCGTTCGGCCAATTGCTGGTCGCGCGCCAGCAGCTGTGCCTTCAGTGCCTGGTCGCTTTCGAAGCGGGCCATGGCCACGTCCAGGATCGGCTGCAGCCGCGCAGCATGGATACCTTCGACGATATAGGCACTGACCCCTGCCTGGATCGCCTGGCG  >  NC_002947/2385335‑2385613
                                                                                                                                               |                                                                                                                                         
ttCGGCAACCTGGATCAGCTTCTGCTGCCGGCTCATGGCCTGGCGACGCATCAGGGTGTAGGCCTGCTCTTCGTTGCAGTCCTTCATTTTCACTCAGCAAGCCCTTGGCCTGCTCGATGCGCTTGCGTTCGGCCAATTGCTGGTcgcgc                                                                                                                                      <  1:209162/149‑1 (MQ=255)
ttCGGCAACCTGGATCAGCTTCTGCTGCCGGCTCATGGCCTGGCGACGCATCAGGGTGTAGGCCTGCTCTTCGTTGCAGTCCTTCATTTTCATCAGCAAGCCCTTGGCCTGCTCGATGCGCTTGCGTTCGGCCAATTGCTGGTcgcgcg                                                                                                                                     >  3:25841/1‑149 (MQ=255)
  cGGCAACCTGGATCAGCTTCTGCTGCCGGCTCATGGCCTGGCGACGCATCAGGGTGTAGGCCTGCTCTTCGTTGCAGTCCTTCATTTTCATCAGCAAGCCCTTGGCCTGCTCGATGCGCTTGCGTTCGGCCAATTGCTGGTCgcgcgc                                                                                                                                   <  1:39836/148‑1 (MQ=255)
            cagcAGCTTCTGCTGCCGGCTCATGGCCTGGCGACGCATCAGGGTGTAGGCCTGCTCTTCGTTGCAGTCCTTCATTTTCATCAGCAAGCCCTTGGCTTGCTCGATGCGCTTGCGTTCGGCCAATTGCTGGTCGCGCGCCAGCAGCtgt                                                                                                                         <  1:67315/145‑1 (MQ=255)
              tCAGCTTCTGCTGCCGGCTCATGGCCTGGCGACGCATCAGGGTGTAGGCCTGCTCTTCGTTGCAGTCCTTCATTTTCATCAGCAAGCCCTTGGCCTGCTCGATGCGCTTGCGTTCGGCCAATTGCTGGTCGCGCGCCAGCAGCTGTGcc                                                                                                                      <  6:3380/149‑1 (MQ=255)
                  cTTCTGCTGCCGGCTCATGGCCTGGCGACGGATCAGGGTGTAGGCCTGCTCTTCGTTGCAGTCCTTCATTTTCATCAGCAAGCCCTTGGCCTGCTCGATGCGCTTGCGTTCGGCCAATTGCTGGTCGCGCGCCAGCAGCTGTGCCtt                                                                                                                    >  1:22780/1‑147 (MQ=255)
                  cTTCTGCTGCCGGCTCATGGCCTGGCGACGCATCAGGGTGTAGGCCTGCTCTTCGTTGCAGTCCTTCATTTTCATCAGCAAGCCCTTGGCCTGCTCGATGCGCTTGCGTTCGGCCGATTGCTGGCGCGCGACcagcaattggc                                                                                                                        >  3:88378/1‑137 (MQ=255)
                  cTTCTGCTGCCGGCTCATGGCCTGGCGACGCATCAGGGTGTAGGCCTGCTCTTCGTTGCAGTCCTTCATTTTCATCAGCAAGCCCTTGGCCTGCTCGATGCGCTTGCGTTCGGCCGATTGCTGGCGCGCGACcagcaattggc                                                                                                                        <  4:88378/143‑7 (MQ=255)
                     ctgctgCCGGCTCATGGCCTGGCGACGCATCAGGGTGTAGGCCTGCTCTTCGTTGCAGTCCTTCATTTTCATCAGCAAGCCCTTGGCCTGCTCGATGCGCTTGCGTTCGGCCAATTGCTGGTCGCGCGCCAGCAGCTGTGCCTTCAGTg                                                                                                               >  8:17671/1‑149 (MQ=255)
                        ctgcCGGCTCATGGCCTGGCGACGCATCAGGGTGTAGGCCTGCTCTTCGTTGCAGTCCTTCATTTTCATCAGCAAGCCCTTGGCCTGCTCGATGCGCTTGCGTTCGGCCAATTGCTGGTCGCGCGCCAGCAGCTGTGCCTTCAGTGcc                                                                                                             <  3:15637/148‑1 (MQ=255)
                        ctgcCGGCTCATGGCCTGGCGACGCATCAGGGTGTAGGCCTGCTCTTCGTTGCAGTCCTTCATTTTCATCAGCAAGCCCTTGGCCTGCTCGATGCGCTTGCGTTCGGCCAATTGCTGGTCGCGCGCCAGCAGCTGTGCCTTCAGTGCCt                                                                                                            <  7:17671/149‑1 (MQ=255)
                             ggCTCATGGCCTGGCGACGCATCAGGGTGTAGGCCTGCTCTTCGTTGCAGTCCTTCATTTTCATCAGCAAGCCCTTGGCCTGCTCGATGCGCTTGCGTTCGGCCAATTGCTGGCGCGCGACcagcaat                                                                                                                            <  8:195425/128‑3 (MQ=255)
                             ggCTCATGGCCTGGCGACGCATCAGGGTGTAGGCCTGCTCTTCGTTGCAGTCCTTCATTTTCATCAGCAAGCCCTTGGCCTGCTCGATGCGCTTGCGTTCGGCCAATTGCTGGCGCGCGACcagcaat                                                                                                                            >  7:195425/1‑126 (MQ=255)
                                          gCGACGCATCAGGGTGTAGGCCTGCTCTTCGTTGCAGTCCTTCATTTTCATCAGCAAGCCCTTGGCCTGCTCGATGCGCTTGCGTTCGGCCAATTGCTGGTCGCGCGCCAGCAGCTGTGCCTTCAGTGCCTGGTCGCTTTCGAAGCgg                                                                                           <  5:32360/148‑1 (MQ=255)
                                            gACGCATCAGGGTGTAGGCCTGCTCTTCGTTGCAGTCCTTCATTTTCATCAGCAAGCCCTTGGCCTGCTCGATGCGCTTGCGTTCGGCCAATTGCTGGTCGCGCGCCAGCAGCTGTGCCTTCAGTGCCTGGTCGCTTTCGAAGCGGGc                                                                                         >  5:172014/1‑148 (MQ=255)
                                               gCATCAGGGTGTAGGCCTGCTCTTCGTTGCAGTCCTTCATTTTCATCAGCAAGCCCTTGGCCTGCTCGATGCGCTTGCGTTCGGCCAATTGCTGGTCGCGCGCCAGCAGCTGTGCCTTCAGTGCCTGGTCGCTTTCGAAGCGGGCCATg                                                                                     <  8:84497/149‑1 (MQ=255)
                                               gCATCAGGGTGTAGGCCTGCTCTTCGTTGCAGTCCTTCATTTTCATCAGCAAGCCCTTGGCCTGCTCGATGCGCTTGCGTTCGGCCAATTGCTGGTCGCGCGCCAGCAGCTGTGCCTTCAGTGCCTGGTCGCTTTCGAAGCGGGCCATg                                                                                     <  6:204665/149‑1 (MQ=255)
                                                              ccTGCTCTTCGTTGCAGTCCTTCATTTTCATCAGCAAGCCCTTGGCCTGCTCGATGCGCTTGCGTTCGGCCAATTGCTGGTCGCGCGCCAGCAGCTGTGCCTTCAGTGCCTGGTCGCTTTCGAAGCGGGCCATGGCCACGTCCAGGATc                                                                      >  4:149845/1‑149 (MQ=255)
                                                                                  ttCATTTTCATCAGCAAGCCCTTGGCCTGCTCGATGCGCTTGCGTTCGGCCAATTGCTGGTCGCGCGCCAGCAGCTGTGCCTTCAGTGCCTGGTCGCTTTCGAAGCGGGCCATGGCCACGTCCAGGATCGGCTGCAGCCGCGCAGCATg                                                  >  6:146842/1‑149 (MQ=255)
                                                                                      ttttCATCAGCAAGCCCTTGGCCTGCTCGATGCGCTTGCGTTCGGCCAATTGCTGGTCGCGCGCCAGCAGCTGTGCCTTCAGTGCCTGGTCGCTTTCGAAGCGGGCCATGGCCACGTCCAGGATCGGCTGCAGCCGCGCAGCATGGATa                                              >  8:222323/1‑149 (MQ=255)
                                                                                         tcatcaGCAAGCCCTTGGCCTGCTCGATGCGCTTGCGTTCGGCCAATTGCTGGTCGCGCGCCAGCAGCTGTGCCTTCAGTGCCTGGTCGCTTTCGAAGCGGGCCATGGCCACGTCCAGGATCGGCTGCAGCCGCGCAGCATGGATACCt                                           >  6:59415/1‑149 (MQ=255)
                                                                                         tcatcaCCAAGCCCTTGGCCAGCTCGATGCGCTTGCGTTCGGCCAATTGCTGGTCGCGCGCCAGCAGCTGAGCCTTCAGTGGCTGGTCGATTTCGAAGCGGTCGATGGCCACGTCCAGGATCGGCTGCAGGCGCGCAGCTTGGATACCt                                           >  5:145000/1‑149 (MQ=255)
                                                                                           atcaGCAAGCCCTTGGCCTGCTCGATGCGCTTGCGTTCGGCCAATTGCTGGTCGCGCGCCAGCAGCTGTGCCTTCAGTGCCTGGTCGCTTTCGAAGCGGGCCATGGCCACGTCCAGGATCGGCTGCAGCCGCGCAGCATGGATACCTTc                                         <  5:146842/149‑1 (MQ=255)
                                                                                               aGCAAGCCCTTGGCCTGCTCGATGCGCTTGCGTTCGGCCAATTGCTGGTCGCGCGCCAGCAGCTGTGCCTTCAGTGCCTGGTCGCTTTCGAAGCGGGCCATGGCCACGTCCAGGATCGGCTGCAGCCGCGCAGCATGGATACCTTcga                                       <  7:211462/148‑1 (MQ=255)
                                                                                               aGCAAGCCCTTGGCCTGCTCGATGCGCTTGCGTTCGGCCAATTGCTGGTCGCGCGCCAGCAGCTGTGCCTTCAGTGCCTGGTCGCTTTCGAAGCGGGCCATGGCCACGTCCAGGATCGGCTGCAGCCGCGCAGCATGGATACCTTcga                                       <  7:222323/148‑1 (MQ=255)
                                                                                                gCAAGCCCTTGGCCTGCTCGATGCGCTTGCGTTCGGCCAATTGCTGGTCGCGCGCCAGCAGCTGTGCCTTCAGTGCCTGGTCGCTTTCGAAGCGGGCCATGGCCACGTCCAGGATCGGCTGCAGCCGCGCAGCATGGATACCTTcgac                                      >  2:251990/1‑148 (MQ=255)
                                                                                                      ccTTGGCCTGCTCGATGCGCTTGCGTTCGGCCAATTGCTGGTCGCGCGCCAGCAGCTGTGCCTTCAGTGCCTGGTCGCTTTCGAAGCGGGCCATGGCAACGTCCAGGATCGGCTGCAGCCGCGCAGCATGGATACCTTCGACGATATAg                               >  6:51505/1‑149 (MQ=255)
                                                                                                           gCCTGCTCGATGCGCTTGCGTTCGGCCAATTGCTGGTCGCGCGCCAGCAGCTGTGCCtt                                                                                                                    >  1:137873/1‑59 (MQ=255)
                                                                                                           gCCTGCTCGATGCGCTTGCGTTCGGCCAATTGCTGGTCGCGCGCCAGCAGCTGTGCCtt                                                                                                                    <  2:137873/59‑1 (MQ=255)
                                                                                                                        gcTTGCGTTCGGCCAATTGCTGGTCGCGCGCCAGCAGCTGTGCCTTCAGTGCCTGGTCGCTTTCGAAGCGGGCCATGGCCACGTCCAGGATCGGCTGCAGCCGCGCAGCATGGATACCTTCGACGATATAGGCACTGACCcctgcctg              >  7:174010/1‑148 (MQ=255)
                                                                                                                            gCGTTCGGCCAATTGCTGGTCGCGCGCCAGCAGCTGTGCCTTCAGTGCCTGGTCGCTTTCGAAGCGGGCCATGGCCACGTCCAGGATCGGCTGCAGCCGCGCAGCATGGATACCTTCGACGATATAGGCACTGACCCCTGCCTGGATCg         >  4:215785/1‑149 (MQ=255)
                                                                                                                             cGTTCGGCCAATTGCTGGTCGCGCGCCAGCAGCTGTGCCTTCAGTGCCTGGTCGCTTTCGAAGCGGGCCATGGCCACGTCCAGGATCGGCTGCAGCCGCGCAGCATGGATACCTTCGACGATATAGGCACTGACCCCTGCCTGGATCGc        <  1:251990/149‑1 (MQ=255)
                                                                                                                                  ggCCATTTGCTGGACGCGCGCGAGCAGGTGTGCCTTCAGTGCCTGGTGTCTTTCGAAGCGGGCCATGGCCACGTCCAGGATCCGCTGCAGCCGCGCAGCATGGATACCTTCGAGGATATAGCCACTGACCCCTGCCTGGATCGCCTGgc   <  6:145000/149‑1 (MQ=255)
                                                                                                                                  ggCCAATTGCTGGTCGCGCGCCAGCAGCTGTGCCTTCAGTGCCTGGTCGCTTTCGAAGCGGGCCATGGCCACGTCCAGGATCGGCTGCAGCCGCGCAGCATGGATACCTTCGACGATATAGGCACTGACCCCTGCCTGGATCGCCTGgc   <  3:149845/149‑1 (MQ=255)
                                                                                                                                   gCCAATTGCTGGTCGCGCGCCAGCAGCTGTGCCTTCAGTGCCTGGTCGCTTTCGAAGCGGGCCATGGCCACGTCCAGGATCGGCTGCAGCCGCGCAGCATGGATACCTTCGACGATATAGGCACTGACCCCTGCCTGGATCGCCTGgcg  >  4:221035/1‑149 (MQ=255)
                                                                                                                                   gCCAATTGCTGGTCGCGCGCCAGCAGCTGTGCCTTCAGTGCCTGGTCGCTTTCGAAGCGGGCCATGGCCACGTCCAGGATCGGCTGCAGCCGCGCAGCATGGATACCTTCGACGATATAGGCACTGACCCCTGCCTCGATCGCCTGgcg  >  8:125258/1‑149 (MQ=255)
                                                                                                                                      agctGCTGGCGCGCGACCAGCAGCTGTGCCTTCAGTGCCTGGTCGCTTTCGAAGCGGGCCATGGCCACGTCCAGGATCGGCTGCAGCCGCGCAGCATGGATACCTTCGACGATATAGGCACTGACCCCTGCCTGGATCGCCTGgc   >  6:106587/4‑145 (MQ=255)
                                                                                                                                               |                                                                                                                                         
TTCGGCAACCTGGATCAGCTTCTGCTGCCGGCTCATGGCCTGGCGACGCATCAGGGTGTAGGCCTGCTCTTCGTTGCAGTCCTTCATTTTCATCAGCAAGCCCTTGGCCTGCTCGATGCGCTTGCGTTCGGCCAATTGCTGGTCGCGCGCCAGCAGCTGTGCCTTCAGTGCCTGGTCGCTTTCGAAGCGGGCCATGGCCACGTCCAGGATCGGCTGCAGCCGCGCAGCATGGATACCTTCGACGATATAGGCACTGACCCCTGCCTGGATCGCCTGGCG  >  NC_002947/2385335‑2385613

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: