Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 508,593 A→G 10.0% intergenic (+34/+282) 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,5930AG10.0% 188.4 / 10.0 60intergenic (+34/+282)ybaQ/copAputative DNA‑binding transcriptional regulator/copper transporter
Reads supporting (aligned to +/- strand):  ref base A (27/27);  new base G (2/4);  total (29/31)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.44e-02

GTTGATCTTTGGGAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTAC  >  NC_000913/508446‑508704
                                                                                                                                                   |                                                                                                               
gTTGATCTTTGGGAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATaa                                                                                                                >  1:1446290/1‑149 (MQ=255)
         tGGGAAGTTGAAAACAACATGCGCGCCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATACGGCGTTCAc                                                                                                       >  3:1236824/1‑149 (MQ=255)
         tGGGAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCAc                                                                                                       >  3:361420/1‑149 (MQ=255)
                    aaaCAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCgcatccgg                                                                                            >  4:435833/1‑149 (MQ=255)
                                      ggAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAg                                                                                                               >  3:1075601/1‑112 (MQ=255)
                                      ggAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAg                                                                                                               <  4:1075601/112‑1 (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)
                                                                  gCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGt                                                                                                           >  2:978307/1‑88 (MQ=255)
                                                                  gCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGt                                                                                                           <  1:978307/88‑1 (MQ=255)
                                                                  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                                                                                            <  3:469182/89‑1 (MQ=255)
                                                                                aCGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCgcatccgg                                                                                            >  4:469182/1‑89 (MQ=255)
                                                                                      tGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTg                            >  4:132144/1‑147 (MQ=255)
                                                                                      tGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTg                            <  3:132144/147‑1 (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                                      >  2:1320926/1‑128 (MQ=255)
                                                                                               tGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACaaa                                      <  1:1320926/128‑1 (MQ=255)
                                                                                               tGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGcc                 <  4:361420/149‑1 (MQ=255)
                                                                                                  aaaaaaGTTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGa              >  2:801986/1‑149 (MQ=255)
                                                                                                  aaaaaaGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGc                                                                                   <  1:1005655/80‑1 (MQ=255)
                                                                                                  aaaaaaGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGc                                                                                   >  2:1005655/1‑80 (MQ=255)
                                                                                                  aaaaaaGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGa              >  4:1438130/1‑149 (MQ=255)
                                                                                                   aaaaaGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGc                                             <  1:163071/117‑1 (MQ=255)
                                                                                                   aaaaaGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGc                                             >  2:163071/1‑117 (MQ=255)
                                                                                                    aaaaGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATa            >  1:825316/1‑149 (MQ=255)
                                                                                                     aaaGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATaa           <  3:1438130/149‑1 (MQ=255)
                                                                                                     aaaGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATaa           <  1:1290412/149‑1 (MQ=255)
                                                                                                     aaaGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATaa           <  2:825316/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                                                      <  1:1111295/91‑1 (MQ=255)
                                                                                                                    ccAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTc                                                      >  2:1111295/1‑91 (MQ=255)
                                                                                                                     cAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGt      <  1:618871/138‑1 (MQ=255)
                                                                                                                     cAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGt      >  2:618871/1‑138 (MQ=255)
                                                                                                                           aTCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCAACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGt      >  1:238176/1‑132 (MQ=255)
                                                                                                                           aTCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCAACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGt      <  2:238176/132‑1 (MQ=255)
                                                                                                                                   aTCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCa                                                                               >  4:34522/1‑51 (MQ=255)
                                                                                                                                   aTCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCa                                                                               <  3:34522/51‑1 (MQ=255)
                                                                                                                                    tctcGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCCGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAAc                         >  2:99437/1‑104 (MQ=255)
                                                                                                                                    tctcGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCCGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAAc                         <  1:99437/104‑1 (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)
                                                                                                                                                   |                                                                                                               
GTTGATCTTTGGGAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTAC  >  NC_000913/508446‑508704

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