Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,153,794 C→A 10.6% intergenic (+15/‑222) sibB → / → mdtA sRNA antisense regulator of toxic IbsB protein/multidrug efflux system, subunit A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,153,7940CA10.6% 149.3 / 8.3 47intergenic (+15/‑222)sibB/mdtAsRNA antisense regulator of toxic IbsB protein/multidrug efflux system, subunit A
Reads supporting (aligned to +/- strand):  ref base C (20/22);  new base A (2/3);  total (22/25)
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.88e-01

GTAGAGCGGGGTTTCCCCCGCCCTGGTAGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCAAATGTGCCCGTCATTCAGACGATTCCAGACAGTGTTTCATAAT  >  NC_000913/2153650‑2153937
                                                                                                                                                |                                                                                                                                               
gTAGAGCGGGGTTTCCCCCGCCCTGGTAGTCTTAGTAAGCGGGGAATCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCt                                                                                                                                             <  1:140864/149‑1 (MQ=255)
gTAGAGCGGGGTTTCCCCCGCCCTGGTAGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCt                                                                                                                                             <  4:293698/149‑1 (MQ=255)
    agCGGGGTTTCCCCCGCCCTGGTAGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGGGACGCTGaaa                                                                                                                                         <  2:252269/149‑1 (MQ=255)
     gCGGGGTTTCCCCCGCCCTGGTAGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAAt                                                                                                                                        >  1:291666/1‑149 (MQ=255)
                       tGGTAGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTCCCTCCttt                                                                                                                      >  3:3340/1‑149 (MQ=255)
                         gTAGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCg                                                                                                                    >  4:447059/1‑149 (MQ=255)
                               ttAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGc                                                                                                              >  3:189905/1‑149 (MQ=255)
                                  gTAAGCGGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAg                                                                                                           <  4:466356/149‑1 (MQ=255)
                                     aGCGGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAg                                                                                                        <  4:3340/149‑1 (MQ=255)
                                          ggAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTAt                                                                                                   <  3:266976/149‑1 (MQ=255)
                                           gAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGACATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATc                                                                                                  <  4:276726/149‑1 (MQ=255)
                                                ttATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTaa                                                                                             >  4:82196/1‑149 (MQ=255)
                                                ttATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTaa                                                                                             >  4:219501/1‑149 (MQ=255)
                                                ttATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGGGACGCTGaaacg                                                                                                                                       <  4:231693/107‑3 (MQ=25)
                                                ttATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGGGACGCTGaaacg                                                                                                                                       >  3:231693/1‑105 (MQ=25)
                                                       tAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAg                                                                                      >  2:242241/1‑149 (MQ=255)
                                                           agCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATaa                                                                                  >  3:31896/1‑149 (MQ=255)
                                                            gCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTACCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATCCCCGCAGGCTATCGCTAACGCGtcgctacg                                                                                 >  1:341695/1‑142 (MQ=21)
                                                             caccacGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATaaaa                                                                                <  3:82196/149‑1 (MQ=255)
                                                                cacGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCTGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAg                                                                             >  3:307054/1‑149 (MQ=255)
                                                                 acGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGt                                                                            <  4:189905/149‑1 (MQ=255)
                                                                    atgatgAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTc                                                                         <  1:242241/149‑1 (MQ=255)
                                                                      gatgaGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCt                                                                       <  2:291666/149‑1 (MQ=255)
                                                                      gatgaGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCt                                                                       <  3:447059/149‑1 (MQ=255)
                                                                        tgaGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAAc                                                                                            <  1:38570/126‑1 (MQ=255)
                                                                        tgaGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAAc                                                                                            >  2:38570/1‑126 (MQ=255)
                                                                          aGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTgtgt                                                                   <  4:443295/149‑1 (MQ=255)
                                                                           gTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGtt                                                                  >  2:246956/1‑149 (MQ=255)
                                                                                       tGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGc                                                      <  4:31896/149‑1 (MQ=255)
                                                                                           ccTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGc                                                  >  3:5168/1‑149 (MQ=255)
                                                                                            cTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTc                                                                         <  2:300486/125‑1 (MQ=255)
                                                                                            cTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTc                                                                         >  1:300486/1‑125 (MQ=255)
                                                                                             tttCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTc                                                <  2:332606/149‑1 (MQ=255)
                                                                                                  ttatttatGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTc                                                                         <  4:196533/119‑1 (MQ=255)
                                                                                                  ttatttatGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTc                                                                         >  3:196533/1‑119 (MQ=255)
                                                                                                  ttatttatGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCaaa                                           >  3:17276/1‑149 (MQ=255)
                                                                                                        atGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGttatta                                                              <  1:52872/124‑1 (MQ=255)
                                                                                                        atGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGttatta                                                              >  2:52872/1‑124 (MQ=255)
                                                                                                             cccTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCAAATGTGCCCGTCa                                >  1:424164/1‑149 (MQ=255)
                                                                                                                    tCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCAAATGTGCCCGTCATTCAGAc                         <  4:307054/149‑1 (MQ=255)
                                                                                                                        gAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCAAATGTGCCCGTCATTCAGACGAtt                     <  2:30723/149‑1 (MQ=255)
                                                                                                                             gCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCAAATGTGCCCGTCATTCAGACGATTCCAGa                <  3:219501/149‑1 (MQ=255)
                                                                                                                               tttCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCAAATGTGCCCGTCATTCAGACGATTCCAGACa              >  4:307578/1‑149 (MQ=255)
                                                                                                                                      tttCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCAAATGTGCCCGTCATTCAGACGATTCCAGACAGTGTTTc       <  4:17276/149‑1 (MQ=255)
                                                                                                                                      tttCAGCGGGACGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGAt                                                        >  3:353928/1‑100 (MQ=21)
                                                                                                                                      tttCAGCGGGACGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGAt                                                        <  4:353928/100‑1 (MQ=21)
                                                                                                                                           gCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGCCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCAAATGTGCCCGTCATTCAGACGATTCCAGACAGTGTTTCATAAt  >  1:39649/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                               
GTAGAGCGGGGTTTCCCCCGCCCTGGTAGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCAAATGTGCCCGTCATTCAGACGATTCCAGACAGTGTTTCATAAT  >  NC_000913/2153650‑2153937

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