Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,187,323 (G)5→6 intergenic (+25/+57) rcnB → / ← yehA periplasmic protein involved in nickel/cobalt export/putative fimbrial adhesin YehA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,187,3181.G95.2% 132.3 / ‑0.2 42intergenic (+20/+62)rcnB/yehAperiplasmic protein involved in nickel/cobalt export/putative fimbrial adhesin YehA
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base G (23/17);  total (24/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 = 1.00e+00

TGCCCGCGCCTGATGCCGGGACGCACTGGACCTATATGGGTGGCGCGTACGTGTTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAACTTTGTTTGTATCGTAGA  >  NC_000913/2187180‑2187450
                                                                                                                                            |                                                                                                                                      
tGCCCGCGCCTGATGCCGGGACGCACTGGACCTATATGGGTGGCGCGTACGTGTTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAg                                                                                                                                       <  1:812048/140‑1 (MQ=255)
tGCCCGCGCCTGATGCCCGGACGCACTGGACGTATGTGGATGGCGCGTTCGTGTTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAg                                                                                                                                       <  2:315775/140‑1 (MQ=255)
  cccGCGCCTGATGCCGGGACGCACTGGACCTATATGGGTGGCGCGTACGTGTTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGTCGGTGAGATTTTTTATAATCGCTAAAAAAAACCCCCTCATCATGAggg                                                                                                                                    >  2:330638/1‑140 (MQ=255)
    cgcgCCTGATGCCGGGACGCACTGGACCTATATGGGTGGCGCGTACGTGTTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAggggg                                                                                                                                  <  1:156845/140‑1 (MQ=255)
                        aCTGGACCTATATGGGTGGCGCGTACGTGTTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGGAAATGCAGACACCTTGTTa                                                                                                              <  2:26360/140‑1 (MQ=255)
                          tGGACCTATATGGGTGGCGCGTACGTGTTGATAAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGGAAATGCAGACACCTTGTTAtt                                                                                                            >  2:599126/1‑140 (MQ=255)
                                  tatGGGTGGCGCGTACGTGTTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGGAAATGCAGACACCTTGTTATTTTttatta                                                                                                    >  2:560808/1‑140 (MQ=255)
                                    tGGGTGGCGCGTACGTGTTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGGAAATGCAGACACCTTGTTATTTTttattatt                                                                                                  >  1:199024/1‑140 (MQ=255)
                                        tGGCGCGTACGTGTTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGcc                                                                                              >  1:728693/1‑140 (MQ=255)
                                         ggCGCGTACGTGTTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCa                                                                                             >  2:786719/1‑140 (MQ=255)
                                         ggCGCGTACGTGTTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCa                                                                                             >  1:770577/1‑140 (MQ=255)
                                          gcgcGTACGTGTTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAAATGCAAACAGCTTGTTATTGTTTATTATTAGCCAc                                                                                            >  1:364239/1‑140 (MQ=255)
                                                    gtTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTct                                                                                  >  2:576792/1‑140 (MQ=255)
                                                       gATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTcttgc                                                                               <  2:770577/140‑1 (MQ=255)
                                                       gATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTcttg                                                                                <  2:22254/139‑1 (MQ=255)
                                                       gATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTcttg                                                                                >  1:22254/1‑139 (MQ=255)
                                                       gATCAGCGACACCGACGGTAAAATCAGTAAATCCTACGACGGTGAGAGTTGTTATCATCGCTAAAAAAAGTCCCCTCATCATGGGGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTcttgc                                                                               <  2:193923/140‑1 (MQ=255)
                                                         tCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTc                                                                                   <  1:609490/134‑1 (MQ=255)
                                                         tCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTc                                                                                   >  2:609490/1‑134 (MQ=255)
                                                         gcAGCGACACCGACGTTAAAATCATAAAATCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTcttgctt                                                                             <  1:233994/139‑1 (MQ=255)
                                                             cGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGt                                                                                    <  2:49129/129‑1 (MQ=255)
                                                             cGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGt                                                                                    >  1:49129/1‑129 (MQ=255)
                                                              gACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGttat                                                                        <  2:728693/140‑1 (MQ=255)
                                                              gACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGttat                                                                        <  2:199024/140‑1 (MQ=255)
                                                                   cGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTcttg                                                                                >  1:611589/1‑127 (MQ=255)
                                                                   cGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTcttg                                                                                <  2:611589/127‑1 (MQ=255)
                                                                   cGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTc                                                                   <  1:846930/140‑1 (MQ=255)
                                                                     aCGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGGAAATGCAGACACCTTCTTATTTTTGATTGTTAGCCACTTGCTCGTCGTGCTTGTTATTAGTCgc                                                                 >  2:233994/1‑139 (MQ=255)
                                                                                 ttAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGGAAATGCAGACCCCTTGTTATTTTTTATTATTAGCCAGTTGGTCGGCTTGCTTGGTATTGGTCGTATTTCACGTTGa                                                     >  2:254895/1‑140 (MQ=255)
                                                                                 ttAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGa                                                     >  2:794421/1‑140 (MQ=255)
                                                                                 ttAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGa                                                     >  1:724435/1‑140 (MQ=255)
                                                                                    aaGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTa                                                  <  1:560808/140‑1 (MQ=255)
                                                                                         tACGACGGTGAGATTGTTTATCCTCCCTACAAAAAGCTCCCTCATCATGAGGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCg                                             <  1:254895/140‑1 (MQ=255)
                                                                                          acgacgGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCAATTGCTCGTATTGGTTGCTATTAGTCTTATTTCACGTTGATTAATGCgg                                            >  2:270087/1‑140 (MQ=255)
                                                                                               ggTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGcc                                       >  1:780780/1‑140 (MQ=255)
                                                                                               ggTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGcc                                       <  1:599126/140‑1 (MQ=255)
                                                                                                 tGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTCGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTc                                     >  1:725933/1‑140 (MQ=255)
                                                                                                   agaTTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCa                                   >  2:446202/1‑140 (MQ=255)
                                                                                                             tcatcGCTAAAAAAAGCCCCCTCATCATGAGGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATATCACGTTCGTTAATGCGGTTGACTCCAGTGCGCCAGa                         >  1:418220/1‑140 (MQ=255)
                                                                                                                   cTAAAAAAAGCCCCCTCATCATGAGGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAAc                   <  2:780780/140‑1 (MQ=255)
                                                                                                                          aaGCCCCCTCATCATTAGGGGGGAACTGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAACtttgttt            <  1:270087/140‑1 (MQ=255)
                                                                                                                              ccccTCATCATGAGGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCTTCCAGTGCGCCAGATTTAACtttgttt            >  1:444260/1‑136 (MQ=255)
                                                                                                                              ccccTCATCATGAGGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAACtttgttt            <  2:444260/136‑1 (MQ=255)
                                                                                                                                 cTCATCATGAGGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGACTCGAGTGCGCCGGATTGAACTTTGTTTGTATCGt     >  2:231664/1‑140 (MQ=255)
                                                                                                                                  tcatcaTGAGGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAACTTTGTTTGTATCGTa    <  1:554706/140‑1 (MQ=255)
                                                                                                                                    atcatGAGGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATACGGTTGCCTCCAGTGCGCCAGATTTAACTTTGTTTGTATCGTAGa  >  1:92480/1‑140 (MQ=255)
                                                                                                                                    atcatGAGTGGGGAAAGCCAGACACCTTGTGTTTTGTTATTATTAGCCTCTTGCTCGTCTTGCTTGTTATTCGTCGTATTTCGCGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAACTTTGTTTGTATCGTAGa  <  2:418220/140‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                      
TGCCCGCGCCTGATGCCGGGACGCACTGGACCTATATGGGTGGCGCGTACGTGTTGATCAGCGACACCGACGGTAAAATCATTAAAGCCTACGACGGTGAGATTTTTTATCATCGCTAAAAAAAGCCCCCTCATCATGAGGGGGAAATGCAGACACCTTGTTATTTTTTATTATTAGCCACTTGCTCGTCTTGCTTGTTATTAGTCGTATTTCACGTTGATTAATGCGGTTGCCTCCAGTGCGCCAGATTTAACTTTGTTTGTATCGTAGA  >  NC_000913/2187180‑2187450

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