Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 508,600 T→G 7.8% intergenic (+41/+275) ybaQ → / ← copA putative DNA‑binding transcriptional regulator/copper transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913508,6000TG 7.8% 220.3 / 10.1 64intergenic (+41/+275)ybaQ/copAputative DNA‑binding transcriptional regulator/copper transporter
Reads supporting (aligned to +/- strand):  ref base T (29/30);  new base G (2/3);  total (31/33)
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.47e-01

TGGGAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATT  >  NC_000913/508455‑508730
                                                                                                                                                 |                                                                                                                                  
tGGGAAGTTGAAAACAACATGCGCGCCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATACGGCGTTCAc                                                                                                                                 >  3:1236824/1‑149 (MQ=255)
tGGGAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCAc                                                                                                                                 >  3:361420/1‑149 (MQ=255)
           aaaCAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCgcatccgg                                                                                                                      >  4:435833/1‑149 (MQ=255)
                                       ggACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCa                                                                                                            >  4:33161/1‑131 (MQ=255)
                                       ggACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCa                                                                                                            <  3:33161/131‑1 (MQ=255)
                                                 aaaCAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGa                                                                                                               >  2:974155/1‑118 (MQ=255)
                                                 aaaCAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGa                                                                                                               <  1:974155/118‑1 (MQ=255)
                                                       tGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCATAACATCGCTTCATCTCGGATGCCGGATGCGGCGTTAACGCCttatccggcctac                                                                                                                 <  4:1009697/110‑5 (MQ=21)
                                                       tGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGGATGCCGGATGCGGCGTGAACGCCttatccggcctac                                                                                                                 >  3:1009697/1‑106 (MQ=21)
                                                         gCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGc                                                                                                     >  3:857824/1‑120 (MQ=255)
                                                         gCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGc                                                                                                     <  4:857824/120‑1 (MQ=255)
                                                           tGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGcc                                                                      >  2:334132/1‑149 (MQ=255)
                                                           tGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGGATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAga                                                                                                             <  1:295900/110‑2 (MQ=18)
                                                                   tGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACt                                                              <  3:435833/149‑1 (MQ=255)
                                                                    ggCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACaaa                                                                <  3:88656/146‑1 (MQ=255)
                                                                    ggCACGCCGTGAAGAGCGCGCAAATAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACaaa                                                                >  4:88656/1‑146 (MQ=255)
                                                                       aCGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATGCGGCGTGAAcgccttat                                                                                                                          >  4:629925/1‑81 (MQ=37)
                                                                       aCGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATGCGGCGTGAAcgccttat                                                                                                                          <  3:629925/85‑5 (MQ=37)
                                                                       aCGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCgcatccgg                                                                                                                      >  4:469182/1‑89 (MQ=255)
                                                                       aCGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCgcatccgg                                                                                                                      <  3:469182/89‑1 (MQ=255)
                                                                             tGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTg                                                      <  3:132144/147‑1 (MQ=255)
                                                                             tGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTg                                                      >  4:132144/1‑147 (MQ=255)
                                                                              gAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAAc                                                   <  1:334132/149‑1 (MQ=255)
                                                                                 gagCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATGCGGCGTGAACGCCttatccg                                                                                                                       <  3:1086161/78‑1 (MQ=25)
                                                                                 gagCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGt                                                >  2:1290412/1‑149 (MQ=255)
                                                                                      tGCAAAAAAGGTCGCGTGAGGCCAAAACATCGCTTCATCTCGTAGGCCGTATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCATATCAGGCCTACAAACTTGTGCCTGAACCGTAGGcc                                           <  4:1236824/149‑1 (MQ=255)
                                                                                      tGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACaaa                                                                <  1:1320926/128‑1 (MQ=255)
                                                                                      tGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACaaa                                                                >  2:1320926/1‑128 (MQ=255)
                                                                                      tGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGcc                                           <  4:361420/149‑1 (MQ=255)
                                                                                         aaaaaaGTTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGa                                        >  2:801986/1‑149 (MQ=255)
                                                                                         aaaaaaGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGc                                                                                                             >  2:1005655/1‑80 (MQ=255)
                                                                                         aaaaaaGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGc                                                                                                             <  1:1005655/80‑1 (MQ=255)
                                                                                         aaaaaaGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGa                                        >  4:1438130/1‑149 (MQ=255)
                                                                                          aaaaaGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGc                                                                       >  2:163071/1‑117 (MQ=255)
                                                                                          aaaaaGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGc                                                                       <  1:163071/117‑1 (MQ=255)
                                                                                           aaaaGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATa                                      >  1:825316/1‑149 (MQ=255)
                                                                                            aaaGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATaa                                     <  2:825316/149‑1 (MQ=255)
                                                                                            aaaGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATaa                                     <  1:1290412/149‑1 (MQ=255)
                                                                                            aaaGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATaa                                     <  3:1438130/149‑1 (MQ=255)
                                                                                                gTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTg                                                                                                      <  2:157283/80‑1 (MQ=255)
                                                                                                gTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTg                                                                                                      >  1:157283/1‑80 (MQ=255)
                                                                                                     gTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTAc                            <  1:801986/149‑1 (MQ=255)
                                                                                                           ccAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTc                                                                                >  2:1111295/1‑91 (MQ=255)
                                                                                                           ccAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTc                                                                                <  1:1111295/91‑1 (MQ=255)
                                                                                                            cAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGt                                >  2:618871/1‑138 (MQ=255)
                                                                                                            cAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGt                                <  1:618871/138‑1 (MQ=255)
                                                                                                                  aTCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCAACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGt                                <  2:238176/132‑1 (MQ=255)
                                                                                                                  aTCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCAACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGt                                >  1:238176/1‑132 (MQ=255)
                                                                                                                          aTCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCa                                                                                                         <  3:34522/51‑1 (MQ=255)
                                                                                                                          aTCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCa                                                                                                         >  4:34522/1‑51 (MQ=255)
                                                                                                                           tctcGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCCGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAAc                                                   <  1:99437/104‑1 (MQ=255)
                                                                                                                           tctcGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCCGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAAc                                                   >  2:99437/1‑104 (MQ=255)
                                                                                                                              cGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGt                                >  2:620623/1‑120 (MQ=40)
                                                                                                                              cGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGt                                <  1:620623/120‑1 (MQ=39)
                                                                                                                                   gCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTa                                               >  3:826066/1‑100 (MQ=255)
                                                                                                                                   gCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTa                                               <  4:826066/100‑1 (MQ=255)
                                                                                                                                           aGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCAtt  >  2:270539/1‑137 (MQ=255)
                                                                                                                                           aGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCAtt  <  1:270539/137‑1 (MQ=255)
                                                                                                                                            ggCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTAt                                                                            <  2:1465392/62‑1 (MQ=255)
                                                                                                                                            ggCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTAt                                                                            >  1:1465392/1‑62 (MQ=255)
                                                                                                                                             gCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTtgt                                                           >  2:341859/1‑78 (MQ=255)
                                                                                                                                             gCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTtgt                                                           <  1:341859/78‑1 (MQ=255)
                                                                                                                                             gCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCgcatccg                  >  2:1088923/1‑119 (MQ=40)
                                                                                                                                             gCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCgcatccg                  <  1:1088923/119‑1 (MQ=40)
                                                                                                                                                 |                                                                                                                                  
TGGGAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATT  >  NC_000913/508455‑508730

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑