Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 964,152 G→A 9.1% intergenic (+40/‑168) ihfB → / → ycaI integration host factor (IHF), DNA‑binding protein, beta subunit/ComEC family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913964,1520GA 9.1% 128.6 / 4.1 44intergenic (+40/‑168)ihfB/ycaIintegration host factor (IHF), DNA‑binding protein, beta subunit/ComEC family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (25/15);  new base A (2/2);  total (27/17)
Fisher's exact test for biased strand distribution p-value = 6.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.35e-01

ACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAATTCCCACCAGACAACCGCTCAACAAAGTTGCA  >  NC_000913/964005‑964293
                                                                                                                                                   |                                                                                                                                             
aCCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCgg                                                                                                                                              <  1:382681/149‑1 (MQ=255)
   ggACGTAATCCGAAGACGGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTg                                                                                                                                           <  1:19895/149‑1 (MQ=255)
      cGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCCGACAGGTGCtt                                                                                                                                        <  4:243868/149‑1 (MQ=255)
       gTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCttt                                                                                                                                       >  1:422283/1‑149 (MQ=255)
         aaTCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCCGACAGGTGCttttc                                                                                                                                     <  4:126330/149‑2 (MQ=255)
           tCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTc                                                                                                                                   >  1:256331/1‑149 (MQ=255)
               aaGACTGGCGATAAAGTAGAACTGGAAGGAAAATACATTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCAttt                                                                                                                               <  2:406134/149‑1 (MQ=255)
                 gACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTtct                                                                                                                             <  4:46772/149‑1 (MQ=255)
                    tGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTa                                                                                                                          >  3:256056/1‑149 (MQ=255)
                      gCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGc                                                                                                                                          >  1:343873/1‑131 (MQ=255)
                      gCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGc                                                                                                                                          <  2:343873/131‑1 (MQ=255)
                       cGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATc                                                                                                                       <  3:464872/149‑1 (MQ=255)
                        gATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCtt                                                                                                                                        <  2:524235/131‑1 (MQ=255)
                        gATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCtt                                                                                                                                        >  1:524235/1‑131 (MQ=255)
                                 gAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATctggaactgga                                                                                                             >  1:158775/1‑149 (MQ=255)
                                  aaCTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATctggaactggaa                                                                                                            <  4:54758/149‑1 (MQ=255)
                                       gaaggaaAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGc                                                                                                       >  1:283373/1‑149 (MQ=255)
                                                    ttCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGtttt                                                                                          <  2:543752/149‑1 (MQ=255)
                                                     tCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTg                                                                                         <  4:256056/149‑1 (MQ=255)
                                                        tCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAAc                                                                                      >  3:399233/1‑149 (MQ=255)
                                                        tCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGAGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAAc                                                                                      >  1:449790/1‑149 (MQ=255)
                                                          aCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAg                                                                                    >  1:69698/1‑149 (MQ=255)
                                                            tttAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGtt                                                                                  >  4:533111/1‑149 (MQ=255)
                                                            tttAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGtt                                                                                  >  2:323381/1‑149 (MQ=255)
                                                                 aCCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCAc                                                                             >  2:57294/1‑149 (MQ=255)
                                                                    tGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCt                                                                          >  1:500952/1‑149 (MQ=255)
                                                                          aGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGtt                                                                    >  1:129842/1‑149 (MQ=255)
                                                                                    gATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATtcttct                                                          >  4:368048/1‑149 (MQ=255)
                                                                                     aTCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTg                                                         >  4:234790/1‑149 (MQ=255)
                                                                                            cAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGc                                                  >  1:365841/1‑149 (MQ=255)
                                                                                                     cGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGaa                                         <  2:500952/149‑1 (MQ=255)
                                                                                                      ggTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAAc                                        <  2:69698/149‑1 (MQ=255)
                                                                                                       gTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCTGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGaaag                                       >  4:502983/1‑147 (MQ=255)
                                                                                                             ttttttACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCaaaa                                  >  4:415285/1‑148 (MQ=255)
                                                                                                             ttttttACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCaaaa                                  <  3:415285/148‑1 (MQ=255)
                                                                                                                     tCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAATTCccacca                         >  3:571775/1‑149 (MQ=255)
                                                                                                                     tCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAATTCccacca                         <  3:234790/149‑1 (MQ=255)
                                                                                                                     tCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAAATCTTCTGAATATGCCTCGGGGAACGCAAAATTCctacca                         >  4:365360/1‑149 (MQ=255)
                                                                                                                       aaaCTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAATTCCCACCAGa                       >  1:541176/1‑149 (MQ=255)
                                                                                                                       aaaCTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAATTCCCACCAGa                       >  1:56124/1‑149 (MQ=255)
                                                                                                                                     aaaaaGCACCCGACAGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAATTCCCACCAGACAACCGCTcaacaa         >  1:263180/2‑149 (MQ=255)
                                                                                                                                      aaaaGCACCCGACAGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAATTCCCACCAGACAACCGCTCAACaaa        >  2:110181/1‑149 (MQ=255)
                                                                                                                                            aCCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAATTCCCACCAGACAACCGCTCAACAAAGTTGca  <  2:158775/149‑1 (MQ=255)
                                                                                                                                                 tCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGttt                                                                                           <  4:272111/55‑1 (MQ=255)
                                                                                                                                                 tCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGttt                                                                                           >  3:272111/1‑55 (MQ=255)
                                                                                                                                                   |                                                                                                                                             
ACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAATTCCCACCAGACAACCGCTCAACAAAGTTGCA  >  NC_000913/964005‑964293

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