Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,842,438 A→C 9.2% intergenic (+24/+135) ascB → / ← hycI cryptic 6‑phospho‑beta‑glucosidase/protease involved in processing C‑terminal end of HycE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,842,4380AC 9.2% 180.3 / 8.7 54intergenic (+24/+135)ascB/hycIcryptic 6‑phospho‑beta‑glucosidase/protease involved in processing C‑terminal end of HycE
Reads supporting (aligned to +/- strand):  ref base A (23/26);  new base C (3/2);  total (26/28)
Fisher's exact test for biased strand distribution p-value = 6.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

ACGCTACGGCTTTGTCTTTGTTGACCGTGACGACGCAGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTT  >  NC_000913/2842291‑2842583
                                                                                                                                                   |                                                                                                                                                 
aCGCTACGGCTTTGTCTTTGTTGACCGTGACGACGCAGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGaa                                                                                                                                                  <  1:379202/149‑1 (MQ=255)
 cGCTACGGCTTTGTCTTTGTTGACCGTGACGACGCAGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAAc                                                                                                                                                 <  4:597545/149‑1 (MQ=255)
 cGCTACGGCTTTGTCTTTGTTGACCGTGACGACGCAGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAAc                                                                                                                                                 <  4:491210/149‑1 (MQ=255)
    tACGGCTTTGTCTTTGTTGACCGTGACGACGCAGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGcc                                                                                                                                              >  1:1262304/1‑149 (MQ=255)
    tACGGCTTTGTCTTTGTTGACCGTGACGACGCAGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGcc                                                                                                                                              <  2:791157/149‑1 (MQ=255)
          tttgtctttgtTGACCGTGACGACGCAGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATcc                                                                                                                                        <  1:389544/149‑1 (MQ=255)
          tttgtctttgtTGACCGTGACGACGCAGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATcc                                                                                                                                        >  2:389544/1‑149 (MQ=255)
          tttgtctttgtTGACCGTGACGACGCAGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATAAGGCGTTCACGccgcat                                                                                                                                          >  3:792859/1‑143 (MQ=255)
              tctttgtTGACCGTGACGACGCAGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATc                                                                                                                                         >  2:971362/1‑144 (MQ=255)
              tctttgtTGACCGTGACGACGCAGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATc                                                                                                                                         <  1:971362/144‑1 (MQ=255)
                tttgtTGACCGTGACGACGCAGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCt                                                                                                                                             <  1:209257/138‑1 (MQ=255)
                tttgtTGACCGTGACGACGCAGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCt                                                                                                                                             >  2:209257/1‑138 (MQ=255)
                tttgtTGACCGTGACGACGCAGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                                    >  2:72469/1‑147 (MQ=255)
                tttgtTGACCGTGACGACGCAGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                                    <  1:72469/147‑1 (MQ=255)
                tttgtTGACCGTGACGACGCAGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTa                                                                                                                                  <  2:1262304/149‑1 (MQ=255)
                   gttgACCGTGACGACGCAGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAt                                                                                                                               <  2:155693/149‑1 (MQ=255)
                      gACCGTGACGACGCAGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTc                                                                                                                            <  2:271012/149‑1 (MQ=255)
                             acgacgCAGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGaa                                                                                                                     >  1:616575/1‑149 (MQ=255)
                                      gCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAgg                                                                                                            >  2:1400118/1‑149 (MQ=255)
                                        aaCGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGcc                                                                                                          >  3:542720/1‑149 (MQ=255)
                                        aaCGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGcc                                                                                                          <  2:616575/149‑1 (MQ=255)
                                            gTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGAt                                                                                                      >  3:1243863/1‑149 (MQ=255)
                                                   gACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAagcgca                                                                                               <  4:542720/149‑1 (MQ=255)
                                                   gACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAagcgca                                                                                               <  4:1243863/149‑1 (MQ=255)
                                                          aCGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATAAGGCGTTCACGccgcatc                                                                                                                                         <  4:769065/100‑6 (MQ=38)
                                                          aCGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATAAGGCGTTCACGccgcatc                                                                                                                                         >  3:769065/1‑95 (MQ=38)
                                                               tAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAagc                                                                                                  <  4:844882/134‑1 (MQ=255)
                                                               tAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAagc                                                                                                  >  3:844882/1‑134 (MQ=255)
                                                                   aaaTCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGa                                                                                                               <  1:1205808/117‑1 (MQ=255)
                                                                   aaaTCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGa                                                                                                               >  2:1205808/1‑117 (MQ=255)
                                                                    aaTCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCAt                                                                                         <  1:308614/138‑1 (MQ=255)
                                                                    aaTCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCAt                                                                                         >  2:308614/1‑138 (MQ=255)
                                                                                        aaaGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAagcgca                                                                                               <  3:527149/112‑1 (MQ=255)
                                                                                        aaaGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAagcgca                                                                                               >  4:527149/1‑112 (MQ=255)
                                                                                        aaaGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTgg                                                          <  3:254293/149‑1 (MQ=255)
                                                                                               ttGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAagc                                                                                                  <  3:297093/102‑1 (MQ=255)
                                                                                               ttGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAagc                                                                                                  >  4:297093/1‑102 (MQ=255)
                                                                                                        aTGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGaa                                          <  1:1400118/149‑1 (MQ=255)
                                                                                                            ggAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATc                                                                                        <  2:429581/99‑1 (MQ=255)
                                                                                                            ggAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATc                                                                                        >  1:429581/1‑99 (MQ=255)
                                                                                                                 aTTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGaa                                                                                                              <  1:1384549/72‑1 (MQ=255)
                                                                                                                 aTTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGaa                                                                                                              >  2:1384549/1‑72 (MQ=255)
                                                                                                                            gTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCTAGAACa                                        >  2:955406/1‑131 (MQ=255)
                                                                                                                            gTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACa                                        <  1:955406/131‑1 (MQ=255)
                                                                                                                             tAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTc                                >  4:622025/1‑138 (MQ=255)
                                                                                                                             tAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTc                                <  3:622025/138‑1 (MQ=255)
                                                                                                                              aaCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTaa                    >  2:1074459/1‑149 (MQ=255)
                                                                                                                                         ataaggcgTTCACGCCGCATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCa                                                 <  2:730090/105‑1 (MQ=38)
                                                                                                                                         ataaggcgTTCACGCCGCATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCa                                                 >  1:730090/5‑109 (MQ=255)
                                                                                                                                           gcggcgTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTActc       >  4:725096/1‑149 (MQ=255)
                                                                                                                                             ggcgTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCGTTTGGCATCAGAGCTGGTTCAcc                                                    >  4:1052613/1‑102 (MQ=255)
                                                                                                                                             ggcgTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCAcc                                                    <  3:1052613/102‑1 (MQ=255)
                                                                                                                                              gcgTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCaa                                             >  2:977174/1‑108 (MQ=255)
                                                                                                                                              gcgTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCaa                                             <  1:977174/108‑1 (MQ=255)
                                                                                                                                                gTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTcttctt  >  4:1156961/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
ACGCTACGGCTTTGTCTTTGTTGACCGTGACGACGCAGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTT  >  NC_000913/2842291‑2842583

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