Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,498,896 A→T 8.3% intergenic (+21/‑42) ydcH → / → rimL protein YdcH/ribosomal‑protein‑L12‑serine acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,498,8960AT 8.3% 137.9 / 5.0 49intergenic (+21/‑42)ydcH/rimLprotein YdcH/ribosomal‑protein‑L12‑serine acetyltransferase
Reads supporting (aligned to +/- strand):  ref base A (24/20);  new base T (2/2);  total (27/22)
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.88e-01

AAAGGAAGGTTCCGACGGTCGAGGGTACAATGCGGAAGTGGTCCGCATGAAAAAACAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACATGCTGTTGCAGAAAATCACGTCAAACCTCTTTATCAGTTAATCTGTAAAAATAAAACCTG  >  NC_000913/1498755‑1499038
                                                                                                                                             |                                                                                                                                              
aaaGGAAGGTTCCGACGGTCGAGGGTACAATGCGGAAGTGGTCCGCATGAAAAAACAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGCGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGccagccag                                                                                                                                               >  1:529560/1‑143 (MQ=255)
    gaaggTTCCGACGGTCGAGGGTACAATGCGGAAGTGGTCCGCATGAAAAAACAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAggcgg                                                                                                                                           >  2:230233/1‑143 (MQ=255)
    gaaggTTCCCGCAGCCGGGGGTACAATGCGGAAGTGAGCCGCATCAAAAAACAAAAGCTACAGATAAAAGATGAGACGCTCAAACACCGGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAggcgg                                                                                                                                           <  2:399635/143‑1 (MQ=255)
        gTTCCGACGGTCGAGTGTACAATGCGGAAGTGGTCCGCATGAAAAAACAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCTGGCTGGCGGCtta                                                                                                                                       <  1:55744/143‑2 (MQ=255)
        gTTCCGACGGGCGAGGGTACAATGCGGAAGTGGTCCGCATGAAAAAACAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCTGGCTGGCGGCtta                                                                                                                                       <  1:610931/143‑2 (MQ=255)
              aCGGTCGAGGGTACAATGCGGAAGTGGTCCGCATGAAAAAACAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCGAGCCCGGCGGCTTTTTTaac                                                                                                                                 >  2:457861/1‑143 (MQ=255)
                agtgaggggTACAGCGCGGCAGGGGGTCGGATGAAAAGAGAAAAGCTACAGTTAAAAGGTGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTaacaa                                                                                                                               <  2:373084/137‑1 (MQ=255)
                        gccacaTGCGGAAGTGGTCCGCATGAAAAAACAAAAGCTACAGTTAAAAGATGAGATACTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATgg                                                                                                                       <  1:61248/138‑1 (MQ=255)
                              cgCGGAAGTGGTGAGCATGAAAAAACAAAAGCTAGAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGAc                                                                                                                 <  2:529224/142‑1 (MQ=255)
                                 ggAAGTGGTCCGCATGAAAAAACAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCAGACAGCCAGGCGGCTTATTTAACAACTGCATGGAtcaactag                                                                                                              >  2:149794/1‑136 (MQ=255)
                                       ggTCCGCATGAAAAAACAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATaa                                                                                                        <  2:423115/143‑1 (MQ=255)
                                         tCCGCAAGAAAAAACAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCTGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGa                                                                                                      <  2:889001/143‑1 (MQ=255)
                                          ccGCATGAAAAAACAAAAGCTACAGTTAAAAGATGAGATGCTCAAAAACCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATACGAt                                                                                                     >  2:45138/1‑143 (MQ=255)
                                            gCATGAAAAAACAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCATGCAGCAGGAGAGCGTCAAACAGCTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACAGGAGATAAGATGa                                                                                                   >  2:12967/1‑143 (MQ=255)
                                             cATGAAAAAACAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGAc                                                                                                  >  2:310975/1‑143 (MQ=255)
                                              aTGAAAAAACAAAAGCTACAGTTAAAAGATGAGATGCGCAAAATCCTGCAGCAGGAGAGCGTCAAAGGGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACt                                                                                                 <  1:830353/143‑1 (MQ=255)
                                                 aaaaaaCAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGaa                                                                                              >  1:589113/1‑143 (MQ=255)
                                                   aaaaCAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGTGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAAc                                                                                            >  1:512249/1‑143 (MQ=255)
                                                   aaaaCAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAAc                                                                                            >  1:672962/1‑143 (MQ=255)
                                                   aaaaCAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTCTTTTAACAACTCCCTGGGGCGCCTGGCGATAAGATGACTGTAAc                                                                                            >  2:408444/1‑143 (MQ=255)
                                                      aCAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGAt                                                                                         <  2:589113/143‑1 (MQ=255)
                                                        aaaaGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATaa                                                                                       <  1:193760/143‑1 (MQ=255)
                                                         aaaGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCATAAGCCGCGAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGACATAAGATGACTGAAACGATaaa                                                                                      >  2:415173/1‑143 (MQ=255)
                                                               aCAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGCAACGATAAAAGTAAg                                                                                >  2:541027/1‑143 (MQ=255)
                                                                  gTTAAAAGATAAGACGCTCAAAATGCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGa                                                                             <  1:913491/143‑1 (MQ=255)
                                                                   ttAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATGAGATGACTGAAACGATAAGAGTCAGCGaa                                                                            >  1:475970/1‑143 (MQ=255)
                                                                   ttAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGAGAAGATGGCTGAAACGATAAAAGTAGGGGaa                                                                            >  2:855419/1‑143 (MQ=255)
                                                                    gaaaaGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGGAAGCGAAt                                                                           <  1:537560/142‑1 (MQ=255)
                                                                       aaGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCAc                                                                        <  1:665630/143‑1 (MQ=255)
                                                                       aaGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCAc                                                                        <  1:508815/143‑1 (MQ=255)
                                                                        agatgagatgCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGACTCACt                                                                       >  2:396376/1‑143 (MQ=255)
                                                                         gttgagatgCTCAAACTCCTGCAGCAAGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACtt                                                                      <  2:143887/141‑1 (MQ=255)
                                                                           tgagatgCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGa                                                                    <  2:894219/143‑1 (MQ=255)
                                                                            gagatgCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGaa                                                                   >  1:660946/1‑143 (MQ=255)
                                                                                 gCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATAGAATGAATTACa                                                              >  2:320535/1‑143 (MQ=255)
                                                                                 gCTCAAAAGCCTGGAGCAGGAGGGCGTCAAAGAGGTGTAAACTTGCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACa                                                              <  2:449884/143‑1 (MQ=255)
                                                                                                   ggAGATCGTCAAAGAGGTGTACACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACATGCTGTTGCAGAAAATCa                                            <  1:58018/143‑1 (MQ=255)
                                                                                                     agagCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACATGCTGTTGCAGAAAATCACg                                          <  1:873954/143‑1 (MQ=255)
                                                                                                         cGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACATGCTGGTGCAGAAAATCCCGGCaa                                      >  2:424557/1‑143 (MQ=255)
                                                                                                          gTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACATGCTGTTGCAGAAAATCACGTCaaa                                     <  1:545112/143‑1 (MQ=255)
                                                                                                           tCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACATGCTGTTGCAGAAAATCACGTCAAAc                                    >  1:594753/1‑143 (MQ=255)
                                                                                                              aaGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACATGCTGTTGCAGAAAAACACGTCAAACctc                                 >  2:870310/1‑143 (MQ=255)
                                                                                                               agagGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGCGATAAGATGACTGAAACGATAATAGTAAGCGAATCACTTGAATTACATGCTGTTGCAGAAAATCACGTCAAACctct                                <  1:689975/143‑1 (MQ=255)
                                                                                                                      tAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGGTAAAAGTAAGCGAATCACTTGAATTACATGCTGTTGCCGCAAATCACGTCAAACCCCATTATccg                         >  2:704202/1‑141 (MQ=255)
                                                                                                                        aaCTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGAGAAGATGACTGAAACGATAAAAGTAAGCGAATCACGTGAATTACATGCTGTTGCAGGAAATTCCGTCAACCCTCtttcgccgct                       >  2:779089/1‑136 (MQ=255)
                                                                                                                               cTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAAGTACATGCTGTTGCAGAAAATCACGTCAAACCTCTTTATCCGTTAATCGGt                >  2:479524/1‑143 (MQ=255)
                                                                                                                                aaaGCCGCCTGGCTGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACATGCTGTTGCAGAAAATCACGTCAAACCTCTTTATCAGTTAATCTGTa               >  1:411096/2‑143 (MQ=255)
                                                                                                                                aaaGCCGCCTGGCTGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACATGCTGTTGCAGAAAATCACGTCAAACCTCTCTATCAGTTAATCTGTa               >  2:503464/2‑143 (MQ=255)
                                                                                                                                 aaGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACATGCTGTTGCAGAAAATCACGTCAAACCTCTTTATCAGTTACTCTGTaa              >  1:362545/1‑143 (MQ=255)
                                                                                                                                    ccgccAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACATGCTGTTGCAGAAAATCACGTCAAACCTCTTTATCAGTTAATCTGTaaaaa           >  2:719176/1‑143 (MQ=255)
                                                                                                                                          gccagGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACATGCTGTTGCAGAAAATCACGTCAAACCTCTTTATCAGTTAATCTGTAAAAATaaacc     <  1:674170/143‑3 (MQ=255)
                                                                                                                                             agGTGGTTTTTTTAATAACTGAATGGATTGACTGAAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACATGCTGTTGCAGAAAATCACGTCAAACCTCTTTATCAGTTAATCTGTAAAAATAAAACCTg  <  2:481793/143‑1 (MQ=255)
                                                                                                                                             agGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACATGCTGTTGCAGAAAATCACGTCAAACCTCTTTATCAGTTAATCTGTAAAAATAAAACCTg  <  1:546587/143‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                              
AAAGGAAGGTTCCGACGGTCGAGGGTACAATGCGGAAGTGGTCCGCATGAAAAAACAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACATGCTGTTGCAGAAAATCACGTCAAACCTCTTTATCAGTTAATCTGTAAAAATAAAACCTG  >  NC_000913/1498755‑1499038

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: