Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,498,896 A→T 8.2% 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.2% 136.4 / 4.5 49intergenic (+21/‑42)ydcH/rimLprotein YdcH/ribosomal‑protein‑L12‑serine acetyltransferase
Reads supporting (aligned to +/- strand):  ref base A (25/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 = 3.95e-01

GAAAGGAAGGTTCCGACGGTCGAGGGTACAATGCGGAAGTGGTCCGCATGAAAAAACAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACATGCTGTTGCAGAAAATCACGTCAAACCTCTTTATCAGTTAATCTGTAAAAATAAAACCT  >  NC_000913/1498754‑1499037
                                                                                                                                              |                                                                                                                                             
gAAAGGAAGGTTCCGACGGTCGAGGGTACAATGCGGAAGTGGTCCGCATGAAAAAACAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCgccagcca                                                                                                                                               >  1:147702/1‑143 (MQ=255)
 aaaGGAAGGTTCCGACGGTCGAGGGTACAATGCGGAAGTGGTCCGCATGAAAAAACAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGccagccag                                                                                                                                              >  2:371945/1‑143 (MQ=255)
  aaggaaggTGCCGACGGTCGAGGGTACAATGCGGAAGTGGGCCGCATGAAAAAACAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAgg                                                                                                                                             <  1:466919/143‑1 (MQ=255)
         gTTCCGACGAGCGAGGTTACAACGCGGAAGTAGTCCGCATGAAAAAACAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCTGGCTGGCGGCtta                                                                                                                                      <  2:309228/143‑2 (MQ=255)
           tacGACGGTCGAGGGTACAAGGCGGAAGTGGTCCGCATGAAAAAACAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCttttt                                                                                                                                    <  2:373359/141‑1 (MQ=255)
           tCCGACGGTCGAGGGTACAATGCGGAAGTGGTCCGCATGAAAAAACAAAAGCTACAGATAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCttttt                                                                                                                                    <  1:546177/143‑1 (MQ=255)
           tCCGACGGTCGAGGGTACAATGCGGAAGGAGGCCGCATGAAAAAACAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCttttt                                                                                                                                    <  1:96407/143‑1 (MQ=255)
              gACGGTCGAGGGTACAATGCGGAAGTGGTCCGCATGAAAAAACAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCTGGCTGGCGGCtt                                                                                                                                       <  1:49267/137‑1 (MQ=255)
              gACGGTCGAGGGTACAATGCGGAAGTGGTCCGCATGAAAAAACAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCTGGCTGGCGGCtt                                                                                                                                       >  2:49267/1‑137 (MQ=255)
                             gaTGCGGAAGTGGGCCGCATGAAAAAACAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTg                                                                                                                  <  2:350414/142‑1 (MQ=255)
                                 cGGAAGTGGTCAGTATGATGAAACAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTg                                                                                                              <  1:97130/143‑1 (MQ=255)
                                  ggAAGTGGTCCGCATGAAAAAACAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTgg                                                                                                             >  1:153040/1‑143 (MQ=255)
                                       tGGTCCGCATGAAAAAACAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACCACTGCATGGATTGACTGGAGATa                                                                                                        >  1:261350/1‑143 (MQ=255)
                                          ttcGCATGAATAAACAAAATCTACAGTTATAAGATGAGATGCTCAAATTCCTGCAGCAGGAGAGCGTCAATGAGGTGTTTACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGa                                                                                                     <  1:204765/141‑1 (MQ=255)
                                              cATGAAAAAACAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGAGGAc                                                                                                 >  1:32547/1‑143 (MQ=255)
                                               aTGAAAAAACAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACt                                                                                                >  1:265891/1‑143 (MQ=255)
                                                   aaaaaCAAAAGCTACAGTTAAAAGATGAGAGGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGaaa                                                                                            <  2:124413/143‑1 (MQ=255)
                                                    aaaaCAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAAc                                                                                           >  2:366267/1‑143 (MQ=255)
                                                      aaCAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGAGAAGATGACTGAAACGa                                                                                         >  2:54184/1‑143 (MQ=255)
                                                                 cAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGc                                                                              >  2:220180/1‑143 (MQ=255)
                                                                  aGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCg                                                                             >  1:270681/1‑143 (MQ=255)
                                                                          gatgagatgCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATGACtc                                                                     >  1:129908/1‑142 (MQ=255)
                                                                          gatgagatgCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACtt                                                                     <  2:99068/143‑1 (MQ=255)
                                                                            tgagatgCTTAAAATCCTGTAGTAGGAGAGCTTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGa                                                                   <  2:213729/143‑1 (MQ=255)
                                                                                     cAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACAATCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACATGc                                                          >  2:139052/1‑143 (MQ=255)
                                                                                      aaaaTCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGAGTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACATGca                                                         >  2:405174/1‑142 (MQ=255)
                                                                                         aTCCTGCAGCAGGAGGGCGTCAAAGAGGTGTAAACTTGCCTAAGCCGCCAGACAGGTGACGTTGTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACATGCTGtt                                                      <  2:190265/143‑1 (MQ=255)
                                                                                          tCCTGCAGCAGGCGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACATGCTGTTg                                                     <  2:196267/143‑1 (MQ=255)
                                                                                          tCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACATGCTGTTg                                                     >  2:81675/1‑143 (MQ=255)
                                                                                           ccTGCAGCAGAACAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGTGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACATGCTGTTGc                                                    <  2:193223/143‑1 (MQ=255)
                                                                                              gcagcaGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACATGCTGTTGCAGa                                                 >  2:436828/1‑143 (MQ=255)
                                                                                                agcagGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACATGCTGTTGCAGaaa                                               <  1:371945/143‑1 (MQ=255)
                                                                                                agcagGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATCACATGCTGTTGACGaaa                                               >  2:287390/1‑143 (MQ=255)
                                                                                                 gcagGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAGTCACTTGAATTACATGCTGTTGCAGaaaa                                              >  1:149779/1‑143 (MQ=255)
                                                                                                   agGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACATGCTGTTGCAGAAAATc                                            <  1:462487/143‑1 (MQ=255)
                                                                                                     gagaGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACATGCTGTTGCAGAAAATCAc                                          <  1:85066/143‑1 (MQ=255)
                                                                                                     gagaGCGTCAAAGAGGTGTAAACTGTCCTAAGCCGCGAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAAGCACTTGAATTACATGCTGTTGCAGAAAATAAc                                          >  2:16515/1‑143 (MQ=255)
                                                                                                      agagCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACATGCTGTTGCAGAAAAGCACg                                         >  2:252974/1‑143 (MQ=255)
                                                                                                        agCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACATGCTGTTGCAGAAAATCACGtg                                       >  2:256632/1‑142 (MQ=255)
                                                                                                           gTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACATGCTGTTGCAGAAAATCACGTCaaa                                    <  1:474573/143‑1 (MQ=255)
                                                                                                           gTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACACAACCGATAATTGTTTGCGGATCACTTGAATTACATGCTGTTGCAGAAAATCACGTCaaa                                    >  1:123580/1‑143 (MQ=255)
                                                                                                            tcgaaGATGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGATATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGCATTACATGCTGTTGCAGAAAATCACGTCAAAc                                   <  2:331699/140‑1 (MQ=255)
                                                                                                              aaaGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACATGCTGTTGCAGAAAATCACGTCAAACct                                 >  2:142982/1‑143 (MQ=255)
                                                                                                               aaGAGGTGTAAACTTTCCTAAGCAGCCAGCCAGGCGGCTTTTTTAACACCTACATAGATTGACTCAACATAAGATTACTGAAACGATAACAATAAGCGAATCACTTTAATTACATGCTCTTTCACAAAAACAAATCAAAcgcc                                >  2:361362/1‑140 (MQ=255)
                                                                                                                agagGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACATGCTGTTGCAGAAAATCACGTCAAACctct                               <  2:393153/143‑1 (MQ=255)
                                                                                                                      gtAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACATGCTGTTGCAGAAAATCACGTCAAACCTCTTTATCa                         <  1:130567/143‑1 (MQ=255)
                                                                                                                       tAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGGCTGAAACGATAGAAGTAAGCGAATCAGTTGAATTACATGCTGTTGCAGTCAATCACGTCTAACCTCTTTATCAg                        >  1:105793/1‑143 (MQ=255)
                                                                                                                                 tAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACATGCTGTTGCAGAAAATCACGTCAAACCTCTTTATCAGTTAATCTGTa              >  1:243287/1‑143 (MQ=255)
                                                                                                                                  aaGCCGCCTGGCTGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACATGCTGTTGCAGAAAATCACGTCAAACCTCTTTATCAGTTAATCTGTaa             >  1:415744/1‑143 (MQ=255)
                                                                                                                                     ccgccAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACATGCTGTTGCAGAAAATCACGTCAAACCTCTTTATCAGTTAATCTGTaaaaa          >  2:392892/1‑143 (MQ=255)
                                                                                                                                     ccgccAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACATGCTGTTGCAGAAAATCACGTCAAACCTCTTTATCAGTTAATCTGTaaaaa          >  1:270888/1‑143 (MQ=255)
                                                                                                                                          agccagGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACATGCTGTTGCAGAAAATCACGTCAAACCTCTTTATCAGTTAATCTGTAAAAATaaaa     <  1:279569/143‑1 (MQ=255)
                                                                                                                                             cagGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACATGCTGTTGCAGAAAATCACGTCAAACCTCTTTATCAGTTAATCTGTAAAAATAAAACCt  >  1:356070/1‑143 (MQ=255)
                                                                                                                                              |                                                                                                                                             
GAAAGGAAGGTTCCGACGGTCGAGGGTACAATGCGGAAGTGGTCCGCATGAAAAAACAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACATGCTGTTGCAGAAAATCACGTCAAACCTCTTTATCAGTTAATCTGTAAAAATAAAACCT  >  NC_000913/1498754‑1499037

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: