Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 508,652 C→T intergenic (+93/+223) ybaQ → / ← copA DNA‑binding transcriptional regulator YbaQ/Cu(+) exporting P‑type ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913508,6520CT95.0% 51.5 / ‑3.8 20intergenic (+93/+223)ybaQ/copADNA‑binding transcriptional regulator YbaQ/Cu(+) exporting P‑type ATPase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (10/9);  minor base . (0/1);  total (10/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

CAGTCGATTTTTGGCTAAACCTCCAGGCGGCGGTTGATCTTTGGGAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGC  >  NC_000913/508414‑508856
                                                                                                                                                                                                                                               |                                                                                                                                                                                                            
caGTCGATTTTTGGCTAAACCGCCAGGCGGCGGTTGATCTTTGGGAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCGTAATATGTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTTTTATCAGGCCTa                                                                                                                                                                                                  <  1:168418/251‑1 (MQ=255)
             gCTAAACCTCCAGGCGGCGGGTGATCTTTGGGAAGTTGCAAACAACCTGCGCACCCAGGAAGACTTGGGACGGATTGAAACAGGGGCTGAATATTTGGCCCGCCGTTAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGCTAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTTTTCTCAAGCCTCCCAACTTGTGCCt                                                                                                                                                                                     >  1:73137/1‑251 (MQ=38)
                                  tGATCTTTGGGAAGTTGAAAACAACAGGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTTTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAgg                                                                                                                                                                 <  1:311667/250‑1 (MQ=255)
                                      cTTTGGGAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGGGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTTTTATCAGGCCTACAAACTtgtn                                                                                                                                                                                        <  1:49247/223‑2 (MQ=255)
                                                                                ggccGGATTGAAACATTGGCTGAATATTTGGCACGCCGTAAAGAGCGTGCAAAAAAGTTCGCGTAAGGCCAAAACATCGCTGCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTTTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGc                                                                                                                  <  1:98029/248‑1 (MQ=255)
                                                                                    ggATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTTTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGAAAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTAc                                                                                                              >  1:402720/1‑251 (MQ=255)
                                                                                            aCAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCTTCCGGCATCTGAGCACCAGTGCCGGATGCTACGCGTACCCGTTTTAGCAGGCCTACAAACTTGTGCCTGAACCGTAGTCTGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCCCCATTGCCTGATGCTACGCTTACGCGTCTTa                                                                                                      <  1:448542/251‑1 (MQ=37)
                                                                                             cAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATACGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTTTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAgg                                                                                                                                                                 >  1:227674/1‑191 (MQ=255)
                                                                                                               cACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTTTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGa                                                                                                                     >  1:144532/1‑217 (MQ=255)
                                                                                                                                 naaaaaaGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTTTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGcc                                                                                                                                                                         >  1:89221/2‑147 (MQ=255)
                                                                                                                                     aaaGGTCGNGNNNGNNCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCNCATCCGGCATCCGAGNACNATTGNCTGNTGCGACGCTTACCCGTTTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCgcatccgg                                                                                                                                               >  2:316466/1‑169 (MQ=255)
                                                                                                                                     aaaGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCTGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTTTTAGCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCgcatccgg                                                                                                                                               <  1:316466/169‑1 (MQ=255)
                                                                                                                                            gcgTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTTTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGNCCGNNNNNggnn                                                                                                                                                               >  2:272910/1‑144 (MQ=255)
                                                                                                                                            gcgTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCAAGCACCATTGCCTAATGCGACGCTTACCCGTTTTATCAGGCCTACAAACTTGTGACTGAACCGTCGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGc                                                                                                                                      >  1:245240/1‑171 (MQ=18)
                                                                                                                                             cgTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATACGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTTTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATCCGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTa                                                                                             >  1:72371/1‑211 (MQ=255)
                                                                                                                                                                                gATAAGGCGGTCACGCCGCGTCCGTCATCCGAGCACCATTGCCTGATGCGACTCTTACCCAGTTTATCAGGCCTACAAAGTTGTGCCTGAACCGTAGGCCGGATAAGGCTTTTACGTCGCATCCGGCATCCGAGCACCACGGCCTGATGCGGCGCTTACTCGTCTTATCAGGCCTGCAAGCCTCAGCCTGAACCGTTGGCCAGATAAGGTGGTGACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATgc                  <  1:384942/251‑1 (MQ=21)
                                                                                                                                                                                   aaGGCGTTCACTCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTGACCCGTTTTATCAGGCCTACAAACTTGTGCCTGATCCGTATGCCGGATAAGGCTTTGACGCCGCATCCGGCAGCCGATCACCATTGCCTGATGCGACTCTTACGCGTCTTATCAGGCCTACAAACCGGTGCCTGAACCTTAGGCGGGATAAGGCGTTGACGCCGCATCCGGCAGTCATGCTTCGATGCCAAATgcgc                <  1:65033/250‑1 (MQ=25)
                                                                                                                                                                                              gccgcATCAGGCATCCGAGCACCATTGACTGATGCGAGGCTTACCCGTTTTAACAGGCCTACAAACTTGTACCTGACCCGGAGGCCGGNNAAGGCGnnt                                                                                                                                                            >  2:34214/1‑99 (MQ=14)
                                                                                                                                                                                               ccgcATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTTTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGAt                           <  1:152210/227‑1 (MQ=255)
                                                                                                                                                                                                cgcATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTTTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAagcagc  <  1:34214/251‑1 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                            
CAGTCGATTTTTGGCTAAACCTCCAGGCGGCGGTTGATCTTTGGGAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGC  >  NC_000913/508414‑508856

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