Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,498,894 C→G 8.4% intergenic (+19/‑44) 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,8940CG 8.4% 132.0 / 6.2 49intergenic (+19/‑44)ydcH/rimLprotein YdcH/ribosomal‑protein‑L12‑serine acetyltransferase
Reads supporting (aligned to +/- strand):  ref base C (25/19);  new base G (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

GAAAGGAAGGTTCCGACGGTCGAGGGTACAATGCGGAAGTGGTCCGCATGAAAAAACAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACATGCTGTTGCAGAAAATCACGTCAAACCTCTTTATCAGTTAATCTGTAAAAATAAAA  >  NC_000913/1498754‑1499034
                                                                                                                                            |                                                                                                                                            
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       >  1:270888/1‑143 (MQ=255)
                                                                                                                                     ccgccAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACATGCTGTTGCAGAAAATCACGTCAAACCTCTTTATCAGTTAATCTGTaaaaa       >  2:392892/1‑143 (MQ=255)
                                                                                                                                          agccagGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACATGCTGTTGCAGAAAATCACGTCAAACCTCTTTATCAGTTAATCTGTAAAAATaaaa  <  1:279569/143‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                            
GAAAGGAAGGTTCCGACGGTCGAGGGTACAATGCGGAAGTGGTCCGCATGAAAAAACAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACATGCTGTTGCAGAAAATCACGTCAAACCTCTTTATCAGTTAATCTGTAAAAATAAAA  >  NC_000913/1498754‑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: