Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,498,892 A→T 8.2% intergenic (+17/‑46) 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,8920AT 8.2% 142.6 / 4.3 49intergenic (+17/‑46)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 = 7.57e-01

CCAGAAAGGAAGGTTCCGACGGTCGAGGGTACAATGCGGAAGTGGTCCGCATGAAAAAACAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACATGCTGTTGCAGAAAATCACGTCAAACCTCTTTATCAGTTAATCTGTAAAAATAAAA  >  NC_000913/1498751‑1499034
                                                                                                                                             |                                                                                                                                              
ccAGCAAGGAAGGTGCCGACGGTCGATGGTACAATGCGGAAGTGGTCCGCATGAAAAAACAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCgccag                                                                                                                                               <  2:434598/143‑1 (MQ=255)
   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                                                                                                                                    >  2:49267/1‑137 (MQ=255)
                 gACGGTCGAGGGTACAATGCGGAAGTGGTCCGCATGAAAAAACAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCTGGCTGGCGGCtt                                                                                                                                    <  1:49267/137‑1 (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)
                                                                                                                                             |                                                                                                                                              
CCAGAAAGGAAGGTTCCGACGGTCGAGGGTACAATGCGGAAGTGGTCCGCATGAAAAAACAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACATGCTGTTGCAGAAAATCACGTCAAACCTCTTTATCAGTTAATCTGTAAAAATAAAA  >  NC_000913/1498751‑1499034

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: